Estudo numérico de escoamento turbulento em padrão anular gás-líquido em dutos verticais

Multiphase flows in ducts can adopt several morphologies depending on the mass fluxes and the fluids properties. Annular flow is one of the most frequently encountered flow patterns in industrial applications. For gas liquid systems, it consists of a liquid film flowing adjacent to the wall and a gas core flowing in the center of the duct. This work presents a numerical study of this flow pattern in gas liquid systems in vertical ducts. For this, a solution algorithm was developed and implemented in FORTRAN 90 to numerically solve the governing transport equations. The mass and momentum conservation equations are solved simultaneously from the wall to the center of the duct, using the Finite Volumes Technique. Momentum conservation in the gas liquid interface is enforced using an equivalent effective viscosity, which also allows for the solution of both velocity fields in a single system of equations. In this way, the velocity distributions across the gas core and the liquid film are obtained iteratively, together with the global pressure gradient and the liquid film thickness. Convergence criteria are based upon satisfaction of mass balance within the liquid film and the gas core. For system closure, two different approaches are presented for the calculation of the radial turbulent viscosity distribution within the liquid film and the gas core. The first one combines a k- Ɛ one-equation model and a low Reynolds k-Ɛ model. The second one uses a low Reynolds k- Ɛ model to compute the eddy viscosity profile from the center of the duct right to the wall. Appropriate interfacial values for k e Ɛ are proposed, based on concepts and ideas previously used, with success, in stratified gas liquid flow. The proposed approaches are compared with an algebraic model found in the literature, specifically devised for annular gas liquid flow, using available experimental results. This also serves as a validation of the solution algorithm

Modelagem numérica de escoamento em padrão anular com transferência de calor e massa

Annular flow is the prevailing pattern in transport and energy conversion systems and therefore, one of the most important patterns in multiphase flow in ducts. The correct prediction of the pressure gradient and heat transfer coefficient is essential for optimizing the system s capacity. The objective of this work is to develop and implement a numerical algorithm capable of predicting hydrodynamic and thermal characteristics for upflow, vertical, annular flow. The numerical algorithm is then complemented with the physical modeling of phenomena that occurs in this flow pattern. These are, turbulence, entrainment and deposition and phase change. For the development of the numerical model, axial diffusion of heat and momentum is neglected. In this way the time-averaged equations are solved in their parabolic form obtaining the velocity and temperature profiles for each axial step at a time, together with the global parameters, namely, pressure gradient, mean film thickness and heat transfer coefficient, as well as their variation in the axial direction. The model is validated for the following conditions: fully-developed laminar flow with no entrainment; fully developed laminar flow with heat transfer, fully-developed turbulent flow with entrained drops, developing turbulent annular flow with entrained drops, and turbulent flow with heat transfer and phase change

Influência da vermiculita expandida na resistência à compressão e reologia de pastas cimentantes de poços de petróleo

Several problems related to the loss of hydraulic seal in oilwells, causing gas migration and/or contamination of the production zone by water, have been reported. The loss of the hydraulic seal is a consequence of cracks which can be occasioned either by the invasion of gas during the wait on cement or by the expansion of the casing causing the fracture of the cement sheath. In case of the pressure of the formation is higher than the pressure in the annulus, gas can migrate into the slurry and form microannulus, which are channels where gas migrates after the cement is set. Cracks can be also occasioned by the fracture of the cement sheath when it does not withstand the thermal and dynamic loads. In reservoirs where the oil is heavy, steam water injection operation is required in order to get the oil flowing. This operation increases the temperature of the casing, and then it expands and causes the fracture of the cement sheath in the annulus. When the failures on the cement are detected, remedial cementing is required, which raise costs caused by the interventions. Once the use of cement in the construction civil sector is older than its use in the petroleum sector, it is common to bring technologies and solutions from the civil construction and apply them on the petroleum area. In this context, vermiculite, a mineral-clay widely encountered in Brazil, has been used, on its exfoliated form, in the civil construction, especially on the manufacture of lights and fireproof concretes with excellent thermal and acoustical properties. It has already been reported in scientific journals, studies of the addition of exfoliated vermiculite in Portland cements revealing good properties related to oilwell cementing operations. Thus, this study aimed to study the rheological behavior, thickening time, stability and compressive strength of the slurries made of Portland cement and exfoliated vermiculite in 5 different compositions, at room temperature and heated. The results showed that the compressive strength decreased with the addition of exfoliated vermiculite, however the values are still allowed for oiwell cementing operations. The thickening time of the slurry with no exfoliated vermiculite was 120 min and the thickening time of the slurry with 12 % of exfoliated vermiculite was 98 min. The stability and the rheological behavior of the slurries revealed that the exfoliated vermiculite absorbed water and therefore increased the viscosity of the slurries, even though increasing the factor cement-water. The stability experiment carried out at 133 ºF showed that, there was neither sedimentation nor reduction of the volume of the cement for the slurry with 12 % of exfoliated vermiculite. Thus, the addition of exfoliated vermiculite accelerates the set time of the cement and gives it a small shrinkage during the wait on cement, which are important to prevent gas migration

Avaliação do comportamento reológico de pastas de cimento para poços de petróleo com adição de plastificantes

The isolation of adjacent zones encountered during oilwell drilling is carried out by Portland-based cement slurries. The slurries are pumped into the annular positions between the well and the casing. Their rheological behavior is a very important component for the cementing process. Nowadays, several alternative materials are used in oilwell cementing, with goal the modification and the improvement of their properties, mainly the increase of the fluidity. And this can be reached by using plasticizers additives able to account for different oilwell conditions, yielding compatible cement slurries and allowing enough time for the complete cementing operation. If the rheological properties of the slurry are properly characterized, the load loss and flow regime can be correctly predicted. However, this experimental characterization is difficult. Rheological models capable of describing the cement slurry behavior must be capable of predicting the slurry cement deformation within reasonable accuracy. The aim of this study was to characterize rheologically the slurries prepared with a especial class of Portland cement, water and plasticizers based on lignosulfonate, melamine and polycarboxylate at temperatures varying from 27°C to 72°C. The tests were carried out according to the practical recommendations of the API RP 10B guidelines. The results revealed a great efficiency and the dispersive power of the polycarboxylate, for all temperatures tested. This additive promoted high fluidity of the slurries, with no sedimentation. High lignosulfonate and melamine concentrations did not reduce the rheological parameters (plastic viscosity and yield stress) of the slurries. It was verified that these additives were not compatible with the type of cement used. The evaluated rheological models were capable of describing the behavior of the slurries only within concentration and temperature ranges specific for each type of additive

Simulador de escoamento multifásico em poços de petróleo (SEMPP)

The multiphase flow occurrence in the oil and gas industry is common throughout fluid path, production, transportation and refining. The multiphase flow is defined as flow simultaneously composed of two or more phases with different properties and immiscible. An important computational tool for the design, planning and optimization production systems is multiphase flow simulation in pipelines and porous media, usually made by multiphase flow commercial simulators. The main purpose of the multiphase flow simulators is predicting pressure and temperature at any point at the production system. This work proposes the development of a multiphase flow simulator able to predict the dynamic pressure and temperature gradient in vertical, directional and horizontal wells. The prediction of pressure and temperature profiles was made by numerical integration using marching algorithm with empirical correlations and mechanistic model to predict pressure gradient. The development of this tool involved set of routines implemented through software programming Embarcadero C++ Builder® 2010 version, which allowed the creation of executable file compatible with Microsoft Windows® operating systems. The simulator validation was conduct by computational experiments and comparison the results with the PIPESIM®. In general, the developed simulator achieved excellent results compared with those obtained by PIPESIM and can be used as a tool to assist production systems development

Sistema inteligente para apoio à decisão na operação de uma malha de escoamento de petróleo

The petroleum production pipeline networks are inherently complex, usually decentralized systems. Strict operational constraints are applied in order to prevent serious problems like environmental disasters or production losses. This paper describes an intelligent system to support decisions in the operation of these networks, proposing a staggering for the pumps of transfer stations that compose them. The intelligent system is formed by blocks which interconnect to process the information and generate the suggestions to the operator. The main block of the system uses fuzzy logic to provide a control based on rules, which incorporate knowledge from experts. Tests performed in the simulation environment provided good results, indicating the applicability of the system in a real oil production environment. The use of the stagger proposed by the system allows a prioritization of the transfer in the network and a flow programming

Simulação do escoamento monofásico em um estágio de uma bomba centrífuga utilizando técnicas de fluidodinâmica computacional

Oil production and exploration techniques have evolved in the last decades in order to increase fluid flows and optimize how the required equipment are used. The base functioning of Electric Submersible Pumping (ESP) lift method is the use of an electric downhole motor to move a centrifugal pump and transport the fluids to the surface. The Electric Submersible Pumping is an option that has been gaining ground among the methods of Artificial Lift due to the ability to handle a large flow of liquid in onshore and offshore environments. The performance of a well equipped with ESP systems is intrinsically related to the centrifugal pump operation. It is the pump that has the function to turn the motor power into Head. In this present work, a computer model to analyze the three-dimensional flow in a centrifugal pump used in Electric Submersible Pumping has been developed. Through the commercial program, ANSYS® CFX®, initially using water as fluid flow, the geometry and simulation parameters have been defined in order to obtain an approximation of what occurs inside the channels of the impeller and diffuser pump in terms of flow. Three different geometry conditions were initially tested to determine which is most suitable to solving the problem. After choosing the most appropriate geometry, three mesh conditions were analyzed and the obtained values were compared to the experimental characteristic curve of Head provided by the manufacturer. The results have approached the experimental curve, the simulation time and the model convergence were satisfactory if it is considered that the studied problem involves numerical analysis. After the tests with water, oil was used in the simulations. The results were compared to a methodology used in the petroleum industry to correct viscosity. In general, for models with water and oil, the results with single-phase fluids were coherent with the experimental curves and, through three-dimensional computer models, they are a preliminary evaluation for the analysis of the two-phase flow inside the channels of centrifugal pump used in ESP systems

Desenvolvimento de um simulador de bombeio por cavidades progressivas

The method of artificial lift of progressing cavity pump is very efficient in the production of oils with high viscosity and oils that carry a great amount of sand. This characteristic converted this lift method into the second most useful one in oil fields production. As it grows the number of its applications it also increases the necessity to dominate its work in a way to define it the best operational set point. To contribute to the knowledge of the operational method of artificial lift of progressing cavity pump, this work intends to develop a computational simulator for oil wells equipped with an artificial lift system. The computational simulator of the system will be able to represent its dynamic behavior when submitted to the various operational conditions. The system was divided into five subsystems: induction motor, multiphase flows into production tubing, rod string, progressing cavity pump and annular tubing-casing. The modeling and simulation of each subsystem permitted to evaluate the dynamic characteristics that defined the criteria connections. With the connections of the subsystems it was possible to obtain the dynamic characteristics of the most important arrays belonging to the system, such as: pressure discharge, pressure intake, pumping rate, rod string rotation and torque applied to polish string. The shown results added to a friendly graphical interface converted the PCP simulator in a great potential tool with a didactic characteristic in serving the technical capability for the system operators and also permitting the production engineering to achieve a more detail analysis of the dynamic operational oil wells equipped with the progressing cavity pump

Dimensionamento de sistemas de aterramento em linhas de transmissão para escoamento eficiente de correntes de surtos atmosféricos

This work presents in a simulated environment, to analyze the length of cable needed counterweight connected to ground rod, able to avoid the phenomenon of flashover return, back flashover, the insulator chains of transmission lines consisting of concrete structures when they are subjected to lightning standardized regarding certain resistivity values of some kinds of soil and geometric arrangements of disposal of grounding systems structures

Simulação computacional do escoamento em bombas de cavidades progressivas

The use of Progressing Cavity Pumps (PCPs) in artificial lift applications in low deep wells is becoming more common in the oil industry, mainly, due to its ability to pump heavy oils, produce oil with large concentrations of sand, besides present high efficiency when compared to other artificial lift methods. Although this system has been widely used as an oil lift method, few investigations about its hydrodynamic behavior are presented, either experimental or numeric. Therefore, in order to increase the knowledge about the BCP operational behavior, this work presents a novel computational model for the 3-D transient flow in progressing cavity pumps, which includes the relative motion between rotor and stator, using an element based finite volume method. The model developed is able to accurately predict the volumetric efficiency and viscous looses as well as to provide detailed information of pressure and velocity fields inside the pump. In order to predict PCP performance for low viscosity fluids, advanced turbulence models were used to treat, accurately, the turbulent effects on the flow, which allowed for obtaining results consistent with experimental values encountered in literature. In addition to the 3D computational model, a simplified model was developed, based on mass balance within cavities and on simplification on the momentum equations for fully developed flow along the seal region between cavities. This simplified model, based on previous approaches encountered in literature, has the ability to predict flow rate for a given differential pressure, presenting exactness and low CPU requirements, becoming an engineering tool for quick calculations and providing adequate results, almost real-time time. The results presented in this work consider a rigid stator PCP and the models developed were validated against experimental results from open literature. The results for the 3-D model showed to be sensitive to the mesh size, such that a numerical mesh refinement study is also presented. Regarding to the simplified model, some improvements were introduced in the calculation of the friction factor, allowing the application fo the model for low viscosity fluids, which was unsuccessful in models using similar approaches, presented in previous works

Simulação do escoamento multifásico no interior de bombas de cavidades progressivas metálicas

The progressing cavity pumping (PCP) is one of the most applied oil lift methods nowadays in oil extraction due to its ability to pump heavy and high gas fraction flows. The computational modeling of PCPs appears as a tool to help experiments with the pump and therefore, obtain precisely the pump operational variables, contributing to pump s project and field operation otimization in the respectively situation. A computational model for multiphase flow inside a metallic stator PCP which consider the relative motion between rotor and stator was developed in the present work. In such model, the gas-liquid bubbly flow pattern was considered, which is a very common situation in practice. The Eulerian-Eulerian approach, considering the homogeneous and inhomogeneous models, was employed and gas was treated taking into account an ideal gas state. The effects of the different gas volume fractions in pump volumetric eficiency, pressure distribution, power, slippage flow rate and volumetric flow rate were analyzed. The results shown that the developed model is capable of reproducing pump dynamic behaviour under the multiphase flow conditions early performed in experimental works

Modelagem dinâmica do escoamento de um sistema de elevação por plunger lift

A critical problem in mature gas wells is the liquid loading. As the reservoir pressure decreases, gas superficial velocities decreases and the drag exerted on the liquid phase may become insufficient to bring all the liquid to the surface. Liquid starts to drain downward, flooding the well and increasing the backpressure which decreases the gas superficial velocity and so on. A popular method to remedy this problem is the Plunger Lift. This method consists of dropping the "plunger"to the bottom of the tubing well with the main production valve closed. When the plunger reaches the well bottom the production valve is opened and the plunger carry the liquid to the surface. However, models presented in literature for predicting the behavior in plunger lift are simplistic, in many cases static (not considering the transient effects). Therefore work presents the development and validation of a numerical algorithm to solve one-dimensional compressible in gas wells using the Finite Volume Method and PRIME techniques for treating coupling of pressure and velocity fields. The code will be then used to develop a dynamic model for the plunger lift which includes the transient compressible flow within the well

Degradação do fungicida tiofanato metílico pelo processo foto-Fenton em reatores anular e solar

Effluents from pesticide industries have great difficulty to decontaminate the environment and, moreover, are characterized by high organic charge and toxicity. The research group Center for Chemical Systems Engineering (CESQ) at the Department of Chemical Engineering of Polytechnical School of University of São Paulo and Department of Chemical Engineering, Federal University of Rio Grande do Norte have been applying the Advanced Oxidation Processes (AOP's) for the degradation of various types of pollutants. These processes are based on the generation of hydroxyl radicals, highly reactive substances. Thus, this dissertation aims to explore this process, since it has been proven to be quite effective in removing organic charge. Therefore, it was decided by photo-Fenton process applied to the degradation of the fungicide Thiophanate methyl in aqueous system using annular reactor (with lamp Philips HPLN 125W) and solar. The samples were collected during the experiment and analyzed for dissolved organic carbon (TOC) using a Shimadzu TOC (Shimadzu 5050A e VCP). The Doehlert experimental design has been used to evaluate the influence of ultraviolet radiation, the concentrations of methyl thiophanate (C12H14N4O4S2), hydrogen peroxide (H2O2) and iron ions (Fe2+), among these parameters, was considered the best experimental conditions, [Fe2+] = 0.6 mmol/L and [H2O2] = 0.038 mol/L in EXP 5 experiment and in SOL 5 experiment, obtaining a percentage of TOC removal of 60% in the annular reactor and 75% in the solar reactor

Modelagem e simulação da deposição de parafinas em escoamento turbulento

The flow assurance has become one of the topics of greatest interest in the oil industry, mainly due to production and transportation of oil in regions with extreme temperature and pressure. In these operations the wax deposition is a commonly problem in flow of paraffinic oils, causing the rising costs of the process, due to increased energy cost of pumping, decreased production, increased pressure on the line and risk of blockage of the pipeline. In order to describe the behavior of the wax deposition phenomena in turbulent flow of paraffinic oils, under different operations conditions, in this work we developed a simulator with easy interface. For that we divided de work in four steps: (i) properties estimation (physical, thermals, of transport and thermodynamics) of n-alkanes and paraffinic mixtures by using correlations; (ii) obtainment of the solubility curve and determination the wax appearance temperature, by calculating the solid-liquid equilibrium of parafinnic systems; (iii) modelling wax deposition process, comprising momentum, mass and heat transfer; (iv) development of graphic interface in MATLAB® environment for to allow the understanding of simulation in different flow conditions as well as understand the matter of the variables (inlet temperature, external temperature, wax appearance temperature, oil composition, and time) on the behavior of the deposition process. The results showed that the simulator developed, called DepoSim, is able to calculate the profile of temperature, thickness of the deposit, and the amount of wax deposited in a simple and fast way, and also with consistent results and applicable to the operation

Variação da qualidade da água do escoamento superficial de duas bacias de drenagem de Natal/RN - Brasil

The urban drainage is one of the powers of environmental sanitation and its scope is the quantitative and qualitative aspects. In decision making of managers and the engineering aspects of design are almost always taken into account only the quantitative aspects. However, the waters of the runoff have the highest concentrations of pollutants at the beginning of precipitation. Thus, if the plot pollution removed, the remaining portion can be used for other purposes. This work has aimed to present the variation of water quality of two drainage basins in the city of Natal / RN-Brazil to support the implementation of drainage to consider the qualitative aspect, and identify potential for the use of water. The basins (M and C) are analyzed closed-type, are in the urban area, are predominantly residential occupation and its waters are used for detention ponds and infiltration. The samples were divided into three phases, the first two direct to final points in a basin and the third in traps distributed over the surface drainage. The parameters had been analyzed were pH, conductivity, dissolved oxygen, Color, Turbidity, COD, Ammonia, nitrite, nitrate, total phosphorus, orthophosphate, Sediments solids, total solids, chloride, sulfate, alkalinity, calcium, magnesium, sodium, potassium, Heavy Metals (Chromium, Cadmium, Lead, Zinc and Copper), Eschichia coli and total coliforms. The parameters studied showed high initial pollution load, events and located in different proportions, except nitrite, heavy metals and biological indicators. The size of the surface drainage and topographic its features influence the quality of water. However, the form of sampling is crucial in the qualitative study in the basin. The samplers developed at work, were generated economic and representative results. The urban rainwater presents organic faecal indicators. The runoff of water from both basins shows no risk of salinity and sodicity for use in irrigation, should be noted the content of chloride in the choice of method of irrigation