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

Projeto de um sistema emulador de escoamentos e vaso de separação primária

The purpose of this study was to develop a pilot plant which the main goal is to emulate a flow peak pressure in a separation vessel. Effect similar that is caused by the production in a slug flow in production wells equipped with the artificial lift method plunger lift. The motivation for its development was the need to test in a plant on a smaller scale, a new technique developed to estimate the gas flow in production wells equipped with plunger lift. To develop it, studies about multiphase flow effects, operation methods of artificial lift in plunger lift wells, industrial instrumentation elements, control valves, vessel sizing separators and measurement systems were done. The methodology used was the definition of process flowcharts, its parameters and how the effects needed would be generated for the success of the experiments. Therefore, control valves, the design and construction of vessels and the acquisition of other equipment used were defined. One of the vessels works as a tank of compressed air that is connected to the separation vessel and generates pulses of gas controlled by a on/off valve. With the emulator system ready, several control experiments were made, being the control of peak flow pressure generation and the flow meter the main experiments, this way, it was confirmed the efficiency of the plant usage in the problem that motivated it. It was concluded that the system is capable of generate effects of flow with peak pressure in a primary separation vessel. Studies such as the estimation of gas flow at the exit of the vessel and several academic studies can be done and tested on a smaller scale and then applied in real plants, avoiding waste of time and money.

Sistema inteligente para detecção de vazamentos em dutos de petróleo usando transformada Wavelet e redes neurais

This work consists in the use of techniques of signals processing and artificial neural networks to identify leaks in pipes with multiphase flow. In the traditional methods of leak detection exists a great difficulty to mount a profile, that is adjusted to the found in real conditions of the oil transport. These difficult conditions go since the unevenly soil that cause columns or vacuum throughout pipelines until the presence of multiphases like water, gas and oil; plus other components as sand, which use to produce discontinuous flow off and diverse variations. To attenuate these difficulties, the transform wavelet was used to map the signal pressure in different resolution plan allowing the extraction of descriptors that identify leaks patterns and with then to provide training for the neural network to learning of how to classify this pattern and report whenever this characterize leaks. During the tests were used transient and regime signals and pipelines with punctures with size variations from ½' to 1' of diameter to simulate leaks and between Upanema and Estreito B, of the UN-RNCE of the Petrobras, where it was possible to detect leaks. The results show that the proposed descriptors considered, based in statistical methods applied in domain transform, are sufficient to identify leaks patterns and make it possible to train the neural classifier to indicate the occurrence of pipeline leaks

Controlador preditivo não linear aplicado ao controle de golfadas em processos de produção de petróleo

Slugging is a well-known slugging phenomenon in multiphase flow, which may cause problems such as vibration in pipeline and high liquid level in the separator. It can be classified according to the place of its occurrence. The most severe, known as slugging in the riser, occurs in the vertical pipe which feeds the platform. Also known as severe slugging, it is capable of causing severe pressure fluctuations in the flow of the process, excessive vibration, flooding in separator tanks, limited production, nonscheduled stop of production, among other negative aspects that motivated the production of this work . A feasible solution to deal with this problem would be to design an effective method for the removal or reduction of the system, a controller. According to the literature, a conventional PID controller did not produce good results due to the high degree of nonlinearity of the process, fueling the development of advanced control techniques. Among these, the model predictive controller (MPC), where the control action results from the solution of an optimization problem, it is robust, can incorporate physical and /or security constraints. The objective of this work is to apply a non-conventional non-linear model predictive control technique to severe slugging, where the amount of liquid mass in the riser is controlled by the production valve and, indirectly, the oscillation of flow and pressure is suppressed, while looking for environmental and economic benefits. The proposed strategy is based on the use of the model linear approximations and repeatedly solving of a quadratic optimization problem, providing solutions that improve at each iteration. In the event where the convergence of this algorithm is satisfied, the predicted values of the process variables are the same as to those obtained by the original nonlinear model, ensuring that the constraints are satisfied for them along the prediction horizon. A mathematical model recently published in the literature, capable of representing characteristics of severe slugging in a real oil well, is used both for simulation and for the project of the proposed controller, whose performance is compared to a linear MPC

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

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

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

Uso de técnicas de otimização baseadas em derivadas como suporte do planejamento operacional de redes de distribuição de energia elétrica

The usual programs for load flow calculation were in general developped aiming the simulation of electric energy transmission, subtransmission and distribution systems. However, the mathematical methods and algorithms used by the formulations were based, in majority, just on the characteristics of the transmittion systems, which were the main concern focus of engineers and researchers. Though, the physical characteristics of these systems are quite different from the distribution ones. In the transmission systems, the voltage levels are high and the lines are generally very long. These aspects contribute the capacitive and inductive effects that appear in the system to have a considerable influence in the values of the interest quantities, reason why they should be taken into consideration. Still in the transmission systems, the loads have a macro nature, as for example, cities, neiborhoods, or big industries. These loads are, generally, practically balanced, what reduces the necessity of utilization of three-phase methodology for the load flow calculation. Distribution systems, on the other hand, present different characteristics: the voltage levels are small in comparison to the transmission ones. This almost annul the capacitive effects of the lines. The loads are, in this case, transformers, in whose secondaries are connected small consumers, in a sort of times, mono-phase ones, so that the probability of finding an unbalanced circuit is high. This way, the utilization of three-phase methodologies assumes an important dimension. Besides, equipments like voltage regulators, that use simultaneously the concepts of phase and line voltage in their functioning, need a three-phase methodology, in order to allow the simulation of their real behavior. For the exposed reasons, initially was developped, in the scope of this work, a method for three-phase load flow calculation in order to simulate the steady-state behaviour of distribution systems. Aiming to achieve this goal, the Power Summation Algorithm was used, as a base for developing the three phase method. This algorithm was already widely tested and approved by researchers and engineers in the simulation of radial electric energy distribution systems, mainly for single-phase representation. By our formulation, lines are modeled in three-phase circuits, considering the magnetic coupling between the phases; but the earth effect is considered through the Carson reduction. It s important to point out that, in spite of the loads being normally connected to the transformer s secondaries, was considered the hypothesis of existence of star or delta loads connected to the primary circuit. To perform the simulation of voltage regulators, a new model was utilized, allowing the simulation of various types of configurations, according to their real functioning. Finally, was considered the possibility of representation of switches with current measuring in various points of the feeder. The loads are adjusted during the iteractive process, in order to match the current in each switch, converging to the measured value specified by the input data. In a second stage of the work, sensibility parameters were derived taking as base the described load flow, with the objective of suporting further optimization processes. This parameters are found by calculating of the partial derivatives of a variable in respect to another, in general, voltages, losses and reactive powers. After describing the calculation of the sensibility parameters, the Gradient Method was presented, using these parameters to optimize an objective function, that will be defined for each type of study. The first one refers to the reduction of technical losses in a medium voltage feeder, through the installation of capacitor banks; the second one refers to the problem of correction of voltage profile, through the instalation of capacitor banks or voltage regulators. In case of the losses reduction will be considered, as objective function, the sum of the losses in all the parts of the system. To the correction of the voltage profile, the objective function will be the sum of the square voltage deviations in each node, in respect to the rated voltage. In the end of the work, results of application of the described methods in some feeders are presented, aiming to give insight about their performance and acuity

Uso de técnicas de otimização baseadas em derivadas como suporte do planejamento operacional de redes de distribuição de energia elétrica

The usual programs for load flow calculation were in general developped aiming the simulation of electric energy transmission, subtransmission and distribution systems. However, the mathematical methods and algorithms used by the formulations were based, in majority, just on the characteristics of the transmittion systems, which were the main concern focus of engineers and researchers. Though, the physical characteristics of these systems are quite different from the distribution ones. In the transmission systems, the voltage levels are high and the lines are generally very long. These aspects contribute the capacitive and inductive effects that appear in the system to have a considerable influence in the values of the interest quantities, reason why they should be taken into consideration. Still in the transmission systems, the loads have a macro nature, as for example, cities, neiborhoods, or big industries. These loads are, generally, practically balanced, what reduces the necessity of utilization of three-phase methodology for the load flow calculation. Distribution systems, on the other hand, present different characteristics: the voltage levels are small in comparison to the transmission ones. This almost annul the capacitive effects of the lines. The loads are, in this case, transformers, in whose secondaries are connected small consumers, in a sort of times, mono-phase ones, so that the probability of finding an unbalanced circuit is high. This way, the utilization of three-phase methodologies assumes an important dimension. Besides, equipments like voltage regulators, that use simultaneously the concepts of phase and line voltage in their functioning, need a three-phase methodology, in order to allow the simulation of their real behavior. For the exposed reasons, initially was developped, in the scope of this work, a method for three-phase load flow calculation in order to simulate the steady-state behaviour of distribution systems. Aiming to achieve this goal, the Power Summation Algorithm was used, as a base for developping the three phase method. This algorithm was already widely tested and approved by researchers and engineers in the simulation of radial electric energy distribution systems, mainly for single-phase representation. By our formulation, lines are modeled in three-phase circuits, considering the magnetic coupling between the phases; but the earth effect is considered through the Carson reduction. Its important to point out that, in spite of the loads being normally connected to the transformers secondaries, was considered the hypothesis of existence of star or delta loads connected to the primary circuit. To perform the simulation of voltage regulators, a new model was utilized, allowing the simulation of various types of configurations, according to their real functioning. Finally, was considered the possibility of representation of switches with current measuring in various points of the feeder. The loads are adjusted during the iteractive process, in order to match the current in each switch, converging to the measured value specified by the input data. In a second stage of the work, sensibility parameters were derived taking as base the described load flow, with the objective of suporting further optimization processes. This parameters are found by calculating of the partial derivatives of a variable in respect to another, in general, voltages, losses and reactive powers. After describing the calculation of the sensibility parameters, the Gradient Method was presented, using these parameters to optimize an objective function, that will be defined for each type of study. The first one refers to the reduction of technical losses in a medium voltage feeder, through the installation of capacitor banks; the second one refers to the problem of correction of voltage profile, through the instalation of capacitor banks or voltage regulators. In case of the losses reduction will be considered, as objective function, the sum of the losses in all the parts of the system. To the correction of the voltage profile, the objective function will be the sum of the square voltage deviations in each node, in respect to the rated voltage. In the end of the work, results of application of the described methods in some feeders are presented, aiming to give insight about their performance and acuity

Análise da eficiência de aquecedor solar espiral com e sem efeito estufa para piscinas

The technical and economic viability of solar heating for swimming pools is unquestionable, besides there it replaces the high costs and environmental impacts of conventional supply of energy, and it improves an optimization in the pool heating uses. This work applies the principles of the greenhouse effect: advanced thermodynamics, heat retention and equalization of temperature, to optimize the solar heating equipment, reducing the area required by collectors as much as 40% (still estimated value) for commercial collectors, with minor architectural and aesthetic impacts on the environment. It features a solar heating alternative in pools, whose main characteristics: low cost, simplicity in manufacturing and assembly and a faster heating. The system consists of two collectors spiral hoses made of polyethylene with a hundred meters each, and working on a forced flow, with only one pass of the working fluid inside the coils, and is used to pump itself treatment of pool water to obtain the desired flow. One of the collectors will be exposed to direct solar radiation, and the other will be covered by a glass slide and closed laterally, so providing the greenhouse effect. The equipment will be installed in parallel and simultaneously exposed to the sun in order to obtain comparative data on their effectiveness. Will be presented results of thermal tests for this the two cases, with and without transparent cover. Will be demonstrated, by comparison, the thermal, economic and material feasibility of these systems for heating swimming pools.

Análise da magnetohidrodinâmica com transferência de calor em canais de placas paralelas via transformação integral

The main goal of the present work is related to the dynamics of the steady state, incompressible, laminar flow with heat transfer, of an electrically conducting and Newtonian fluid inside a flat parallel-plate channel under the action of an external and uniform magnetic field. For solution of the governing equations, written in the parabolic boundary layer and stream-function formulation, it was employed the hybrid, numericalanalytical, approach known as Generalized Integral Transform Technique (GITT). The flow is sustained by a pressure gradient and the magnetic field is applied in the direction normal to the flow and is assumed that normal magnetic field is kept uniform, remaining larger than any other fields generated in other directions. In order to evaluate the influence of the applied magnetic field on both entrance regions, thermal and hydrodynamic, for this forced convection problem, as well as for validating purposes of the adopted solution methodology, two kinds of channel entry conditions for the velocity field were used: an uniform and an non-MHD parabolic profile. On the other hand, for the thermal problem only an uniform temperature profile at the channel inlet was employed as boundary condition. Along the channel wall, plates are maintained at constant temperature, either equal to or different from each other. Results for the velocity and temperature fields as well as for the main related potentials are produced and compared, for validation purposes, to results reported on literature as function of the main dimensionless governing parameters as Reynolds and Hartman numbers, for typical situations. Finally, in order to illustrate the consistency of the integral transform method, convergence analyses are also effectuated and presented

Utilização de uma coluna de flotação para remoção de compostos orgânicos da água de produção

Innovative technologies using surfactant materials have applicability in several industrial fields, including petroleum and gas areas. This study seeks to investigate the use of a surfactant derived from coconut oil (SCO saponified coconut oil) in the recovery process of organic compounds that are present in oily effluents from petroleum industry. For this end, experiments were accomplished in a column of small dimension objectifying to verify the influence of the surfactant SCO in the efficiency of oil removal. This way, they were prepared emulsions with amount it fastens of oil (50, 100, 200 and 400 ppm), being determined the great concentrations of surfactant for each one of them. Some rehearsals were still accomplished with produced water of the industry of the petroleum to compare the result with the one of the emulsions. According to the experiments, it was verified that an increase of the surfactant concentration does not implicate in a greater oil removal. The separation process use gaseous bubbles formed when a gas stream pass a liquid column, when low surfactant concentrations are used, it occurs the coalescence of the dispersed oil droplets and their transport to the top of the column, forming a new continuous phase. Such surfactants lead to a gas-liquid interface saturation, depending on the used surfactant concentration, affecting the flotation process and influencing in the removal capacity of the oily dispersed phase. A porous plate filter, with pore size varying from 40 to 250 mm, was placed at the base of the column to allow a hydrodynamic stable operation. During the experimental procedures, the operating volume of phase liquid was held constant and the rate of air flow varied in each experiment. The resulting experimental of the study hydrodynamic demonstrated what the capturing of the oil was influenced by diameter of the bubbles and air flow. With the increase flow of 300 about to 900 cm3.min-1, occurred an increase in the removal of oil phase of 44% about to 66% and the removal kinetic of oil was defined as a reaction of 1° order

Evolução da biodegradabilidade da matéria orgânica em um sistema de lagoas de estabilização

Waste stabilization ponds (WSP) have been widely used for sewage treatment in hot climate regions because they are economic and environmentally sustainable. In the present study a WSP complex comprising a primary facultative pond (PFP) followed by two maturation ponds (MP-1 and MP-2) was studied, in the city of Natal-RN. The main objective was to study the bio-degradability of organic matter through the determination of the kinetic constant k throughout the system. The work was carried out in two phases. In the first, the variability in BOD, COD and TOC concentrations and an analysis of the relations between these parameters, in the influent raw sewage, pond effluents and in specific areas inside the ponds was studied. In the second stage, the decay rate for organic matter (k) was determined throughout the system based on BOD tests on the influent sewage, pond effluents and water column samples taken from fixed locations within the ponds, using the mathematical methods of Least Squares and the Thomas equation. Subsequently k was estimated as a function of a hydrodynamic model determined from the dispersion number (d), using empirical methods and a Partial Hydrodynamic Evaluation (PHE), obtained from tracer studies in a section of the primary facultative pond corresponding to 10% of its total length. The concentrations of biodegradable organic matter, measured as BOD and COD, gradually reduced through the series of ponds, giving overall removal efficiencies of 71.95% for BOD and of 52.45% for COD. Determining the values for k, in the influent and effluent samples of the ponds using the mathematical method of Least Squares, gave the following values respectively: primary facultative pond (0,23 day-1 and 0,09 day-1), maturation 1 (0,04 day-1 and 0,03 day-1) and maturation 2 (0,03 day-1 and 0,08 day-1). When using the Thomas method, the values of k in the influents and effluents of the ponds were: primary facultative pond (0,17 day-1 and 0,07 day-1), maturation 1 (0,02 day-1 and 0,01 day-1) and maturation 2 (0,01 day-1 and 0,02 day-1). From the Partial Hydrodynamic Evaluation, in the first section of the facultative pond corresponding to 10% of its total length, it can be concluded from the dispersion number obtained of d = 0.04, that the hydraulic regime is one of dispersed flow with a kinetic constant value of 0.20 day-1

Utilização de uma coluna de flotação para remoção de compostos orgânicos da água de produção

Innovative technologies using surfactant materials have applicability in several industrial fields, including petroleum and gas areas. This study seeks to investigate the use of a surfactant derived from coconut oil (SCO saponified coconut oil) in the recovery process of organic compounds that are present in oily effluents from petroleum industry. For this end, experiments were accomplished in a column of small dimension objectifying to verify the influence of the surfactant SCO in the efficiency of oil removal. This way, they were prepared emulsions with amount it fastens of oil (50, 100, 200 and 400 ppm), being determined the great concentrations of surfactant for each one of them. Some rehearsals were still accomplished with produced water of the industry of the petroleum to compare the result with the one of the emulsions. According to the experiments, it was verified that an increase of the surfactant concentration does not implicate in a greater oil removal. The separation process use gaseous bubbles formed when a gas stream pass a liquid column, when low surfactant concentrations are used, it occurs the coalescence of the dispersed oil droplets and their transport to the top of the column, forming a new continuous phase. Such surfactants lead to a gas-liquid interface saturation, depending on the used surfactant concentration, affecting the flotation process and influencing in the removal capacity of the oily dispersed phase. A porous plate filter, with pore size varying from 40 to 250 mm, was placed at the base of the column to allow a hydrodynamic stable operation. During the experimental procedures, the operating volume of phase liquid was held constant and the rate of air flow varied in each experiment. The resulting experimental of the study hydrodynamic demonstrated what the capturing of the oil was influenced by diameter of the bubbles and air flow. With the increase flow of 300 about to 900 cm3.min-1, occurred an increase in the removal of oil phase of 44% about to 66% and the removal kinetic of oil was defined as a reaction of 1° order.

Análise hidrodinâmica e morfodinâmica do complexo estuarino do Rio Piranhas-Açu/RN, Nordeste do Brasil

The aim of this study was the seasonal characterization of the morphology, sedimentology and hydrodynamic of the Açu, Cavalos and Conchas estuaries. These estuaries are inserted in a semi-arid climate area and form the mouth of the hydrographic basin of the Piranhas-Açu river, that represent the discharge of the largest watershed in the state. They are embedded in an environment consisting of a fluvial-marine floodplain, mangrove ecosystem, sandbanks, fields of dunes, spits and sandy beaches. Adjacent to the natural units are the main local socioeconomic activities (oil industry, salt industry, shrimp farming, fishing and tourism) are dependent on this river and its conservation. The environmental monitoring is necessary because it is an area under constant action of coastal processes and at high risk of oil spill. The acquisition and interpretation of hydrodynamic, sonographic and sediment data was conducted in two campaigns, dry season (2010) and rainy season (2011), using respectively the current profiler ADCP Doppler effect, the side-scan sonar and Van Veen sampler. In these estuaries: Açu, Cavalos and Conchas were identified the following types of bedforms: flatbed and Dunes 2-D and 3-D (small to medium size), generated at lower flow regime (Froude number <1). Structures such as ripples were observed in the Açu estuary mouth. The higher values of flow discharge and velocity were recorded in the Açu estuary (434,992 m³.s-¹ and 0,554 m.s-¹). In rainy season, despite the record of highest values of discharge and flow velocities at the mouth, the energy rates upstream did not differ much from the data of the dry season. However, in all estuaries were recorded an increase in speed and flow, with reservation to the flow in the Açu estuary and flow at the mouth of the Conchas estuary. Sediment grain sizes tend to increase towards the mouth of the estuary and these ranged from very fine sand to very coarse sand, medium sand fraction being the most recurrent. Based on the data acquired and analyzed, the estuaries Açu, Cavalos and Conchas are classified as mixed , dominated by waves and tides. According to their morphology, they are classified as estuaries constructed by bar and according to the classification by salinity, estuaries Conchas and Cavalos were ranked as hypersaline estuaries, and Açu as hypersaline and vertically well mixed type C