Estudo paramétrico da recuperação de óleo no processo de drenagem gravitacional com injeção de CO2

The gas injection has become the most important IOR process in the United States. Furthermore, the year 2006 marks the first time the gas injection IOR production has surpassed that of steam injection. In Brazil, the installation of a petrochemical complex in the Northeast of Brazil (Bahia State) offers opportunities for the injection of gases in the fields located in the Recôncavo Basin. Field-scale gas injection applications have almost always been associated with design and operational difficulties. The mobility ratio, which controls the volumetric sweep, between the injected gas and displaced oil bank in gas processes, is typically unfavorable due to the relatively low viscosity of the injected gas. Furthermore, the difference between their densities results in severe gravity segregation of fluids in the reservoirs, consequently leading to poor control in the volumetric sweep. Nowadays, from the above applications of gas injection, the WAG process is most popular. However, in attempting to solve the mobility problems, the WAG process gives rise to other problems associated with increased water saturation in the reservoir including diminished gas injectivity and increased competition to the flow of oil. The low field performance of WAG floods with oil recoveries in the range of 5-10% is a clear indication of these problems. In order to find na effective alternative to WAG, the Gas Assisted Gravity Drainage (GAGD) was developed. This process is designed to take advantage of gravity force to allow vertical segregation between the injected CO2 and reservoir crude oil due to their density difference. This process consists of placing horizontal producers near the bottom of the pay zone and injecting gás through existing vertical wells in field. Homogeneous models were used in this work which can be extrapolated to commercial application for fields located in the Northeast of Brazil. The simulations were performed in a CMG simulator, the STARS 2007.11, where some parameters and their interactions were analyzed. The results have shown that the CO2 injection in GAGD process increased significantly the rate and the final recovery of oil

Influência da adição de gres porcelanato no comportamento mecânico e microestrutural em pastas de cimento portland para cimentação de poços de petróleo

Cementation operation consists in an extremely important work for the phases of perforation and completion of oil wells, causing a great impact on the well productivity. Several problems can occur with the cement during the primary cementation, as well as throughout the productive period. The corrective operations are frequent, but they are expensive and demands production time. Besides the direct cost, prejudices from the interruption of oil and gas production till the implementation of a corrective operation must be also taken into account. The purpose of this work is the development of an alternative cement paste constituted of Portland cement and porcelainized stoneware residue produced by ceramic industry in order to achieve characteristics as low permeability, high tenacity, and high mechanical resistance, capable of supporting various operations as production or oil wells recuperation. Four different concentration measures of hydrated paste were evaluated: a reference paste, and three additional ones with ceramic residue in concentrations of the order of 10%, 20% and 30% in relation to cement dough. High resistance and low permeability were found in high concentration of residues, as well as it was proved the pozolanic reactivity of the residue in relation to Portland cement, which was characterized through x-ray and thermogravimetry assays. It was evident the decrease of calcium hydroxide content, once it was substituted by formation of new hydrated products as it was added ceramic residue

Nova técnica para obtenção de fraturas com altíssima condutividade em poços de petróleo

Increase hydrocarbons production is the main goal of the oilwell industry worldwide. Hydraulic fracturing is often applied to achieve this goal due to a combination of attractive aspects including easiness and low operational costs associated with fast and highly economical response. Conventional fracturing usually involves high-flowing high-pressure pumping of a viscous fluid responsible for opening the fracture in the hydrocarbon producing rock. The thickness of the fracture should be enough to assure the penetration of the particles of a solid proppant into the rock. The proppant is driven into the target formation by a carrier fluid. After pumping, all fluids are filtered through the faces of the fracture and penetrate the rock. The proppant remains in the fracture holding it open and assuring high hydraulic conductivity. The present study proposes a different approach for hydraulic fracturing. Fractures with infinity conductivity are formed and used to further improve the production of highly permeable formations as well as to produce long fractures in naturally fractured formations. Naturally open fractures with infinite conductivity are usually encountered. They can be observed in rock outcrops and core plugs, or noticed by the total loss of circulation during drilling (even with low density fluids), image profiles, pumping tests (Mini-Frac and Mini Fall Off), and injection tests below fracturing pressure, whose flow is higher than expected for radial Darcian ones. Naturally occurring fractures are kept open by randomly shaped and placed supporting points, able to hold the faces of the fracture separate even under typical closing pressures. The approach presented herein generates infinite conductivity canal held open by artificially created parallel supporting areas positioned both horizontally and vertically. The size of these areas is designed to hold the permeable zones open supported by the impermeable areas. The England & Green equation was used to theoretically prove that the fracture can be held open by such artificially created set of horizontal parallel supporting areas. To assess the benefits of fractures characterized by infinite conductivity, an overall comparison with finite conductivity fractures was carried out using a series of parameters including fracture pressure loss and dimensionless conductivity as a function of flow production, FOI folds of increase, flow production and cumulative production as a function of time, and finally plots of net present value and productivity index

ECONOMIC PROFITABILITY OF THE BAKKEN, NORTH DAKOTA UNCONVENTIONAL OIL PLAYS BASED ON A TYPICAL WELL PERFORMANCE WITH CURRENT MARKET CONDITIONS

Increases in oil prices after the economic recession have been surprising for domestic oil production in the United States since the beginning of 2009. Not only did the conventional oil extraction increase, but unconventional oil production and exploration also improved greatly with the favorable economic conditions. This favorable economy encourages companies to invest in new reservoirs and technological developments. Recently, enhanced drilling techniques including hydraulic fracturing and horizontal drilling have been supporting the domestic economy by way of unconventional shale and tight oil from various U.S. locations. One of the main contributors to this oil boom is the unconventional oil production from the North Dakota Bakken field. Horizontal drilling has increased oil production in the Bakken field, but the economic issues of unconventional oil extraction are still debatable due to volatile oil prices, high decline rates of production, a limited production period, high production costs, and lack of transportation. The economic profitability and viability of the unconventional oil play in the North Dakota Bakken was tested with an economic analysis of average Bakken unconventional well features. Scenario analysis demonstrated that a typical North Dakota Bakken unconventional oil well is profitable and viable as shown by three financial metrics; net present value, internal rate of return, and break-even prices.

La exploración y producción de hidrocarburos en yacimientos no convencionales, frente a los principios de prevención, precaución y desarrollo sostenible

El presente trabajo de grado busca exponer el panorama actual de la exploración y producción de hidrocarburos en Yacimientos No Convencionales, realizada utilizando el fracturamiento hidráulico – Fracking – cómo método para lograr mejores condiciones físicas en el reservorio que permitan la extracción del recurso. El método mencionado es estudiado a la luz de los principios de prevención, precaución y desarrollo sostenible, que rigen la política ambiental en Colombia, con el objetivo de analizar los posibles impactos ambientales y sociales que se puedan generar por el desarrollo de la actividad extractiva de hidrocarburos en Yacimientos No Convencionales. Para finalmente mostrar que el Fracking como actividad está legalmente permitida en Colombia, y la legislación vigente responde a los principios mencionados previamente.

Flow of Complex Fluids in Geological Fractures

In this study, the lubrication theory is used to model flow in geological fractures and analyse the compound effect of medium heterogeneity and complex fluid rheology. Such studies are warranted as the Newtonian rheology is adopted in most numerical models because of its ease of use, despite non-Newtonian fluids being ubiquitous in subsurface applications. Past studies on Newtonian and non-Newtonian flow in single rock fractures are summarized in Chapter 1. Chapter 2 presents analytical and semi-analytical conceptual models for flow of a shear-thinning fluid in rock fractures having a simplified geometry, providing a first insight on their permeability. in Chapter 3, a lubrication-based 2-D numerical model is first implemented to solve flow of an Ellis fluid in rough fractures; the finite-volumes model developed is more computationally effective than conducting full 3-D simulations, and introduces an acceptable approximation as long as the flow is laminar and the fracture walls relatively smooth. The compound effect of shear-thinning fluid nature and fracture heterogeneity promotes flow localization, which in turn affects the performance of industrial activities and remediation techniques. In Chapter 4, a Monte Carlo framework is adopted to produce multiple realizations of synthetic fractures, and analyze their ensemble statistics pertaining flow for a variety of real non-Newtonian fluids; the Newtonian case is used as a benchmark. In Chapter 5 and Chapter 6, a conceptual model of the hydro-mechanical aspects of backflow occurring in the last phase of hydraulic fracturing is proposed and experimentally validated, quantifying the effects of the relaxation induced by the flow.

Coupled wellbore-reservoir simulation of thermal effects during multiphase CO2 injection

T2Well-ECO2M is a coupled wellbore reservoir simulator still under development at Lawrence Berkeley National Laboratory (USA) with the ability to deal with a mixture of H2O-CO2-NaCl and includes the simulation of CO2 phase transition and multiphase flow. The code was originally developed for the simulation of CO2 injection into deep saline aquifers and the modelling of enhanced geothermal systems; however, the focus of this research was to modify and test T2Well-ECO2M to simulate CO2 injection into depleted gas reservoirs. To this end, the original code was properly changed in a few parts and a dedicated injection case was developed to study CO2 phase transition inside of a wellbore and the corresponding thermal effects. In the first scenario, the injection case was run applying the fully numerical approach of wellbore to formation heat exchange calculation. Results were analysed in terms of wellbore pressure and temperature vertical profiles, wellhead and bottomhole conditions, and characteristic reservoir displacement fronts. Special attention was given to the thorough analysis of bottomhole temperature as the critical parameter for hydrate formation. Besides the expected direct effect of wellbore temperature changes on reservoir conditions, the simulation results indicated also the effect of CO2 phase change in the near wellbore zone on BH pressure distribution. To test the implemented software changes, in a second scenario, the same injection case was reproduced using the improved semi-analytical time-convolution approach for wellbore to formation heat exchange calculation. The comparison of the two scenarios showed that the simulation of wellbore and reservoir parameters after one year of continuous CO2 injection are in good agreement with the computation time to solve the time-convolution semi-analytical reduced. The new updated T2Well-ECO2M version has shown to be a robust and performing wellbore-reservoir simulator that can be also used to simulate the CO2 injection into depleted gas reservoirs.

Scavenging processes of atmospheric particulate matter: a numerical modeling of case studies

Below cloud scavenging processes have been investigated considering a numerical simulation, local atmospheric conditions and particulate matter (PM) concentrations, at different sites in Germany. The below cloud scavenging model has been coupled with bulk particulate matter counter TSI (Trust Portacounter dataset, consisting of the variability prediction of the particulate air concentrations during chosen rain events. The TSI samples and meteorological parameters were obtained during three winter Campaigns: at Deuselbach, March 1994, consisting in three different events; Sylt, April 1994 and; Freiburg, March 1995. The results show a good agreement between modeled and observed air concentrations, emphasizing the quality of the conceptual model used in the below cloud scavenging numerical modeling. The results between modeled and observed data have also presented high square Pearson coefficient correlations over 0.7 and significant, except the Freiburg Campaign event. The differences between numerical simulations and observed dataset are explained by the wind direction changes and, perhaps, the absence of advection mass terms inside the modeling. These results validate previous works based on the same conceptual model.

Kinematic analysis and membership status of TWA22 AB

Context. TWA22 was initially regarded as a member of the TW Hydrae association (TWA). In addition to being one of the youngest (approximate to 8 Myr) and nearest (approximate to 20 pc) stars to Earth, TWA22 has proven to be very interesting after being resolved as a tight, very low-mass binary. This binary can serve as a very useful dynamical calibrator for pre-main sequence evolutionary models. However, its membership in the TWA has been recently questioned despite due to the lack of accurate kinematic measurements. Aims. Based on proper motion, radial velocity, and trigonometric parallax measurements, we aim here to re-analyze the membership of TWA22 to young, nearby associations. Methods. Using the ESO NTT/SUSI2 telescope, we observed TWA22 AB during 5 different observing runs over 1.2 years to measure its trigonometric parallax and proper motion. This is a part of a larger project measuring trigonometric parallaxes and proper motions of most known TWA members at a sub-milliarcsec level. HARPS at the ESO 3.6 m telescope was also used to measure the system's radial velocity over 2 years. Results. We report an absolute trigonometric parallax of TWA22 AB, pi = 57.0 +/- 0.7 mas, corresponding to a distance 17.5 +/- 0.2 pc from Earth. Measured proper motions of TWA 22AB are mu(alpha) cos(delta) = -175.8 +/- 0.8 mas/yr and mu delta = -21.3 +/- 0.8 mas/yr. Finally, from HARPS measurements, we obtain a radial velocity V(rad) = 14.8 +/- 2.1 km s(-1). Conclusions. A kinematic analysis of TWA22 AB space motion and position implies that a membership of TWA22 AB to known young, nearby associations can be excluded except for the beta Pictoris and TW Hydrae associations. Membership probabilities based on the system's Galactic space motion and/or the trace-back technique support a higher chance of being a member to the beta Pictoris association. Membership of TWA22 in the TWA cannot be fully excluded because of large uncertainties in parallax measurements and radial velocities and to the uncertain internal velocity dispersion of its members.

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.

Cluster identification and characterization in the riser of a circulating fluidized bed from numerical simulation results

A methodology of identification and characterization of coherent structures mostly known as clusters is applied to hydrodynamic results of numerical simulation generated for the riser of a circulating fluidized bed. The numerical simulation is performed using the MICEFLOW code, which includes the two-fluids IIT`s hydrodynamic model B. The methodology for cluster characterization that is used is based in the determination of four characteristics, related to average life time, average volumetric fraction of solid, existing time fraction and frequency of occurrence. The identification of clusters is performed by applying a criterion related to the time average value of the volumetric solid fraction. A qualitative rather than quantitative analysis is performed mainly owing to the unavailability of operational data used in the considered experiments. Concerning qualitative analysis, the simulation results are in good agreement with literature. Some quantitative comparisons between predictions and experiment were also presented to emphasize the capability of the modeling procedure regarding the analysis of macroscopic scale coherent structures. (c) 2007 Elsevier Inc. All rights reserved.

Finite element computational modeling of externally bonded CFRP composites flexural behavior in RC beams

This paper focuses on the flexural behavior of RC beams externally strengthened with Carbon Fiber Reinforced Polymers (CFRP) fabric. A non-linear finite element (FE) analysis strategy is proposed to support the beam flexural behavior experimental analysis. A development system (QUEBRA2D/FEMOOP programs) has been used to accomplish the numerical simulation. Appropriate constitutive models for concrete, rebars, CFRP and bond-slip interfaces have been implemented and adjusted to represent the composite system behavior. Interface and truss finite elements have been implemented (discrete and embedded approaches) for the numerical representation of rebars, interfaces and composites.

On the finite element modeling of fatigue crack growth in pressurized cylindrical shells

A methodology for the computational modeling of the fatigue crack growth in pressurized shell structures, based on the finite element method and concepts of Linear Elastic Fracture Mechanics, is presented. This methodology is based on that developed by Potyondy [Potyondy D, Wawrzynek PA, Ingraffea, AR. Discrete crack growth analysis methodology for through crack in pressurized fuselage structures. Int J Numer Methods Eng 1995;38:1633-1644], which consists of using four stress intensity factors, computed from the modified crack integral method, to predict the fatigue propagation life as well as the crack trajectory, which is computed as part of the numerical simulation. Some issues not presented in the study of Potyondy are investigated herein such as the influence of the crack increment size and the number of nodes per element (4 or 9 nodes) on the simulation results by means of a fatigue crack propagation simulation of a Boeing 737 airplane fuselage. The results of this simulation are compared with experimental results and those obtained by Potyondy [1]. (C) 2008 Elsevier Ltd. All rights reserved.

Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.

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.