Aplicação da técnica de elementos finitos na estratificação de armazenadores térmicos com barreiras de meios porosos

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Simulação numérica de fluxo de águas subterrâneas na fazenda Cara Preta - Santa Rita do Passa Quatro - SP

The term model refers to any representation of a real system. The use of models in Hydrogeology can be valuable predictive tools for management of groundwater resources. The numeric models of groundwater flow, object of this study, consist on a set of differential equations that describe the water flow in the porous medium. In this context, numeric simulations were made for a sub-basin located at Cara Preta farm – Santa Rita do Passa Quatro – SP. The aquifer at the local is composed by rocks of Pirambóia Formation, which is part of Guarani Aquifer System. It was developed a conceptual model from previous studies in the area, and from that, simulations were made through the software Visual Modflow®. The conceptual model established previously was considered consistent through the results of simulation.

Comportamento das curvas de saturação da água e LNAPL em solo da formação Rio Claro

Pós-graduação em Geociências e Meio Ambiente - IGCE

Método multigrid algébrico: reutilização das estruturas multigrid no transporte de contaminantes

A necessidade de obter solução de grandes sistemas lineares resultantes de processos de discretização de equações diferenciais parciais provenientes da modelagem de diferentes fenômenos físicos conduz à busca de técnicas numéricas escaláveis. Métodos multigrid são classificados como algoritmos escaláveis.Um estimador de erros deve estar associado à solução numérica do problema discreto de modo a propiciar a adequada avaliação da solução obtida pelo processo de aproximação. Nesse contexto, a presente tese caracteriza-se pela proposta de reutilização das estruturas matriciais hierárquicas de operadores de transferência e restrição dos métodos multigrid algébricos para acelerar o tempo de solução dos sistemas lineares associados à equação do transporte de contaminantes em meio poroso saturado. Adicionalmente, caracteriza-se pela implementação das estimativas residuais para os problemas que envolvem dados constantes ou não constantes, os regimes de pequena ou grande advecção e pela proposta de utilização das estimativas residuais associadas ao termo de fonte e à condição inicial para construir procedimentos adaptativos para os dados do problema. O desenvolvimento dos códigos do método de elementos finitos, do estimador residual e dos procedimentos adaptativos foram baseados no projeto FEniCS, utilizando a linguagem de programação PYTHONR e desenvolvidos na plataforma Eclipse. A implementação dos métodos multigrid algébricos com reutilização considera a biblioteca PyAMG. Baseado na reutilização das estruturas hierárquicas, os métodos multigrid com reutilização com parâmetro fixo e automática são propostos, e esses conceitos são estendidos para os métodos iterativos não-estacionários tais como GMRES e BICGSTAB. Os resultados numéricos mostraram que o estimador residual captura o comportamento do erro real da solução numérica, e fornece algoritmos adaptativos para os dados cuja malha retornada produz uma solução numérica similar à uma malha uniforme com mais elementos. Adicionalmente, os métodos com reutilização são mais rápidos que os métodos que não empregam o processo de reutilização de estruturas. Além disso, a eficiência dos métodos com reutilização também pode ser observada na solução do problema auxiliar, o qual é necessário para obtenção das estimativas residuais para o regime de grande advecção. Esses resultados englobam tanto os métodos multigrid algébricos do tipo SA quanto os métodos pré-condicionados por métodos multigrid algébrico SA, e envolvem o transporte de contaminantes em regime de pequena e grande advecção, malhas estruturadas e não estruturadas, problemas bidimensionais, problemas tridimensionais e domínios com diferentes escalas.

Avaliação da influência de cinza da casca de arroz em cimentos para poços de petróleo submetidos à injeção de vapor

Reservoirs that present highly viscous oils require methods to aid in their recovery to the surface. The elev ated oil viscosity hinders its flow through porous media and conventional recovery methods have not obtained significant efficiency. As such, the injection of steam into the reservoir through an injection well has been the most widely used method of therma l recovery, for it allows elevated volumes of recovery due to the viscosity reduction of the oil, facilitating the oil’s mobility within the rock formation and consequently into the production well where it will be exploited. On the other hand, the injecti on of vapor not only affects the fluids found in the rock pores, but the entire structure that composes the well where it is injected due to the high temperatures used in the process. This temperature increment is conducted to the cement, found in the annu lus, responsible for the isolation of the well and the well casing. Temperatures above 110 ̊C create new fazes rich in calcium in the cement matrix, resulting in the reduction of its permeability and the consequential phenomenon of mechanical resistance ret rogression. These alterations generate faults in the cement, reducing the well’s hydraulic isolation, creating insecurity in the operations in which the well will be submitted as well as the reduction of its economic life span. As a way of reducing this re trograde effect, this study has the objective of evaluating the incorporation of rice husk ash as a mineral additive substitute of silica flour , commercially utilized as a source of silica to reduce the CaO/SiO 2 ratio in the cement pastes submitted to high temperatures in thermal recovery. Cement pastes were formulated containing 20 and 30% levels of ash, apart from the basic paste (water + cement) and a reference paste (water + cement + 40% silica flour) for comparison purposes. The tests were executed th rough compression resistance tests, X - Ray diffraction (XRD) techniques, thermogravimetry (TG), scanning electron microscopy (SEM) and chemical anal ysis BY X - ray fluorescence (EDS) on the pastes submitted to cure at low temperatures (45 ̊C) for 28 days following a cure at 280 ̊C and a pressure of 2,000 PSI for 3 days, simulating vapor injection. The results obtained show that the paste containing 30% r ice shell ash is satisfactory, obtaining mechanical resistance desired and equivalent to that of the paste containing 40% silica flour, since the products obtained were hydrated with low CaO/SiO 2 ratio, like the Tobermorita and Xonotlita fases, proving its applicability in well subject to vapor injection.

Estudo paramétrico da injeção de polímeros em reservatórios de petróleo

Developing an efficient methodology for oil recovery is extremely important . Within the range of enh anced oil recovery, known as EOR, the injection of polymer solutions becomes effective in controlling the mobility of displacing fluid . This method consists of adding polymers to the injection water to increase its viscosity, so that more water diffuses in to the porous medium and increasing the sweep efficiency in the reservoir. This work is studied by numerical simulation , application of the injection polymer solution in a homogeneous reservoir , semisynthetic with similar characteristics to the reservoirs of the Brazilian Northeast , numerical simulations were performed using thermal simulator STARS from CMG (Computer Modelling Group ). The study aimed to analyze the influence of some parameters on the behavior of reservoir oil production, with the response to cumulative production. Simulations were performed to analyze the influence of water injection, polymer solution and alternating injection of water banks and polymer solution, comparing the results for each simulated condition. The primary outcomes were: oil viscosity, percentage of injected polymer, polymer viscosity and flow rate of water injection. The evaluation of the influence of variables consisted of a complete experimental design followed a Pareto analysis for the purpose of pointing out which va riables would be most influential on the response represented b y the cumulative oil production . It was found that all variables significantly influenced the recovery of oil and the injection of polymer solution on an ongoing basis is more efficient for the cumulative production compared to oil recovery by continuous water injection. The primary recovery show ed low levels of oil production , water injection significantly improves the pro duction of oil in the reservoir , but the injection of polymer solution em erges as a new methodology to increase the production of oil, increasing the life of the well and possible reduction of water produced.

Estudo das forças capilar, viscosa e gravitacional através de grupos adimensionais na simulação da drenagem gravitacional assistida por gás (GAGD)

The distribution and mobilization of fluid in a porous medium depend on the capillary, gravity, and viscous forces. In oil field, the processes of enhanced oil recovery involve change and importance of these forces to increase the oil recovery factor. In the case of gas assisted gravity drainage (GAGD) process is important to understand the physical mechanisms to mobilize oil through the interaction of these forces. For this reason, several authors have developed physical models in laboratory and core floods of GAGD to study the performance of these forces through dimensionless groups. These models showed conclusive results. However, numerical simulation models have not been used for this type of study. Therefore, the objective of this work is to study the performance of capillary, viscous and gravity forces on GAGD process and its influence on the oil recovery factor through a 2D numerical simulation model. To analyze the interplay of these forces, dimensionless groups reported in the literature have been used such as Capillary Number (Nc), Bond number (Nb) and Gravity Number (Ng). This was done to determine the effectiveness of each force related to the other one. A comparison of the results obtained from the numerical simulation was also carried out with the results reported in the literature. The results showed that before breakthrough time, the lower is the injection flow rate, oil recovery is increased by capillary force, and after breakthrough time, the higher is the injection flow rate, oil recovery is increased by gravity force. A good relationship was found between the results obtained in this research with those published in the literature. The simulation results indicated that before the gas breakthrough, higher oil recoveries were obtained at lower Nc and Nb and, after the gas breakthrough, higher oil recoveries were obtained at lower Ng. The numerical models are consistent with the reported results in the literature

Modelagem multiescala de fenômenos eletrocinéticos em meios porosos carregados eletricamente: aplicação a meios porosos argilosos

In this master thesis, we propose a multiscale mathematical and computational model for electrokinetic phenomena in porous media electrically charged. We consider a porous medium rigid and incompressible saturated by an electrolyte solution containing four monovalent ionic solutes completely diluted in the aqueous solvent. Initially we developed the modeling electrical double layer how objective to compute the electrical potential, surface density of electrical charges and considering two chemical reactions, we propose a 2-pK model for calculating the chemical adsorption occurring in the domain of electrical double layer. Having the nanoscopic model, we deduce a model in the microscale, where the electrochemical adsorption of ions, protonation/ deprotonation reactions and zeta potential obtained in the nanoscale, are incorporated through the conditions of interface uid/solid of the Stokes problem and transportation of ions, modeled by equations of Nernst-Planck. Using the homogenization technique of periodic structures, we develop a model in macroscopic scale with respective cells problems for the e ective macroscopic parameters of equations. Finally, we propose several numerical simulations of the multiscale model for uid ow and transport of reactive ionic solute in a saturated aqueous solution of kaolinite. Using nanoscopic model we propose some numerical simulations of electrochemical adsorption phenomena in the electrical double layer. Making use of the nite element method discretize the macroscopic model and propose some numerical simulations in basic and acid system aiming to quantify the transport of ionic solutes in porous media electrically charged.

Modelagem multiescala de fenômenos eletrocinéticos em meios porosos carregados eletricamente: aplicação a meios porosos argilosos

In this master thesis, we propose a multiscale mathematical and computational model for electrokinetic phenomena in porous media electrically charged. We consider a porous medium rigid and incompressible saturated by an electrolyte solution containing four monovalent ionic solutes completely diluted in the aqueous solvent. Initially we developed the modeling electrical double layer how objective to compute the electrical potential, surface density of electrical charges and considering two chemical reactions, we propose a 2-pK model for calculating the chemical adsorption occurring in the domain of electrical double layer. Having the nanoscopic model, we deduce a model in the microscale, where the electrochemical adsorption of ions, protonation/ deprotonation reactions and zeta potential obtained in the nanoscale, are incorporated through the conditions of interface uid/solid of the Stokes problem and transportation of ions, modeled by equations of Nernst-Planck. Using the homogenization technique of periodic structures, we develop a model in macroscopic scale with respective cells problems for the e ective macroscopic parameters of equations. Finally, we propose several numerical simulations of the multiscale model for uid ow and transport of reactive ionic solute in a saturated aqueous solution of kaolinite. Using nanoscopic model we propose some numerical simulations of electrochemical adsorption phenomena in the electrical double layer. Making use of the nite element method discretize the macroscopic model and propose some numerical simulations in basic and acid system aiming to quantify the transport of ionic solutes in porous media electrically charged.

Adição de superplastificante e anti-segregante à base de CMC em pastas de cimento para Poços de Petróleo

Oil well cementing materials consist of slurries of Special class Portland cement dispersed in water. Admixtures can be used to provide the necessary fluidity, so the material can be efficiently pumped down as well as penetrate porous rocks with controlled filter loss. Construction admixtures can be used to modify the properties of oil well cements provided they can withstand and hold their properties at the higher than ambient temperatures usually encountered in oil fields. In civil construction, superplasticizer play the role of dispersants that reduce the facto r of water cement improve mechanical properties and fluidity of the cement, whereas anti-segregation agents improve the workability of the slurry. In the present study, oil well cement slurries were produced adding both a dispersant and an anti-segregation agent conventionally used in Portland CPII-Z-32 RS cement aiming at materials for primary cementing and squeeze operations. Three basic aspects were evaluated: fluidity, filter loss and the synergetic effect of the admixtures at two temperatures, i.e., 27°C and 56°C, following API RP 10B practical recommendations. The slurries were prepared using admixture concentrations varying from 2.60 Kgf/m3 (0.02 gallft3) to 5.82 Kgf/m3 (0.045 galJft3) BWOC. The density of the slurries was set to 1.89 g/cm3 (15.8 Ib/gal). 0.30 to 0.60% BWOC of a CMC-based anti-segregation agent was added to the cement to control the filter loss. The results showed that the addition of anti-segregation at concentrations above 0.55% by weight of cement resulted in the increased viscosity of the folders in temperatures evaluated. The increasing the temperature of the tests led to a reduction in the performance of anti-segregation. At concentrations of 5.20 kgf/m3 (0,040 gallft3) and 5.82 Kgf/m3 (0,045 gal/ft 3) observed a better performance of the properties evaluated in the proposed system. At low temperature was observed instability in the readings of rheology for all concentrations of anti-segregation. Contents that increasing the concentration of anti¬-segregation is limited concentrations greater than 0.55 % BWOC of the CMC in temperature analyzed. The use of the system with CMC promoted a good performance against the properties evaluated. The principal function of anti¬-segregation was optimized with increasing concentration of superplasticizer, at temperatures above the 2rC. The study of the behaviour of systemic additives, resulting in slurries of cement, which can be optimized face studies of other intrinsic properties in oil fields

Efeito da incorporação de diatomita na estabilidade e permeabilidade de compósitos cimentíceos espumados aplicados a poços de petróleo

Lightweight oilwell cement slurries have been recently studied as a mean to improve zonal isolation and sheath-porous formation adherence. Foamed slurries consisting of Portland cement and air-entraining admixtures have become an interesting option for this application. The loss in hydrostatic pressure as a consequence of cement hydration results in the expansion of the air bubbles entrapped in the cement matrix, thus improving the sheath-porous formation contact. Consequently, slurries are able to better retain their water to complete the hydration process. The main objective of the present study was to evaluate the effect of the addition of an air-entraining admixture on the density, stability and permeability of composite slurries containing Portland cement and diatomite as light mineral load. Successful formulations are potential cementing materials for low fracture gradient oilwells. The experimental procedures used for slurry preparation and characterization were based on the American Petroleum Institute and ABNT guidelines Slurries containing a pre-established concentration of the air-entraining admixture and different contents of diatomite were prepared aiming at final densities of 13 to 15 lb/gal. The results revealed that the reduction of 15 to 25% of the density of the slurries did not significantly affect their strength. The addition of both diatomite and the air-entraining admixture increased the viscosity of the slurry providing better air-bubble retention in the volume of the slurry. Stable slurries depicted bottom to top density variation of less than 1.0 lb/gal and length reduction of the stability sample of 5.86 mm. Finally, permeability coefficient values between 0.617 and 0.406 mD were obtained. Therefore, lightweight oilwell cement slurries depicting a satisfactory set of physicochemical and mechanical properties can be formulated using a combination of diatomite and air-entraining admixtures for low fracture gradient oilwells

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

Modelagem analítica e experimental da filtração em meios porosos

Deep bed filtration occurs in several industrial and environmental processes like water filtration and soil contamination. In petroleum industry, deep bed filtration occurs near to injection wells during water injection, causing injectivity reduction. It also takes place during well drilling, sand production control, produced water disposal in aquifers, etc. The particle capture in porous media can be caused by different physical mechanisms (size exclusion, electrical forces, bridging, gravity, etc). A statistical model for filtration in porous media is proposed and analytical solutions for suspended and retained particles are derived. The model, which incorporates particle retention probability, is compared with the classical deep bed filtration model allowing a physical interpretation of the filtration coefficients. Comparison of the obtained analytical solutions for the proposed model with the classical model solutions allows concluding that the larger the particle capture probability, the larger the discrepancy between the proposed and the classical models

Modelagem e previsão da perda de injetividade em poços canhoneados

Waterflooding is a technique largely applied in the oil industry. The injected water displaces oil to the producer wells and avoid reservoir pressure decline. However, suspended particles in the injected water may cause plugging of pore throats causing formation damage (permeability reduction) and injectivity decline during waterflooding. When injectivity decline occurs it is necessary to increase the injection pressure in order to maintain water flow injection. Therefore, a reliable prediction of injectivity decline is essential in waterflooding projects. In this dissertation, a simulator based on the traditional porous medium filtration model (including deep bed filtration and external filter cake formation) was developed and applied to predict injectivity decline in perforated wells (this prediction was made from history data). Experimental modeling and injectivity decline in open-hole wells is also discussed. The injectivity of modeling showed good agreement with field data, which can be used to support plan stimulation injection wells

Simulação de fluxo de fluidos em meios porosos desordenados uma análise de efeito de escala na estimativa da permeabilidade e do coeficiente de arrasto

The present study provides a methodology that gives a predictive character the computer simulations based on detailed models of the geometry of a porous medium. We using the software FLUENT to investigate the flow of a viscous Newtonian fluid through a random fractal medium which simplifies a two-dimensional disordered porous medium representing a petroleum reservoir. This fractal model is formed by obstacles of various sizes, whose size distribution function follows a power law where exponent is defined as the fractal dimension of fractionation Dff of the model characterizing the process of fragmentation these obstacles. They are randomly disposed in a rectangular channel. The modeling process incorporates modern concepts, scaling laws, to analyze the influence of heterogeneity found in the fields of the porosity and of the permeability in such a way as to characterize the medium in terms of their fractal properties. This procedure allows numerically analyze the measurements of permeability k and the drag coefficient Cd proposed relationships, like power law, for these properties on various modeling schemes. The purpose of this research is to study the variability provided by these heterogeneities where the velocity field and other details of viscous fluid dynamics are obtained by solving numerically the continuity and Navier-Stokes equations at pore level and observe how the fractal dimension of fractionation of the model can affect their hydrodynamic properties. This study were considered two classes of models, models with constant porosity, MPC, and models with varying porosity, MPV. The results have allowed us to find numerical relationship between the permeability, drag coefficient and the fractal dimension of fractionation of the medium. Based on these numerical results we have proposed scaling relations and algebraic expressions involving the relevant parameters of the phenomenon. In this study analytical equations were determined for Dff depending on the geometrical parameters of the models. We also found a relation between the permeability and the drag coefficient which is inversely proportional to one another. As for the difference in behavior it is most striking in the classes of models MPV. That is, the fact that the porosity vary in these models is an additional factor that plays a significant role in flow analysis. Finally, the results proved satisfactory and consistent, which demonstrates the effectiveness of the referred methodology for all applications analyzed in this study.