SGP - TR - 30 PROCEEDINGS FOURTH WORKSHOP GEOTHERMAL RESERVOIR ENGINEERING. Editors. December13-15, , 1978 SGP - TR - 30 CONF

Similar documents
AN EVALUATION OF JAMES' EMPIRICAL FORMULAE FOR THE DETERMINATION OF TWO-PHASE FLOW CHARACTERISTICS IN GEOTHERMAL WELLS

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

Compressibility Effects

Appendix I: Derivation of the Toy Model

Assume that the water in the nozzle is accelerated at a rate such that the frictional effect can be neglected.

Lecture 12: Chemical reaction equilibria

NAME TEMPERATURE AND HUMIDITY. I. Introduction

Theoretical study of third virial coefficient with Kihara potential

Numerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment

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

Methods for Determination of Mean Speckle Size in Simulated Speckle Pattern

5 th grade Common Core Standards

COASTAL ENGINEERING Chapter 2

OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION

Thermodynamics Partial Outline of Topics

Comparing Several Means: ANOVA. Group Means and Grand Mean

ChE 471: LECTURE 4 Fall 2003

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

General Chemistry II, Unit I: Study Guide (part I)

Process Engineering Thermodynamics E (4 sp) Exam

More Tutorial at

Module 4: General Formulation of Electric Circuit Theory

Fill in your name and ID No. in the space above. There should be 11 pages (including this page and the last page which is a formula page).

Calculating the optimum pressure and temperature for vacancy minimization from theory; Niobium is an example. Jozsef Garai

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

Verification of Quality Parameters of a Solar Panel and Modification in Formulae of its Series Resistance

Drought damaged area

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 11 (3/11/04) Neutron Diffusion

Synchronous Motor V-Curves

Lecture 24: Flory-Huggins Theory

CHAPTER 8 APPLICATIONS OF EMPIRICAL FETCH-LIMITED SPECTRAL. FORMULAS to GREAT LAKES WAVES 1. Paul C. Liu, Mi ASCE

Chapter 4. Unsteady State Conduction

NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION

Modeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function

Problem 1 Known: Dimensions and materials of the composition wall, 10 studs each with 2.5m high

7 TH GRADE MATH STANDARDS

Determination of Static Orientation Using IMU Data Revision 1

Sections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.

Chapter 17: Thermodynamics: Spontaneous and Nonspontaneous Reactions and Processes

Measurement of Radial Loss and Lifetime. of Microwave Plasma in the Octupo1e. J. C. Sprott PLP 165. Plasma Studies. University of Wisconsin DEC 1967

bulk velocity through orifice,

Kinematic transformation of mechanical behavior Neville Hogan

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

Surface and Contact Stress

Analysis on the Stability of Reservoir Soil Slope Based on Fuzzy Artificial Neural Network

BOUNDED UNCERTAINTY AND CLIMATE CHANGE ECONOMICS. Christopher Costello, Andrew Solow, Michael Neubert, and Stephen Polasky

Entropy. Chapter The Clausius Inequality and Entropy

(Communicated at the meeting of January )

Use a lens holder fabricated from SiC. SiC has a larger CTE than C-C, i.e. it is better matched to the SFL6.

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

39th International Physics Olympiad - Hanoi - Vietnam Theoretical Problem No. 1 /Solution. Solution

TOPPER SAMPLE PAPER 2 Class XII- Physics

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

Power Formulas for Various Energy Resources and Their Application

Ray tracing equations in transversely isotropic media Cosmin Macesanu and Faruq Akbar, Seimax Technologies, Inc.

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

, which yields. where z1. and z2

Thermodynamics EAS 204 Spring 2004 Class Month Day Chapter Topic Reading Due 1 January 12 M Introduction 2 14 W Chapter 1 Concepts Chapter 1 19 M MLK

Simulation of the Coating Process

Department of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets

INTERNATIONAL CENTRE FOR THEORETICAL PHYSICS

Inflow Control on Expressway Considering Traffic Equilibria

A mathematical model for complete stress-strain curve prediction of permeable concrete

MODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b

Feasibility Testing Report Muscle Optimization Stage I

Electric Current and Resistance

ANALYTICAL SOLUTIONS TO THE PROBLEM OF EDDY CURRENT PROBES

Comparison of two variable parameter Muskingum methods

Instructions: Show all work for complete credit. Work in symbols first, plugging in numbers and performing calculations last. / 26.

Review of the Roll-Damping, Measurements in the T-38 Wind Tunnel

1/2 and e0 e s ' 1+ imm w 4 M s 3 πρ0 r 3 m. n 0 ktr. .Also,since n 0 ktr 1,wehave. 4 3 M sπρ 0 r 3. ktr. 3 M sπρ 0

Applying Kirchoff s law on the primary circuit. V = - e1 V+ e1 = 0 V.D. e.m.f. From the secondary circuit e2 = v2. K e. Equivalent circuit :

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

Fundamental Concepts in Structural Plasticity

BASIC DIRECT-CURRENT MEASUREMENTS

MXNSSTRY OF SUPPLY AERONAUTICAL RESEARCH COUNCIL REPORTS AND MEMORANDA. of 7 a Pitot ~ntake the Air Stream. X. PLACE and R.

7.0 Heat Transfer in an External Laminar Boundary Layer

Lead/Lag Compensator Frequency Domain Properties and Design Methods

Computational modeling techniques

CHEM Thermodynamics. Change in Gibbs Free Energy, G. Review. Gibbs Free Energy, G. Review

Chapter 39. A GUIDE TO THE DESIGN OP AIR BUBBLERS FOR MELTING ICE Simon Ince Hydraulics Section, National Research Council Ottawa, Canada

Unit code: H/ QCF level: 5 Credit value: 15 OUTCOME 3 - STATIC AND DYNAMIC FLUID SYSTEMS TUTORIAL 3 - VISCOSITY

Course Stabilty of Structures

Thermodynamics and Equilibrium

Biplots in Practice MICHAEL GREENACRE. Professor of Statistics at the Pompeu Fabra University. Chapter 13 Offprint

ACOUSTIC LOGGING GUIDED WAVES IN TRANSVERSELY ISOTROPIC FORMATIONS

Heat Effects of Chemical Reactions

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.

Design and Analysis of Gas Turbine Blade by Potential Flow Approach

An Eigenvalue Based Acoustic Impedance Measurement Technique

On Huntsberger Type Shrinkage Estimator for the Mean of Normal Distribution ABSTRACT INTRODUCTION

I. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is

Interference is when two (or more) sets of waves meet and combine to produce a new pattern.

THE TOPOLOGY OF SURFACE SKIN FRICTION AND VORTICITY FIELDS IN WALL-BOUNDED FLOWS

Chapter 23 Electromagnetic Waves Lecture 14

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

CALCULATION OF BRAKING FORCE IN EDDY CURRENT BRAKES

3. Design of Channels General Definition of some terms CHAPTER THREE

A study on GPS PDOP and its impact on position error

A.P. CHEMISTRY. SOLUTIONS AND ACID BASE CHEMISTRY. p 1

Transcription:

SGP - TR - 30 SGP - TR - 30 CON-781222-26 PROCEEDINGS OURTH WORKSHOP GEOTHERMAL RESERVOIR ENGINEERING Paul Paul Krugerand and Henry.. Ramey, Ramey., r. r. Editrs December13-15, 13-15., 1978 DISTRIBUTION O THiS [;(~l~1e1h IS UNLlrl1lTED

AN EVALUATION O AMES' EMPIRICAL ORMULAE OR THE DETERMINATION O TWO-PHASE LOW CHARACTERISTICS IN GEOTHERMAL WELLS P. Cheng and M. Karmarkar University f Hawaii Hnlulu, Hawaii 96822 Intrductin One f the mst ecnmical and simple methds f determinatin f tw-phase flw parameters in gethermal well testing is the s-called ames' methd [1,2]. The methd cnsists f the measurements f lip pressure (p), and the flw rate f water (w) by a cnventinal weir. The stagnatin enthalpy (h ) is then determined frm a plt shwing versus w/po.96 which is gmpirically determined by ames [1,2]. The 0 mass flw rate is then determined frm the fllwing empirical frmula 0.96 G = 11,400 ~1.102 (1) where G is the ttal mass flw rate in 1b m /sec-ft 2, p is the lip pressure in psia, and h O is the specific enthalpy in BTU/lb m. The abve relatin is empirically determined fr discharge pressure up t 64 psia and pipe diamters up t 8". r pipe diamters smaller than 0.2", it has been suggested that the value f 11,400 be replaced by 12,800. In view f the widespread use f the ames' methd, it is imprtant t assess its accuracy and range f applicability. h Tw-Phase Critical lw Thery In this paper we shall cmpare the wellbre discharge characteristics btained frm ames' empirical frmulae t thse predicted by auske's tw-phase critical flw thery [3]. auske suggested that in tw-phase flw the maximum discharge rate may nt necessarily be accmplished by a shck frnt. He prpsed that at the critical flw cnditin the abslute value f the pressure gradient at a given lcatin is maximum but finite fr a given flw rate r quality, i.e., (dp/dz)g = maximum and finite,,x (2) where z is the crdinate alng the streamwise directin, and x the quality f the saturated mixture. Under the assumptins f (i) annular flw pattern, (ii) tw-phases -207-

in thermal equilibrium, (iii) negligible frictinal lss, and (iv) nedimensinal steady flw, auske (3] btained the fllwing analytical expressin fr the critical flw rate f a saturated mixture: (3) - 2xk 2 +k 2» dx/dp], g = 32.2 c 1b -ft m -=-----=-2 ' and 1b f -sec k v g and denting the specific vlume f the saturated vapr and liquid respectively. Thus, the critical flw rate can be calculated frm Eq.(3) if the steam quality and the lip pressure are knwn. The crrespnding stagnatin enthalpy can be determined frm 2 3 2 2 h f (l-x) + h x + (G /2g )«x v /Rg ) g c g 3 2 2 +(l-x) v f /(l-rg) )/, (4) where R is the gas vid fractin which is related t steam quality [4]. In cmpa~isn with experimental data, Levy [4] fund, hwever, that using Eq.(4) fr the cmputatin f h O wuld lead t under-predictin f the mass flw rate. r this reasn we shall cmpute the stagnatin enthalpy n the basis f a hmgeneus mdel, i.e., 2 2 h O = h f (l-x) + hgx + G v h /2g c, (5) where v h = v f (l-x) + vgx and = 778 ft-ibm/btu The weir flw rate is then determined frm w = G(l - x) (6) Results and Discussin r a given set f values f lip pressure and steam quality and with the data f saturated steam-water prperties [5], Eq.(3) can be used fr the cmputatin f ttal mass flw rate G. The stagnatin where Q [(l-x+kx) x (dv/dp) + (v g (1+2kx-2x) +v f (2xk-2k -208-

enthalpy and the weir flw rate are then determined frm Eqs.(5) and (6). The results f the cmputatins fr the lip pressure frm 14.7 psia t 200 psia fr gethermal well testing applicatins are pltted in igs. 1 and 2. When the lip pressure and the weir flw rate are measured in a gethermal well test. the stagnatin enthalpy f the reservir. the steam quality at the well head. and the ttal mass flw rate can easily be determined frm these plts. T assess the accuracy f the ames' methd. calculatins were carried ut fr five different sets f lip pressure and weir flw rate using ames' empirical frmulae and auske's theretical predictin (i.e igs. 1 and 2). The results fr ttal mass flw rate. the stagnatin enthalpy. and the steam quality are tabulated in Table 1 fr cmparisn. It is shwn that the results based n the tw methds differ within 8%. References L ames. R "Measurement f Steam Water Mixtures Discharging at the Speed f Sund t the Atmsphere. New Zealand Engineering. pp. 437-441 (1966). 2. ames. R. "Steam-Water Critical lw Thrugh Pipes." Prc. Inst. Mech. Engrs. v. 176. pp. 741-748 (1962). 3. auske. R "Cntributin t the Thery f Tw-Phase. One Cmpnent Critical lw." Argnne Natinal Labratry. Rept. N. ANL-6633 (1962). 4. Levy. S "Predictin f Tw-Phase Critical lw Rate." A.S.M.E. Paper 64-RT-8. 5. Keenan.. and Kays. Keys,. G., Steam Tables, hn Wiley and Sns. Inc New Yrk. N.Y. (1969). -209-

600 500, u."... 400 E -0, 0:: 300 ~ LL Lip Pressure (psia) 3= 200 100 400 500 600 700 800 900 1000 1100 1200 1300 Stagnatin Enthalpy, h (Btu/Ibm) ig. 1. Weir lw Rate vs. Stagnatin Enthalpy at Selected Values f Lip Pressure Accrding t auske's Thery --l-~,...---,----,----.---r---r-...--.--'--:~~~~~-- -210-

600 500 I U." '- 400 E...0 -. ō ly: 300 ~ LL Lip Pre ssure (psia) ::; 200 100 O-+---r-.---~---,---,------.-~..,.----,----r-.-.,...--,-----r=~-t- 0.2 0.3 0.4 0.5 06 0.7 Steam Quality, X ig. 2. Weir lw Rate vs. Steam Quality at Selected Values f Lip Pressure Accrding t auske1s Thery. -211-0.8 0.9 1.0

Table 1. COMPARISON O RESULTS BASED ON THE AMES' METHOD AND AUSKE'S ANALYTICAL MODEL Case p w h G x Methd (psia) (1 b m /sec-ft2) BTU71bm lb m /sec-ft 2 14.7 40 730.90 88.44.54 auske () 800.00 95.13.58 ames () 2 25.0 85.5 698.78 750.00 164.59 170.06.48.50 3 60.0 226.0 697.79 403.39 715.00 415.42.44.46 4 100.0 105.0 1004.52 419.84 985.00 476.60.75.78 5 150.0 53.0 1148.87 1130.00 523.41 590.00.90.90-212-