Análise de viabilidade de injeção de fluidos alternativos ao vapor em reservatórios de óleo pesado

Many of hydrocarbon reserves existing in the world are formed by heavy oils (°API between 10 and 20). Moreover, several heavy oil fields are mature and, thus, offer great challenges for oil industry. Among the thermal methods used to recover these resources, steamflooding has been the main economically viable alternative. Latent heat carried by steam heats the reservoir, reducing oil viscosity and facilitating the production. This method has many variations and has been studied both theoretically and experimentally (in pilot projects and in full field applications). In order to increase oil recovery and reduce steam injection costs, the injection of alternative fluid has been used on three main ways: alternately, co-injected with steam and after steam injection interruption. The main objective of these injection systems is to reduce the amount of heat supplied to the reservoir, using cheaper fluids and maintaining the same oil production levels. This works discusses the use of carbon dioxide, nitrogen, methane and water as an alternative fluid to the steam. The analyzed parameters were oil recoveries and net cumulative oil productions. The reservoir simulation model corresponds to an oil reservoir of 100 m x 100 m x 28 m size, on a Cartesian coordinates system (x, y and z directions). It is a semi synthetic model with some reservoir data similar to those found in Brazilian Potiguar Basin. All studied cases were done using the simulator STARS from CMG (Computer Modelling Group, version 2009.10). It was found that waterflood after steam injection interruption achieved the highest net cumulative oil compared to other fluids injection. Moreover, it was observed that steam and alternative fluids, co-injected and alternately, did not present increase on profitability project compared with steamflooding

Análise de viabilidade de injeção de fluidos alternativos ao vapor em reservatórios de óleo pesado

Many of hydrocarbon reserves existing in the world are formed by heavy oils (°API between 10 and 20). Moreover, several heavy oil fields are mature and, thus, offer great challenges for oil industry. Among the thermal methods used to recover these resources, steamflooding has been the main economically viable alternative. Latent heat carried by steam heats the reservoir, reducing oil viscosity and facilitating the production. This method has many variations and has been studied both theoretically and experimentally (in pilot projects and in full field applications). In order to increase oil recovery and reduce steam injection costs, the injection of alternative fluid has been used on three main ways: alternately, co-injected with steam and after steam injection interruption. The main objective of these injection systems is to reduce the amount of heat supplied to the reservoir, using cheaper fluids and maintaining the same oil production levels. This works discusses the use of carbon dioxide, nitrogen, methane and water as an alternative fluid to the steam. The analyzed parameters were oil recoveries and net cumulative oil productions. The reservoir simulation model corresponds to an oil reservoir of 100 m x 100 m x 28 m size, on a Cartesian coordinates system (x, y and z directions). It is a semi synthetic model with some reservoir data similar to those found in Brazilian Potiguar Basin. All studied cases were done using the simulator STARS from CMG (Computer Modelling Group, version 2009.10). It was found that waterflood after steam injection interruption achieved the highest net cumulative oil compared to other fluids injection. Moreover, it was observed that steam and alternative fluids, co-injected and alternately, did not present increase on profitability project compared with steamflooding

Estudo paramétrico da segregação gravitacional na injeção contínua de vapor

Steam injection is the most used method of additional recovery for the extraction of heavy oil. In this type procedure is common to happen gravitational segregation and this phenomenon can affect the production of oil and therefore, it shoulds be considered in the projects of continuous steam injection. For many years, the gravitational segregation was not adequately considered in the calculation procedures in Reservoir Engineering. The effect of the gravity causes the segregation of fluids inside the porous media according to their densities. The results of simulation arising from reservoirs could provide the ability to deal with the gravity, and it became apparent that the effects of the gravity could significantly affect the performance of the reservoir. It know that the gravitational segregation can happen in almost every case where there is injection of light fluid, specially the steam, and occurs with greater intensity for viscous oil reservoirs. This work discusses the influence of some parameters of the rock-reservoir in segregation as viscosity, permeability, thickness, cover gas, porosity. From a model that shows the phenomenon with greater intensity, optimized some operational parameters as the rate flow rate steam, distance between the wells injector-producer, and interval of completion which contributed to the reduction in gravity override, thus increasing the oil recovery. It was shown a greater technical-economic viability for the model of distance between the wells 100 m. The analysis was performed using the simulator of CMG (Computer Modeling Group-Stars 2007.11, in which was observed by iterating between studied variables in heavy oil reservoirs with similar characteristics to Brazilian Northeast

Estudo paramétrico da segregação gravitacional na injeção contínua de vapor

Steam injection is the most used method of additional recovery for the extraction of heavy oil. In this type procedure is common to happen gravitational segregation and this phenomenon can affect the production of oil and therefore, it shoulds be considered in the projects of continuous steam injection. For many years, the gravitational segregation was not adequately considered in the calculation procedures in Reservoir Engineering. The effect of the gravity causes the segregation of fluids inside the porous media according to their densities. The results of simulation arising from reservoirs could provide the ability to deal with the gravity, and it became apparent that the effects of the gravity could significantly affect the performance of the reservoir. It know that the gravitational segregation can happen in almost every case where there is injection of light fluid, specially the steam, and occurs with greater intensity for viscous oil reservoirs. This work discusses the influence of some parameters of the rock-reservoir in segregation as viscosity, permeability, thickness, cover gas, porosity. From a model that shows the phenomenon with greater intensity, optimized some operational parameters as the rate flow rate steam, distance between the wells injector-producer, and interval of completion which contributed to the reduction in gravity override, thus increasing the oil recovery. It was shown a greater technical-economic viability for the model of distance between the wells 100 m. The analysis was performed using the simulator of CMG (Computer Modeling Group-Stars 2007.11, in which was observed by iterating between studied variables in heavy oil reservoirs with similar characteristics to Brazilian Northeast

Hybrid modeling of convective laminar flow in a permeable tube associated with the cross-flow process

The confined flows in tubes with permeable surfaces arc associated to tangential filtration processes (microfiltration or ultrafiltration). The complexity of the phenomena do not allow for the development of exact analytical solutions, however, approximate solutions are of great interest for the calculation of the transmembrane outflow and estimate of the concentration, polarization phenomenon. In the present work, the generalized integral transform technique (GITT) was employed in solving the laminar and permanent flow in permeable tubes of Newtonian and incompressible fluid. The mathematical formulation employed the parabolic differential equation of chemical species conservation (convective-diffusive equation). The velocity profiles for the entrance region flow, which are found in the connective terms of the equation, were assessed by solutions obtained from literature. The velocity at the permeable wall was considered uniform, with the concentration at the tube wall regarded as variable with an axial position. A computational methodology using global error control was applied to determine the concentration in the wall and concentration boundary layer thickness. The results obtained for the local transmembrane flux and the concentration boundary layer thickness were compared against others in literature. (C) 2007 Elsevier B.V. All rights reserved.

Desenvolvimento de pastas de cimento flexíveis para poços de petróleo submetidos à injeção de vapor utilizando látex do tipo SBR

Novel cementing materials formulations containing flexible polymeric admixtures have been studied aiming at improving the mechanical behavior of oil well cement slurries submitted to steam injection. However, research activities in this sector are still under development. The steam injected directly into the well causes casing dilation, which after a reduction in temperature, tends to return to its original dimensions, resulting in crack formation and hydraulic isolation loss of the well, which will result in shortening of well life. In this scenario, the objective of the present study was to evaluate the mechanical behavior of Portland-based slurries containing SBR latex, applied in oil well cementing of wells submitted to steam injection. Were formulated slurries with densities of 1.797 g/cm3 (15.0 lb/Gal) and 1.869 g/cm3 (15.6 lb/Gal), containing admixtures with a latex concentration of 0; 66.88; 133.76; 200.64 and 267.52 L/m3 (0, 0.5, 1.0, 1.5 and 2.0 gpc). Tests including rheology, fluid loss control, thickening time, API compressive strength and splitting tensile strength, beyond steam injection simulation. Microstrutural characteristics of the slurries were also performed (XRD, TG, FTIR and SEM). The results showed that increasing the polymer concentration increased in the rheological properties and fluid loss, and a decrease in the elasticity modulus of the cement slurries. The results obtained showed that the slurries can be applied in cementing operations of oil wells submitted to steam injection.

Desenvolvimento de pastas de cimento flexíveis para poços de petróleo submetidos à injeção de vapor utilizando látex do tipo SBR

Novel cementing materials formulations containing flexible polymeric admixtures have been studied aiming at improving the mechanical behavior of oil well cement slurries submitted to steam injection. However, research activities in this sector are still under development. The steam injected directly into the well causes casing dilation, which after a reduction in temperature, tends to return to its original dimensions, resulting in crack formation and hydraulic isolation loss of the well, which will result in shortening of well life. In this scenario, the objective of the present study was to evaluate the mechanical behavior of Portland-based slurries containing SBR latex, applied in oil well cementing of wells submitted to steam injection. Were formulated slurries with densities of 1.797 g/cm3 (15.0 lb/Gal) and 1.869 g/cm3 (15.6 lb/Gal), containing admixtures with a latex concentration of 0; 66.88; 133.76; 200.64 and 267.52 L/m3 (0, 0.5, 1.0, 1.5 and 2.0 gpc). Tests including rheology, fluid loss control, thickening time, API compressive strength and splitting tensile strength, beyond steam injection simulation. Microstrutural characteristics of the slurries were also performed (XRD, TG, FTIR and SEM). The results showed that increasing the polymer concentration increased in the rheological properties and fluid loss, and a decrease in the elasticity modulus of the cement slurries. The results obtained showed that the slurries can be applied in cementing operations of oil wells submitted to steam injection.

Parallel modeling of the electric field distribution in the brain

The term "Brain Imaging" identi�es a set of techniques to analyze the structure and/or functional behavior of the brain in normal and/or pathological situations. These techniques are largely used in the study of brain activity. In addition to clinical usage, analysis of brain activity is gaining popularity in others recent �fields, i.e. Brain Computer Interfaces (BCI) and the study of cognitive processes. In this context, usage of classical solutions (e.g. f MRI, PET-CT) could be unfeasible, due to their low temporal resolution, high cost and limited portability. For these reasons alternative low cost techniques are object of research, typically based on simple recording hardware and on intensive data elaboration process. Typical examples are ElectroEncephaloGraphy (EEG) and Electrical Impedance Tomography (EIT), where electric potential at the patient's scalp is recorded by high impedance electrodes. In EEG potentials are directly generated from neuronal activity, while in EIT by the injection of small currents at the scalp. To retrieve meaningful insights on brain activity from measurements, EIT and EEG relies on detailed knowledge of the underlying electrical properties of the body. This is obtained from numerical models of the electric �field distribution therein. The inhomogeneous and anisotropic electric properties of human tissues make accurate modeling and simulation very challenging, leading to a tradeo�ff between physical accuracy and technical feasibility, which currently severely limits the capabilities of these techniques. Moreover elaboration of data recorded requires usage of regularization techniques computationally intensive, which influences the application with heavy temporal constraints (such as BCI). This work focuses on the parallel implementation of a work-flow for EEG and EIT data processing. The resulting software is accelerated using multi-core GPUs, in order to provide solution in reasonable times and address requirements of real-time BCI systems, without over-simplifying the complexity and accuracy of the head models.

Simulation by discrete mass modeling of offshore wind turbine system with DC link

This paper presents an integrated model for an offshore wind turbine taking into consideration a contribution for the marine wave and wind speed with perturbations influences on the power quality of current injected into the electric grid. The paper deals with the simulation of one floating offshore wind turbine equipped with a permanent magnet synchronous generator, and a two-level converter connected to an onshore electric grid. The use of discrete mass modeling is accessed in order to reveal by computing the total harmonic distortion on how the perturbations of the captured energy are attenuated at the electric grid injection point. Two torque actions are considered for the three-mass modeling, the aerodynamic on the flexible part and on the rigid part of the blades. Also, a torque due to the influence of marine waves in deep water is considered. Proportional integral fractional-order control supports the control strategy. A comparison between the drive train models is presented.

The modeling of turbulent reactive flows based on multiple mapping conditioning

A new modeling approach-multiple mapping conditioning (MMC)-is introduced to treat mixing and reaction in turbulent flows. The model combines the advantages of the probability density function and the conditional moment closure methods and is based on a certain generalization of the mapping closure concept. An equivalent stochastic formulation of the MMC model is given. The validity of the closuring hypothesis of the model is demonstrated by a comparison with direct numerical simulation results for the three-stream mixing problem. (C) 2003 American Institute of Physics.

Cell-based computational modeling of vascular morphogenesis using Tissue Simulation Toolkit.

Computational modeling has become a widely used tool for unraveling the mechanisms of higher level cooperative cell behavior during vascular morphogenesis. However, experimenting with published simulation models or adding new assumptions to those models can be daunting for novice and even for experienced computational scientists. Here, we present a step-by-step, practical tutorial for building cell-based simulations of vascular morphogenesis using the Tissue Simulation Toolkit (TST). The TST is a freely available, open-source C++ library for developing simulations with the two-dimensional cellular Potts model, a stochastic, agent-based framework to simulate collective cell behavior. We will show the basic use of the TST to simulate and experiment with published simulations of vascular network formation. Then, we will present step-by-step instructions and explanations for building a recent simulation model of tumor angiogenesis. Demonstrated mechanisms include cell-cell adhesion, chemotaxis, cell elongation, haptotaxis, and haptokinesis.

Analysis and spatial modeling of the indoor radon Swiss data

The present research deals with an important public health threat, which is the pollution created by radon gas accumulation inside dwellings. The spatial modeling of indoor radon in Switzerland is particularly complex and challenging because of many influencing factors that should be taken into account. Indoor radon data analysis must be addressed from both a statistical and a spatial point of view. As a multivariate process, it was important at first to define the influence of each factor. In particular, it was important to define the influence of geology as being closely associated to indoor radon. This association was indeed observed for the Swiss data but not probed to be the sole determinant for the spatial modeling. The statistical analysis of data, both at univariate and multivariate level, was followed by an exploratory spatial analysis. Many tools proposed in the literature were tested and adapted, including fractality, declustering and moving windows methods. The use of Quan-tité Morisita Index (QMI) as a procedure to evaluate data clustering in function of the radon level was proposed. The existing methods of declustering were revised and applied in an attempt to approach the global histogram parameters. The exploratory phase comes along with the definition of multiple scales of interest for indoor radon mapping in Switzerland. The analysis was done with a top-to-down resolution approach, from regional to local lev¬els in order to find the appropriate scales for modeling. In this sense, data partition was optimized in order to cope with stationary conditions of geostatistical models. Common methods of spatial modeling such as Κ Nearest Neighbors (KNN), variography and General Regression Neural Networks (GRNN) were proposed as exploratory tools. In the following section, different spatial interpolation methods were applied for a par-ticular dataset. A bottom to top method complexity approach was adopted and the results were analyzed together in order to find common definitions of continuity and neighborhood parameters. Additionally, a data filter based on cross-validation was tested with the purpose of reducing noise at local scale (the CVMF). At the end of the chapter, a series of test for data consistency and methods robustness were performed. This lead to conclude about the importance of data splitting and the limitation of generalization methods for reproducing statistical distributions. The last section was dedicated to modeling methods with probabilistic interpretations. Data transformation and simulations thus allowed the use of multigaussian models and helped take the indoor radon pollution data uncertainty into consideration. The catego-rization transform was presented as a solution for extreme values modeling through clas-sification. Simulation scenarios were proposed, including an alternative proposal for the reproduction of the global histogram based on the sampling domain. The sequential Gaussian simulation (SGS) was presented as the method giving the most complete information, while classification performed in a more robust way. An error measure was defined in relation to the decision function for data classification hardening. Within the classification methods, probabilistic neural networks (PNN) show to be better adapted for modeling of high threshold categorization and for automation. Support vector machines (SVM) on the contrary performed well under balanced category conditions. In general, it was concluded that a particular prediction or estimation method is not better under all conditions of scale and neighborhood definitions. Simulations should be the basis, while other methods can provide complementary information to accomplish an efficient indoor radon decision making.

Simulation and Analysis of Arterial Traffic Operations Along the U.S. 61 Corridor in Burlington, Iowa, 1998

The Center for Transportation Research and Education (CTRE) used the traffic simulation model CORSIM to access proposed capacity and safety improvement strategies for the U.S. 61 corridor through Burlington, Iowa. The comparison between the base and alternative models allow for evaluation of the traffic flow performance under the existing conditions as well as other design scenarios. The models also provide visualization of performance for interpretation by technical staff, public policy makers, and the public. The objectives of this project are to evaluate the use of traffic simulation models for future use by the Iowa Department of Transportation (DOT) and to develop procedures for employing simulation modeling to conduct the analysis of alternative designs. This report presents both the findings of the U.S. 61 evaluation and an overview of model development procedures. The first part of the report includes the simulation modeling development procedures. The simulation analysis is illustrated through the Burlington U.S. 61 corridor case study application. Part I is not intended to be a user manual but simply introductory guidelines for traffic simulation modeling. Part II of the report evaluates the proposed improvement concepts in a side by side comparison of the base and alternative models.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Kinetic modeling of the biosorption of Cd2+ ions from aqueous solutions onto Eichhornia crassipes roots using potentiometry: low-cost alternative to conventional methods

This work presents the use of potentiometric measurements for kinetic studies of biosorption of Cd2+ ions from aqueous solutions on Eichhornia crassipes roots. The open circuit potential of the Cd/Cd2+ electrode of the first kind was measured during the bioadsorption process. The amount of Cd2+ ions accumulated was determined in real time. The data were fit to different models, with the pseudo-second-order model proving to be the best in describing the data. The advantages and limitations of the methodology proposed relative to the traditional method are discussed.