CARBONATED WATER INJECTION FOR EOR IN ONE DIMENSIONAL FLOW WITH CONSTANT PRESSURE BOUNDARIES

Similar documents
rcrit (r C + t m ) 2 ] crit + t o crit The critical radius is evaluated at a given axial location z from the equation + (1 , and D = 4D = 555.

A Study of Water and Carbonated Water Injection with Constant Pressure. Boundaries. Huan Yang

Question 2-1. Solution 2-1 CHAPTER 2 HYDROSTATICS

The calculation method of small-scale water injection multiple in water drive reservoirs

Material Balance Equations

NEW COREFLOOD INTERPRETATION METHOD BASED ON DIRECT PROCESSING OF IN-SITU SATURATION DATA

Monitoring and Predicting CO 2. Flooding Using Material Balance Equations

Y.J. Cho **, Hazim Awbi** & Taghi Karimipanah* *) Fresh AB, SWEDEN **) University of Reading, UK

Harmonic Motion (HM) Oscillation with Laminar Damping

Main Menu. SEG Houston 2009 International Exposition and Annual Meeting. Summary

th th th The air-fuel ratio is determined by taking the ratio of the mass of the air to the mass of the fuel,

1) p represents the number of holes present. We know that,

TWO-PHASE SIMULATION AND CORRELATION OF NMR MAGNETIC DECAY IN EQUILATERAL TRIANGULAR PORES

20th International Congress of Mechanical Engineering THREE-DIMENSIONAL MULTIPHASE BLACK-OIL SIMULATOR BASED ON GLOBAL MASS FRACTION

Inertial Mass of Charged Elementary Particles

Two-Phase Simulation and Correlation of NMR Magnetic Decay in Equilateral Triangular Pores.

B g. B w. C w Water compressibility (pressure -1 ) G i G p m N N p. P Pressure R p Cumulative producing gas-oil ratio (st.vol./st.

A MODEL FOR THERMOPHYSICAL PROPERTIES OF CO 2 -BRINE MIXTURES AT ELEVATED TEMPERATURES AND PRESSURES

h = 0 km for your table in (a)], and for T0 = 316 K. Then obtain a Homework 1 AERE355 Fall 2017 Due 9/1(F)

The Effect of Salt Concentration on the Conductivity and Viscosity of Binary Mixed Electrolyte Solutions

Heterogeneous versus Homogeneous electron transfer reactions at liquid-liquid interfaces: The wrong question?

Chapter 5 Generalized Material Conservation Equations and Special Cases

3. Mass Transfer with Chemical Reaction

Chemical Engineering 170 October 16, 2003 Midterm Exam Closed Book and Closed Notes One 8.5 x 11 in. page of notes allowed

Unit 3. Electrochemistry

Simulation of Push-pull Multi-output Quasi-resonant Converter

A Study to Determine Reservoir Parameters and Remaining Oil Saturation by using Inter-well Tracer Test

Lecture 2: Single-particle Motion

Scheduling of a pipeless multi-product batch plant using mixed-integer programming combined with heuristics

SIMPLE NUMERICAL METHOD FOR KINETICAL INVESTIGATION OF PLANAR MECHANICAL SYSTEMS WITH TWO DEGREES OF FREEDOM

SOFT MASSIVE SPRING Objectives: Apparatus: Introduction:

Lecture 11 DAMPED AND DRIVEN HARMONIC OSCILLATIONS. Composition of harmonic oscillations (1) Harmonic motion diff. equation is: -linear -uniform

Supporting information

Temperature Regime of Boreholes: Cementing of Production Liners

SECTION I (Multiple Choice Questions)

Monitoring and Characterization Based Upon Geophysical Onset Times*

[Abdurahman, 2(11): November, 2013] ISSN: Impact Factor: 1.852

Inverse Problems for Great Solar Energetic Particle Event: Monitoring by Network Stations and Forecasting

Study Group Report: Plate-fin Heat Exchangers: AEA Technology

Modeling Crystallization in CMSMPR: Results of Validation Study on Population Balance Modeling in FLUENT

Examiner: Dr. Mohamed Elsharnoby Time: 180 min. Attempt all the following questions Solve the following five questions, and assume any missing data

FE Analysis of a Vibrating Rigid Circular Piston in Water

Lecture 17: Free Energy of Multi-phase Solutions at Equilibrium

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA

Advanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell

Physics 321 Solutions for Final Exam

HEAT CONDUCTION IN CONVECTIVELY COOLED ECCENTRIC SPHERICAL ANNULI A Boundary Integral Moment Method

Solute permeation through lipid membrane: Statistical Rate Theory approach

INTRODUCTION TO ENZYME KINETICS

Pressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects

7.0 Heat Transfer in an External Laminar Boundary Layer

EE247B/ME218: Introduction to MEMS Design Lecture 7m1: Lithography, Etching, & Doping CTN 2/6/18

Vane geometry effect on lubrication conditions between vane tip and cam-ring in hydraulic vane machines

2 Principles of Heat Transfer and Thermodynamics N 2. and H 2. r fr. 1.4 kj/kgk for carbohydrates. 1.6 kj/kgk for proteins. c p. 1.

ON THE EFFECTIVENESS OF POROSITY ON UNSTEADY COUETTE FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT

Materials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals of Diffusion

Nanoscale Heat Transfer using Phonon Boltzmann Transport Equation

Example

Lecture XXX. Approximation Solutions to Boltzmann Equation: Relaxation Time Approximation. Readings: Brennan Chapter 6.2 & Notes. Prepared By: Hua Fan

PHASE IDENTIFICATION AND SATURATION DETERMINATION IN CARBON DIOXIDE FLOODING OF WATER FLOODED CHALK USING X-RAY COMPUTED TOMOGRAPHY

Advances in Environmental Biology

( ) kt. Solution. From kinetic theory (visualized in Figure 1Q9-1), 1 2 rms = 2. = 1368 m/s

2015 Regional Physics Exam Solution Set

ChE 471: LECTURE 4 Fall 2003

Unit 11 Solutions- Guided Notes. What are alloys? What is the difference between heterogeneous and homogeneous mixtures?

Two-Dimensional COMSOL Simulation of Heavy-Oil Recovery by Electromagnetic Heating

ANALYTICAL MODEL FOR PREDICTING STRESS-STRAIN BEHAVIOUR OF BACTERIAL CONCRETE

DIFFERENT PHYSICAL PROPERTIES OF FEW AMINO ACIDS FOR FIVE DIFFERENT TEMPERATURES IN AQUEOUS SODIUM ACETATE SOLUTION

Flow-Net-Computed Uplift Pressures Along Concrete Monolith/Rock Foundation Interface

Appendix I: Derivation of the Toy Model

Effects of piezo-viscous dependency on squeeze film between circular plates: Couple Stress fluid model

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture

PHYSICS 151 Notes for Online Lecture #23

Semester 2 AP Chemistry Unit 12

Homework 1 AERE355 Fall 2017 Due 9/1(F) NOTE: If your solution does not adhere to the format described in the syllabus, it will be grade as zero.

Analysis of Surface Flaw in Pressure Vessels by a New 3-Dimensional Alternating Method

Estimation of Thermodynamic Properties and Ionic Equilibria of Cobalt Chloride Solution at 298 K

Time varying fields and Maxwell's equations Chapter 9

Summary. Introduction. Basic equations. Simulation of gas and oil flow in hydrocarbon reservoirs. P. Cole & P. A. Collins

More Tutorial at

The Journal of Supercritical Fluids

2.5 - Case studies. Outline. Materials for oars Materials for flywheel Materials for spark plug insulators Air cylinders for trucks

ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS

Module 4: General Formulation of Electric Circuit Theory

Analysis of the Effect of Die Angle Variation on the Behavior of Temperature Profile in Extrusion Process by using ANSYS Poly Flow Software

Chapter 9 Chemical Reactions NOTES

Panel Analysis of Well Production History in the Barnett Shale

Aircraft Performance - Drag

Least Squares Optimal Filtering with Multirate Observations

Mass transport with varying diffusion- and solubility coefficient through a catalytic membrane layer

Hollow block s structural impact on the thermal behaviour of PCM wall

Displacement and Deflection Sensitivity of Gas-coupled Laser Acoustic. Detector

Kinetic Model Completeness

UC A-G Chemistry. Gorman Learning Center (052344) Basic Course Information

COUNTER FLOW HEAT EXCHANGERS' IRREVERSIBILITY MINIMIZATION

The nature acids: taste sour bases: taste bitter, feel slippery

CHEMICAL EQUILIBRIUM

Chapter 5: Diffusion (2)

Considering Cable Stretch in Logging Applications

Latz, p. Technjscpe. De lit. Observations made on board of. Dutch ships. June 1967

Transcription:

SCA013-090 1/6 CARBONATED WATER INJECTION FOR EOR IN ONE DIMENSIONAL FLOW WITH CONSTANT PRESSURE BOUNDARIES H. Yang, L.A. Jaes, T.E. Jhansen Merial University f Nefundland, St. Jhn s, NL, Canada This paper as prepared fr presentatin at the Internatinal Sypsiu f the Sciety f Cre Analysts held in Napa Valley, Califrnia,, USA, 16-19 Septeber, 013 ABSTRACT A three cpnent cpsitinal nuerical del has been develped t study the perfrance f carbnated ater injectin at cnstant pressure bundaries. This calculatin ethd includes the effects f il viscsity reductin, il vlue change and interfacial tensin variatin. T cases have been studied: High and l interfacial tensin cnditins using carbnated ater injectin. The results are cpared ith siple ater injectin. Additinal il recvery is attributed t the decrease f il viscsity and il-ater interfacial tensin. INTRODUCTION Enhanced il recvery ethds include carbn dixide (CO ) flding here the CO is iscible r iiscible ith the il at reservir cnditins. Hever, pr seep efficiency has been reprted [1] due t high bility f gas and gravity driven gas verride leading t preature CO breathrugh. Carbnated ater injectin (CWI) ay help alleviate the l seep efficiency during CO injectin, here CO is disslved in ater prir t injectin and the CO partitins t the il phase upn ater-il cntact. CWI has t ain advantages; 1) the CO disslved in the il phase changes the il viscsity and hence the bility rati, and ) experientally the interfacial tensin (IFT) beteen the ater and il phases is reduced [] resulting in an iprved verall perfrance f CWI cpared t ater injectin [3]. Early n, carbn dixide slvent flding prcess as investigated in a labratry flding experient here additinal il as recvered by the carbn dixide slutin drive [4]. A calculatin ethd has been develped fr CWI [5]. The slutin using this ethd revealed the crucial effect f viscsity reductin and il selling using CWI. Results fr high pressure icrdel experients indicate better il recvery cpared t siple ater flding [6]. Even under the iiscible CO -assisted ater flding, the results shed a significant iprveent f il recvery ver siple ater flding [6]. In this paper, a ne-diensinal three cpnent set f equatins is develped t study the reservir behavir and perfrance f CWI. A free CO gas phase is nt present. The effect f viscsity and interfacial tensin as a functin f changing pressure and CO cncentratin is cnsidered, as bth are iprtant factrs in enhanced il recvery (EOR).

SCA013-090 /6 MATHEMATICAL MODEL The atheatical del as cnstructed using the flling assuptins: 1) Fl is ne diensinal and hrizntal; ) There is n surce r sin ter beteen the injectr and prducer; 3) The reservir pressure is sufficiently high that n free CO r hydrcarbn gases exist in the reservir; 4) Initially CO is disslved in the il phase and bth phases fl siultaneusly; 5) CO partitins beteen the ater and il phases but there is n ass transfer f ater and il cpnents; 6) CO diffusin ithin the phases is ignred; 7) Equilibriu beteen il and ater saturated ith CO is reached instantaneusly; 8) The reservir fratin has cnstant prsity and pereability; and 9) The injectin and prductin pressures are ept cnstant. The cpsitinal del is based n the ass cnservatin fr each cpnent ater, CO and il. The syste f equatins representing three cpnents, t phase * siultaneus fl (here is the prsity, S is the saturatin, is the density rati beteen stc tan cnditin and reservir cnditin, c is the ass cncentratin, is the density, u is the Darcy velcity, K and r are the pereability and relative pereability respectively, ith subscripts and superscripts CO, = ater and = il) can be ritten as: u [ S (1 )] [ (1 )] * cc c * c (1) t x here u [ S (1 c )] [ (1 c )] t * c * c x u u [ S c S c ] [ c c ] t ST ST ST ST * c * c * c * c x u p,,,. ST r * K RC x The su f the t phase saturatins ust equal ne and the cncentratin f the CO in the il ( c ) is a functin f the cncentratin f CO in the ater phase and the partitin c cefficient ( ). S c () (3) S 1 (4) c. (5) c c The five equatins and five unnns: S, S,, p are slved using a blccentred IMPES ethd by finite differences. c, c c c CO transfers fr the carbnated ater t il thereby changing the il prperties during the flding. We assue that the CO in ater and il are at equilibriu based n the

SCA013-090 3/6 pressure dependent slubility in il and ater. The flling relatinship describes the partitining f the CO if ater and il ere in equilibriu ith a CO gas phase: here y c, x c, respectively, and c c x x y (6) c c c c c x are the ass fractin f CO in the gas, ater and il phases, are the le equilibriu cefficients in ater-co, il-co c systes. In this study, t-phase fl (ater and il) is assued. We used the crrelatins given by Eera and Saa [7] and Duan and Sun [8] t predict the CO slubility in crude il and ater. The relatinship f CO ass cncentratins in liquid phase is given as: c c (7) c c c c The partitin cefficient ( ) f CO in three-cpnent t-phase fluid syste is defined as the rati beteen ass equilibriu cefficients c c and can be calculated fr CO slubility dels f ater-co and il-co systes. The partitin cefficient is a functin f ttal pressure in istheral carbnated ater flding. Oil viscsity decreases ith increasing CO, thus enhancing the il bility. In the istheral reservir cnditin, il viscsity varies ainly due t the change f CO ass cncentratin and ttal pressure. We used the crrelatin presented by Eera and Saa fr viscsity as a functin f CO cncentratin [7]. The ater viscsity is inially affected by the disslved CO. Lieise, the ater viscsity is a functin f teperature, pressure and salt cncentratin but nt CO cncentratin [10]. We assue that the carbnated ater viscsity reains cnstant. The il density increases ith increased CO slubility hich results fr higher pressures [9]. A density crrelatin [7] has been develped hich is calculated fr pressure, teperature, il specific gravity and the il density at bubble pint pressure. The effect f disslved CO n ater density is assued t be negligible. The pressure dependent ater density as generated fr a PVT sftare pacage. The ater-il interfacial tensin decreases due t the utual slubility f CO. Hever, due t liited literature available n the effect f interfacial tensin in a three-cpnent liquid-liquid syste, e assue a siple linear crrelatin based n reservir teperature and CO ass cncentratin. It shs an inverse relatinship beteen interfacial tensin and teperature: 0.05 c +0. 0 0. 0043 ( T 80). (8) c Fr siplificatin, the flling crrelatin is assued beteen IFT and residual il saturatin at a l interfacial tensin regin (0.00015 N/ < < 0.00 N/) [11]:

SCA013-090 4/6 S r 0 0.039678 ( N / ) Sr. (9) RESULTS & DISCUSSION T cases are presented based n t different reservir teperature cnditins. The ai f this study is t evaluate the effects f CWI under bth l (case 1) and high (case ) IFT cnditins. The cases ere develped under cnstant pressure bundary cnditins and cpared t WI. The paraeters used in the case study are list in Table. 1. The relative pereabilities used reflect a ater et syste. The saturatin prfiles f CWI in case 1 (l teperature, high IFT), case (high teperature, l IFT), and WI at breathrugh tie (38 days) are shn by Fig. 1. Because the aunt f residual il ill nt be affected in the high IFT (> 0.00 N/), the axiu ater saturatin (1-S r ) in case 1 is the sae as the ne in WI. In the high teperature syste, IFT decreases ith increasing CO cncentratin; hence, re il has been recvered under ler IFT cnditins. This can be verified by ater saturatin prfile f case here the ater saturatin in case is uch higher clse t the injectin pint. Hever, after a sharp reductin the curve flls the sae saturatin prfile as it des in case 1 due t a ler CO slubility. The verlapping ater saturatin prfile f the three prcesses in Fig. 1 iplies the sae fluid behavir. This cnfirs the calculatin results fr previus r [5] hich stated that due t the cntact ith il the initial carbnated ater injected lses its CO and then prceeds as plain ater. Thus, the CO ves behind the pure ater in CWI. Althugh the residual il saturatin stays cnstant ith bth CWI and WI under the high IFT, an additinal il recvery is bserved ith CWI. Fig. shs a cparisn beteen ater saturatin prfiles f case 1 and plain ater injectin after 00 days at injectin. The reservir rc is saturated ith re ater using CWI cpared t WI, as tie prgresses. This is ainly because f the il viscsity reductin ith increasing CO cncentratin hich can be seen fr Fig. 3. The teperature in case is uch higher than the ne in case 1 decreasing the interfacial tensin accrding t eqtn. 8. Fig. 4 shs the change f residual il ith CO slubility in the il phase. Once the CO ass cncentratin reaches 56%, l IFT is established leading t ler residual il saturatin. The cuulative aunt f il prduced in case 1, case and ater flding are pltted in Fig. 5. The CWI, ith a l IFT, has the best result flled by CWI in a high IFT cnditin. Cpared t CWI, less il can be recvered by ater flding. CONCLUSION In this paper a cpsitinal del as develped t study the perfrance f carbnated ater injectin ith cnstant pressure bundaries. A cparisn as ade beteen carbnated ater flding and plain ater flding. The results sh that

SCA013-090 5/6 additinal il recvery using carbnated ater injectin is due t the il viscsity reductin by disslved CO. In additin, if the IFT can be reduced dn t a certain regin, significant increase in il recvery ill be bserved due t reductin in residual il saturatin. Hever, since this certain regin f IFT can hardly be achieved by regular CWI, the il viscsity reductin plays the ain rle in il recvery enhanceent. Figure 1. Water saturatin after 38 days Figure. Water saturatin after 00 days Figure 3. CO ass cncentratin & il viscsity after 00 days fr case 1 Figure 4. CO ass cncentratin & residual il saturatin after 00 days fr case Figure 5. Oil prductin after 00 days

SCA013-090 6/6 Table 1. Paraeters used in case study Data CWI CWI WI Data Case 1 Case Case 1 Case T ( C) 80 50 80 ST (g/ 3 ) 874. in c 0.0386 _ ST (g/ 3 ) 999 c c 0.578 _ in c c 0.0018718 _ res c c 1.9 10-4 _ r res c r 1S Sr 0.8( ) 1S S c r S r 0 0.3 S c 0.5 L () 100 r r WI S Sr 0.( ) 1S S p in (MPa) 3. 3.3 ϕ 0.18 p res (MPa) 3.1 K ( ) 1 10-1 p ut (MPa) 3.0 μ i (Pa.s) 0.009 μ (Pa.s) 0.001 REFERENCES 1. Patel, P., R. Christan, and J. Gardner, Investigatin f Unexpectedly L Field-Observed Fluid Mbilities during Se CO Tertiary Flds, SPE, Nv 1987, Vl., N. 4, 507-513.. Mungan, Necttin, Rle f Wettability and Interfacial Tensin in Water Flding, SPE, Jun 1964, Vl. 4, N., 115-13. 3. Dng, Y., B. Dindru, and C. Ishizaa, An Experiental Investigatin f Carbnated Water Flding, SPE Internatinal, 011. 4. Hl, L., Carbn Dixide Slvent Flding fr Increased Oil Recvery, SPE, 1959, Vl. 16, 5-31. 5. Nevers, D., Nle, A Calculatin Methd fr Carbnated Water Flding, SPE, Mar 1964, Vl. 4, N. 1, 9-0. 6. Shrabi, M., M. Riazi, and M. Jailahady, Enhanced Oil Recvery and CO Strage by Carbnated Water Injectin, IPTC, Dec 009. 7. Eera, M., H. Sara, A Genetic Algrith-Based Mdel t Predict CO-il Physical Prperties fr Dead and Live Oil, JCPT, Feb 008, Vl. 47, N.. 8. Duan, Z., R. Sun, An Iprved Mdel Calculating CO Slubility in Pure Water and Aqueus NaCl Slutins fr 73 t 533K and fr 0 t 000 bar, Che. Gel., 003, Vl. 193, 57 71. 9. Deruiter, R., L. Nash, and M. Singletary, Slubility and Displaceent Behavir f a Viscus Crude ith CO and Hydrcarbn Gases, SPERE, May 1994, Vl. 9, N., 101-106. 10. Kestin, J., C. Liffrd, and H. Khalifa, Viscsity f Aqueus Ptassiu Chlride Slutins in the Teperature Range 5-150 and the Pressure Range 0-30 MPa, J. f Cheical and Engineering Data, 1978, Vl. 4,N., 11-16. 11. Kuar, S., S. Trabzadeh, and L. Handy, Relative Pereability Functins fr High and L Tensin Syste at Elevated Teperatures, SPE Califrnia Reginal Meeting, Mar 1985. c r