Leaver, Jonathan D. (1986

Download as pdf or txt
Download as pdf or txt
You are on page 1of 199

SGP-TR-95

A Technical Review of Interference


Testing with Application in the
Ohaaki Geothermal Field
Jonathan D. Leaver

March 1986

Financial support was provided through the Stanford


Geothermal Program under Department of Energy Contract
No. DE-AT03-80SF11459 and by the Department of Petroleum
Engineering, Stanford University

Stanford Geothermal Program


Interdisciplinary Research in
Engineering and Earth Sciences
STANFORD UNIVERSITY
Stanford, California
ACKNOWLEDGEMENTS

This project would not have been possible without the professional advice and support of As-

SOC. Prof. Abraham Sageev and Prof. Henry J. Ramey Jr. of Stanford University Petroleum

Engineering Department.

The assistance of Mr. Paul F. Bixley, Ministry of Works and Development, Wairakei, New

Zealand, and Dr. Malcolm Grant, Department of Scientific and Industrial Reseadch, Wellington,

New Zealand, in obtaining the data and background reports was also greatly appreciated.

Test data for the "C" series of tests was collected by Mr. David M. Wilson, M.W.D., Wairak-

ei, and kindly filtered and supplied by Dr. Mark J. McGuinness, Applied Math. Division,

D.S.I.R., Wellington.

Financial support was provided by M.W.D., New Zealand Energy Department, Stanford Geoth-

ermal Program, DOE Contract No. DE-AT02-80SF11459, and by Stanford University.


To
Cynthia
and
Trudie, Rosanne and Luke.

...
-111-
ABSTRACT

The chronology and mathematical development of interference testing and relatdd earth tide arid

compressibility effects are reviewed. Hydrological information from field tests performed in

the Ohaaki geothermal field since first development in 1965 is summarized.

Data from 12 interference tests (four on the same doublet) performed in the de&pOhaaki reset1

voir between 1979 and 1983 are analyzed using conventional log-log and semf-log type cunf

matching techniques along with recently developed semi-log type curve matdhing technique$
e
I
which allow linear boundary detection without the necessity to develop two semi-log straiglnr

lines or the need to know reservoir parameters.

Data in nine tests were recorded using water level chart recording devices &le data in thq

remaining three were recorded using quartz crystal pressure gauges. Both instruments have a

resolution of about 100Pa.

Some tests show the presence of a no-flow boundary for which the inference ellipses haw$

been located. Interpretation is made on the probable location of the no-flow boundary. A study

of the significance of early time data and other factors shows that vastly different reservoir

parameters can be obtained unless careful interpretation is made.

A hydrological model of the Ohaaki geothermal field is presented.


TABLE OF CONTENTS

Page

ABSTRACT ..................................................................................................................... ii

TABLE OF CONTENTS ............................................................................................... ...


111

LIST OF FIGURES ........................................................................................................ vi

LIST OF TABLES ........................................................................................................... xii

.
1 INTRODUCTION ............................................................................................. ..........
+ 1

2.MATHEMATICAL DEVELOPMENT ................................................ ..........,.......... 5

2.1 LOG-LOG TYPE CURVES .................................................................................. 5

2.1.1 Line Source (Theis 1935) ...................................................................... 5

2.1.2 Linear Boundary Drawdown (Stallman 1952) ......................................... 6

2.1.3 Linear Boundary Drawdown and Buildup (Eipper 1985) ........................ 8

2.2 SEMI-LOG TYPE CURVES ................................................................................ 8

2.2.1 Linear Boundary Drawdown (Sageev et al . 1985)................................... 8


2.2.2 Linear Boundary Buildup (Fox 1984) ...................................................... 10

2.3 INFERENCE ELLIPSE .......................................................................................... 11

2.4 TRANSMISSIVITY AND STORATIVITY........................................................... 12

.
3 BAROMETRIC AND EARTH TIDE EFFECTS ..................................................... 13

3.1 INTRODUCTION................................................................................................... 13

3.2 BAROMETRIC EFFECTS ..................................................................................... 14

3.3 EARTH TIDES (Sageev et a2. 1986) .................................................................... 14

3.4 PRESSURE SIGNIFICANCE LEVEL................................................................... 18

4. COMPRESSIBILITY .................................................................................................. 20
5. EARLY TIME DATA MATCHING (Sageev et al . 1986) ....................................... 22

5.1 THEORY ................................................................................................................ 22

5.2 MISSING DATA .................................................................................................... 27

6.HYDROLOGY OF THE OHAAKI GEOTHERMAL FIELD ............................... 31

6.1 BACKGROUND .................................................................................................... 31

6.2 GENERAL ............................................................................................................. 33

6.3 EAST BANK ......................................................................................................... 34

6.4 WEST BANK ........................................................................................................ 35

6.5 BROADLANDS DACITE .................................................................................... 36

6.6 OHAAKI RHYOLITE ................................................................................ ........... 38


6.7 COMMENT ........................................................................................................... 40

.
7 ANALYSIS OF INTERFERENCE TESTS DATA IN THE OHAAKI

WEST BANK ........................................................................................................ 41

7.1 MEASUREMENT ACCURACY .............................................................. .,........... 41


7.2 WELL LOCATION .................................................................................... + ........... 42 .
~

7.3 DATA FROM WATER LEVEL CHART RECORDERS ........................ 4 ........... 43 I


I
7.3.1 Test B1: BR13 Response to BR23 Discharge (Leaver et a1. 198$)......... 43 i
7.3.2 Test B2: BR23 Response to BR13 Discharge .......................................... 53 I

7.3.3 Test B3: BR23 Response to BR19 Discharge .............................. ........... 58


7.3.4 Test B4: BR23 Response to BR13 Injection ............................................ 63

7.3.5 Test B6: BR23 Response to BR20 Discharge .......................................... 69

7.3.6 Test B7: BR34 Response to BR31 Discharge .......................................... 75

7.3.7 Test B8: BR34 Response to BR23 Discharge .......................................... 80

7.3.8 Test B9: BR34 Response to BR19 Discharge .......................................... 84

7.3.9 Test B10: BR23 Response to BR13 Discharge ........................................ 91

7.4 DATA FROM QUARTZ CRYSTAL GAUGES .................................................. 97

7.4.1 Test C1: BR13 Response to BR20 Discharge .......................................... 97

-iv-
7.4.2 Test C2: BR23 Response to BR20 Discharge .......................................... 103

7.4.3 Test C3: BR34 Response to BR20 Discharge .......................................... 109

8.RESULTS .................................................................................... ................................116

8.1 TRANSMISSIVITIES AND STORATIVITIES ........................ ........................ ...116


8.2 LINEAR BOUNDARY LOCATION ..................................................................... 118

9.DISCUSSION .............................................................................................................. 120

9.1 GENERAL (Sugeev et a1. 1986) ........................................................................... 120

9.2 ANALYSIS TECHNIQUES ................................................................................... 120

9.3 DATA MATCHING ............................................................................................... 121

9.4 RESERVOIR PARAMETERS .................................................................. .,........... 122


9.5 OFTIMIZING TEST PARAMETERS ........................................................ ,........... 122
.
10 CONCLUSIONS ....................................................................................................... 124

NOTATION .............................................................................................................. 125

REFERENCES ........................................................................................ ....,........... 127


APPENDICES ................................................................................................. i. ...........135
A. COMPUTER PROGRAMME FOR THEIS LINE SOURCE SOLUTION

(Eipper 1985) ......................................................................................................... 135


B . SPECIFIC PRESSURE VERSES TIME DATA ..................................................... 136

-V-
LIST OF FIGURES

Page

2.1 Log-log Drawdown Response of a Line Source Well in an Infinite Reservbir

(after Theis 1935) ...................................................................................... ........... 6


2.2 Log-log Drawdown Response for Reservoir with Linear Boundary

(after Stallman 1952) ............................................................................................. 7

2.3 Log-log Drawdown and Buildup Response for Reservoir with Linear

Boundary (after Eipper 1985) ................................................................. 4 ........... 9 1


I

2.4 Semi-log Drawdown Response for Reservoir in a Semi-infinite System.....,........... 9 ~

I
2.5 Semi-log Drawdown Response for Reservoir with Linear Boundary

(after Sageev 1985) ............................................................................................... 10

2.6 Semi-log Buildup Response for Reservoir with Linear No-Flow Boundary

(after Fox 1984) .................................................................................................... 11

2.7 Schematic Showing Notation for Linear Boundary Detection .............................. 12

3.1 Atmospheric Pressure Variations During the "C"Tests ......................................... 15

3.2 Barometric and Earth Tide Fluctuations During the "C"Tests .............................. 15

3.3 Modelling of the Earth Tide Response Using Two Superposed

Sine Waves On a Linear Barometric Trend (after Sageev et al. 1986) ..,........... 16

3.4 Match of the Synthetic and Real Earth Tide Responses

(after Sageev et al. 1986) .......................................................................... 1........... 17

3.5 Superposition of the Line Source With a Single Sine Wave

(after Sageev et al. 1986) ...................................................................................... 18

-vi-
3.6 Test C3:Log-log Matches of Data to Stallman Type Curves ................................ 19

4.1 Variation in Water Compressibility with Temperature and Pressure

(after Home 1980).................................................................................................. 21

4.2 Variation in Steam Compressibility with Pressure (after Home 1980) ....,............21

5.1 Log-log type curve for a Semi-infinite system (after Stallman 1952 ).................. 22

5.2 Semi-log type curve for a Semi-infinite system ....................................................... 23

5.3 Log-log Match of the Line Source To the Curve for r2/rl = 1.1

(after Sageev et al . 1986) ...................................................................................... 24

5.4 Log-log Response of the Line Source and a Semi-infinite System

With r2hl = 2.0 (after Sageev et al . 1986) ............................................. ............ 25


5.5 Log-log Match of the Line Source and a Semi-infinite System With

r2/rl = 2.0 (after Sageev et a1. 1986) ................................................................... 26

5.6 Log-log Type Curve for Dimensionless Pressures and Semi-log

Dimensionless Derivatives for a Semi-infinite System

(after Sageev et al . 1986) ......................................................................... ............ 27


5.7 Log-log Graph of Data From Test C3 ....................................................... i ............ 28

5.8 Hypothetical Cases for Missing Pressure Data for r2/rl = 2.0

(after Sageev et al. 1986) ...................................................................................... 29

5.9 Four Possible Log-log matches of Test C3 data to the Stallman Type Curves

(after Sageev et al . 1986) ...................................................................................... 30

6.1 Locality Map for the Ohaaki Geothermal Field ....................................................... 31

6.2 Well Layout of the Ohaaki Geothermal Field ........................................... ........... 32


6.3 North-West South-East Geologic Cross Section Showing Postulated

Flows in the Ohaaki Reservoir .................................................................. ........... 33

-vii-
6.4 West Bank Well Pressure Profiles (after Bixley 1982) ......................................... 36

6.5 West Bank Well Grouping From Pressure Gradient Characteristics

(after Bixley 1982) ................................................................................................. 37

6.6 Areal View of the Ohaaki Field Showing West Bank Gravity Anomaly ............... 38
6.7 Areal View of the Ohaaki Field Showing Significant Geology and Well

Layout for West Bank Rhyolite Interference Tests

(after McGuinness 1985) ....................................................................................... 39

7.1 Source and Observation Well Location for Interference Tests ............................... 42

7.2 Test B1: Cartesian Graph of the Interference Data

(after Leaver et a1. 1985) ...................................................................................... 43

7.3 Test B1: Log-log Match of Data to Eipper Type Curves

(after Leaver et a1. 1985) ...................................................................................... 45

7.4 Test B1: Semi-log Match of Data to Semi-infinite Type Curves

(after Leaver et al . 1985) ...................................................................................... 47

7.5 Test B1: Match of Buildup Data on a Horner Graph

(after Leaver et a1. 1985) .......................................................................... b ........... 47


7.6 Test B1: Semi-log Match of Drawdown Data to Sageev et a1. Type Curves

(after Leaver et a1. 1985) ...................................................................................... 49

7.7 Test B1: Incorrect Semi-log Match of Buildup Data to Fox Type Curve

(after Leaver et a1. 1985) ...................................................................................... 50

7.8 Test B1: Correct Semi-log Match of Buildup Data to Fox Type Curve

(after Leaver et a1. 1985) ...................................................................................... 51

7.9 Test B1: Inference Ellipse Location (after Leaver et a1. 1985) ............................ 52

7.10 Test B2: Cartesian Graph of the Interference Data ............................................... 53

7.1 1 Test B2: Log-log Match of Data to Stallman Type Curves ...................... ........... 56
7.12 Test B2: Semi-log Match of Data to Semi-infinite Type Curves ......................... 57

7.13 Test B3: Cartesian Graph of the Interference Data ............................................... 58

7.14 Test B3: Log-log Match of Data to Stallman Type Curves ................................... 60

7.15 Test B3: Semi-log Match of Data to Semi-infinite Type Curves ......................... 61

7.16 Test B3: Inference Ellipse Location ...................................................................... 62

7.17 Test B4: Cartesian Graph of the Interference Data ............................................... 63

7.18 Test B4: Log-log Match of Data to Stallman Type Curves .................................. 65

7.19 Test B4: Semi-log Match of Data to Semi-infinite Type Curves ......................... 67

7.20 Test B6: Cartesian Graph of the Interference Data ............................................... 69

7.21 Test B6: Log-log Match of Data to Eipper Type Curves for r2/r1-4.0 ....i........... 71
7.22 Test B6: Log-log Match of Data to Line Source Solution ................................... 72

7.23 Test B6: Semi-log Match of Data to Semi-infinite Type Curves for r2/r144.0 .... 73
7.24 Test B6: Semi-log Match of Data to Line Source Solution ................................... 74

7.25 Test B6: Match of Buildup Data on a Horner Graph .............................. ........... 74

7.26 ...............................................
Test B7: Cartesian Graph of the Interference Data 75

7.27 Test B7: Log-log Match of Data to Stallman Type Curves ................................... 77

7.28 Test B7: Semi-log Match of Data to Semi-infinite Type Curves ......................... 78

7.29 Test B7: Inference Ellipse Location ...................................................................... 79

7.30 Test B8: Cartesian Graph of the Interference Data ............................................... 80

7.31 Test B8: Log-log Match of Data to Line Source Solution..................................... 82

7.32 Test B8: Semi-log Match of Data to Line Source Solution ................................... 83

7.33 ............................................... 84
Test B9: Cartesian Graph of the Interference Data

7.34 Test B9: Log-log Match of Data to Eipper Type Curves for rz/rl -1.5 ............... 86

7.35 Test B9: Log-log Match of Data to Line Source Solution......................... ........... 87

-ix-
7.36 Test B9: Semi-log Match of Data to Semi-infinite Type Curves for r2/r1t1.5 .... 88
7.37 Test B9: Semi-log Match of Data to Line Source Solution ................................... 89

7.38 Test B9: Match of Buildup Data on a Horner Graph for rdr1=1.5 ....................... 90

7.39 Test B9: Match of Buildup Data on a Horner Graph to Line Source ................... 90

7.40 Test B10: Cartesian Graph of the Interference Data ............................................. 91

7.41 Test B10: Log-log Match of Data to Line Source Solution................................... 93

7.42 Test B10: Semi-log Match of Data to Line Source Solution................................. 94

7.43 Test B10: Semi-log Match of Drawdown Data to Sageev et al . Type C w e s ..... 95
7.44 Test B10: Inference Ellipse Location .................................................................... 96

7.45 ...............................................
Test C1: Cartesian Graph of the Interference Data 97

7.46 Test C1: Log-log Match of Data to Eipper Type Curves ..................................... 99

7.47 Test C1: Semi-log Match of Data to Semi-infinite Type Curves ......................... 100

7.48 Test C1: Match of Buildup Data on a Horner Graph .......................................... 101

7.49 Test C 1: Inference Ellipse Location ...................................................................... 102

7.50 ................................... 1...........103


Test C2: Cartesian Graph of the Interference Data

7.5 1 Test C2: Log-log Match of Data to Eipper Type Curves ..................................... 105

7.52 Test C2: Semi-log Match of Data to Semi-infinite Type Curves ............. ...........106
b

7.53 Test C2: Match of Buildup Data on a Horner Graph .......................................... 106

7.54 Test C2: Inference Ellipse Location ...................................................................... 108

7.55 Test C3: Cartesian Graph of the Interference Data ............................................... 109

7.56 Test C3: Log-log Match of Data to Eipper Type Curves for rz/rl =1.5 ...............111
7.57 Test C3: Log-log Match of Data to Line Source Solution..................................... 111

7.58 Test C3: Semi-log Match of Data to Semi-infinite Type Curves for rz/rl c 1.5 ..113
7.59 Test C3: Semi-log Match of Data to Line Source Solution................................... 113

7.60 Test C3: Match of Buildup Data on a Horner Graph for rdrl = 1.5 .................... 114

-X-
7.61 Test C3: Match of Buildup Data on a Homer Graph to Line Source ................... 115

8.1 Inference Ellipse Locations for Tests B1, B3, B7, B10, C1, C2 ..............,...........119

-xi-
LIST OF TABLES

Page

7.1 Test B1: Test Specifications ....................................................................... ........... 44


7.2 Test B2: Test Specifications .................................................................................. 54

7.3 Test B2: Flow Rate Data ........................................................................... .......... 55


7.4 Test B3: Test Specifications ..................................................................... i. .......... 59
7.5 Test B4: Test Specifications ................................................................................. 64
I

7.6 Test B6: Test Specifications . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70


7.7 Test B7: Test Specifications .................................................................................76
7.8 Test B7: Flow Rate Data .................................................................................. 77

7.9 Test B8: Test Specifications.................................................................................. 81

7.10 Test B8: Flow Rate Data ...................................................................... ........... 82


7.11 Test B9: Test Specificatio............................................................................i ........ 85
~

7.12 Test B9: Flow Rate Data ......................................................................... i’


i .......... 86
I
7.13 Test B10: Test Specifications................................................................................. ~

92

7.14 Test C1: Test Specifications................................................................................... 98

7.15 Test C2: Test Specificatio............................................................................I ........ 104

7.16 Test C3: Test Specificatio..................................................................................... 110

8.1 Transmissivities and Storativities................................................................ 1...........1161

B1 Test B1: Specific Pressure vs Time Da................................................................ 136

B2 Test B2: Specific Pressure vs Time Data ............................................................. 138

B3 Test B3: Specific Pressure vs Time Da............................................................... 140

-xii-
B4 Test B4: Specific Pressure vs Time Data ............................................................. 142
B5 Test B6: Specific Pressure vs Time Da............................................................... 144

B6 Test B7: Specific Pressure vs Time Da.................................................. ,............145


B7 Test B8: Specific Pressure vs Time Data ............................................................. 146

B8 Test B9: Specific Pressure vs Time Data ............................................................. 147


B9 Test B10: Specific Pressure vs Time Data.............................................................. 149 ,

B10 Test C1: Specific Pressure vs Time Data.............................................................. 151

B11 Test C2: Specific Pressure vs Time Data .............................................................. 166


B12 Test C3: Specific Pressure vs Time Data.............................................................. 175 I

1
I,

-1

1. INTRODUCTION

Interference testing in geothermal fields is one of the principal tools used to determine bulk

transmissivities (kh)and storativities ($c,h) and to locate hydrological heterogeneities such as

no-flow and pressure support boundaries. Changes in the hydrology of geotheml reservoirs in

geologic time can be caused through mineral deposition, changes in the heat flux into or out o f ,

the system, crustal movements, changes in fluid components (brine, gas) moving into or out of

the reservoir and exploitation. Exploitation increases by orders of magnitude the rate of natural

processes in the reservoir, with mineral deposition and reduction in the heat of the system be-~

ing the most significant. Interference tests are both an essential and economical tool in assess-'

ment of the extractable heat capacity of a geothermal field and in monitoring changes in reser-

voir characteristics as exploitation proceeds. The past 50 years has seen steakiy progress in1

developing basic fluid flow theory to the practical techniques available in 1985, A brief chro-
I

nology of the developments in interference testing is described below.


I

Theis [1935] published the line source solution based on heat transfer analogy and used time1

superposition to simulate a finite flow period. The stated assumptions were:

"(i) Entirely homogeneous bodies.


(ii) Infinite areal extent.

(iii) Well penetrates the entire thickness of the water body.

(iv) Coeficient of transmissibility is constant at all times and in all places.

(v) Pumped well has infinitesimal diameter.

(vi) Applicable only to unconfined water bodies in which the water in the volume of

sediments through which the water table has fallen is discharged instantaneously

with the fall of the water table."

Theis [1935] also noted that:

"The eflect of boundaries can be considered by more elaborate analyses lonce they
-2-

are located."

However the formal publication of these analyses awaited a further 17 years.

Jacob [1941] circumspectly used the term "interference testing" to describe the1 change in pres1

sure in a shut-in well caused by changes in flow rate at other wells:

"The interference of two wells, in the technical sense, is determined by a cornparis- 1

on of the pa0 of one of two wells when both are $owing freely, with He pa0 of ~

each well when the other is shut o#?."

Jacob also makes mention of the match between data and a theoretical type 'curve providin$

the parameters required to determine permeability and porosity thereby possiblp being the fir6t
to use type curve matching techniques. ~

Muskat [1937] used space superposition to sum pressure changes from each ifldividual well

any point in a reservoir at steady state thereby finding the total pressure chanke at that

Kazmann [1946] used space superposition of source and image wells to detemine the

to a recharge boundary. Van Everdingen and Hursr [1949] showed that

could be described by the Diffusivity equation. The use of Laplace transforms was shown

simplify the solution process for many problems. Stallman [1952]

included the effect of a no-flow boundary during drawdown. H u m [1960] 4s a corollary

work on interference between oil fields showed how to determine pressure m a t i o n s due

two wells in a single layer heterogeneous reservoir. Barenblatt

blurt et al. [1960] introduced the concept of a double porosity model contaikng matrix an
6
fissures. Fluid was assumed to flow from the matrix to the fissures under psebdo steady sta

conditions. Mathematical development utilized the average matrix and fissure pressures meas-

ured at the same reservoir location. Hantush [1960] analytically derived type curves for th
pressure response in a reservoir bounded above by a leaky aquifer. The difference betweelp
c
these curves and the line source solution was small for compressible systems which generalljl

meant the properties of the caprock could not be determined. Witherspoon [1962] resolved thi

problem by using well measurements in the aquifer overlying the main reservipir to determidl
P
-3-

the leaky caprock properties. Ferris et af. [1962] presented methods for modelling linear fault

combinations in a single homogeneous isotropic medium. Warren and Root [I9631 simplified

the analysis of double prosity systems by describing the pressure distribution io terms of onl4

two parameters - the interporosity flow parameter (1)and a storativity ratio (a). Davis an
Hawkins [1963] used the semi-log doubling of slope characteristic for wells of constant skin id
d
a semi-infinite homogeneous isotropic medium to locate the position of a noiflow boundary.

The technique used a semi-log dimensionless pressure-time graph on which the intersection

point of the two straight lines occurred at the same dimensionless time regardlless of the dis-
I

tance to the barrier. It was applicable to both drawdown and buildup tests. Papctdopufos [19651
I
produced a solution for the analysis of interference tests in homogeneous anisotxtopic reservoirs.~

Ramey [1970] developed the solution of Papadopufos for application to petrolebm engineering/

problems. Kazerni et af. E19691 extended the Warren and Root double porosity model for a

single well test to include the pressure response at an observation well. System characteristics~
I

were determined from the difference in early time behavior between the doul$e porosity and

isotropic homogeneous models. Earfougher and R m e y [19731 published the bolution for in

terference effects in bounded systems. Tables were presented of dimensionless pessure and di]
1
mensionless time for various reservoir shapes and pressure observation points. Koefoed [1974

continued the work of Witherspoon by publishing type curves for determining the leakage fac
1
tor of a reservoir with a caprock. Jargon [1976] used a finite difference numeriical simulator ld
generate interference test data with wellbore storage and skin at the sourcd well. Resul

showed that neglecting wellbore storage causes underestimation of transmissivib and overed

mation of storativity. Increased skin damage was shown to prolong wellbore 'storage effects

Vela [1977] described the location of a linear boundary detected from a single interference tes

as a tangent to an inference ellipse. Najun'eta [1979] continued work on doubile porosity sysi
1I
tems and included the effect of transient interporosity flow in the analysis. Deruyck et al~

[1982] presented two type curves for double porosity reservoirs. These could b~ used to deter4
1
mine whether the flow from the less permeable to the more permeable medium was transient o
4
-4-

pseudo steady state and also to find the nature of the double porosity system. The curves were

equally applicable to either single layered fractured or multilayered homogeneous reservoirs.

Sugeev et ul. [1985] mathematically collapsed a semi-log version of Stallman’s type curves to

produce a single semi-log type curve applicable for drawdown tests. The technique allowed

the inference ellipse of a constant pressure or a no-flow boundary to be determiped without re-1

quiring any knowledge of reservoir parameters. The method is applicable for b6undaries locat-

ed more than 5rl from either the source or observation well. Fox [1984] used lple unpublished~
I
technique of Sogeev et al. [ 19851 to produce a single semi-log type curve which was applicable1

for buildup tests. Graphs were developed using both total time and Hornek time. Eippep
I

[1985] extended Stallman’s work and produced a set of log-log type curves shawing the effect/

of either a no-flow or constant pressure boundary on both drawdown and buildup. Sage&

E19851 examined the effects of a steam cap near a doublet and produced a method for detei

mining the probable location of the steam cap based on the transmissivity and storativity value

obtained from conventional techniques.


-5-

2. MATHEMATICAL DEVELOPMENT

The development of fundamental equations for type curves relies upon the assumption thdt

fluid flow in the reservoir can be described by the diffusivity equation. In radial coordinateb

this is: I

The equation is derived using conservation of mass, Darcy's Law, and an ehuation of stad,

and is a second order linear homogeneous partial differential equation. Due t0 the linear pr
(si-
perty, multiples of any solution to the diffusivity equation and its associated boundary con

1
tions will also be solutions. Hence space superposition of constant rate line sources is possibl

allowing the generation of single and multiple linear boundaries. Important as$umptions in th

development of the diffusivity equation are:

Horizontal radial flow,

Negligible gravity effects.

Homogeneous isotropic porous medium.

A single fluid of small and constant compressibility.

Darcy flow in the formation.

Permeability, compressibility, viscosity and porosity are independent of pressure.

LOG-LOG TYPE CURVES

2.1.1. Line Source (Theis)

The solution of the pressure response at any point in an infinite reservoir due 'to the dischxg

of a well of infinitesimal diameter at constant rate was given by Theis [19351:


I
I
-6-

The exponential integral (Ei) is defined as:

A dimensionless pressure-time graph is shown in Figure 2.1.

Q,
c
2
u)

e 1
Q
u)

-
u)
8
C
0
-0i
*1
t
0,
.-E
0

0*01
0.1 1 10 100 1oao
Dimensionless Time
Fig. 2.1 Log-log drawdown response of a Line Source Well in an infinite reservoir (after Theis 1935).

2.1.2. Linear Boundary Drawdown (Srullman)

The solution of the pressure response at any point in a semi-infinite reservoir was given by

Stallman [1952] superposing an image source in space to create the effects of a b e a r no-flow

boundary. The dimensionless pressure at any spatial point in the reservoir will be the sum of

the pressures due to the source and the image wells.

where:
-7-

Substituting the exponential integral gives: Figure 2.2 is a portion of the figure presented by

Stallman [19521.

A dimensionless pressure-time graph which is a portion of the figure presented by Stallman

[1952], is shown in Figure 2.2.

0
c
3
u)
u)

t 1
n
(I)

-
u)
0
C
.-
0
u)0*1
C
Q
.-
0
E

0.01
0.1 1 10 100 1000
Dimensionless Time I
Fig. 2.2 Log-log drawdown response for reservoir with linear boundary (after Stallman 1952).

The lower most curve is the line source curve for a r2/rl ratio of infinity, while the upper most

w e is for a pressure point at the boundary, r2/r1 = 1. Figure 2.2 is a portion of the figure

presented by Stallman [1952].


-8-

2.1.3. Linear Boundary Drawdown and Buildup (Eipper)

Buildup behavior is incorporated into the response by using superposition in time. Zero flow

rate is reflected mathematically by adding at the time of shut-in (tp) the press* response of

source and image wells which have equal and opposite flow rates to the doublet used during

drawdown:

Substituting the exponential integral gives:

+ E i [ z -4
A dimensionless pressure-time log-log graph is shown in Figure 2.3.

2.2. SEMI-LOG TYPE CURVES

2.2.1. Linear Boundary Drawdown (Sugeev et al.)

Semi-log type curves enable better matching of late time data as the higher pressure and late

time scales are expanded from the log-log scale. Sageev et al. [1985] noted the si@larity in the

shape of semi-log graphs of the Stallman type curves (Figure 2.4) for r2/rl ratios greater than

10. These curves and those for a constant pressure boundary were mathematickilly collapsed

giving a single semi-log type curve which could be used to determine the infereim ellipse for

any linear boundary. The semi-log curves were arbitrarily collapsed onto the cutve for rzlrl =

100. The new curves used modified dimensionless pressure-time parameters which are defined

as follows:
-9-

10 I I I

Q)
L
3
(I)
(b

t 1
L
(b
(I)
-
Q)

C
.-0 0.1
(I)
C
0)
.-
0
E

0.01
0.1 1 10 100 1Ob0
Dimensioniess Time
Fig. 2.3 Log-log drawdown and buildup response for reservoir with linear boundary (after Eipper
1985).

0)
L
3
cn 3
(b

Q)

L
c
(I)
cn 2
- 0,
C
0

L
1

0
0* I 1 10 100 1Ob0
Dimensionless Time

Fig. 2.4 Semi-log drawdown response for reservoir in a semi-infinite system.


- 10-

A dimensionless pressure-time semi-log graph is shown in Figure 2.5.

Q)
L
8
3
u)
u)
Q)
L
a. 6
u)
u)
-C
0,

.-0 4
u)
C
0,
.-
E
0
0 2
.-
@
c
.-
W
0 I I 1 I I1111 I I I I I1111 I I I I IIIII I 1 I II Ill1
t 0'
10 100 1000 10000 1e+O5
Modified Dimensionless Time
Fig. 2.5 Semi-log drawdown response for reservoir with linear boundary (after Sugeev 1985).

2.2.2. Linear Boundary Buildup (Fox)

Fox [19841 used the unpublished technique of Sugeev et ul. [1985] to produce a single type

curve to analyze buildup data. Mathematical development used the log approxbation to the

line source solution to simplify the exact solutions. The behavior for buildup data in a semi-

infinite system is characterized by two semi-log straight lines. Fox equated &e amount by

which each of the two straight lines was required to be shifted in pressure and Homer time at

the point of intersection of the two. The semi-log curves were arbitrarily collapsed onto the

curve for r2/rl = 100. This gave two equations from which the modified dimensilonless param-

eters were determined :


- 11 -

(2.10)

(2.11)

A dimensionless pressure-time graph is shown in Figure 2.6.

e *
3
u)
0
Q
L
& 6
(I)
u)
-QC
0
--
u)
4
C
Q
.-E
P
P Z
.-0 )
.-
c
0
P o
1 10 100 1000 101
Modified Horner Time
Fig. 2.6 Semi-log buildup response for reservoir with linear no-flow boundary (after Fox 1984).

23. INFERENCE ELLIPSE

Vela [1977] proved that in a homogeneous isotropic semi-infinite reservoir the lidear boundary '

is located tangent to an inference ellipse. The wells are the focii of the ellipse. Vela's equation

is redefined below using the x-y origin at the source well (Figure 2.7).

(2.12)

A schematic layout for notation used in linear boundary detection with the infemnce ellipse is

shown in Figure 2.7.


I 11

- 12-

I Image Wen

........
........

........
........
........
........
........
........
........ ..........
..........
..........
..........
..........
..........
..........
...............
...............
...............

Infer&a Ellipse

1.01325~10'~pDmaeh
kh= d-m
2x PSmatch

where:

1 d = 1.01325~10'~durcy
- 13-
3. BAROMETRIC AND EARTH TIDE EFFECTS

3.1. INTRODUCTION

The increased use of sensitive quartz crystal gauges for pressure measuremedts has meant

greater emphasis on filtering out earth tide and barometric effects to obtain a clean interference

pressure response. The response of reservoirs to barometric pressure, earth and oceanic tides

and rainfall has long been observed with the magnitude of the pressure response depending on

the fluid viscosity ,permeability, porosity and total compressibility.

Barometric variations can produce reservoir responses of more than 8 kPa while earth tide

responses in the reservoir are usually less than 1 kpa.

Studies of barometric and earth tide effects have been camed out since the beginning of the

century. Young [1913] examined the responses of fluid levels in bore holes with a tidal varia- ~

tion. LaRocque [1941] observed fluctuations of water levels in wells during pedods of strong

earthquake activity. Jucob [1944] correlated ground water levels with rainfall. Bredehoeji

E19671 and Bodvarsson [1970] studied the relation between rock characteristics &Id amplitude

of the pressure response in an open aquifer. Thorsreinsson and Eliusson [1900] correlated

pressures in the Laugarnes geothermal field in Iceland with oceanic tides. Khuraha [1976] and

Srrobel et al. [1976] showed that tidal phenomena also affect closed systems. dlrdirty [1978]

developed the following expression for the amplitude of the pressure response of a closed sys-

tem to earth tides. The equation has been adapted to petroleum notation:

Sensitivity of the response to the forcing frequency and total compressibility was also exam-

ined.
- 14-
33. BAROMETRIC PRESSURE CORRECTION

To calculate pressure changes in the reservoir from water level measurements aorrection must

be made for changes in barometric pressure. A force balance in the wellbore gives:

where:

and:

t\b = barometric efficiency

(3.4)

The xometric effic.acy is determined by compaug changes in reservoir ptessure due to

barometric fluctuations under steady state conditions with the barometric changps recorded at

the surface. The barometric efficiency used for these tests was 0.85 after Grunt [1980].

An important assumption in determining the reservoir pressure from water level measurements

is that the density (pw) of the fluid column in the wellbore remains constant. Regular tempera-
ture or pressure profiles are recommended during pressure monitoring in the observation well

to monitor drift which may lead to incorrect interpretation. Pressure measurements during

drawdown in Test B2 may have been affected by changing wellbore density ttlthough other

factors cannot be ruled out (refer 0 7.2.2).

33. EARTH TIDES (Sageev et al. 1986)


I
The effects of earth tides may be significant during the early time response of an observation
I

well. On a dimensionless match to the line source, there are three parameters that determine I

the importance of earth tides: the amplitude, the phase with respect to the start df the test, and 1 ~

the frequency. Figure 3.1 presents a record of the barometric pressure at the Ohaaki field dur- I

ing the the time that the flow test was carried out. The first ten days of this recofd are present-

ed in Figure 3.2. In the five day period between 24 and 144 hours a linear ovetall decline of
1 I /

- 15 -

96 I

0 500 loo0
TIME (Hours)

Fig. 3.1 Atmospheric pressure variations during the " C Tests.

101

IO0

99

98

97

96

95
0 24 48 72 96 120 144 168 192 216 240
TIHE [Hours)

Fig. 3.2 Barometric and earth tide fluctuations during the "C" Tests.
- 16-
the barometric pressure is evident. In addition, the daily oscillations are distinct, and are a

combination of two sine waves with different amplitudes.

The first sine wave has a daily cycle, and the second sine wave has two cycles @r day. This

second sine wave with the higher frequency is responsible for the small depressi4ns at the top

of each pressure cycle that occurs at noon time in each of the five days. Hence, there is one

linear overall pressure decline that represents the regional barometric trend, and tdvo oscillatory

pressure functions, superimposed on the linear pressure trend, caused by the sun &d the moon.

Figure 3.3 presents a simulated pressure response of these three pressure

-a

-4

-6

-8

-10

0 1 .2 3 4 5
TIME (Doysl
Fig. 3.3 Modelling of the earth tide response using two superposed sine waves on a laear barometric
trend (after Sugeev et uf. 1986).

functions. The thick oscillating curve that displays similar characteristics tu the Ohaaki

barometric record, is the sum of the three other curves on the figure. The amlplitude of the
- 17 -

high frequency wave is 0.7 the amplitude of the low frequency wave. Also, tliere is a phase

shift of 0 . h assigned to the high frequency wave.

Figure 3.4 presents a match of the simulated barometric pressure of

101

loo

CI

0 W
n
25
w 98
a
3
v1
v)
w
a 97
n

96

85
0 24 48 72 96 120 144 168 182 216 240
TIME (Hours1
Fig. 3.4 Match of the synthetic and real earth tide responses (after Sageev et al. 1986).

Figure 3.3 to the five day record presented in Figure 3.2. The match is quite gwd. The inten-

tion is to demonstrate that knowing the barometric pressure during the flow test,, and matching

it to simulated responses may help in the future to filter out their effects from feservoir pres-

sure data. Also, it may be possible after deconvolving the barometric pressure^ into its com-

ponents, to evaluate reservoir properties such as storativity and transmissivity. Figure 3.5

presents three hypothetical earth tide effects superimposed on the line source curVe. The curve

denoted by A has a dimensionless pressure (Po) amplitude of 0.05, a zero phase,, and a dimen-

sionless time (rD) frequency of 0.2. The large amplitude causes the oscillations @fthe pressure

to extend up to a dimensionless time of 1. The other two curves denoted with B'and C have a

dimensionless pressure amplitude of 0.01, and there are no oscillations presenb Both these

curves have a zero phase, but have different frequencies. Curve B has double th$ frequency of

curves A and C. With combinations of different phases and frequencies it is possbble to get os-
- 18-

10

w
m
3 1
v)
v)
W
a
0
v)
gJ 0.1
J
z
0
u
v)
z
W
I: 0.01
u
0

0.01 0.1 1 10 100


DIMENSIONLESS TIME
Fig. 3.5 Superposition of the Line Source with a single sine wave (after Sageev et al. 1986).

cillations at early time even with low amplitudes.

The first data of the log-log matches from Test C3 (Figure 3.6) indicate a dimer@ionlesspres-

sure amplitude of 0.01 and a dimensionless time frequency of 0.2. This is one bf the reasons

for deviation of early time data below a specific pressure level of about 40 kPf-s/m' (6 3.4)

from the line source. The linear barometric pressure decline was not considered since it had a
smaller effect than the earth tide effects. Some processing methods may call for' smoothing of

pressure data or clipping of anomalies where the pressure declines unexpectedlyl The magni-

tude of interference pressure changes in the early time flow period may be of the kame order of

magnitude as the earth tide and barometric pressure changes. Under such conditiobs, smoothing

or clipping may lead to loss of information.

3.4. PRESSURE SIGNIFICANCE LEVEL

The pressure level at which the data can be considered significant can be d e t e q n e d from the

magnitude of the earth tide response as the data have been corrected for baronjetric
I
pressure

variations only. Figures 3.3 and 3.4 show that the amplitude of the earth tide flubtuations dur-
10

W
a
3 1
tn
tn
W
0:
a
v)
v) 0.1
W
-I
z
0
c(

tn
2:
y 0.01
c)
0

0.001
0.01 0.1 1 10 loo ~

DIMENSIONLESS TIME
Fig. 3.6 Test C3:Log-log matches of data to SruZZman type curves. I I

ing the "C"series of tests was about 500 Pa. Fluctuations of more than about a~tenthof a log

cycle cause problems in obtaining a unique match to the line source solution.lTherefore the

magnitude of the measured pressure fluctuations is required to be approximately 5 times the

earth tide variation in order to be significant. That is the pressure change in the bservoir must
I
be a minimum of 3 kPu. Since test data is normalized against flow rate befole plotting the

4
specific pressure significance level for each test must be determined by dividin the minimum

absolute pressure change of 3 kPa by the flow rate. For most tests the flow ratd was about 70

Ns which means the specific pressure significance level is 40 kPu-s/m3.


- 20 -
4. COMPRESSIBILITY

Determination of the effective compressibility of both the reservoir fluid and thb formation is

au important factor is ensuring accurate data analysis. Several publications have jcontributed to

the more accurate determination of total compressibility for multiphase systems.

Perrine [1956] used the following expression for two phase total system compreskibility:

Ct = S d W+ sgcg+ Cf (4.1)

Martin [19591 developed a more exact expression for total compressibility:

R m e y [1964] used the correlations of Standing [1952] to obtain approximatiod to the partial

derivatives reducing the amount of work required to calculate compressibilities. ~Horne[19801

developed graphs showing the variation in water and steam compressibilities wi@ temperature
l

and pressure. These graphs are presented in Figures 4.1 and 4.2 . Grant andiSorey [1979]

used conservation of mass and energy to derive the total compressibility for bhase change
I

alone in a two phase geothermal system:

A useful approximation to this equation using the Cluusius-Clupeyron equation whs also given.

+ [ ( l + ) P ~ f + WWPWCW]
4.O1x1r2
cph = P-l.& (4.4)

Grant et al. [1982] give a modification to this equation for non-condensible gases.

-
1
= 9.97~10-' 4 pi$ + 2.13~10-~
pg p2F1 (4.5)
- 21 -

~ kPa"
WATER COMPRESSIBILITY c , 10'
2500 k h

moo LP.

1500 tp.

loo00 kP.
P500 kP.
15000 W a

180 200 220 240 260 280 300


TEYPERATilRE *C

Fig. 4.1 Variation in water compressibility with temperature and pressure (after Home la80).

STEAM COMPRESSXI3LITY c,x lo3 kPa"


2.0

1.5

1.o

0.5

2000 4000 6000 8000 10000 12000


PRESSURE 19.

Fig. 4.2 Variation in steam compressibility with pressure (after Home 1980).
- 22 -
5. EARLY TIME DATA MATCHING
(Sugeev er al. 1986)

5.1. THEORY
,
I

The dimensionless pressure responses of semi-infinite systems bounded by arb impermeable


linear boundary are presented in Figures 5.1 and 5.2.

0.01 0.1 1 10 100 lo00


DIMENSIONLESS TIME

Fig. 5.1 Log-log type curve for a semi-infinite system (after Stallman, 1952 ).

Figure 5.1 is similar to the log-log Stallman [19521 type curve. The ratio r2/rl vdies between 1

and 10. The lowermost curve in Figure 5.1 represents the Theis line source solution. The up-

permost curve represents a semi-infinite system where the observation point is idjacent to the
I

boundary, yielding an r2/r1of unity. In view of Equation 2.4, this curve is a sum~oftwo identi-

cal exponential integrals, and is shifted by a factor of 2 in the vertical directiol from the line
I

source curve. As the value of the ratio r2/r1 increases, the dimensionless pressure response

departs from the line source solution at later times, as the effects of the imperme/able boundary
I
become significant at the observation point. I
7
- 23 -

a*o

6.0
10.0

0.01 0.1 1 10 100 1-


DIMENSIONLESS TIME

Fig. 5.2 Semi-log type c w e for a semi-infinite system.

Figure 5.3 is a semi-log presentation of the same data as in Figure 5.1. The curves having r2/rl

ratios between 1 and 1.5 are closely spaced in both Figures 5.1 and 5.2, makih log-log and I
I
I
semi-log type curve matching difficult,

The similarity of the curves for small values of the distance ratio is presented in Figure 5.3. In

this figure, the line source curve is translated along the time and pressure axes~by factors of

1.075 and 1.96 respectively, to match the curve for an r2/rl of 1.1. The two cur& match well.

The square of the difference between these two curves is presented in the lower1 thin curve in

Figure 5.3. The two minima in the error curve represent the fact that the two

cross each other twice.


mlatched curves

The difference between the two curves is smaller than b e resolution


I

expected even from very sensitive pressure recording devices. Hence, transient pressure data
from wells near an impermeable linear boundary with an r2/r1 ratio of 1.1 can beimatched suc-

cessfully to the line source solution, Based on the translation of the line sourk curve, this

match would yield almost the correct storativity but a transmissivity which is la factor of 2

lower. No indication of the presence of an impermeable boundary would be prese t.


1I
- 24 -

la

W
[L
=) 0.1
u)
v)
W
E
n
v)
E 0.01
1
z
E!
v)
z
W
rO.001
c(

o.Ooo1

1.-OE
a I 0.1 1 10 100 lo00
DIMENSIONLESS TIME

Fig. 5.3 Log-log match of the Line Source to the curve for r2h1 = 1.1 (after Sugeev et 41. 1986).

A standard procedure in analyzing interference pressure data that match the line sburce solution 1
is to evaluate the minimum reservoir area that is free of any boundaries. In the cCse of a linear

boundary, this area is elliptical in shape, as described by Vela [1977]. This anilysis with in-

terference data coming from a configuration with a distance ratio less than aboht 2 will over

estimate the reservoir area free of boundaries. The magnitude of this over estidation depends
~

on the duration of the test. For example, if a pressure response with a distan? ratio of 1.1
I ,
matches the line source up to a dimensionless time of 20, the minimum distanb ratio deter-

mined from this match would be 10. This match would locate the linear bound* at a distance

ten times further than the actual boundary location.

The results from an interference test with a distance ratio smaller than 2 should be compared

with results from other responses in observation wells or the source well. In d e case of the

source well, the distance ratio is large and the distance ratio and transmissivity mby be estimat-

ed. The difference between the transmissivities of the source well and the observdtion well may
- 25 -
indicate the correct match of the interference data.

Although the curves for r2/rl ratios of 1.1 and 1.5 are very similar to the line sobrce, the curve

for an r2/r1 ratio of 2 is unique, and cannot be matched to the line source curire. Figure 5.4

presents the two curves without any translation This yields

w
a
2 0.
v)
W
K
a.
u)
E 0.01
-1
B
CI
v)
z
W
ro.wi
Y
0

o.Ooo1

10-05

DIMENSIONLESS TIME
Fig. 5.4 Log-log response of the Line Source and a semi-infinite system with r2/rl = 2.pI (after Sugeev
et al. 1986).

an early time match to the infinite acting line source response, but after a dim(nsion1ess time
I

of 1, the two curves are different. This can also be seen in the error curve. Thd error is small
at early time, and increases rapidly with time. Figure 5.5 presents a late time mz/tch of the line

source to the curve for an rdrl ratio of 2. Here, the early time line source +havior is not

matched, as described by the error curve. The minimum in the error curve is caused by the two

curve crossing one another. Hence, the line source curve cannot match simultanebusly the early

time and the late time responses of an observation well with a r2/rl ratio of 2, an4 the detection

of a linear boundary is possible.


1(

W
a
=I 0.1
v)
v)
W
a
n
v)
E 0.01
-1
z
E!
v)
z
W
YO.001
u
0

o.oO01

I .-os
0.01 0.1 1 10 la, lo00
DlHENSIONLESS TIME

Fig. 5.5 Log-log match of the Line Source and a semi-infinite system with r2/rl = 2.0 (#fer Sugeev et
al. 1986).
,
The log-log and the semi-log responses are summarized in Figure 5.6. The up@ family of

curves represents the log-log pressure responses, as presented in Figure 5.1. The lower family
of curves, translated by a factor of 5 for display purposes, represents the semitlog pressure
I

derivatives. From the pressure derivative curves, it is clear that for r2/rl ratios greabr than 5 the

semi-log type curve matching technique is applicable because there is a distihct transition

between the infinite and the semi-infinite behaviors. Pressure responses from cakes with r2/rl

ratios between 2 and 5 can be analyzed using the type curves presented her% since these ,
responses are significantly different from the line source. Also, pressure responses~with r2/r1 ra-

tios smaller than 2 can be matched to the line source curve, and may lead to erro$eous estima-
tions of some reservoir parameters.
- 27 -

10

W
(Y
3 1
v)
v)
W
nz
Q
tn
0.1
-l
z
0
c(

v)
z
W
r 0.01
Y
0

0.001
0.01 0.1 1 10 100 loo0 1c D

DIMENSIONLESS TIME
Fig. 5.6 Log-log type curve for dimensionless pressures and semi-log dimensionless +rivatives for a
semi-infinite system (after Sageev et al. 1986).

5.2. MISSING DATA


The detection of linear boundaries from transient pressure tests is based upon1 the transition

between the infinite acting and semi-infinite acting pressure responses. If this trbnsition is not

present little can be done to detect the linear boundary. Some tests may have thei required time

span for the detection of a linear boundary, but, for human or mechanical reasonp, some of the

pressure data are missing. Local sampling problems that arise from discretizing/ the time and

pressure domains are not considered. Local sampling is assumed adequate. Coqcern is with a

time period of missing pressure data that is significantly longer than the samplirig intervals re-

quired for transient pressure analysis.

Figure 5.7 presents an example of missing pressure data from Test C3 in the Obaaki geother- 1

mal field. Test C3 was started on May 1, 1984 and the drawdown portion of it hasted for 339

hours. The active well, BR20, produced at a constant rate of 84 Zls, and pressures were meas-

ured with a quartz crystal gauge at the observation well, BR34. The distance between the two

wells is 1145 meters (Table 7.16). The pressure scale in Figure 5.7 is normalizqd by the flow
I
1 10 100 loo0 10000 1.*05
TIME (Hours)

Fig. 5.7 Log-log plot of data From Test C3.

rate measured during the test. Pressures were not recorded from the 12th to thd 54th hour of

the test. The early time data of this test below a specific pressure level of la kPu-s/m3 are

affected by earth tides and baromemc pressure changes.

To determine the significance of the missing data in Test C3 consider a hypoqetical case of

the interference pressure response from a semi-infinite reservoir with a distance rbtio of 2. Fig-

ure 5.8 presents four cases of pressure responses, three of which have missing pressure data.

Figure 5.8A is the complete response that spans the transition flow period, and c@ be success-

fully analyzed for the detection of the linear boundary. In Figure 5.8B, a log cydle of the data

during the transition is missing. Yet, the early time infinite acting response andl the late time

semi-infinite response are well defined. The two portions of the pressure ddta cannot be

matched to the line source simultaneously, and the data can still be successfullj analyzed for

detecting the linear boundary.

The pressure response presented if Figure 5.8C has a longer period of missingldata than the

response of Figure 5.8B. Yet this pressure response can be successfully analyked since the

infinite and the semi-infinite responses are well defined. In the pressure respons$ presented in
I

Figure 5.8D, all the infinite acting pressure response is missing. The late timd semi-infinite
- 29 -

A B

w
[tl
3
cn
cn
w
[tl
a
v)
cn
w
-1
z
0
U
v)
z
w
x
Y
0

D
DIMENSIONLESS TIME
I
Fig. 5.8 Hypothetical cases for missing pressure data for rdrl = 2.0 (after Sugeev et al. 1986).

I ~

response can erroneously be matched to the line source curve, hence, yielding i b r r e c t reser- 1

voir parameters, and providing no indication of the presence of a linear boundary. ~

The first ten hours of the test from Ohaaki are influenced by earth tides and bdometric pres- ~

sures, making the condition of data in this test similar to the response presented in Figure

5.8D.Various matches of the drawdown data of this test are presented in FigureI5.9. The dis-
tance ratio varies between 1.2 and infinity (for the line source). The storativities lderived from

the log-log match vary by a factor of less than 2, but the transmissivities vary within an order

of magnitude. Without information from other flow test, the analysis of the prdssure data is

ambiguous. In this case, the pressure data are missing in the time period that is required to es-

tablish the infinite acting response of the reservoir, and the late time data repredent the semi-

infinite behavior. Hence, a unique analysis of this tests is not possible.


10

W
Q:
3 1
tn
v)
W
a
0
tn
v) 0.1
W
.
l
0
z
Y
v)
z
y 0.01
U
0

0.001
0.01 0.1 1 10 loo
DIMENSIONLESS TIME
Fig. 5.9 Four possible log-log matches of Test C3 data to the Stallman type curves.
- 31 -
6. HYDROLOGY OF THE OHAAKI GEOTHERMAL F 4 L D

6.1. BACKGROUND

The Ohaaki geothermal field is located in the Taupo volcanic zone, 30 km nod-east of Lake

Taupo, near the center of the Nonh Island of New Zealand (Figure 6.1).

Fig. 6.1 Locality map for the Ohaaki Geothermal Field.

The Waikato River, the largest in the North Island, bisects the field in a north-e+st south-west

direction. A total of 44 wells have currently been drilled (Figure 6.2). These rbnge in depth
- 32 -

Fig. 6.2 Well layout of the Ohaaki Geothermal Field.

from 358 meters to 2587 meters with an average of 1256 meters (Hedenquist 1985). The east
bank land area contains 13 productive wells and is known as Broadlands. The &st bank land i

area which contains 17 productive wells and includes the power station site, is known as

Ohaaki. The term Ohaaki is used in this report to designate the total geotherm4 system as a

consequence of the siting of the power station. I

Drilling began in the west bank of Ohaaki at well BR1 in 1965. The first well to $e successful-

ly discharged was BR2 in 1966. Between 1967 and 1971 large scale discharge besting of the

field was undertaken with 35 Mr of fluid being withdrawn from the reservoir (Hint and Hicks
- 33 -

1975). Reservoir pressures in the producing area dropped during this period by 1.2-1.8 MPa.

From 1971 to 1973 the field was at near shutdown as other energy projects tookldevelopmental

priority. Then, between 1974 and 1983, discharge testing continued with 19 Mrl of fluid being

produced of which 5 Mr was reinjected. Reservoir pressures during this period recovered by

600-800 WU (Grunt et al. 1982).

6.2. GENERAL
Hydrologically, the Ohaaki field consists of three reservoirs (Grunt et ul. 1983, McGuinness

1985). The main reservoir is located in the Rangitaiki ignimbrite and Rautawiri breccia forma-

tions below 500 meters depth (Figure 6.3). Geologically, the lower limit

Depth (m)
~~

0
I ~ Ohaaki ~ Waikato
I
l I

750

1500

NW SE
Fig. 6.3 North-west south-east geologic cross section showing postulated flows in the Oqaaki reservoir. 1

of the reservoir could be defined by the upper limit of the basement greywacke which extends 1

from 1000 meters depth in the east at BR7 to over 1700 meters depth at BR115 in the west

(Grindley and Browne 1968, Browne 1973). During discharge and subsequent rdcovexy of the

field between 1967 and 1979 noticeable changes in pressure and temperature bccurred only

above 1000 meters depth. This may represent the exploitable limit of the field, albough a tem-

perature saddle below 1000 meters between the east and west banks may be a
I
- 34 -
flow below that depth (Grunt 1983). Pressure gradients in the reservoir indicate that vertical

permeability is much lower than horizontal permeability. Vertical permeability wSs roughly es-

timated at about 3 millidarcies (Grunt 1983).

63. EASTBANK

On the east bank, upflows of hot fluid from depth may occur near wells BR36 and BR43 as

temperatures show a continued increase with depth in this area (Grunt 1983). Hi& temperature

chloride waters are saturated with silica and calcium carbonate with the average qarbon dioxide

concentration in the water at depths of 400 - 800 meters being 0.6% by weigljt. The partial
pressure of carbon dioxide in the reservoir varies from 0.8 - 3 MPa, but at d e p q below 1000

meters partial pressures are probably over 10 MPa. The high gas concentrationg appear to be

related to the basement greywackes and argillites. This is confirmed by higher gas concentra-

tions found in those wells which penetrate the basement rocks. (Ellisand Mahon, 1977). Most

rising hot fluid on the east bank probably travels west because only there, in the1 Ohaaki pool,

does any significant surface discharge of reservoir fluid occur. Steeper pressure gladients in the

region near BR24 and BR26 indicate a possible flow barrier to this area. Shallo+ ground wa-

ters in the Broadlands dacite above the postulated east bank upflow zone (Fidre 6.3) show

steam dilution effects but no significant addition of deep chloride water. m e estimated

minimum fluid temperature in the upflow zone is 310°C. Geochemical analysis iof fluid from

east bank production wells show evidence of both boiling and dilution with bidarbonate rich

condensates (Hedenquist 1983). Condensates can be formed by the deep 310°C fluid flowing

upward to the base of the sealing mudstone and siltstone at 500 meters depth. *ere the pres-

sure reduction would cause about 10% of the total mass to flash off into the oveislying Broad-

lands dacite (Grunt 1983) while the other 90% is assumed to travel to the west bink within the

hot permeable region of the field. Lumped parameter modelling of the east baqk (Grunt and

Iles 1982) showed the flow from the east to the west was about 3.5 kt/yr-Pa. qodelling also

showed that total recharge to the east bank was about 4 kt/yr-Pa while recharde to the west

bank not including that from the east bank was about 1.5 ktlyr-Po. Hence the kast bank has
~ ___ ~ _ _ _ ~ _ _
I - 1

- 35 -
better recharge than the west bank (Grunt 1983).

6.4. WEST BANK

On the west bank total heat discharge from 1965 to 1979 was 5x1Ol6J. Change$ in well tem-

perature profiles indicate that only 20% of the enthalpy of the discharged flui4 was not re-

placed (Grunt 1983). Pressure and temperature profiles measured in the central west bank show

no evidence of recharge of brine from directly below. Profiles measured in well$ BR2, 17, 18

have shown some signs of the existence of upflow in the past but these have 'not persisted.

Geochemical analysis of production fluid from the west bank indicated that boilihg with minor

dilution is taking place in the production zone (Hedenquist 1983).

Outside the field pressure gradients were 400 f f u below reservoir pressures befofle exploitation

ensuring that during early drawdown of the reservoir no significant recharge of +ld water was

likely (Grunt 1982). Detailed analysis of pressure measurements in nearly all wdst bank wells

(Figure 6.4) both before and after exploitation showed that the wells fell into thre/: groups each ~

with a characteristic linear formation pressure profile (Bixley 1982). Wells wid the smallest
~

pressure gradients were contained within the 270°C contour at 900 meters depth1 and this was

considered to define the productive area of the field (Figure 6.5). Wells included lwere BR5, 8,

I
11, 12, 13, 14, 19, 21, 24, 28. Wells offset more than 500 meters from this conto r showed the

greatest formation pressure gradients and were considered very poorly connectep to the deep

reservoir. Wells included in this group were BR6, 32, 33, 37, 39. Some wells showed forma-

tion pressure gradients which were intermediate between the previous two cases b d were con-

sidered to be poorly connected to the deep reservoir. Wells included were BRIO, 16, 29, 34,

38. Wells BR5, 12, 34 were exceptions to the rule showing fair connection to &e deep reser-

voir possibly as a consequence of north-east trending faults in this region.

Gravity changes indicate that total recharge to 1983 to both east and west banks ib about 25 Mt

or about half the total mass discharge. Gravity changes affected an area grebter than that

covered by the 5 ohm-m resistivity boundary and were consistent with a reservoid of 1.5 km ra-

dius and thickness of 500 meters (Hunt 1985). A persistent negative trend in grqvity measure-
- 36 -

400 -100
OITOlnrt pmduclivr

-
n
E
4
v!
-
n
E Pmductivr Ohrakl I
c

n 800
d
QI

-500
-t
5
a
w

1200 I I I -900
1 3 5 7 Q 11

Pressure (MPa)
Fig. 6.4 West Bank well pressure profiles (after Sixley 1982).

ments extending up to 2 km west of a line bounded by wells BR38 and BR39 map indicate that

the main reservoir is hydrologically connected to this part of the field (Figure 6.6). The fact

that gravity surveys show a continued loss of mass in this area could be due tu the buffering '
effect of a deep two phase zone and/or low permeability which has delayed respOnse to earlier

reservoir drawdown. Measurements in the east of the field between BR14 and dR32 reflect a

similar trend.

6.5. BROADLANDS DACITE (McGuinness 1985)

Overlying the Ohaaki field south of a line joining wells BR 3, 14, M6 is the seqond reservoir

unit, the Broadlands dacite. This formation contains liquid brine at about 140°C, has transmis-

sivity about 370 d-m and storativity of about 3.0xlO-'mlWa. A 35 day interference test in

January 1984 in which source well BR40 was produced at 44 kgls and pressureis observed at
- 37 -

Fig. 6.5 West Bank well grouping from pressure gradient characteristics (after SirZey 1982).

well BR5, both with permeable zones in the Broadlands dacite showed no observdble boundary

effects. The region of influence covered a radius of about 2.5 km around BR5 arid BR40. No

response in wells BR33, M10, M11, 910 with permeable zones in the Ohaaki rhholite implies

near isolation of this geological unit form the Broadlands dacite. The upper bound permeability

estimate in the barrier was about 10 md.

The Broadlands dacite shows no connection with the deep reservoir as discharge of well

BRM9 with permeable connection in the Broadlands dacite produced no observable response in

wells BR13, 23 having connection to the deep reservoir. Due to its hydrological ibolation from
- 38 -

Fig. 6.6 Areal view of the Ohaaki Field showing West Bank gravity anomaly.

the deep reservoir and apparent large extent the Broadlands dacite has been moo@d for reinjec-
I

tion of waste fluid from the power station.

6.6. OHAAKI RHYOLITE

Overlying the west bank is the third reservoir unit, the Ohaaki rhyolite. Interferebce tests were

carried out from November 1983 - May 1984 in which the response of wells BR12, 33, 37,

M7, M10, M11, 9/0 was monitored to the discharge of wells BRll, 22, 33, 40, M9 (Figure

6.7). Wells BR12, 37, M7 showed no response to production confirming their &or hydrologi- '
- 39 -

BRMl

Fig. 6.7 Areal view of the Ohaaki Field showing significant geology and Well Layout for West Bank
Rhyolite interference tests (after McGuinness 1985).

cal connection to the deep west bank reservoir. A response to drawdown in the Ohaaki rhyol-

ite was measured in shallow ground water well BR9/0 with slow pressure reicovery during

buildup. The linear drawdown pressure response of observation wells BR33, MlO, M11, 910

and the slow buildup, indicated that the Ohaaki rhyolite was an open tank reserkoir with poor

recharge (McGuinness 1985). Analysis of interference test data gave the porosiQ area (+A) of

the open tank as 0.23 k d assuming that the box is filled with liquid water at 140" C. Ground

water well BR9/0 is known to have the same response as the Ohaaki Pool, the single major
-40-

chloride spring in the Ohaaki field. However ground water wells further than abo+t 500 meters

radius from BR9/0 on the west bank have shown little response and east bank sldallow ground

water wells show no response to discharge of the deep reservoir.

Between 1970 and 1974 pressure rises were measured at the base of the Ohaaki rhyolite at 450

- 540 meters depth indicating the existence of a two phase zone in this region (Grunt 1983).
This zone disappeared after 1974. Pressure changes at the surface in the Ohaaki @ool were 5 -
10 times less than those occurring at the base of the Ohaaki rhyolite. The Ohaaki rhyolite and

the deep west bank reservoir were assumed to be separated by an aquitard codisting of low

permeability siltstones and Huka mudstones (Grunt1983). Pressures in the Ohaa# rhyolite ap-

pear to be buffered either by a lateral connection to another near surface reservoit or by a two
phase zone within the rhyolite. Lack of evidence of large lateral flows in the rhydlite make the

two phase zone option more credible.

6.7. COMMENT
At present, a 110 MW(e) power station is under construction on the west bank a( Ohaaki with

generation expected to commence 1989. Despite over 20 years of researchi during field

development, enigmas remain. Two phase zones and isolated areas of low permgability in the

main hydrological units complicate interference test interpretation. The role of fa9lting in con-

necting and distributing colder recharge fluid to the deep reservoir has not been fully esta-

blished. The gravity anomaly west of the 5 o h - m resistivity boundary on the des, bank has
I

not been explained by independent measurements. The lower limit of the deep rejervoir is ten-

tatively assumed as 1000 meters however evidence of pressure drawdown in BRb8 extends to

nearly 1200 meters and cores from BR34 show hydrothermal mineral depositio$ in low per-

meability veins in the basement greywacke at 2590 meters (Brairhwuite 1979). “$e extent and

type of connections of the Ohaaki rhyolite and Broadlands dacite reservoir unid with forma-
tions outside the primary production area have not been established. The answers to these enig-
mas may have to await the next decade when large scale production will inevibbly provide

enough data to obtain solutions to the complex problems involved.


-41 -
7. ANALYSIS OF INTERFERENCE TESTS DATA IN THE OHAkKI WEST BAN

7.1. MEASUREMENT ACCURACY I

Two types of instruments were used to measure pressures in the interference test$. These were:

(i) Water level chart recording devices.

(ii) Quartz crystal gauges.

The water level chart recorder has an accuracy of about 100 Pa. and was used iin the observa-

tion wells for the "B" series of tests. Barometric pressures for this series of tedts were meas-

ured using a barograph which also has an accuracy of about 100 Pa.

Quartz crystal gauges were used in the "C*series of tests. The accuracy of the quartz crystal
gauges varies depending on the pressure range, but is typically lo-' of the maxiqum range. At-

mospheric pressures were monitored with an accuracy of 1 Pa, while reservoir prbssures had an

accuracy of 100 Pa.


1 I1
I

- 42 -

7.2. WELL LOCATION

The locations of observation and source wells for the interference tests are sh@wnin Figure

7.1.

Fig. 7.1 Source and observation well location for interference tests .
- 43 -

7 3 . DATA FROM WATER LEVEL CHART RECORDERS

73.1. TEST B1: BR13 RESPONSE TO BR23 (Leaver et al. 1985)

Specific pressure verses time data for this test are presented in Table B1 and plQtted in carte-

sian form in Figure 7.2. Test specifications are presented in Table 7.1.

600
c)

E
3
0
\
(I)
I 400
0
0
y.
Y

u)
u)
Q,
L 200
[L

0
0)
P
v)

I
0 25 50 75 100 125 150
Time (hr)
Fig. 7.2 Cartesian graph of the interference data (after haver et al. 1985).

For a minimum pressure significance level of 3 kPu (refer $3.4)and an average ebrly time flow

rate of 84 Ns the minimum specific pressure significance level is 60 kPu-s/m3.

73.1.1. Log-log Analysis I

A log-log pressure-time graph of the data is drawn and matched on the Eipper /dimensionless
,
pressure-time graph of the same scale (Figure 7.3). Dimensionless pressure and !dimensionless '
time are defined as:
-44-

TABLE 7.1
TEST B1: BR13 RESPONSE TO BR23 DISCHARGE

TEST SPECIFICATIONS

Observation Well Source Well


1
Well No. BR13 BR23
I

Coordinates S 79705.90fr W 45687.30 fi S 78923.12fr W 45215127fr


(Broadlands Grid)

R.L. (C.H.F.) 293.5 m 291.8 m


(Motwiki Datum)

Permeable Depth 915 f 30 m 1015 - 1055


Drilled Depth 1081 m 1097 m
Open Hole D i m 197 mm 194 mm
1I
Interwell Distance 279 m
Discharge Rate 52.4 11s
Discharge Temp. 245°C
Discharge Enthalpy 1050 kllkg

Recording Meter Water Level Chart Recorder I

Flow Start Time 1230 Dec. 12 1979


I
~

Drawdown Period 49.5 h


Total Test Time 142.5 h

The graph contains three curves. Curve 1 represents the line source solution, while curves 2 1

and 3 represent type curves for drawdown and buildup which include the presenc): of a no-flow

boundary (Stallman 1952) with r2/tl=9.25.The early time data is matched to $e line source

solution giving a match point of: I


~
- 45 -

10

s
0
$ 1
h
0

-:
.- 0.1
C
e
.-
0
E

I I I
0.01 Y

0.1 1 10 100 loo0


Dimensionless Time
Fig. 7.3 Log-log match of data to Eipper type curves (after Leaver et al. 1985).

t = lOhr -tD-- 6.7


4

Values of transmissivity and storativity derived from the match point are presented in Table

8.1.
,
In conventional log-log analysis the late time data is matched to the appropriate Ftullrnun type
I

curve and the distance to the no-flow boundary or the ratio rzlrl determined. Froni Figure 7.3 it

can be seen that there is insufficient late time data to make an accurate match to one of the

Stallman curves, yet the early time infinite acting response matches reasonably wtll to the line

source solution

The buildup portion of the test initially follows the log-log buildup curve up to about 70 hours.

Then, the pressure builds up faster than expected for a system with one no-flow linear boun-

dary. This rapid buildup indicates the presence of pressure support in the system. This could
I

- 46 -
be caused by:

(i) Recharge to the system.

(ii) A two phase region in the reservoir.

(iii) Injection into the reservoir at the time of testing.

(iv) Instrument error.

The effect of this pressure support is not apparent in the drawdown portion of thei test, suggest-

ing that the distance between the source well and the pressure support is greated than the dis-

tance between the source well and the no-flow boundary. This boundary effect d+es not appear

in Test C2 and Test C1 shows a possible further barrier effect at late time. Further considera-

tion of this trend is therefore not considered warranted.

73.1.2. Semi-log Analysis

Conventional analysis involves matching the data to a semi-log Stallman plot which has the

effect of expanding the log-log graph in the higher pressure and time ranges. Fi+e 7.4 shows

that definition of the rz/rl ratio is still difficult.

Better definition can be obtained using the type of Sageev et al. [1985]. a s analysis in-

volves using the semi-log type curve formed by mathematically collapsing the s e ~ - l o gplots of

Stallman's type curves. The method gives accurate boundary definition providingj the r2/r1ratio

is greater than 10 and definition to 20% of the true ratio for values greater thaq 5. The graph

uses modified dimensionless coordinate axes for pressure and time (Figure 715). Modified

pressure and time are defined respectively as:

Figure 7.5 presents a semi-log match of the dimensionless pressure data to +e type curve
1 I I

- 47 -

Oimeneionioee T i n e
Fig. 7.4 Semi-log match of data to semi-infinite type curves (after &aver et al. 1985).

10

8 c
3
00
m
0
L
6 a

4 .-
OUFV. no.3

0
10 100 1000 10000
flodified Dimensionless Time
Fig. 7.5 Semi-log match of drawdown data to Sugeev et al. type curves (after Leaver el al. 1985).
- 48 -
presented by Sageev et al. [ 19851. The pressure match point is:

Substituting the semi-log pressure match point into Equation 7.3 and solving for the ratio rz/rl

yields:

rdrl = exp (O.Cb2.37 + h99) = 9.25 (7.6)

From the semi-log graph it can be noted that in order to get a unique match, only one semi-log

straight line and the transition period are required. The data show that during *e drawdown

period the first semi-log straight line and the transition developed, but the secbnd semi-log

straight line did not.

73.13. Horner Analysis

Figure 7.6 presents a conventional dimensionless pressure Horner plot of the builbup data with

the buildup curve for a linear boundary at rzhl = 9.25 shown on the same graph.

An alternative Horner plot can also be used which enables the boundary distance to be deter-

mined from a single buildup type curve. This analysis involves using the semi-lbg type curve

for long producing times (Ramey et al. 1973) formed by mathematically collapsihg the family

of Horner curves describing buildup behavior affected by a linear no-flow bbundary (Fox

1984). This graph uses modified dimensionless coordinate axes for pressure and Horner time.

The modified pressure and modified Homer time are defined respectively as:
I I
I

- 49 -

1 10 100
Hocner fine
Fig. 7.6 Match of buildup data on a Horner plot (after Leuver et al. 1985).

and are arbitrarily defined for do = 100 (do = -)


r3 and tpo = lo6.
rw

Type curve matching the data requires careful interpretation. A simple best fit /of the data on

the type curve would give a match point of (Figure 7.7):

from Figure 7.3:

tpo = (49.5)(0.67) = 33.17


from Equation 7.8:

solving for r2/rl gives:

-
e(3.4-o.o)

100 I” = 6.3 (7.10)


I I

1 10 100 1000
Modified Horner Time
Fig. 7.7 Incorrect Semi-log match of buildup data to Fox type curve (after Leuver et al. 1985).

1
This analysis is incorrect. The log-log analysis (Figure 7.3) shows that( the data aft
io = 43 (t = 64hr) are influenced by the pressure support and that only the buildpp data recorde

before this time can be expected to match the type curve. Inspection of the se#u-log plot (Fi
ure 7.4) shows that only the first semi-log straight line is present and that d e transition an

second semi-log straight lines are masked by the effect of the pressure support, Analysis of tht
i
I

drawdown data gives rz/rl = 9.25. The correct match point for the buildup ddta can be deteli-
I
mined by substituting this value in Equation 6. This gives: I

The correct match can now be made (Figure 7.8).

For this test the semi-log and log-log graphs alone confirm the location of no-flow bow-

dary. The distance between the source well and the observation well, rl, is 2719 meters. Hence
,
the distance between the observation well and the image well, r2, is:
I I

- 51 -

81 .

/ Horner Time

4 . .
Rorruro l o l c h point0
0 pi 4.2

1 10 100 IO00 10000 1e+&


Modified Horner T i m e
Fig. 7.8 Correct Semi-log match of buildup data to Fox type curve (after Leaver et a/. 1985).

r2 = (279)(9.25) = 2582 meters (7.121

The corresponding inference ellipse to which the postulated linear boundab is tangent i

presented in Figure 7.9.


s
I ,

- 52 -

Fig. 7.9 Inference ellipse location (after Leaver et al. 1985).


- 53 -
73.2. TEST B2: BR23 RESPONSE TO BR13 DISCHARGE

Specific pressure verses time data for this test are presented in Table B2 and plotted in carte-

sian form in Figure 7.10. Test specifications are shown in Table 7.2 and

1500 '
n
E
3
0 0
\
cb 0.0
1000
a
L,
Y
Y

0 100 200 300 4 0


Time (hr)
Fig. 7.10 Cartesian graph of the interference data

variations in flow rate are detailed in Table 7.3. The flow rate variation for th& test period u

to 190 hours fits the following correlation:


I
I

= 0.0815 - 0.00361 log -


l3hJ (7.13)

,
For a minimum pressure significance level of 3 kPa (refer $3.4) and an average( early time f l o i

rate of 81 Us the minimum specific pressure significance level is 37 kPu-s/m3.

73.2.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched on the Stullmtkn dimensionleds

pressure-time graph (Figure 7.11). On this graph, the lower curve is the Theis line sour+
I

solution while the upper curve is that for a no-flow boundary with r2/rl= 4.5. Early time dab
TABLE 7.2
TEST B2: B E 3 RESPONSE TO BR13 DISCHARGE

TEST SPECIFICATIONS

Observation Well

Well No. BR23 BR13

Coordinates S 78923.128 W 45215.278 S 79705.90fr W 456B7.30ft


(Broadlands Grid)

R.L. (C.H.F.) 291.8 m 293.5 m


(Moturiki Datum)

Permeable Depth 1015 - 1055 915 k 30 m

Drilled Depth 1097 m 1081 m

Open Hole D i m 194 mm 197 mm

Interwell Distance 279 m

Discharge Rate refer Table 7.3

Discharge Temp. 240°C

Discharge Enthalpy 1040 Wlkg


~~ ~

Recording Meter Water Level Chart Recorder

Flow Start Time 1400 Mar. 19 1980

Drawdown Period 382 h

Total Test Time 382 h

fit well to the line source solution for times up to 40 hours and yield a match h i n t of

ps = 10 kPa-slm3 p~ = 0.038

t=10h -tD-
6 - l.l0
- 55 -

TABLE 7.3
TEST B 2 BR23 RESPONSE TO BR13 DISCHARGE

DRAWDOWN FLOW RATE DATA

Time How Rate


(h) (Us)
0.0 81.8
2.0 82.1
23.O 76.3
46.0 74.9
190.0 73.4
214.0 72.7
238.0 72.7
262.0 70.9
331.0 70.9
334.0 63.3
382.0 59.0

After 40 hours the data follow the pressure response for a boundary of r$rl 9 4.5. After 120

hours the data show an upward trend above the response for this boundary. h e trend in the

data, if real, would indicate a further no-flow boundary. However it does not appear in any

(i) Temperature variations in the wellbore.


I
I
(ii) Discharge of another well in the field for which no record was kept.
~

(iii) Malfunction of the recording equipment.


I
I

73.22. Semi-log Analysis

A semi-log graph of the data is shown in Figure 7.12.The lower curve is the Itheis line sour4I

solution while the upper curve is the solution for the linear no-flow boundary kit rdrl = 4.5. 14

the semi-log plot the definition of the data match with the theoretical solutions lis improved foi

higher values of pressure and time. The graph shows that the match of th+ data with thq
- 56 -

10

5
k
-3 0.1

.-li
.-E
0

0.01 I
I . 10k Io. 1.10

0.001 I ...I I I b

0.01 0.1 1 10 100


Dimmioniemm Tima

Fig. 7.11 Log-log match of data to Stallman type curves.

chosen boundary position is good for I I 120 hours.

While the semi-log graph confirms the position of the postulated no-flow bound*, the results

from this analysis for transmissivity and the r2hl ratio are considered anor/ialous as the

matched ratio was not in agreement with the value of 9.3 obtained from TestsiB1 and B10.

The results of this test are presumed to be affected by by one or more of the fackors discussed

in the previous section. The storativity value which is relatively insensitive to dhanges in the

matched pressure level is retained for analysis.


0.01 0.1 1 10 100
Dimencioniees Time
Fig. 7.12 Semi-log match of data to semi-infinite type curves.
- 58 -
73.3. TEST B3: BR23 RESPONSE TO BR19 DISCHARGE

Specific pressure verses time data for this test are presented in Table B3 and plotted in carte-

sian form in Figure 7.13. Test specifications are shown in Table 7.4.

500
-E
; 400
0
\
(Ib

4
Y
300
e

0 50 100 150 20a


Time (hr)
Fig. 7.13 Cartesian graph of the interference data.

For a minimum pressure significance level of 3 f f u (refer §3.4) and an average etrly time flow

rate of 83 Ns the minimum specific pressure significance level is 36 kPa-s/m3.

7.3.3.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched on the Stallman dimensionless

pressure-time graph (Figure 7.14). On this graph, the lower curve is the Theis line source

solution while the upper curve is that for a no-flow boundary with r2/rl= 2.8. E&rly time data

show some deviation from the line source solution for t I 8 hours. Data deviate from the line

source solution for specific pressures of up to 60 ffu-s/m3 which is above the minimum

significance level of 36 kPu-s/m3 and therefore the deviation cannot be attributable to earth tide

effects. One explanation is that at early times the large pressure drop at the wellb4re may cause
- 59 -

TABLE 7.4
TEST B3: BR23 RESPONSE TO BR19 DISCHARGE

TEST SPECIFICATIONS

Observation Well Source Well

Well No. BR23 BR19

Coordinates S 78923.12fr W 45215.27fr S 78224.54 fr W 46154.80 fi


(Broadlands Grid)

R.L. (C.H.F.) 291.8 m 294.4 m


(Motuiki Datum)

Permeable Depth 1015-1055 600 - 1074

Drilled Depth 1097 m 1074 m

Open Hole D i m 194 mm 194 mm

Interwell Distance 356 m

Discharge Rate 83.2 11s

Discharge Temp. -
Discharge Enthalpy 1180 W/kg

Recording Meter Water Level Chart Recorder

Flow Start Time 1500 April 5 1980

Drawdown Period 189 h

Total Test Time 189 h

some non-Darcy flow resulting in a small two phase zone near the wellbore whikh would pro-

vide pressure support at early times only. Unfortunately no pressure data is available at the

source well which could substantiate this effect. Whatever the reason the early time data is im-

portant in obtaining the match point to the line source solution. Linear bolundaries with

r21rl 5 10 can only be accurately located from drawdown data once the infinite acting line

source behavior has been observed. In Figure 7.14 the match to the line sour& only occurs
-60-

10

s
(0

(0
1

t
L
m
-
f
C
0.1

-
0
0
E-
E

O 0.01

0.001
0.01 . 0.1 1 10
Dimensionless Time
Fig. 7.14 Log-log match of data to Stallman type curves.

between times of 8 and 27 hours or half a log cycle.

Matching data between times of 8 and 27 hours yields a match point of:

ps= 10 kPa-slm3 pD = 0.072

f = 10 h -
tD

6
= 0.735

Values of permeability and storativity derived from the match point are pr sen :d in Table 8.1.

After 27 hours the data show a good match to the pressure response for a boundary of

r21r1= 2.8.
I 11
,

- 61 -
73.3.2. Semi-log Analysis

A semi-log graph of the data is shown in Figure 7.15. The lower curve is

0.01 0.1 1 10
DimenoionIer6 T i n e

Fig. 7.15 Semi-log match of data to semi-infinite type curves.

the Theis line source solution while the upper curve is the solution for the linear no-flow boun-

dary at r2h1 = 2.8. In the semi-log plot the definition of the data match with the theoretical

solutions is improved for higher values of pressure and time. The graph shows khat the match '

of the late time data with the chosen boundary position is good. I

The semi-log and log-log graphs confirm the location of the no-flow boundary. The distance

between the source well and the observation well, rl, is 357 meters. Hen* the distance

between the observation well and the image well, r2, is:

r2 = (357)(2.8) = loo0 meters

The corresponding inference ellipse to which the postulated linear no-flow boundlary is tangent

is presented in Figure 7.16.

There is an indication of either tidal or residual barometric effects in the data from the faint os-
- 62 -
cillations on the fit to the theoretical curves.

Fig. 7.16 Inference ellipse location.


~ ~ ~~ ~~ ~~ ~~

- 63 -
73.4. TEST B4: BR23 RESPONSE TO BR13 INJECTION

Specific pressure verses time data for this test are presented in Table B4 and plotted in carte-

sian form in Figure 7.17. Test specifications are shown in Table 7.5.

1000
e

800

600 e

400

200

0
I
0 100 EO0 300 400 500
Time (hr)
Fig. 7.17 Cartesian graph of the interference data.

For a minimum pressure significance level of 3 kPa (refer 33.4) and an average early time flow

rate of 53 US the minimum specific pressure significance level is 57 kpu-s/m3.

73.4.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched on the SraZZman ldimensionless

pressure-time graph (Figure 7.18). On this graph, the lower curve is the Thds line source

solution while the upper curve is that for a no-flow boundary with rzlrl = 3.0. The data give a

fair match to the line source for t I 10 hours and yield a match point of

ps = 10 kPa-s/m3 p~ = 0.0375

tD
t=10h - =1.5
d
-64-

TABLE 7.5
TEST B4: BR23 RESPONSE TO BR13 INJECTION

TEST SPECIFICATIONS

Observation Well Source Well


~ ~

Well No. BR23 BR13

Coordinates S 78923.12fr W 45215.27p S 79705.90fr W 45681.30p


(Broadlands Grid)

R.L. (C.HF.) 291.8 m 293.5 m


(Moturiki Datum)

Permeable Depth 1015 - 1055 915 f 30 m

Drilled Depth 1097 m 1081 m

Open Hole D i m 194 mm 197 mm

Intenwell Distance 279 m

Injection Rate 53.1 Us

Injection Temp. 162°C

Injection Enthalpy 690 Wlkg

Recording Meter Water Level Chart Recorder

Flow Start Time 1400May 13 1980

Drawdown Period 402 h

Total Test Time 402 h

Between 10 and 200 hours the data show a good match to the pressure resporlse for a boun-

dary of rz/rl = 3.0. After 200 hours the data assume a gradient very similar to that of the line

source solution.

The match point of the early time data to the infinite acting curve agrees closely with that ob-

tained from Tests B2 and B10. However the appearance of the boundary effect is known to ac-
- 65 -

10

1
0
5
a 100
a
h'a
-3 0.1

.-
E
.-
0
10
10
Tiw (hr)
im

0.01 wtch point@


c. I O Po' 0.-
1 IO k 1
0
= 1.6

0.001 I I *
0.01 0.1 1 10 100
Oinenrimlesr Time
Fig. 7.18 Log-log match of data to Stallman type curves.

tually occur at rz/rl - 9.3 as shown in Tests B 1 and B 10. The significant deviatibn of this data
from the known solution must be due to injection of brine. The difference in~dimensionless

pressure at the end of the test between the curve for rz/q = 3.0 and the line so+ce is A ~ D= 1
or specific pressure difference ApS= 267 kPa-slm3. This pressure difference is about 29% of the

total pressure rise in the reservoir.

The difference may be due to a composite effect of the 162OC injected fluid nloving into the

reservoir predominantly at 270°C. Calculation of the radius of the bank of injected fluid shows

that the amount of fluid injected represents a very small portion of the total flbid within the

sphere of influence of the test:

Q = qt = 7t&h (7.14)
- 66 -

(7.15)

The permeable depth of BR13 is at 915 meters. Colder fluid injected at this degth could con-

servatively be expected to occupy permeable thickness in the reservoir at or belbw this depth.

The total permeable thickness in the Ohaaki reservoir is unknown but a maximpm permeable

depth of 1200 meters will be assumed. The effective permeable thickness seen by the injected

fluid is then:

h = 1200 - 915 = 285 meters (7.16)

A porosity (q) of 20% is arbitrarily assumed.

ri = 4- (7.17)

ri = 21 meters (7.18)

The distance between the source and observation wells (rl) is 279 meters. The injected fluid

front would not appear to have wide reaching thermal effects on the reservoir.

Fluid at 162°C heated in the reservoir to 270°C would undergo a specific volumd change from

0.1104 m3/kg to 0.1302 m3/kg or a change of 18% which would cause pressure changes of a
similar amount.

Other factors affecting the response could include:

(i) Mineral deposition in the formation (Henley and Harper, 1979).

(ii) Permeability changes due to hydro-fracturing at the source well (Bixley and Grant, 1980).

(iii) Permeability changes due to thermal effects of the injected fluid on the forniation.

(iv) Storativity changes due to condensation of any two phase fluid in the regidn of influence

of the injected fluid.

The data appear to show oscillations about the theoretical curves which may be the result of

earth tides or residual barometric fluctuations.


- 67 -
73.4.2. Semi-log Analysis

A semi-log graph of the data is shown in Figure 7.19. The lower curve is

E
.-
0 1 -

0
0.01 0.1 1 10 100
Dimensionlees T i n e
Fig. 7.19 Semi-log match of data to semi-infinite type curves.

the Theis line source solution while the upper curve is the solution for the linear 'no-flow boun-

dary at r2/r1 = 3.0. In the semi-log plot the definition of the data match with the theoretical

solutions is improved for higher values of pressure and time. The graph shows that the match

of the data with the chosen boundary position is good for r I200 hours. Aftet, this time the

data show a flattening off. Reasons for this trend were discussed in the previous section how-

ever if the flattening of were "real" it would be indicative of contact with p&ssure support

somewhere in the reservoir. This could be due to:

(i) A two phase zone.

(ii) Interaction with a "leak" within the reservoir.

(iii) Unknown discharge of another well in the reservoir.

(iv) Instrument error.


- 68 -
The apparent anomalous location of the no-flow boundary when compared to Tests B1 and

B10 performed on the same doublet preclude further use of the results of this test. It is difficult
to rationalize why the results appear to be so affected by the injection of a "small" amount of
separated brine. Further study of this test using a thermal simulator may help resolve this prob-

lem.
- 69 -
73.5. TEST B6: BR23 RESPONSE TO BR20 DISCHARGE

Specific pressure verses time data for this test are presented in Table B5 and plotted in carte-

sian form in Figure 7.20. Test specifications are shown in Table 7.6.

30 I

0
Y
Y

15 -
(R 0
(R
Q,
c -
n 10
0
0
0,
a 5 -
v)

0 e
0e% 5 10 15 20
Time ( h r )
Fig. 7.20 Cartesian graph of the interference data.

For a minimum pressure significance level of 3 kPu (refer 93.4) and an average aarly time flow

rate of 127 Us the minimum specific pressure significance level is 24 kPu-s/m3.

73.5.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched on the Eipper dimensionless

pressure-time graph (Figure 7.21). A feature of this graph is the short test dwation and the

corresponding small pressure drawdown induced in the reservoir. The maximum specific pres-

sure drop was 21.2 kPu-s/m3 (Table B5) which for a flow rate of 127 11s (Table 7.6)

corresponds to a pressure drop of only 2.7 kPu. This is less than the minimunl specific pres-

sure significance level of 3 P a . Despite this the data have been crudely clipped of secondary

effects and analysis proceeded with to demonstrate aspects particular to tests with small pres-
- 70 -

TABLE 7.6
TEST B6: BR23 RESPONSE TO BR20 DISCHARGE

TEST SPECIFICATIONS

Observation Well Source Well

Well No. BR23 BR20

Coordinates S 78923.128 W 45251.278 S 77671.308 W 4471&.90fr


(Broadlands Grid)

R.L. (C.H.F.) 291.8 m 291.7 m


(Moturiki Datum)

Permeable Depth 1015 - 1055 m 815 m, 945 -1085 m

Drilled Depth 1097 m 997 m

Open Hole Dim. 194 mm 194 mm

Intmell Distance 411 m

Discharge Rate 127 Us

Discharge Temp. 270°C

Discharge Enthalpy

Recording Meter Water Level Chart Recorder

Flow Start Time 1430 Oct.20 1980

Drawdown Period 6h

Total Test Time 21.5 h

sure drawdown.

The drawdown at early times is represented by the Theis line source solution 4s at early time

fluid is drawn from close proximity to the wellbore and the existence of a boubdaq has little

effect on the pressure response at the observation well. The magnitude af the pressure

response increases with time.


-71-

0
5
i 0.1
h
0

-
0

.-

\
p.01
.-
0

1 10 100

1 I
0.001
0.01 0.1 1 10 100
Dineneiunleeo Time
Fig. 7.21 Log-log match of data to Eipper type curves for r21rl 4.0.

(7.19)

(7.20)

High porosity, viscosity and compressibility and low permeability contribute to kxtending the

time it takes for a given pressure change to occur in the reservoir and therefore these attributes

help maintain pressures in the system.

Figure 7.21 shows that the effect of a boundary located at r2/tl= 4.0 is felt at the observation

well after tD&= 1. The lower curve represents the Theis line source solution wW1e the upper

curve is the pressure response for the no-flow boundary at r2/rl= 4.0. The data give a fair

match to the theoretical curves and yield a match point of:

ps = 1 kPa-s/m3 pD = 0.006

t= 10h -
tD
tD
=0.90
- 72 -
A feature of the short duration test is that the shape of the theoretical curves demonstrate how

reservoir pressures continue to fall for a short time after the source well has stopped discharg-

ing. This feature is present for all the Eipper curves but is usually masked by thb log-log scal-

ing of the data.

While the match shown in Figure 7.21 appears to be satisfactory a non-uniqu&ess problems

exist with the data. Only three points are shown on the drawdown with which tb match to the

line source solution. The determination of the existence of a boundary is largely reliant on the

match of the buildup data. Figure 7.21 shows that a good match can be obtained for an r2/rl ra-

tio of 4.0. However an alternative match using the line source solution alone ils presented in

Figure 7.22. Again the match looks plausible.

f::0.1
h
*
-3
.-5
p.01
.-
0

0.001
0.01 0.1 I 10 100
Dinensionleso fine

Fig. 7.22 Log-log match of data to line source solution.

Without additional data this problem could not be resolved. In this case the presdure match ob-

tained from Test C2 on the same doublet could be used and the buildup data thbn fitted to the

appropriate boundary curve (Figure 7.21). With the drawdown data matched at the same pres-

sure level as Test C2 the buildup data match fairly to the pressure response for la boundary of
- 73 -
rdrl = 4.0.

73.53. Semi-log Analysis

Semi-log plots of the data are shown in Figures 7.23.7.24.

0.01 0.1 1 10 100


D i m e m i m i n e Tiw

Fig. 7.23 Semi-log match of data to semi-infinite type c w e s for rdr1 4.0.

The figures show the low pressure response level compared to other tests as dl the semi-log
dimensionless pressure-time plots are to the same scale. Both semi-log plots show a good

match to the late time data emphasising the problem of non-uniqueness.

73.53. Horner Analysis

The Horner plot (Figure 7.25) demonstrates the drawback of obtaining a sigplificant match

when data span less than one-fifth of a Horner time log cycle. The data show a p o d match to

the theoretical buildup curve for r2h1 = 4.0 but would also an equally good makch to the line

source solution.
1, 1 1

- 74 -

0.75

0.5

0.S

0
0.01 0.1 1 10 100
Dimenoionieee Tine

Fig. 7.24 Semi-log match of data to line source solution.

1 10
Horner Time

Fig. 7.25 Match of buildup data on a Horner graph.


- 75 -

73.6. TEST B7: BR34 RESPONSE TO BR31 DISCHARGE

Specific pressure verses time data for this test are presented in Table B6 and plotted in carte-

sian form in Figure 7.26.Test specifications are shown in Table 7.7 and

300 I
n
E
3
250 -
e
e
ui e
I
0
200
a e
X e
Y
-
150 e
cn
(b e
e
ili e
c 100 - e
e
a 0
4 e
0 -
e
50 e

-
Q,
a e
v,
red: I I I I w
0
0 50 100 150 200 250 300
Time ( h r )
Fig. 7.26 Cartesian graph of the interference data.

variations in flow rate are detailed in Table 7.8.

For a minimum pressure significance level of 3 &Pa (refer §3.4)and an average early time flow

rate of 79 Us the minimum specific pressure significance level is 38 kPu-s/m3.

73.6.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched on the Stallmum dimensionless

pressure-time graph (Figure 7.27). On this graph, the lower Curve is the Theis line source

solution while the upper curve is that for a no-flow boundary with rz/rl= 1.8. The data show

an early time deviation from the theoretical boundary curve until the specific pressure draw-
down reaches about 10 kPu-s/m3. Non-uniqueness of match is a problem with this data. The

data has a marginal number of two points above the specific pressure significance level of 38
- 76 -

TABLE 7.7
TEST B7: BR34 RESPONSE TO BR31 DISCHARGE

TEST SPECIFICATIONS

Observation Well Source Well


~

Well No. BR34 BR31

Coordinates S 80200.84 ft W 47490.01 ft S 78498.30ft W 46766.50ft


(Broadlands Grid)

R.L. (C.H.F.) 308.8 m 297.7 m


(MoturikiDatum)

Permeable Depth 730~900-1000*

Drilled Depth 2587 m 1252 m

Open Hole D i m 216 mm 168 mm

Interwell Distance 564 m

Discharge Rate refer Table 7.8

Discharge Temp. -
Discharge Enthalpy 1270 W/kg

Recording Meter Water Level Chart Recorder

Flow Start Time 1130 Mar. 25 1980

Drawdown Period 293 h

Total Test Time 293 h

Pa-s/m3 which are matched to the line source. The match to the curve for rz/rl * 1.8 is at the

limit of resolution (refer §5.2). The final match was influenced by the need to bave an excel-

lent fit of the late time drawdown data which showed an upward trend that could not be fitted

to the line source solution. The chosen match point gave an excellent fit to the data above the
critical specific pressure level for a no-flow boundary located at r2/rl= 1.8.
- 77 -

TABLE 7.8
TEST B7: BR34 RESPONSE TO BR31 DISCHARGE
DRAWDOWN FLOW RATE DATA

293.0

10

j
0
h
-g 0.1

.-6
.-i
10

I'
0.01

' * I....,
0.001 I I
I

0.01 0. 1 1 10 100
Dinoncioniocs T i n e
Fig. 7.27 Log-log match of data to Stallman type curves.

ps = 10 kPa-s/m3 p D = 0.031

t=10h -lD-- 0.037


4
- 78 -
Values of permeability and storativity derived from the match point are presented in Table 8.1.

The data show minor oscillations about the theoretical curves which may be the result of earth

tides or residual barometric fluctuations.

73.6.2. Semi-log Analysis

A semi-log graph of the data is shown in Figure 7.28. The lower curve is

0.01 0.1 1 10 100


Oineneioniooo Time

Fig. 7.28 Semi-log match of data to semi-infinite type curves.

the Theis line source solution while the upper curve is the solution for the linear ino-flow boun-

dary at r2/rl = 1.8. In the semi-log plot the definition of the data match with the theoretical

solutions is improved for higher values of pressure and time. The graph shows that the match

of the late time data with the chosen boundary position is excellent.

The semi-log and log-log graphs confirm the location of the no-flow boundary. The distance

between the source well and the observation well, rl, is 564 meters. Hence the distance

between the observation well and the image well, r2, is:

r2 = (564)(1.8) = 1015 meters


- 79 -
The corresponding inference ellipse to which the postulated linear boundary is tangent is

presented in Figure 7.29.

Fig. 7.29 Inference ellipse location.


- 80 -
73.7. TEST B8: BR34 RESPONSE TO BR23 DISCHARGE

Specific pressure verses time data for this test are presented in Table B7 and pbtted in carte-

sian form in Figure 7.30. Test specifications are shown in Table 7.9 and

300
0
0
E 0
3
250 0
0 e
\
u) e
200 e
n
& 0
x
Y
0
150 e
u) e
u) e
Q,
nc 100 e
e
0 e
Q
50
0
v)
e
0 I I I I 1
0
0 50 100 150 200 250 300
Time (hr)

Fig. 7.30 Cartesian graph of the interference data

variations in flow rate are detailed in Table 7.10. A 32% LOP 1.A flow rate occur&( in the final

10 hours of the test.

For a minimum pressure significance level of 3 kPa (refer $3.4)and an average ebrly time flow

rate of 71 Its the minimum specific pressure significance level is 42 kPu-stm3.

73.7.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched to the Theis line source solu-

tion (Figure 7.31). The data show an excellent match to line source solution at $11 times. The

match is non-unique however as a match could also be obtajned to that of a doundary with

r2/rl of less than 2 (refer 95.2). The match point for the line source solution is:
- 81 -

TABLE 7.9
TEST B8: BR34 RESPONSE TO BR23 DISCHARGE

TEST SPECIFICATIONS

Observation Well Source Well

Well No. BR34 BR23

Coordinates S 80200.84 ft W 47490.01 ft S 78923.12ft W 45215.27 ft


[Broadlands Grid)

R.L. (C.H.F.) 308.2 m 291.8 m


[Moturiki Datum)

Permeable Depth 1015-1055

Drilled Depth 2587 m 1097 m

Open Hole D i m 216 mm 194 mm

Interwell Distance 795 m

Discharge Rate refer Table 7.10

Discharge Temp.

Discharge Enthalpy 1209 kllkg

Recording Meter Water Level Chart Recorder

Flow Start Time 1650 Mar. 2 1981

Drawdown Perid 223 h

Total Test Time 223 h

The value of storativity which is relatively insensitive to the matched pressure level is present-

ed in Table 8.1.
I, 1 1

- 82 -

TABLE 7.10
TEST B8: BR34 RESPONSE TO BR23 DISCHARGE
DRAWDOWN FLOW RATE DATA

Time FlowRate

223.0 70.9
233.0 48.5

10

0.1

0.01

0.001
0.01
I
/ . . .I

0.1
I

1
I

10
Y

100
Dimensionless Time
Fig. 731 Log-log match of data to Theis type curve.

73.73. Semi-log Analysis

The data show an excellent fit to the Theis line source solution but this match is non-unique

for the reasons discussed in the previous section (Figure 7.32).


- 83 -

0.01 0.1 1 IO 100


Dinoneionless Time

Fig. 7.32 Semi-log match of data to line source solution.


- 84 -
73.8. TEST B9: BR34 RESPONSE TO BR19 DISCHARGE

Specific pressure verses time data for this test are presented in Table B8 and plotted in carte-

sian form in Figure 7.33. Test specifications are shown in Table 7.l l and

E 3000
4

2000

(I)
(I)
0,
c
1000
0
0,
P
v)

0
0 100 200 300 400
I

Time ( h r )
Fig. 7.33 Cartesian graph of the interference data

variations in flow rate are detailed in Table 7.12.The flow rate variation was fitted to the fol-

lowing correlation for the total discharge time:

(7.21)

For a minimum pressure significance level of 3 kPa (refer §3.4)and an average early time flow

rate of 70 Us the minimum specific pressure significance level is 43 kPu-s/m3.

73.8.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched on the Eipper dimensionless

pressure-time graph (Figure 7.34). On this graph, the lower curve is the Theis line source
solution while the upper curve is for a no-flow boundary with rz/rl = 1.5. The 4lrawdown data
- 85 -

TABLE 7.11
TEST B9: BR34 RESPONSE TO BR19 DISCHARGE

TEST SPECIFICATIONS

Observation Well Some Well

Well No. BR34 BR19

Coordinates S 80200.84 fi W 47490.01 fi S 78224.54fi W 46154.80ft


(Broadlands Grid)

R.L.(C.H.F.) 308.2 m 294.4 m


(Moruriki D a m )

Permeable Depth - 600-1074m

Drilled Depth 2587 m 1074 m

Open Hole Diam. 216 mm 194 mm

Interwell Distance 727 m

Discharge Rate refer Table 7.12

Discharge Temp.

I Discharge Enthalpy 1180 Wlkg

Water Level Chart Recorder

Flow Start Time 1200 Mar. 2 1980

Drawdown Period 143 h

Total Test Time 455 h

show an early time deviation from the chosen match to the line source solution until the

specific pressure drawdown reaches about 20 kPu-s/m3 after which an excellent fit is obtained.

The buildup data match well to the buildup portion of the line source solution until a time of

320 hours when the data drop below the theoretical curve indicating interaction with pressure

support in the system. Possible constant pressure sources are discussed in Test B1 (5 7.3.1).
- 86 -

TABLE 7.12
TEST B9: BR34 RESPONSE TO BR19 DISCHARGE
DRAWDOWN FLOW RATE DATA

35.0 68.0
57.0
85.0 67.3
103.0
132.0 66.2

0.01 0.1 1 10 100


Dinenoionlesr T i n e
Fig. 7.34 Log-log match of data to Eipper type curves for rzlrl -1.5.

The match points were:

ps = 10 kPa-slm3 p~ = 0.048

t=100h -tD-- 1.00


6
- 87 -
While the match shown in Figure 7.34 appears to be satisfactory a non-uniquenass problem ex-

ists with the data. Early time line source behavior is masked by the anomalous bhavior of data

for specific pressure levels of less than 20 kPa-s/d. An alternative match of the data to the

Theis line source solution alone is shown in Figure 7.34. This match appears eqpally as plausi-

ble. The reason for this is that the resolution of the data is insufficient to distinguish between

the line source solution and semi-infinite solutions for r2/rrratios 5 2 (refer 05.3. The buildup

data follow the theoretical behavior for only one-half a log cycle providing no additional assis-

tance in defining the match point. Hence the non-uniqueness of the match.

Figure 7.35 shows the match point to the line source solution as:

f
h
-a
&
.-
5
.-
0.1

0.01
0.01 0.1 1 10 100
DinenrionIeee Time

Fig. 7.35 Log-log match of data to line source solution.

pz = 10 kPa-s/m3 p~ = 0.015

rD
r=100h -= 0.78
a
The permeability ratio between the two matches is three, while the storativity ratio is only 1.25
- 88 -
which demonstrates the greater sensitivity of permeability to the matching process. The stora-

tivity derived from this match point which is relatively insensitive to changes in the matched

pressure level is presented in Table 8.1 for the match to the curve for r$rl = 1.5.

Both the drawdown and buildup data show oscillations about the theoretical curyes which may

be the result of earth tides or residual barometric fluctuations.

73.8.2. Semi-log Analysis

Semi-log plots for the r2/tl = 1.5 match and the line source match are shown in Figures 7.36

and

0.01 0.1 1 10 100


Dinensionless fine

Fig. 7.36 Semi-log match of data to semi-infinite type curves for r2/r1 -1.5.

7.37 respectively. Both plots show the same quality of fit with the characteristics idescribed for

the log-log plots except that the early time data deviation is masked by the semi-lag scale.

Drawdown and buildup data again show oscillations about the ,theoretical curves which may be

the result of earth tides or residual barometric fluctuations.


- 89 -

0.01 0.1 1 10 100


Dineneionierr fine
Fig. 7.37 Semi-log match of data to line source solution.

73.83. Horner Analysis

The Horner plot for rz/rl = 1.5 is presented in Figure 7.38 and confirms an excellent match of

the early time buildup data with the theoretical solution even though a unique match is not pos-

sible. Oscillations in the data are still evident and the effect of pressure sup@rt is seen for

Horner times I2. An alternative graph for the line source match (Figure 7.39) shows similar

characteristics.
I, 1 1

- 90 -

1 10 100
Horner Time

Fig. 7.38 Match of buildup data on a Horner graph for r2/rl=1.5 .

1 10 100
Hornr limo

Fig. 7.39 Match of buildup data on a Horner graph to line source solution.
- 91 -
73.9. TEST B10: BR23 RESPONSE TO BR13 SHUT-IN

Specific pressure verses time data for this test are presented in Table B9 and plotted in carte-

sian form in Figure 7.40. Test specifications are shown in Table 7.13.

600
-
E
3
0
\
(rD

b 400
a
r.
Y

(I)
v,

at 200
0
a,
a
cn
0
0 100 200 300
Time Ihr)
Fig. 7.40 Cartesian graph of the interference data.

This test was originally part of Test B2. Large fluctuations in flow rate at the end of the draw-
down period meant that data obtained was not suitable for analysis. After five d$ys of fluctua-

tions, pressures in well BR23 fortuitously reached steady state. The well was then shut and the

buildup performance monitored.

For a minimum pressure significance level of 3 &Pa (refer $3.4) and an average early time flow

rate of 52 Ns the minimum specific pressure significance level is 58 kPu-s/m3.

73.9.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched on the Theis line source solu-

tion graph (Figure 7.41). The data give a good match to the line source for all times except

for the first data point which is below the minimum specific pressure significance level of 58
- 92 -

TABLE 7.13
TEST B1Q BR23 RESPONSE TO BR13 DISCHARGE

TEST SPECIFICATIONS

Observation Well Source Well

Well No. BR23 BR13

Coordinates S 78923.12ji W 45215.27 j2 S 79705.90ji W 4568t.30 ft


(Broadlands Grid)

R.L. (C.H.F.) 291.8 m 293.5 m


(Morwiki Darwn)

Permeable Depth -
1015 1055 915 * 30 m

Drilled Depth 1097 m 1081 m

Open Hole D i m 194 mm 197 mm

Interwell Distance 279 m

Discharge Rate 51.9 11s

(volumetric average rate during drawdown)

Discharge Temp.

Discharge Enthalpy 1050 kllkg


~ - ~~~~

Recording Meter Water Level Chart Recorder

Buildup Start Time 1600 April 21 1980


Total Test Time 234 h

Wu-slm3 and chose after 100 hours of the test where there is a slight indication of a no-flow

boundary.

The match points chosen were:

pa = 10 k~a-strn~p D = 0.04
- 93 -

10

1
f
0
1)
0
c
a

-:
0
0.1

0
-
i-
= 0.01

0.001
0.01 . 0. I 1 10 100
Dimensionless Tine
Fig. 7.41 Log-log match of data to Theis type curve.

r=10h -t D-
6-
Values of permeability and storativity derived from the match point are presented in Table 8.1.

The data appear to show oscillations about the theoretical curves which may be the result of

earth tides or residual barometric fluctuations.

73.9.2. Semi-log Analysis

A semi-log graph of the data on the Theis line source solution is shown in F i a r e 7.42. The

scale of the semi-log graph gives better definition to the no-flow boundary effect which is

shown by the data dropping below the line source solution at r&, = 18 or t = 170 hours.

The data show oscillations about the line source curve which may be the result of earth tides

or residual barometric fluctuations.


- 94 -

0.01 0.1 1 10 100


Dimensionless Time

Fig. 7.42 Semi-log match of data to the line source solution.

Better definition can be obtained using the type of Sugeev et d. [198S]. Figurer 7.43 presents

a semi-log match of the dimensionless pressure data to the type Curve presented by Sugeev et

al. [1985]. The pressure match point is:

PD = 0.0 p i = 2.40
Substituting the semi-log pressure match point into Equation 7.5 and solving for the ratio r2/rl

yields:

rdrl = exp[ @ D - p ~+) I ~ ( I W I ) ] (7.24)

r2ltl = exp(0.0-2.40 + h99) = 9.0 (7.17)

From the semi-log graph it can be noted that in order to get a unique match, only one semi-log

straight line and the transition period are required. The data show that during the drawdown

period the first semi-log straight line and the transition developed, but the seeond semi-log

straight line did not.


- 95 -

T i r (kl
-1. mtch point0

p, = 0 p; 0 9.40
I , . . . , .,.I , , , , , , . . I , , , , ,

10 100 1000 10000 1e+05


Hodif ied Dimensionless Time

Fig. 7.43 Semi-log match of drawdown data to Sugeev er al. type curves.

The corresponding inference ellipse to which the postulated linear boundary is tangent is

presented in Figure 7.44.


1 1

- 96 -

ower

Fig. 7.44 Inference ellipse location.


- 97 -
7.4. DATA FROM QUARTZ CRYSTAL GAUGES

7.4.1. TEST C1: BR13 RESPONSE TO BR20 DISCHARGE

The data for this test are presented in Table B10 and plotted in Cartesian form in Figure 7.45.

Test specifications are shown in Table 7.14.

n
700
5 600

>
0
9
Y
3
0

I
500

400
P,
f 300
QD

hb
200

x
0

v)
100

0
0 168 336 504 672 840 1008 1176
Time (hr)
Fig. 7.45 Cartesian graph of the interference data.

Equipment malfunction meant that no data was recorded for the periods 12 to 5 1 hours, 55 to

69 hours, 163 to 194 hours, 380 to 407 hours, 545 to 578 hours and 580 to 889 hours.

For a minimum pressure significance level of 3 kPa (refer §3.4) and an average mrly time flow

rate of 84 lls the minimum specific pressure significance level is 36 kpu-s/m3.

7.4.1.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched on the Eipper /dimensionless

pressure-time graph (Figure 7.46). O n this graph the lower curve is the Theis li& source solu-

tion while the upper curve is for a no-flow boundary with r2/r1 = 4.0. Despite the missing data
a half log cycle is available for a match with the Theis line source solution withlthe data after
- 98 -

TABLE 7.14
TEST C1: BR13 RESPONSE TO BR20 DISCHARGE
TEST SPECIFICATIONS

Observation Well Source Well

Well No. BR13 BR20

Coordinates S 79705.90jl W 45687.30jl S 7767 1.30ft W 447 14:Wfr


(Broadlands Grid)

R.L. (C.H.F.) 293.5 m 291.7 m


(Moturiki Datum)

Permeable Depth 915 f 30 m 815 m, 945 m, 1045 m

Drilled Depth 1081 m

Open Hole D i m 197 mm


~~ ~~

Interwell Distan; 687 m

Discharge Rate 84.0 11s

Discharge Temp. 260OC

Discharge Enthalpy

Recording Meter quartz crystal gauge

Flow Start Time 1430 May 1 1984

Drawdown Period 339 h

Total Test Time 1217 h

69 hours showing a good fit to the upper no-flow boundary curve.

The buildup data match to the no-flow boundary curve for r21rl = 4.0 up to 540 hlours at which
I
point the data deviate above this curve. If real, this trend would indicate the presence of a
I

second no-flow boundary in the system. No other test was run for a similar length of time so

there is no data with which to confirm this trend. The pressure and time match points for r2/r1
- 99 -

10
1

Fig. 7.46

are:

pa = 10 wa-slrn3 pD = 0.040

Values of permeability and storativity derived from the match point are presented f n Table 8.1.

The drawdown data show minor oscillations about the theoretical curves which may be the

result of earth tides or residual barometric fluctuations. I

7.4.1.2. Semi-log Analysis

A semi-log graph of the data is shown in Figure 7.47. The lower curve is the Thdis line source ~

solution while the upper curve is the solution for a linear no-flow boundary at r2kl = 4.0. This

graph confirms the good fit of the the late time drawdown data and the initial bdildup data to

the upper no-flow boundary curve. ~


0.1 1 10 100
Dinenoionleoe Time

Fig. 7.47 Semi-log match of data to semi-infinite type curves.

Drawdown data again show oscillations about the theoretical curves which may ble the result of

earth tides or residual barometric fluctuations.

7.4.13. Horner Analysis

The Horner plot for r2/rl = 4.0 (Figure 7.48) shows a good match of the early time buildup

data with the theoretical curve. The periodic nature of fluctuations about the theoretical curve

indicate that this is more likely due to earth tides rather than barometric pressure variations.

Deviation of the data above the linear boundary curve for r& = 4.0 at Horner times less than

three (t 2 54Ohours) show possible contact with a further no-flow boundary. Since no other test

showed a similar trend further analysis of the second no-flow boundary has not been undertak-

en.

The Horner plot confirms the match obtained form the semi-log and log-log graphs. The dis-

tance between the source well and observation well, TI, is 687 meters. Hence the distance

between the observation well and the image well, r2, is:
1, I1

- 101 -

1 10 100
brier T i n e

Fig. 7.48 Match of buildup data on a Horner graph.

r2 = (687)(4.0) = 2748 meters

The corresponding inference ellipse to which the first postulated linear no-flow boundary is

tangent is presented in Figure 7.49.


'ower
itation site

E
a

Fig. 7.49 Inference ellipse location.


~~

- 103 -
7.4.2. TEST C2: BR23 RESPONSE TO BR20 DISCHARGE

The data for this test are presented in Table B11 and plotted in Cartesian form in Figure 7.50.

Test specifications are shown in Table 7.15.

CI 700
600 -

so0

400

300

200 -
&
100 r
1
I
I I I I I I I
0
0 168 336 SO4 672 840 1008 1176
Time (hr)
Fig. 7.50 Cartesian graph of the interference data.

I
Equipment malfunction meant that no data was recorded for the periods 12 to 54 hours, 55 to

69 hours, 163 to 194 hours, 380 to 407 hours and 545 to 578 hours.

For a minimum pressure significance level of 3 Wu (refer $3.4) and an average mrly time flow

rate of 84 11s the minimum specific pressure significance level is 36 Wu-s/d.

7.4.2.1. Log-log Analysis

A log-log pressure-time graph of the data is drawn and matched on the Eipper dimensionless

pressure-time graph (Figure 7.51). On this graph the lower curve is the Theis line source solu-

tion while the upper curve is for a no-flow boundary with &,= 4.0. Despite the missing data

a half log cycle is available for a match with the Theis line source solution. With the data after

69 hours showing an excellent fit to the upper no-flow boundary curve.


I, I I

- 104-

TABLE 7.15
TEST C2: BR23 RESPONSE TO BR20 DISCHARGE

TEST SPECIFICATIONS

Observation Well Source Well

Well No. BR23 BR20

Coordinates S 78923.12fl W 45215.27 fr S 77671.30fr W 44714.9a)ft


(BroadlandsGrid)

R.L. (C.H.F.) 291.8 m 291.7 m


(Motwiki D a m )

Permeable Depth -
1015 1055 m 815 m, 945 m, 1045 rn

Drilled Depth 1097 m -


Open Hole D i m 194 mm -
Interwell Distance 411 m

Discharge Rate 84.0 11s

Discharge Temp. 260OC

Discharge Enthalpy

Recording Meter quartz crystal gauge

Flow Start Time 1430 May 1 1984

Drawdown Period 339 h

Total Test Time 579 h

The buildup data show a fair match with the buildup solution for a boundary located at

r2/rl = 4.0. The pressure and time match points for rz/rl = 4.0 are:

p. = 10 kPa-s/m3 p D = 0.060
~

tD
t=100h - --16.0
to
I I I

- 105 -

10

0
s
a
a
0 1
k
0

-P
C
.-0a 0.1
c
.-
0

0.01
0.01
A
I
0.1
I

1
I

10
J
1M
Dimensionless Time
Fig. 7.51 Log-log match of data to Eipper type curves.

Values of permeability and storativity derived from the match point are presented in Table 8.1.

The drawdown data show minor oscillations about the theoretical curves which may be the

result of earth tides or residual baromemc fluctuations.

7.4.2.2. Semi-log Analysis

A semi-log graph of the data is shown in Figure 7.52. The lower curve is the "%&isline source

solution while the upper curve is the solution for a linear no-flow boundary at rdrl = 4.0. This

graph confirms the good fit of the the late time drawdown data and the initial buildup data to

the upper no-flow boundary curve.

Drawdown data again show oscillations about the theoretical curves which may be the result of

earth tides or residual baromemc fluctuations.

7.4.23. Horner Analysis

-
The Homer plot for r2/rl 4.0 (Figure 7.53) shows a fair match of the early time buildup data

with the theoretical curve.


I
I

- 106 -

0.1

Fig. 7.52 Semi-log match of data to semi-infinite type curves.

1 10 100 lo00
tiornor Time

Fig. 7.53 Match of buildup data on a Homer graph.


- 107 -

The value of using the Homer analysis is demonstrated in identifying the slight mismatch. The

mismatch of the buildup data on the log-log and semi-log plots was barely discernible. The

Homer plot expands the time scale on the initial buildup period making any mismatch of data

much more obvious.

The Homer plot confirms the match obtained form the semi-log and log-log graphs. The dis-

tance between the source well and observation well, rI, is 411 meters. Hence the distance

between the observation well and the image well, r2, is:

r2 = (411)(4.0) = 1644 meters

The corresponding inference ellipse to which the postulated linear no-flow boundary is tangent
is presented in Figure 7.54.
I / I

- 108 -

Fig. 7.54 Inference ellipse location.


- 109 -
7.4.3. TEST C3: BR34 RESPONSE TO BR20 DISCHARGE

The data for this test are presented in Table B12 and plotted in Cartesian form in Figure 7.55.

Test specifications are shown in Table 7.16.

700

3
0
2 500 -
b
Q.
Y
Y
400 -

300 -
0)

ti 200 -

0
g
v)
100 -

O l I I I 1 I I I

0 168 336 504 672 840 1008 1176


Time ( h r )
Fig. 7.55 Cartesian graph of the interference data.

Equipment malfunction meant that no data was recorded for the periods 12 to 54 hours, 55 to

69 hours, 163 to 194 hours, 380 to 407 hours and 545 to 578 hours.

For a minimum pressure significance level of 3 kPu (refer 83.4) and an average early time flow

rate of 84 Us the minimum specific pressure significance level is 36 Wu-s/m3. I

7.4.3.1. Log-log Analysis

Log-log pressure-time plots of the data are drawn and matched on the Eipper dimensionless

pressure-time curves (Figures 7.56, 7.57). Both plots show an excellent match to the data 1

above the critical pressure level of 36 kPu-s/m3 thereby demonstrating that from the log-log

plots alone a unique match is not possible. On Figure 7.56, the lower curve *presents the

Theis line source solution while the upper curve is for a no-flow boundary locdted at r2/r1=
- 110 -

TABLE 7.16
TEST C3: BR34 RESPONSE TO BR20 DISCHARGE
TEST SPECIFICATIONS

Observation Well Source Well

Well No. BR34 BR20

Coordinates S 80200.84ft W 47490.018 S 77611.308 W 44714190ft


(Broadlands Grid)

R.L. (C.H.F.) 308.2 m 291.7 m


(Moturiki Datum)

Permeable Depth 815 m, 945 m, 1045 h

Drilled Depth 2587 m

Open Hole D i m 216 mm

Interwell Distance 1145 m

Discharge Rate 84.0 11s

Discharge Temp. 260°C

Discharge Enthalpy -
Recording Meter ~
quartz crystal gauge

Flow Start Time 1430 May 1 1984

Drawdown Period 339 h

Total Test Time 579 h

1.5.

The pressure and time match points are:

ps = 10 Pa-slm3 p~ = 0.030

tD
t=100h -= 0.78
6
I 1

- 111 -

10

i
(D
/ match poinl b

PL I ?
.
I s
IO Dm 0.-
mow l o - o.n
a
0

-P

0
s
.-
K
.-
E
Ool

0.01
0.01 0.1 1
\
10 100
Dimensionless Time
Fig. 7.56 Log-log match of data to Eipper type curves for r2/rl -1.5.

10

lo00 mtch point8


Pb IO Dm 0.0100
1 0

I im w I 0.u

.-
Q
0.1

0.01
0.01 0.1 1
Dimensionless Time
Fig. 7.57 Log-log match of data to line source solution.
1 10 100
- 112-
Figure 7.57 represents the match of the data to the line source solution. The match points for

this graph were:

ps = 10 Wa-slm3 p~ = 0.015

tD
t=100h -= 0.52
6
This interference test poses singularly difficult problems which make a unique /match to the ,
,
data impossible. The specific issue of non-uniqueness is covered in detail in $5.

The observation well BR34 is in a non-productive part of the field and could &+refore be ex-

pected to be close to a hydrological boundary. The distance between the sourcef well, BR20,

and BR34 is 1145 meters which is nearly twice as large as any of the other tesq. A combina-

tion of both these two factors means that a boundary if present could have an/r2/rl ratio of
i

close to unity. If this is the case then the Theis line source behavior may not be 4tectable dur-

ing the drawdown period as it will occur at pressure levels well below the critihl level of 36

kPu-s/m3. Test B9 (0 7.3.8) illustrated this point. Inspection of Figures 7.56, 7.57 bnfirms:

(i) That there can be no match to the line source since the data that do exist in b e early time

region (r&j I0.5)are below the critical significance level of 36 kPu-s/m3. ,

(ii) That there is less than one-third of a log cycle of buildup data which is insubcient to ob-
I
tain a unique match of the data to a buildup curve for a given boundary locdtion.
I

(iii) The drawdown data above the critical specific pressure level lie entirely wi&n the transi-

tion flow region.

The net result of these observations is that a unique match cannot be made frop the log-log

graph. The value of storativity which is relatively insensitive to the matched pdssure level is

presented in Table 8.1.

7.4.3.2. Semi-log Analysis

Semi-log plots of the data are shown in Figures 7.58, 7.59. On Figure 7.58 thd lower curve
- 113 -

Dimeneionieer T i m e

Fig. 7.58 Semi-log match of data to semi-infinite type curves for r$rl ~ 1 . 5 .

f
a

f
-8 2

.-6
E
.-
0

Dinmrioniecr fine

Fig. 7.59 Semi-log match of data to line source solution.


- 114-
represents the Theis line source solution while the upper curve is the solution far a linear no-

flow boundary at r2/rl = 1.5. On Figure 7.59 the data is matched to the line source solution

alone. The excellent fits of data at all times on both plots adds no additional definition to the

boundary location.

7.4.33. Horner Analysis

Horner plots for r2/r1 = 1.5 a m 13r the Theis line source solution are presented ixb Figures 7.60

and 7.61 respectively.

1 10 100 lo00
Wner T i a s

Fig. 7.60 Match of buildup data on a Horner graph for r2/rl I1.5.

The Horner plot expands the time scale on the initial buildup period making any mismatch of

data much more obvious. Unfortunately the missing data along with the difficulty of obtaining

a unique match for the low rz/rl ratio mean that the location of the inference ellipse cannot be

determined.
I / I

- 115 -

3
fP
h
0

-:2
.-5
0
E
.-
1

0
1 10 100
Horne~T i m e

Fig. 7.61 Match of buildup data on a Homer graph to line source.


- 116 -

8. RESULTS

8.1. TRANSMISSIVITIES AND STORATIVITIES

Transmissivities and storativities are calculated using Equations 2.13 and 2.14. Results are

shown in Table 8.1.

TABLE 8.1

TRANSMISSIVITlES AND STORATMTIES

EST OBS. WELL SOURCE WELL INTERWELL *fh


DISTANCE

(m) (tdkp~io~)
3

B1 23 13 279 9.25 94 0.65

B2 23 13 279 2.50

B3 23 19 357 2.80 115 4.40

B4 23 13 279

B6 23 20 411

B7 34 31 ,564 1.80 49 15.00

B8 34 23 795 40.00

B9 34 19 727 5.20

B 10 23 13 '279 9.25 64 2.80

c1 13 20 687 4.00 64 1 .oo

c2 23 20 411 4.00 96 1.30

c3 34 20 1145 1.70

Transmissivities range from 49 d-rn in Test B7 to 115 d-m in Test B3. The average is 80 d-m.

Test B7 was performed between wells BR34 and BR31 which are closer to the resistivity
- 117 -
boundary of the field than any of the other well pairs for which transmissivities were obtained.

The lower permeability recorded for this test may be a reflection of lower permeQbility in this

part of the field. The remaining well tests show relatively uniform values of transmissivity

values with the average of these being 82 d-m and the range being 140% to 60% of t h i s value.

An approximate value of horizontal permeability can be calculated by assuminlg a reservoir

thickness of 700 meters (Figure 6.3):


~

k=- 8Ox1O3
700
k = 114 md
Storativites show a much greater range than the transmissivities possibly reflectiog the forma-

tion of small two phase zones during drawdown in those regions where the permeability is

likely to be low. Tests B7, B8, B9, C3 for which pressures were measured at Bk34 could be

expected to fall into this category. Tests B7 and B8 show significantly larger storativities than

any of the other tests. Test B9 has a large storativity and Test C3 shows a near average stora-

tivity. Flow rates in these tests ranged from about 66 11s to 84 lls. While there codd be expect-

ed to be some significance in the trend of storativities with the location of the doublet in the

field it does not show conclusively in the results. Excluding the anomalously high storativities

for Tests B7 and B8 the average is 2.4 x lo4 d k P a with the range varying from 158% to 30%

of this value. An estimate of the compressibility of the system can be made assuming that the

porosity in the reservoir is 20% and the reservoir thickness is 700 meters.

+~,h= 2.4~10~
dkPa

From Bixley [1982](Figure 6.4)the average reservoir pressure in the productive region of the

reservoir at the average depth of 800 meters is 7 MPa. The average enthalpy per test of

discharged fluid was 1152 Wlkg which corresponds to a brine temperature of 263°C. The

compressibility of 263°C fluid at a pressure of 7 MPa is 1.7~10~


Pa-'. The fact that this value
- 118 -
is close to that calculated using the assumed reselrvoir properties probably indicgtes that there

are no permanent two phase zones present within the sphere of influence of thdse tests. This

but does cover an area roughly approximating that contained within the 5 o&tn

boundary of the field.


4
sphere of influence is assumed not to include that portion of the field covere by the data

which showed both pressure support and no-flow boundary trends in Tests B1, 2,B9 and C3

resistivity

8.2. LINEAR BOUNDARY LOCATION I

The location of the principal no-flow boundary is shown in Figure 8.1. The

the inference ellipses for Tests B1, B3, B7, B10, (31, C2 can be interpreted to

trending hydrological barrier near the resistivity bcundary in the northern part of lthe field. The

NE-SW trend of faults located from geological interpretation support the locatibn of the no-
flow boundary (Figure 6.5). The location is also supported by the fact that the +-them boun-

dary of those wells defined by Sixley [1982] as having the "lowest" pressure grkdient (Figure

6.5) coincides closely in direction and location with the located no-flow boundary

If Test B3 is excluded from the analysis it is poslsible that the no-flow bound
cated in a N-S direction under the power station site. This location is suppo*d
4
geology which shows a fault scarp at about the same location in this direction./However the
could be lo-

by surface

weight of evidence points strongly to the first located position. Other possible lo+tions of bar-
~

riers to the south and east of the field cannot be convincingly supported by the $cation of the

inference ellipses.
1 I /
I

- 119-

Fig. 8.1 Inference ellipse locations for Tests B1,B3,B7, B10,C1,C2.


- 120 -

9. DISCUSSION

9.1. GENERAL (Sugeev et al. 1986)

The main advantage of interference testing over source well testing for detectin$ reservoir lim-
. its is the increased characteristic length scale. When analyzing source well pressure data, the

characteristic length scale is the diameter of the wellbore, that is typically on the order of a few

inches. In interference testing, the characteristic length scale is the distance between the obser-

vation well and the source well, that is on the order of tens or hundreds of mete@. Also, some

near wellbore effects such as wellbore skin in a source well of a constant rate test are not

significant in interference testing.

The main disadvantage of interference testing is the decreased amplitudes df the pressure

changes. As the observation well is located further form the source well, the space resolution

increases but the magnitude of the pressure changes decreases. These are co*peting effects

that have to be addressed during the design stages of a test. When the early time pressure

changes are small, the effects of earth tides and barometric pressure may be slignificant, and

might lead to large uncertainties in the estimation of reservoir parameters.

In the case of interference testing for linear boundary detection, there is an addefl problem that

is purely geometrical. For distance ratios smaller than 2, it is very difficult to detect the pres-

ence of an impermeable linear boundary regardless of the actual distance betwedn the observa-

tion well and the source well.

9.2. ANALYSIS TECHNIQUES

The use of a variety of type curves assists in obtaining the correct match. The log-log Stall-

mun type curves (Figure 2.2) for drawdown are useful in matching the early time data to the
line source solution. Once a match of the early time data has been obtained semi-log plots of

the Stallman type curves (Figure 2.4) can be used to match the late time drawdown data. The

Eipper type curves (Figure 2.3) are useful for obtaining both the early time data match and
- 121 -
where applicable a preliminary match of the buildup data. Once a preliminary match has been

made a Horner plot can c o n h the match of the drawdown data. For r2/rl 2 10 in drawdown
tests the Sugeev et al. type curves (Figure 2.5) can be used to determine the r2kl ratio with thq

Fox type curve (Figure 2.6) used to confirm this ratio on the buildup data. =termination of

the inference ellipse from the rdrl ratio is performed using Vela’s technique (Equation 2.12).

93. DATA MATCHING

Analysis of the test results showed that early time data up to the minimum sigrpificant pressurd

level of 3 kPa (0 3.4) failed to match the theoretical curves due to the effect of earth tides,

Clipping techniques often exclude data which fall when the overall trend is a rise and hen4
I
earth tide oscillations are either filtered out or occur below the resolution of the gauge. Tesq

B3, B6, B7, B9, B10, C3 exhibited this characteristic at early time.

Interaction with more than one hydrological heterogeneity was indicated in Tests B 1, B2, B9,
I

C1 due to the deviation of data at late dimensionless time from the initial matcbed curve. Tesq

B1 and B9 showed pressure support effects while tests B2 and C1 indicatdd contact w i d

second no-flow boundaries. The lack of consistency in these late time trends pdecluded furthd

analysis. All the trends with the exception of Test B2 appear to be credible. Test B2 is les

credible because the appearance of the two no-flow boundaries was not confiqed in Test B1

which was performed immediately following Test B2. It is possible that the pressure supporq

seen in Tests B1 and B9 was due to a two phase zone which later collapsed add therefore
4
did^
not appear in Test C 1. I
The interference response of Test B4 was affected by gravity segregation, tlhermal, moving
I

front and deposition effects associated with the injection of colder separated brine into thq

reservoir. It may be possible to analyze all these effects with the aid of a thermal simulator bu t
this has not been attempted in this study.

Non-uniqueness problems in Tests B6, B8, B9, C3 meant that values of pemeabilities wera

not obtained for these tests but values of storativity which are less affected by the non-1
- 122 -
uniqueness problems were retained in Tests B8, B9, C3. Test B6 was too shdrt for a match.

Tests B6 ,B9 and C3 were expected to show no-flow boundaries but the boundaries if present

occurred at rdrl ratio's of less than about two making detection improbable (0 52). I

9.4. RESERVOIR PARAMETERS

- The principal results obtained were that the average test transmissivity was 80 d-m and tha

average storativity was 2 . 4 ~ 1 0


m~W a with both values having a range of abaut 350% if tha

two high storativity values from Tests B7 and B8 are excluded. For a porosity of 20% an4

reservoir thickness of 700 meters the average test compressibility is 1.7~1O~kPa-'


and thd

average test permeability is 114 md. The estimated compressibility agrees with that for brine ad
I
the average production temperature 263°C indicating that there are no resident two phase zones1

in the field. Interpretation of the superposed inference ellipses indicates that there is a NE-SW

trending no-flow boundary near the extreme northern end of the productive afea of the field
I
I
(Figure 8.1).

9.5. OPTIMIZING TEST PARAMETERS

Interference tests benefit from the following features:

(i) Atmospheric earth tide and barometric pressure measurements measured to 1 Pa which al.

low resolution and deconvolution of earth tide components.


~

(ii) Downhole gauge resolution better than 100 Pa to allow correction of primsry interference/

test data for earth tide effects and also allow the deconvoluted earth tide data to be used~
~

to determine reservoir parameters.

(iii) Flow rates large enough to produce early time data that can be matched to the line source
solution. Large flow rates may produce two phase effects in the field dule to drawdown,

during production.

(iv) Doublets chosen to produce r2/r1 ratios in the range 3 to 8 on the best inffonnation avail-

able to the designer prior to the start of the test.


- 123 -
(v) Cautious use of interference tests with injected fluid as complex analysi$ is required to
determine in-situ reservoir parameters.
I , 1
I

- 124 -
10. CONCLUSIONS

(i) Careful interpretation is required to obtain the correct reservoir propembs from interfeq-

ence tests.

(ii) Assuming a reservoir thickness of 700m the average permeability foII the productivk
Ohaaki reservoir is 110 md.
I
(iii) The average compressibility of 1 . 7 ~ 1 0 derived ~ ' the interferende tests indicateb
~ ~ ~ from
I

that there are no mobile two phase zones in the Ohaaki reservoir within the sphere OF

I
influence of the test.

(iv) A no-flow boundary was located utilizing the inference ellipses deduced $om 6 of the lq

tests analyzed.
I
(v) Tests B1 and B9 showed evidence of late time pressure support whichimay have bee4
due to a two phase zone in the reservoir which later collapsed, while t e s t C1 showe

evidence of a second no-flow boundary.


- 125 -
NOTATION

formation volume factor for brine

formation volume factor for gas

Cf formation compressibility (Wu-')

vapour compressibility (Ma-')

compressibility due to phase change ( f l u - ' )

total compressibility (Pa-')

CW brine compressibility (Pa-')

d well depth (m)

dD dimensionless distance between source well and boundary

permeable zone depth (m)

depth of water surface (m)

mamx shear modulus

h formation thickness (m)

latent heat (Wlkg)

height of water surface above permeable zone (m)

permeability (durcy)

P pressure ( P u )

Pam atmospheric pressure ( P a )

Pc overburden pressure at average reservoir depth (Wu)

Pi initial pressure (Wu)

P6 specific pressure ( ~ u - s / m ~ >

Pmt saturation pressure ( P a )

PD dimensionless pressure

PL modified dimensionless pressure

PLH modified Homer dimensionless pressure


- 126 -
permeable zone pressure at the well ( P u )

saturation pressure (Pa)

volumetric flow rate ( d s )

distance between pressure point and source well (m)

distance between source well and observation well (m)

distance between observation well and the image well (m)

source well radius (m)

storage coefficient

vapor saturation

water saturation

time (s)

dimensionless time

modified dimensionless time

dimensionless Horner time

modified dimensionless Horner time

dimensionless production time

dynamic viscosity (Pu-s)

earth tide frequency

porosity

formation density ( k g / d )

vapor density (kg/m3)

brine density (kg/m3)


/ I

- 127 -
REFERENCES

Arditty P : "The Earth Tide Effects on Petroleum Reservoirs," M.S.thesis, Stu@ord Universid.
I
Stanford Geothermal Programme Report SGP-TR-34 (1979). 1

Barenblatt G I, Zheltov Y P : "Fundamental Equations of Homogeneous Liquids in Fissure1

Rocks," Dokl. A M . Nuuk. U.S.S.R. (June 1960), pp 545-548.

Barenblatt G I, Zheltov Y P, Kochina I N : "Basic Concepts in the Theory of Seepage of

Homogeneous Liquids in Fissured Rocks (Strata)," J. Appl. Math. Mech., 0.S.S.R. (1960),

pp 1286-1303.

Bixley P F : "Broadlands Geothermal Area: Underground Temperatures and so@e Thoughts 04


I
the Reservoir and Future Drilling," Ministry of Works and Development, Wairakei (1976)

Geothermal Circular PFB 1.

Bixley P F : "Hydrology of the Western Broadlands Reservoir," Ministry of Works and


I
Development, Wairakei (Feb. 1982), file 27/35.

Bixley P F, Grant M A : "Reinjection Testing at Broadlands," Proc. 6th Stunfbrd Geothermal

Workshop (1980).

Blakely M R, O'Sullivan M J, Bodvarsson G S: "A Simple Model of the Ohilaki Geothermal

Reservoir," Proc., 5th N.Z.Geothermal Workshop (1983), pp 11-15.

Bodvarsson G : "Confined Fluids as Strain Meters," J. Geophys. Research (1970), 75(14), pd

27 11-27 18.

I
Braithwaite R L : "Perrological Report: Broadlands Drillhole Br34," Gedlogical Survey4

D.S.Z.R., Lower Hutt, N.Z.(June 1979).


- 128 -

Bredehoeft J D : "Response of Well-Aquifer System to Earth Tides," J. GBophys., 72(12)

(1967), pp 3075-3087.

Wellington (1973).
1
Browne P R L : "The Geology, Mineralogy and Geothermometry of the Broa lands Geother-~

mal Field, Taupo Volcanic Zone, New Zealand," Ph.D. dissertation, Vict ria University,l

1
I

Davis G E, Hawkins M F : "Linear Fluid-Barrier Detection by Well Pressure $leasurements,"I


I
J . Pet. Tech. (Oct. 1963), pp 1077-1079.

Deruyck B G, Bourdet D P, DaPratt G, Ramey H J Jr.: "Interpretation of Inte

1
Reservoirs With Double Porosity Behaviour-Theory and Field Examples,', SPE 57th An-1I

n u l Conference, New Orleans (Sept. 1982), SPE paper 11025.


rence Tests in1
I

Earlougher R C, Ramey H J Jr.: "Interference Analysis in Bounded Systemsp J. Can. Pet.

Tech. (0ct.-Dec. 1973), pp 33-45. I I


I
1

i
I

Eipper M E : "Computer Generation of Type Curves,'' M.S.thesis, Stanford U iversiry. Stan-~


I
ford Geothermal Programme Report SGP-TR-86 (1985). I
I
I
I
Ellis A 1, Mahon W A J : Chemistry and Geothermal System, Academic W s s , New York

(1977).
I

Ferris J G, Knowles D B, Brown R H, Stallman R W :"Theory of Aquifer Tejts", U.S.Geol

Survey, Water Supply, paper 1536-E (1962). 1

Fox G : "Linear Boundary Detection using Pressure Buildup Tests," M.S.thesis, Stanfors
University. Stanford Geothermal Programme Report SGP-TR-83 (1984).

Grant M A : "Broadlands Geothermal Field, Interference Tests at Br13, 19 a

Mathematics Division, D.S.I.R.,Wellington, New Zealand (1980), Report


- 129 -
Grant M A : "Broadlands a Gas Dominated Reservoir," Geothermics, Vol. 6 (1677), pp 9-29.

Grant M A : "Changes in Fluid Distribution at Ohaki," Proc., 2nd NZ. GeotGrmal Workshop

(1980), pp 69-74.
~

Grant M A : "A Heat Balance on the Ohaaki Reservoir," AppZied Mathe&zrics Divisionb

D.S.I.R., Wellington, New Zealand (1982), Report MAG35 I

Grant M A : "The Interaction of the Ohaaki Reservoir with Surrounding / F l u i d AppZie4

Mathematics Division, D.S.I.R., Wellington, New Zealand (1982), Report $AG37. I


I

,
Grant M A : "Broadlands Reservoir Structure," AppZied Mathematics Division1 D.S.I.R., Well

lington, New Zealand (Nov. 1983).

Grant M A, Donaldson I G, Bixley P F : Geothermal Reservoir Engineering, hcademic Press,


I
New York (1982).
I
I
I

Grant M A, Iles D W : "A Simple Model of the Connection Between Easteb and Wester;
I

Sections of Broadlands Geothermal Field," Applied Mathematics Division1 D.S.I.R., We14

lington, New Zealand (1982), Report MAG34.

Grant M A, Sorey M L : "The Compressibility and Hydraulic Diffusivity o f / a Steam-Water/


I

Mixture," Water Resources Research (1979), pp 684-686.


I

Grindley G W, Browne P R L : "The Subsurface Geology of the Broadlahds Geothermal

Field," N.Z.GeoZogicaZ Survey (1968), Report 34.

Hall H N : "Compressibility of Reservoir Rocks," Trans., AIME, 198 (1953), pb 309-311.


I

Henley R W, Harper R T : "Water-Rock Interactions During Injection of Hasbed BR2 Wateq

at BR34 - Preliminary Interpretations and Some Implications," Chedistry Division4

D.S.I.R., New Zealand (1979). I ~


- 130 -
Hitchcock G W, Bixley P F : "Observations of the Effect of a Three Year Shuwown at Broad-

lands Geothermal Field," Proc., 2nd UN. Cogerence on the DeveZopdpr and Use of

Geothermul Resources (1973, pp 1657- 166 1.

Hedenquist J W : "Characteristics of Broadlands-Ohaaki Water Chemistry and changes Subjed

to Initial Production," Proc., 5th N.Z.Geothermal Workshop (1983), pp 13b-144.


I

Hedenquist J W, Stewart M K : "Natural C02-Rich Steam-Heated Waters in b e Broadlandst

Ohaaki Geothermal System, New Zealand: Their Chemistry, Distributioq and Corrosiv$

Nature," Trans., Geothermal Resources Council, Vol. 9, Part 2 (Aug. 19851, pp 245-250. '

Horne R N : "Geothermal Reservoir Engineering Course Notes," Japan Gedthermul Enera


Development Centre, Tokyo (1980). i1

Horner D R : "Pressure Build-Up in Wells," Proc. Third World Pet. Congress1 Section 11, E

Leiden, Netherlands (195 1). j

Hunt T M : "Recharge at Broadlands (Ohaaki) Geothermal Field 1967-1983 determined !ion

dI
Repeat Gravity Measurements," Proc., 6th N.Z. Geothermal Workshop 1984), pp 195

200.

Hunt T My Hicks S R : "Repeat Gravity Measurements at Broadlands Geothe

1974," Geophysics Division, D.S.I.R.,Wellington, New Zealand (1975),

Hurst W : "Interference Between Oilfields," Trans., AZME, 219 (1960), pp 175-c92.

Jacob C E : "Coefficients of Storage and Transmissibility Obtained from Pumpbng Tests in th

Houston District, Texas," Trans., Amer. Geophysical Union (1941), pp 74-756.

Jacob C E : "Correlation of Ground-water Levels and Precipitation on Long Island New York,

Trans., Amer. Geophysical Union, Part 2 (1944), pp 564-573.


- 131 -
Jargon J R : "Effect of Wellbore Storage and Well Damage at the Active Well on Interference

Test Analysis," J. Pet. Tech. (Aug. 1976), pp 851-858.

Kazemi H, Memll L S, Jargon J R : "Problems in Interpretation of Pressure FB11-off Tests in

Reservoirs With and Without Fluid Banks," J. Pet. Tech. (Sept. 1972), pp 1t147-1156.

Kazemi H, Seth M S, Thomas G W : "The Interpretation of Interference T e a in Naturally

Fractured Reservoirs With Uniform Fracture Distribution," SOC. Pet. Eng. 9. (Dec. 1969),

pp 463-472.

Kazmann R G : "Notes on Determining the Effective Distance to a Line of Redharge", Trans.,

Amer. Geophysical Union, Part 2 (Dec. 1946),p 854.

Khurana A K : "Influence of Tidal Phenomena on Interpretation of Pressure Buildup and Pulse'

Tests," Australian Pet. Exploration Assn. J., 16(1) (1976),pp 99-105. ,

Koefoed 0 : "Simultaneous Observation of Drawdown Data from Several ObServation Wells1

by Means of a Modification of the Walton Method," J. Hydrology, 21(4) (April 1974),ppi


,
I

3 15-319.

LaRocque G A Jr.: "Fluctuations of Water Levels in Wells in the Los Angeleb Basin During1
I
Five Strong Earthquakes 1933-1940,'' Trans., American Geophys. Union, Part 2 (1941),~

pp 374-386.

Leaver J D, Sageev A, Ramey H J Jr.: "Linear Boundary Detection In a Sin$le Interference,

Test," Proc., 7th NZ. Geothermal Workshop (1985).

Leaver J D, Sageev A, Ramey H J Jr.: "Multiple Well Interference Testing in the Ohaaki
I

Geothermal Field," SPE Regional Calif. Con), Oakland, April 2-4,paper 19122 (1986).

Martin J C : "Simplified Equations of Flow in Gas Drive Reservoirs and the meoretical Foun-

dation of Multiphase Pressure Buildup Analyses," Trans., AIME 216 (19598,pp 309.
I

- 132 -
McGuinness M J : "Interference Testing of the Rhyolites Overlying the Odaaki Geothermal

Reservoir," Trans., Geothermal Resources Council, Vol. 9, Part 2 (Augj 1985), pp 541-
!

545. I I

I
Mueller T D, Witherspoon P A : "Pressure Interference Effects Within Reservoirs 4
Aquifers," J. Pet. Tech. (April 1965), pp 471-474. I I

(1956), pp 482-509.

Ramey H J Jr.: "Rapid Methods for Estimating Reservoir

(April 1964), pp 447-454. ~

ii i
Ramey H J Jr.: "Approximate Solutions for Unsteady Liquid Flow in Composite Reservoirs,p
I I
J. Canadian Pet. Tech. (Jan. 1970), pp 32-37. ~

! I
I
I I
Ramey H J Jr., Kumar A, Gulati M S : Gas Well Test Analysis Under Water rive Condition ,
Amer. Gas Assn., Arlington, VA.( 1973). 4 5
Robinson T W : "Earth Tides Shown by Fluctuations of Water Levels in Welld in New Mexic

and Iowa," Trans., American Geophys. Union, Part 4 (1939), pp 656-666.


~

,
~
p I

Sageev A : "Interference Testing Near a Steam Cap," Trans., Geothermal Re


- 133 -
Vol. 9, Part 2 (Aug. 1985), pp 569-572.

Sageev A, Home R N, Ramey H J Jr.: "Detection of Linear Boundaries by Drawdown Tests;


I

A Semilog Curve Matching Approach," Water Resources Research.,Vol. 21, no.3 (19834
pp 305-310.

Sageev A, Leaver J D, Ramey H J Jr.: "The Significance of Early Time Data In Interferen
cei
Testing For Linear Boundary Detection," Proc. Zlth Stanford Geothemul Workshoq

(1986).
~

Stallman R W : "Nonequilibrium Type Curves Modified for Two-well System," Geologicaj

Survey Groundwater Note 3, U.S.Dept. of Int, Washington D.C. (1952).

Standing M B : Volumerric and Phase Behavior of Oil Field Hydrocarbon Systems, Reinhold1
,
Pub. Co., New York (1952). I

Strobe1 C J, Gulati M S, Ramey H J Jr.: "Reservoir Limit Tests in a Natdrally Fracturedl

Reservoir--A Field Case Study Using Type Curves," J. Pet. Tech. (Sept. ll976), pp 1097-1

1106.

Theis C V : "Relation Between the Lowering of the Piezometric Surface and the Rate and
I

Duration of Discharge of a Well Using Ground-Water Storage," Eos Truros. AGU, 16(2))

(1935), pp 519-524.

Thorsteinsson T, Eliasson J : "Geohydrology of the Laugarnes Hydrotherinal System in

Reykjavik, Iceland," U.N. Syrnpos. on the Devel. and Use of Geothermal Resources, Vol.

3 (1970), pp 1191-1204. I

van Everdingen A F, Hurst W : "The Application of the Laplace Transformatiod to Flow Prob-

lems in Porous Media," AIME,186 (1949), pp 305-324 I

Vela S : "Effect of a Linear Boundary on Interference and Pulse tests - The Elliptical inference^
_ _ _ _ _ _ _ ~

- 134 -
Area,” J . Pet. Tech. (Aug. 1977), pp 947-950.

Warren J E, Root P J : “The Behaviour of Naturally Fractured Reservoirs,” SPE J . (Sept.

1963), pp 245-255.

I
Young A : “Tidal Phenomena at Inland Boreholes Near Cradock,” Trans.,Ropal SOC.South
I

Africa, Vol. 3, Part 1 (1913), pp 61-105.


~

--------- : “Broadlands Geothermal Field Investigation Report,” Power Divisi4m Ministry 4


I
Works and Development, Wellington, New Zealand (1977). ~
I

- 135 -

APPENDIX A

LINE SOURCE TYPE CURVE

(Eipper 1985)

C
c This program calculates PD vs tD/rD**2 for the line source solution curve.
C
c Variables used:
c mmdei=imsl routine for exponential integral solution
c pd=dimensionless pressure
c td=dimensionless time divided by dimensionless radius squared
C
C
c These loops generate tD/rD**2 values between 0.1 and loo00 and calculate,
c corresponding PD values.
C
implicit real*8(a-h,o-z)
dimension td( lOOO),pd(lOOO)
double precision mmdei
iopt=1
n=O
do 10 i=1,6
do 20 j=1,20
tdlog=-2.+i+(j- 1)/20.
tdd=1O.* *tdlOg
if(tddgt10000.)go to 10
n=n+l
td(n)=tdd
arg=-1J(4.*td(n))
pd(n)=-O.5*mmdei(iopt,arg,ier)
20 continue
10 continue
C
c This loop outputs the values for plotting.
C
write(6,1000)n
1OOO format(i3)
do 30 i=l,n
write(6,2OOO)td(i),pd(i)
2000 format(e10.5,5x,el5.7)
30 continue
stop
end
- 136 -

APPENDIX B

SPECIFIC PRESSURE vs TIME DATA

TABLE B1

TEST B 1: BR13 RESPONSE TO BR23 DISCHARGE

BR13 Reswnse to BR23 Discharee--


Y

Specific Pressure
(k~a-slrn’)
0.28053e+02
0.87023e+02
0.12347e+03
0.14866e+03
0.16832e+03
0.18244e+03
0.1958&+03
0.20649e+03
0.217 18e+03
0.22805e+03
0.23588e+03
0.2437Oe+03
0.24924e+03
0.25763e+03
0.26622e+03
0.2717&+03
0.27729e+03
0.28302e+03
0.28855e+03
0.29427e+03
0.30267e+03
0.30878e+03
0.30935e+03
0.31546e+03
0.321 18e+03
0.32385e+03
0.32958e+03
0.33225e+03
0.33798e+03
0.34084e+03
0.34676e+03
0.350 19e+03
- 137 -
BR13 Response to BR23 Discharge-
Time Specific Pressure
(h) (~a-s/rn~)
0.325OOe+02 0.3563&+03
0.335OOe42 0.35859e+03
0.345OOe+02 0.36126e+03
0.355OOe+02 0.36603e+03
0.365OOe+02 0.36889e+03
0.375OOe42 0.3715&+03
0.385OOe+02 0.37 156e+03
0.395OOe+02 0.37729e+03
0.405OOe+02 0.37996e+03
0.4 15OOe+02 0.38569e+03
0.425OOe+02 0.38798e+03
0.435OOe+02 0.39294e+03
0.445OOe+02 0.3979Oe+03
0.455OOe+02 0.40019e+03
0.465OOe+02 0.40248e+03
0.475OOe+02 0.4049&+03
0.485oOe+02 0.40763e+03
0.495oOe+02 0.3965&+03
0.505OOe+02 0.3 3473e+03
0.5 15oOe+02 0.28989e+03
0.52500e+02 0.26374e+03
0.535OOe+02 0.24504e+03
0.545OOe+02 0.2299&+03
0.555OOe+02 0.21889e+03
0.565oOe+02 0.2 1508e+03
0.575OOe+02 0.20191e+03
0.585OOe+02 0.18893e+03
0.595OOe+02 0.18149e+03
0.605OOe+02 0.17233e+03
0.625OOe+02 0.1601 le+03
0.645OOe+02 0.1477 le+03
0.665&+02 0.14O4&+03
0.685OOe+02 O.l3282e+O3
0.705OOe+02 0.12233e+03
0.725OOe+02 0.11279e+03
0.745OOe+02 O.lO439e+03
0.775OOe+02 0.95992e+02
0.805OOe+02 0.89885e+02
0.835OOe+02 0.85305e+02
0.865OOe+02 0.76908e+02
0.895OOe+02 0.72328e+M
0.925OOe+02 0.66412e+02
0.955OOe+02 0.62023e+02
0.985OOe+02 0.55153e+02
0.1015Oe+03 0.51336e+02
0.1O45Oe+03 0.46756e+02
0.1075Oe+03 0.4351 le+M
0.1 105Oe+03 0.37977e+02
0.11350e+O3 0.35115e+02
0.1 165Oe+03 0.32443e+02
0.1 195Oe+03 0.29962e+02
- 138 -

Time Specific Pressure


(h) (k~a-s~rn’)
0.1225Oei-03 0.23282e+02
0.1255Oei-03 0.17557e+02
0.1285Oe+03 0.17557e+02
0.1315Oei-03 0.1584Oe+O2
0.1345Oe43 0.1355Oe+O2
0.1375Oei-03 0.78244e+Ol
0.1425Oei-03 0.1145Oe+01
- 139 -
TABLE B2
TEST B2: BR23 RESPONSE TO BR13 DISCHARGE

BR23 Response to BR13 Discharge-


Specific Pressure
(kPa-s/m3)
0.67500e+Ol
0.1443Oe42
0.3848Oe42
0.57990e42
0.787&+02
0.9404Oe+02
0.10478e43
0.12244e43
0.13221e43
0.14379e43
0.15346e+03
0.16326e+03
0.171 14e+03
0.17954e+03
0.1882 l e 4 3
0.1987Oe+03
0.2055&+03
0.21229e+03
0.21825e+03
0.22421e+03
0.23017e43
0.23613e43
0.24026e+03
0.2454 3e+03
0.25336e+03
0.25932e+03
0.26732e43
0.27741e43
0.29014e43
0.2963Oe+03
0.30640e+03
0.31255e+03
0.31869e+03
0.32879ei-03
0.33679e43
0.34332e43
0.35277e+03
0.36064ei.03
0.36545e43
0.37278e43
0.37998e+03
0.38731e+0?
0.39624ei-02
0.40317e+0?
0.40704e+0?
- 140 -

BR23 Response 1 3R13 Discharge


Time Specific Pressure
(h) (k~a-slm~)
0.66OOOe+02 0.41678e+03
0.68OOOe+02 0.42185e.r-03
0.70000e+02 0.4261 l e 4 3
0.72OOOe+02 0.4269Oei-03
0.74OOOe+O2 0.42929ei-03
0.76OOOe+02 0.43635ei-03
0.78OOOe+02 0.4382Oe+03
0.80000ei-02 0.449 15e+03
0.82OOOe+02 0.45408e+03
0.86OOOe+02 0.46325e43
0.90000e+02 0.47525ei-03
0.94OOOe+02 0.48012ei-03
0.98OOOe+02 0.48768e+03
0.102OOe+03 0.49886e+03
0.106OOe+03 0.51275e+03
0.1 lOOOe+03 0.52003e+03
0.1 14OOe+03 0.52974e+03
0.118OOe+03 0.53553e+03
0.122oOe+03 0.54294e+03
0.126OOe+O3 0.55386ei-03
0.13OOOe+03 0.56357e+03
0.134We+03 0.57 11le+03
0.13800e+03 0.58001e+03
0.14200e+03 0.5 8687e+03
0.146OOe+03 0.59549e+03
0.15OOOe+03 0.60385e+03
0.166OOe+03 0.63523e+03
0.178OOe+03 0.66016e43
0.19OOOe+03 0.681 14e+03
0.202OOe+03 0.7028 le+03
0.214OOe+03 0.7253Oe+03
0.226OOe+03 0.74628e+03
0.238OOe+03 0.7 6755ei-03
0.25OOOe+03 0.7873Oe+03
0.262OOei-03 0.8323Oe+03
0.274OOe+03 0.85 148e+03
0.286OOe+03 0.86869ei-03
0.298&+03 0.88872e+03
0.3 10OOe+03 0.90705e+03
0.332&+03 0.92623e+03
0.33 1oOe+03 0.94048e43
0.334OOe+03 0.10543e+04
0.34600e+03 0.10589e+04
0.358OOe+03 0.10663e+04
0.382OOe+03 0.1 1512e+04
- 141 -
TABLE B3
TEST B3: BR23 RESPONSE TO BR19 DISCHARGE

BR23 Response to BR19 Discharge-

Time Specific Pressure


(h) (wa-s/rn3)
0.10000e+01 0.34856e+01
0.20000e+01 0.10577e42
0.30000e+01 0.15865e42
0.40000e+01 0.2476Oe+02
0.50000e+01 0.29567e+02
0.60000e+01 0.36538e+02
0.70000e+01 0.4363Oe+02
0.80000e+01 0.47236e+02
0.90000e+01 0.52524e+02
0.10000e+02 0.59615e42
0.1 1OOOe+02 0.64904e+02
0.12OOOe+02 0.68389e42
O.l3OOOe+O2 0.7 3678e+02
0.14OOOe+02 0.77284e+02
O.l50Ooe+O2 0.8 149Oe+02
0.1 6oooe+02 0.86779e+02
0.17OOOe+02 0.89784e+02
0.18OOOe+02 0.92909e+02
0.190OOe+02 0.95 192e+02
0.2m+02 0.96875e+02
0.2 1Oooe+O2 0.99639e+02
0.22000e+02 0.10144e+03
0.23OOOe+02 0.10553e4I3
0.24000e+02 0.10721e4I3
0.25OOOe+02 0.11082e+O3
0.26000e+02 0.1 143Oe43
0.27OOOe+02 0.1 1563e43
0.28OOOe+02 0.1 191le+03
0.29OOOe+02 0.12392e43
0.30000e+02 0.1274Oe+03
0.31ooOe+O2 0.12921e+03
0.32OOOe+02 O.l3269e+03
0.33OOOe+02 0.1363Oe43
0.34OOOe+02 0.138 1Oe-i-03
0.35OOOe+02 0.14159e+03
0.36oooe+02 0.14507e43
0.37OOOe+02 0.14688e43
0.37OOOe+02 0.14688e+03
0.38OOOe+02 0.14844e43
0.39OOOe+02 0.15349e43
0.40000e+02 0.15493e+03
0.4 1OOOe+02 0.15673e+03
0.42OOOe+02 0.16022e+03
0.43OOOe+02 0.16202e+03
0.44OOOe+02 0.16755e+03
- 142 -

BR23 Response to BR19 Discharge-


Time Specific Pressure
(h) (~a-slrn’)
0.45OOOe+02 0.17344e+03
0.46OOOe+02 0.17536e+03
0.47OOOe+02 0.17764e43
0.48OOOe+02 0.17969e43
0.49OOOe+02 0.18173e+03
0.50000e+02 0.18401e43
0.5 1OOOe+02 0.1852243
0.52OOOe+02 0.19026e+03
0.55OOOe+02 0.19303e+03
0.57OOOe+02 0.19639e+03
0.59OOOe+02 0.2021&+03
0.61OOOe+02 0.20625e+03
0.63OOOe+02 0.20877e43
0.65OOOe+02 0.2 1298e+03
0.67OOOe+02 0.21827e+03
0.69OOOe+02 0.21827e+03
0.7 1OOOe+02 0.22055e43
0.73OOOe+02 0.22236e+03
0.75OOOe+02 0.2253&+03
0.77OOOe+02 0.22813e+03
0.79OOOe+02 0.23317e+03
0.81OOOe+02 0.23498ei-03
0.89OOOe+02 0.24736e+03
0.93OOOe+02 0.25132m-03
0.97OOOe+02 0.25656e+03
0.101OOe+03 0.2650&+03
0.105OOe+03 0.26851e+03
0.109OOe+03 0.27308e+03
0.1 13OOe+03 0.28017e+03
0.117OOe+03 0.28245ei-03
0.121OOe+03 0.28678e43
0.12500e+03 0.29435e+03
0.129OOe+03 0.29784e+03
0.137OOe+03 0.30925e+03
0.189OOe+03 0.33425e+03
- 143 -
TABLE B4
TEST B4: BR23 RESPONSE TO BR13 INJECTION

BR23 Response to BR13 Discharge-


Time Specific Pressure
(h) (k~a-slrn~)
0.10000e+01 0.14878e+02
0.2m+01 0.51789e+02
0.30000e+01 0.6742Oe+02
0.40000e+01 0.8305 l e 4 2
0.50000e+Ol 0.10979e+03
0.60000e+01 0.12825e+03
0.70000e+01 0.1467Oe+03
0.80000e+01 0.1642%+03
0.90000e+01 0.1791Oe+03
0.10000e+02 0.19379e+03
0.1 1OOOe+02 0.20866e+03
O.l2OOOe+O2 0.22147e+03
0.13OOOe+02 0.23352~03
0.14OOOe+02 0.24557e+03
0.15OoOe+02 0.2548Oe+03
0.16OOoe+02 0.26478e+03
O.l7OOOe+02 0.27213e+03
0.18OOOe+02 0.28493e+03
0.19OOOe+02 0.29736e+03
0.20000e+02 0.30377e+03
0.2 1OOOe+02 0.31299e+03
0.22oooe+02 0.32185e+03
0.23OOOe+02 0.32467e+03
0.24OOOe+02 0.33672e+03
0.25OOOe+02 0.35 14l e 4 3
0.26OOoe+02 0.36347e+03
0.27OOOe+02 0.37175e+03
0.28OOOe+02 0.378 15e+03
0.29OOOe+02 0.38832e+03
0.30000e+02 0.39473e+03
0.31OOOe+02 0.40113e+03
0.320oOe+02 0.40753e+03
0.33OOOe+02 0.4 143l e 4 3
0.34OOOe+02 0.42072e+03
0.35OOOe+02 0.42712e+03
0.36OOOe+02 0.43635e+03
0.37OOOe+02 0.43992e+03
0.38OOOe+02 0.44633e+03
0.39OOOe+02 0.45273e+03
0.40000e+02 0.4563 le+03
0.4 1OOOe+02 0.46271e+03
0.42OOOe+02 0.46667e+03
0.43OOOe+02 0.47382e+03
0.44OOOe+02 0.4806oe+03
0.45OOOe+02 0.48418e+03
- 144 -

Time Specific Pressure


(h) (wa-s/rn3)
0.46OoOe42 0.48776ei-03
0.47OOOe+02 0.49134e43
0.69OOOe42 0.56798e43
0.10000e43 0.64821e 4 3
0.106OOe+03 0.66252e43
0.118OOe+O3 0.68983e43
0.13OOOe+03 0.7177Oe43
0.142OOe43 0.77985e+03
0.154OOe43 0.79699e43
0.178oOe+03 0.81525e43
O.l92OOe+O3 0.83503e43
0.216OOe+03 0.85669e+03
0.24OOOe+O3 0.87702e43
0264OOe43 0.89134e+03
0.402OOe43 0.93089e43
- 145 -
TABLE B5
TEST B6: BR23 RESPONSE TO BR20 DISCHARGE

BR23 ResDonse to BR20 Discharge


1

Specific Pressure
(k~a-slrn~)
0.27559e+01
0.73228e41
0.12913e+02
0.2 118l e 4 2
0.20079e+02
0.18268e42
O.l6299e+02
0.15827e+02
0.14961e42
0.14331e+02
0.14252e+02
0.1378&+02
0.1338&+02
0.13 15&+02
0.12913e+02
0.1244 le+02
0.1 1654e+02
O.l1024e+02
0.10079e42
- 146 -
TABLE B6
TEST B7: BR34 RESPONSE TO BR31 DISCHARGE

-
BR34 Resmnse to BR31 Discharge
-1

Time Specific Pressure


(h) (kpa-s/rn3)
0.60000e+01 0.25OOOe+00
0.70000e+01 0.80000e+00
0.75OOOe+01 0.135OOe+01
0.175OOe+02 0.3 17OOe+01
0.245OOe+02 0.544OOe+01
0.305OOe+02 0.77 lOOe+Ol
0.335OOe42 0.10000e+02
0.41OOOe+02 0.1395Oe42
0.49o00e+02 0.2019Oe+02
0.6 1OOOe+02 0.2847Oe+02
0.73OOOe+02 0.4268Oe+02
0.85OOOe+02 0.58 18Oe+02
0.97OOOe+02 0.7023Oe42
0.109OOe+03 0.7905Oe+02
0.121OOe+03 0.8923Oe+02
0.133OOe+03 0.9782Oe+02
0.145OOe+03 0.10992e+03
0.157OOe+03 0.1 195Oe+03
0.169OOe+03 0.12843e+03
0.1 81OOe+03 0.1379 le+03
0.205OOe+03 0.15938e43
0.229OOe+03 0.17869e+03
0.253OOe+03 0.19884e+03
0.277OOe+03 0.2 1722e+03
0.293OOe+03 0.2275Oe+03
- 147 -
TABLE B7
TEST B8: BR34 RESPONSE TO BR23 DISCHARGE

BR34 ResDonse to BR23 Discharge


Time Specific Pressure
(h) (k~a-strn~)
0.19OOOei-02 0.452OOei-01
0.31OOOe+02 0.1631Oei-02
0.43OOOei-02 0.3539Oei-02
0.55OOOei-02 0.5284Oei-02
0.67OOOei-02 0.7224Oei-02
0.79OOOe+02 0.9293Oei-02
0.91OOOei-02 0.1137l e 4 3
0.103OOe+03 0.1307 l e 4 3
0.115OOe+O3 0.14873ei-03
0.127OOe+03 0.16591ei-03
0.139OOe+03 0.18527ei-03
0.151oOe+03 0.20056ei-03
0.163OOe+03 0.21139ei-03
0.185OOei-03 0.22544ei-03
0.197OOei-03 0.23996ei-03
0.209oOe+03 0.25244e+03
0.221OOe+03 0.26654ei-03
0.223OOe+03 0.26898ei-03
- 148 -
TABLE B8
TEST B9: BR34 RESPONSE TO BR19 DISCHARGE

BR34 Response to BR19 Discharge-


Time Specific Pressure
(h) (k~a-strn~)
0.10000e+01 0.15487e+02
0.30000e+Ol 0.70099e+02
0.50000e+01 0.92554e42
0.70000e+01 0.11047e43
0.90000e+01 0.1 lOOle+03
0.1 1OOOe+02 0.1321 3e+03
0.13OOOe+02 0.1515le+03
0.17OOOe+02 0.18202e+03
0.19OOOe+02 0.2308&+03
031oOoe+02 0.26784e+03
0.23OOOe+02 0.3246243
0.25OOOe+02 0.40577e+03
0.27oOoe+02 0.49367e43
0.29oOoe+02 0.59208e+03
0.3 1OOOe+02 0.63445e43
0.33OOOe+02 0.69849e43
0.35OOOe+02 0.74465e43
0.37OOOe+02 0.79015e+03
0.41OOOe+02 0.8468 l e 4 3
0.45OOOe+02 0.9533&+03
0.49OOOe+02 0. 1O936e+O4
0.53ooOe+02 0.12374e+04
0.57OOOe+02 0.1401Oe+04
0.61OOOe+02 0.14884e+04
0.65OOOe+02 0.15358e4
0.69OOOe+02 0.16 142e+W
0.73OOOe+02 0.17259e+04
0.79OOOe+02 0.1947 6e+04
0.85OOOe+02 0.21123e+04
0.9 1OOOe+02 0.21975e+04
0.97OOOe+02 0.23399e+04
0.103OOe+03 0.25 1 5 7 e 4
0.109OOe+03 0.26475e+04
0.115OOe+03 0.2749&+04
0.13 1OOe+03 0.31562e+04
0.133OOe+03 0.3 188l e d
0.135OOe+03 0.31929e+04
0.143OOe+03 0.32198e+04
0.151OOe+03 0.3 1626e+04
0.159OOe+03 0.30778e+W
0.167OOe+03 0.29246e+04
0.175OOe+03 0.27724em
0.183OOe+03 0.268 17e+04
0.19 1OOe+03 0.25317e+04
0.203OOe+03 0.2400&+04
- 149 -
BR34 Response to BR19 Discharge
Time Specific Pressure
(h) (~a-slmj)
0.215OOei-03 0.21762e+04
0.227OOei-03 0.20773e+04
0.239OOei-03 O.l8514e+04
0.25 1OOei-03 0.17 165e+04
0.26300ei-03 O.l552Oe+04
0.275OOei-03 0.14784ei-04
0.287OOei-03 0.13522e+04
0.299OOei-03 0.13235eM
0.3 1100ei-03 0.12613e+04
0.32800ei-03 0.1 1515e+04
0.335OOei-03 0.10475e+04
0.347OOei-03 0.94134ei-03
0.359OOei-03 0.78263e+03
0.37 1oOei-03 0.73768ei-03
0.383OOei-03 0.53281ei-03
0.395OOei-03 0.44185ei-03
0.407OOei-03 0.32771ei-03
0.43 1OOei-03 0.29385ei-03
0.455OOei-03 0.23693e43
- 150-
TABLE B9
TEST B10: BR23 RESPONSE TO BR13 SHUT-IN

BR23 Response to BR13 Shut-in


Time Specific Pressure
(h) (k~a-s/cu.m~)
0.10000e+01 O.lWe+02
0.20000e+01 0.21495e+02
0.30000e+01 0.36604e+02
0.40000e+01 0.58567e+02
0.50000e41 0.77882e42
0.60000e+01 0.95327e+02
0.70000e+01 0.10763e+03
0.80000e+01 0.1168%+03
0.90000e+01 0.12679e+03
0.10000e+02 0.13738e+03
0.1 lOOOe+02 0.14657e+03
0.12o0oe+02 0.1557&+03
0.13OOOe+02 O.l6433e+03
0.14OOOe+02 0.16947e+03
0.15OOOe+02 0.17586e+03
0.16OOoe+02 0.18178e+03
0.17OOOe+02 0.18769e+03
0.18OOOe+02 0.19299e+03
O.l9OoOe+O2 0.19829e+03
0.20000e+02 0.20187e+03
0.2 1OOOe+02 0.20514e+03
0.22OOOe+02 0.20888e+03
0.23OOOe+02 0.21464e+03
0.24OOOe+02 0.22103e+03
0.25OOOe+02 0.2271Oe+03
0.26OOoe+02 0.23131e+03
0.27OoOe+02 0.23629e+03
0.280oOe+02 0.24299e+03
0.29OOOe+02 0.24673e43
0.30000e+02 0.25093e+03
0.3 1OOOe+02 025607e+03
0.32000e+02 0.26 121e+03
0.33OOOe+02 0.26526e43
0.34OOOe+02 0.26838e+03
0.35OoOe+02 0.27243e+03
0.36OOoe+02 0.27383e43
0.37OOOe+02 0.28209e+03
0.380oOe+02 0.28474e43
0.39OOOe+02 0.2891Oe+03
0.40000e+02 0.29174e+03
0.42OOOe+02 0.29688e43
0.46OOoe+02 0.30639e43
0.48OOOe+02 0.3 1324e+03
0.50000e+02 0.3 182%+03
0.52OOOe+02 0.32243e+03
1 I /

- 151 -
BR23 Reswnse to BR13 Shut-in
Time Specific Pressure
(h) (k~a-s/cu.rn~)
0.54OOOe+02 0.32695e+03
0.56OOOe+02 0.33162e+03
0.58OOOe+02 0.33598e43
0.60000e+02 0.33801e43
0.62OOOe+02 0.34174e+03
0.64OOOe+02 0.34735e43
0.66OOOe+02 0.35125e43
0.68OOOe+02 0.35639e+03
0.70000e+02 0.36059e+03
0.72OOOe+02 0.36698e43
0,74OoOe+02 0.36994e+03
0.76OOOe+02 0.37U3e43
0.78OOOe+02 0.37492e43
0.80000e+02 0.37804e43
0.84OOOe+02 0.38551eA3
0.88OOOe+02 0.3905Oe+03
0.92OOOe+02 0.39533e43
0.96OOOe+02 0.39984e+03
0.10000e+03 0.4015&+03
0.1O4OOe+03 0.40187e43
0.108OOe+03 0.40888e+03
0.1 12OOe+03 0.4 1402e+03
0.1 16OOe+03 0.41542e+03
0.124OOe+03 0.41916~03
0.14OOOe+03 0.43 1oOe+03
0.148OOe+03 0.43941e43
0.156OOe+03 0.44377e43
0.164OOe+03 0.44579e43
0.172OOe+03 0.45062e+03
0.18OOOe+03 0.4553Oe+03
0.188OOei-03 0.4542le+03
0.196OOe+03 0.46121e43
0.204OOe+03 0.4595OeA3
0.2 12OOe+03 0.45701e+03
0.22OOOe+03 0.46308e+03
0.228OOe+03 0.45981e43
0.234OOe+03 0.46199ei-03
- 152 -
TABLE B10
TEST C1: BR13 RESPONSE TO BR20 DISCHARGE

-
BR13 Response to BR20 Discharge
Specific Pressure
(kpa-s/rn3)
0.154768e+01
0.154768e+01
0.154768e+01
0.14881Oe+02
0.28214&+02
0.282142e+02
0.348809e+02
0.482143e+02
0.54881Oe+02
0.615476e+02
0.682 142e+02
0.748809e+02
0.234762e+03
0.271429e+03
0.274762e+03
0.274762e+03
0.28 1429e+03
0.28 1429e+03
0.281429e+03
0.281429e+03
0.288095e+03
0.288095e+03
0.288095e+03
0.288095e+03
0.288095e+03
0.294762e+03
0.294762e+03
0.298095e+03
0.301429e+03
0.294762e+03
0.301429e+03
0.308095e+03
0.308095e+03
0.314762e+03
0.3 14762e+03
0.3 14762e+03
0.321429e+03
0.321429e+03
0.321429e+03
0.328095e+03
0.334762e+03
0.334762e+03
0.334762e+03
0.334762e+03
0.334762e+03
- 153 -
BR13 Response to BR20 Discharge
-
Time specific Pressure
(k~a-slrn~)
0.341428e+03
0.341428e+03
0.348095e+03
0.348095e+03
0.348095e+03
0.351429e+03
0.35476&+03
0.35476&+03
0.354762e+03
0.35476&+03
0.35476&+03
0.354762e+03
0.361429e+03
0.35476%+03
0.354762e+03
0.361429e+03
0.368095e+03
0.368095e+03
0.368095e+03
0.368095e+03
0.368095e+03
0.37476&+03
0.368095e+03
0.374762e+03
0.378095e+03
0.37476%+03
0.374762e+03
0.381428e+03
0.381428e+03
0.381428e+03
0.381428e+03
0.381428e+03
0.388095e+03
0.388095e+03
0.388095e+03
0.388095e+03
0.388095e+03
0.388095e+03
0.39476&+03
0.394762e+03
0.394762e+03
0.401429e+03
0.401429e+03
0.401429e+03
0.404762e+03
0.408095e+03
0.408095e+03
0.408095e+03
0.408095e+03
0.408095e+03
0.408095e+03
-154-

-
BR13 Resmnse to BR20 Discharee
Specific Pressure
(~a-~lrn~)
0.414762e+03
0.4 14762e+03
0.414762e+03
0.41476&+03
0.4 14762e+03
0.4 1476&+03
0.4 14762e+03
0.421428e+03
0.421428e+03
0.421428e+03
0.461428e+03
0.468095e+03
0.468095e+03
0.448095e+03
0.474762e+03
0.474762e+03
0.474762e+03
0.474762e+03
0.474762e+03
0.47476&+03
0.47476&+03
0.481429e+03
0.478095e+03
0.481429e+03
0.481429e+03
0.481429e+03
0.481429e+03
0.48 1429e+03
0.481429et03
0.488095e+03
0.481429e+03
0.488095e+03
0.488095e+03
0.484762e+03
0.488095e+03
0.488095e+03
0.488095e+03
0.488095e+03
0.488095e+03
0.488095e+03
0.494762e+03
0.494762e+03
0.494762e+03
0.491429e+03
0.494762e+03
0.494762e+03
0.494762e+03
0.494762e+O?
0.494762e+O?
0.488095e+O?
0.494762e+O?
- 155 -
BR13 Reswnse to BR20 Dischame
1

Time specific Pressure


(h) (k~a-slrn~)
0.2355OOe+03 0.494762e+03
0.2365OOe+03 0.501429e+03
0.2375OOe+03 0.501429e+03
0.2385OOe+03 0.50 1429e+03
0.2395OOe+03 0.501429e+03
0.2405OOe+03 0.508095e+03
0.U 15OOe+03 0.508095e+03
0.2425OOe+03 0.508095e+03
O.%435OOe+O3 0.508095e+03
0.2445OOe+03 0.508095e+03
O.U55OOe+O3 0.508095e+03
0.2465OOe+03 0.5080!35e+03
0.2475OOe+03 0.5 1476&+03
0.2485OOe+03 0.51 1429e+03
0.2495OOe+03 0.5 14762e+03
0.2505OOe+03 0.5 14762e+03
0.2515OOe+03 0.508095e+03
0.2525OOe+03 0.5 14762e+03
0.2535OOe+03 0.5 14762e+03
0.2545OOe+03 0.514762e+03
0.2555OOe+03 0.518095e+03
0.2565OOe+03 0.514762e+03
0.2585OOe+03 0.514762e+03
0.2595OOe+03 0.514762e+03
0.2605OOe+03 0.52 1429e+03
0.2615OOe+03 0.521429e+03
0.2625OOe+03 0.528095e+03
0.2635OOe+03 0.528095e+03
0.2645OOe+03 0.528095e+03
0.2655OOe+03 0.528095e+03
0.2665OOe+03 0.534762e+03
0.2675OOe+03 0.528095e+03
0.2685OOe+03 0.534762e+03
0.2695OOe+03 0.528095e+03
0.2705OOe+03 0.53476&+03
0.27 15OOe+03 0.534762e+03
0.2725OOe+03 0.534762e+03
0.2735OOe+03 0.534762e+03
0.2745OOe+03 0.541429e+03
0.2755OOe+03 0.538095e+03
0.2765OOe+03 0.53476k+03
0.2775OOe+03 0.54 1429e+03
0.2785OOe+03 0.534762e+03
0.2795OOe+03 0.54 1429e+03
0.2805OOe+03 0.548095e+03
0.28 15OOe+03 0.544762e+03
0.2825OOe+03 0.541429e+03
0.2835OOe+03 0.548095e+03
0.284 5OOe+03 0.544762e+03
0.2855OOe+03 0.548095e+03
0.2865OOe+03 0.554881e+03
- 156 -
BR13 Resmnse to BR20 Discharge
1

Time Specific Pressure


(h) (k~a-s~rn~)
0.2875OOe+03 0.554881e+03
0.2885OOe+03 0.554881e+03
0.2895OOe+03 0.554881e+03
0.2905OOe+03 0.558214e+03
0.29 15OOe+03 0.55488 le+03
0.2925OOe+03 0.55488le+03
0.2935OOe+03 0.561548e+03
0.2945OOe+03 0.55488 le+03
0.2955OOe+03 0.561548e+03
0.2965OOe+03 0.561548e+03
0.2975OOe+03 0.561548e+03
0.2985OOe+03 0.561548e+03
0.2995OOe+03 0.561548e+03
0.3005OOe+03 0.561548e+03
0.3015OOe+03 0.5682 14e+03
0.3025OOe+03 0.561548e+03
0.3035OOe+03 0.561548e+03
0.3045OOe+03 0.568214e+03
0.3055OOe+03 0.5682 14e+03
0.3065OOe+03 0.568214e+03
0.3075OOe+03 0.5682 14e+03
0.3085OOe+03 0.568214e+03
0.3095OOe+03 0.574881e+03
0.3105OOe+03 0.568214e+03
0.3 115OOe+03 0.574881e+03
0.3125OOe+03 0.58 1548e+03
0.3135OOe+03 0.57488 le+03
0.3145OOe+03 0.57488 le+03
0.3155OOe+03 0.574881e+03
0.3 165OOe+03 0.57488 le+03
0.3175OOe+03 0.57488 le+03
0.3 185OOe+03 0.574881e+03
0.3195OOe+03 0.57488 le+03
0.3205OOe+03 0.574881e+03
0.32 15OOe+03 0.57488 le+03
0.3225OOe+03 0.57488 le+03
0.3235OOe+03 0.57488 le+03
0.3245OOe+03 0.574881e+03
0.3255OOe+03 0.581548e+03
0.3265OOe+03 0.57488 le+03
0.3275OOe+03 0.58 1548e+03
0.3285OOe+03 0.57488 le+03
0.3295OOe+03 0.57488 le+03
0.3305OOe+03 0.574881e+03
0.3315OOe+03 0.574881e+03
0.3325OOe+03 0.581548e+03
0.3335OOe+03 0.588214e+03
0.3345OOe+03 0.588214e+03
0.3355OOe+03 0.584881e+03
0.3365OOe+03 0.588214e+03
0.3375OOe+03 0.588214e+03
I
1. 1 1

- 157 -

Time specific pressure


(h) (k~a-strn~)
0.3385OOe+03 0.588214e+03
0.3395OOe+03 0.588214e+03
0.3405OOe+03 0.588214e+03
0.34 15OOe+03 0.578214e+03
0.3425OOe+03 0.568214e+03
0.3435OOe+03 0.561548e+03
0.3445OOe+03 0.554881e+03
0.3455OOe+03 0.541429e+03
0.3465OOe+03 0.534762e+03
0.3475OOe+03 0.528095e+03
0.3485OOe+03 0.5 14762e+03
0.3495OOe+03 0.508095e+03
0.3505OOe+03 0.501429e+03
0.35 15OOe+03 0.49476%+03
0.3525OOe+03 0.488095e+03
0.3535OOe+03 0.484762e+03
0.3545&+03 0.474762e+03
0.3555OOe+03 0.474762e+03
0.3565&+03 0.474762e+03
0.3575&+03 0.474762e+03
0.3585OOe+03 0.468095e+03
0.3595OOe+03 0.468095e+03
0.3605OOe+03 0.461428e+03
0.36 15OOe+03 0.461428e+03
0.3625OOe+03 0.454762e+03
0.3635OOe+03 0.454762e+03
0.3645OOe+03 0.448095e+03
0.3655OOe+03 0.448095e+03
0.3665OOe+03 0.44 1429e+03
0.3675OOe+03 0.434762e+03
0.3685OOe+03 0.4347 62e+03
0.3695OOe+03 0.428095e+03
0.3705OOe+03 0.42 1428e+03
0.37 15OOe+03 0.421428e+03
0.3725OOe+03 0.421428e+03
0.37350Oec03 0.414762e+03
0.3745OOe+03 0.4 11429e+03
0.3755OOe+03 0.408095e+03
0.376500e+03 0.401429e+03
0.3775OOe+03 0.401429e+03
0.3785OOe+03 0.401429e+03
0.3795OOe+03 0.394762e+03
0.4075OOe+03 0.348095e+03
0.4085OOe+03 0.348095e+03
0.4095OOe+03 0.34 1428e+03
0.4 105OOe+03 0.341428e+03
0.41 15OOe+03 0.341428e+03
0.4 125OOe+03 0.341428e+03
0.4135OOe+03 0.341428e+03
0.4145OOe+03 0.334762e+03
0.4 155oOe+03 0.334762e+03
- 158 -

-
BR13 Response to BR20 Discharge
Time Specac Pressure
(h) (k~a-strn~)
0.4165OOe+03 0.334762e+03
0.4 175OOe+03 0.328095e+03
0.4185OOe+03 0.328095e+03
0.4195OOe+03 0.328095e+03
0.4205OOe+03 0.321429e+03
0.4215OOe+03 0.328095e+03
0.4225OOe+03 0.321429e+03
0.4235OOe+03 0.321429e+03
0.4245OOe+03 0.314762+03
0.4255OOe+03 0.314762e+03
0.4265OOe+03 0.3 14762e+03
0.4275OOe+03 0.314762e+03
0.4285OOe+03 0.314762e+03
0.4295OOe+03 0.3 14762e+03
0.430500e+03 0.314762e+03
0.43 15OOe+03 0.314762e+03
0.4325OOe+03 0.3 14762e+03
0.4335OOe+03 0.314762e+03
0.4345OOe+03 0.308095e+03
0.4355OOe+03 0.308095e+03
0.4365OOe+03 0.308095e+03
0.4375OOe+03 0.308095e+03
0.4385OOe+03 0.301429e+03
0.4395OOe+03 0.30 1429e+03
0.4405OOe+03 0.308095e+03
0.4415OOe+03 0.308095e+03
0.4425OOe+03 0.308095e+03
0.4435OOe+03 0.301429e+03
0.4445OOe+03 0.301429e+03
0.4455OOe+03 0.301429e+03
0.4465OOe+03 0.301429e+03
0.4475OOe+03 0.301429e+03
0.4485OOe+03 0.301429e+03
0.4495OOe+03 0.294762e+03
0.4505OOe+03 0.294762e+03
0.45 15OOe+03 0.294762e+03
0.4525OOe+03 0.288095e+03
0.4535OOe+03 0.288095e+03
0.4545OOe+03 0.294762e+03
0.4555OOe+03 0.294762e+03
0.4565OOe+03 0.294762e+03
0.4575OOe+03 0.288095e+03
0.4585&+03 0.288095e+03
0.4605OOe+03 0.2947 62e+03
0.46 15OOe+03 0.288095e+03
0.4625OOe+03 0.288095e+03
0.4635OOe+03 0.288095e+03
0.4645OOe+03 0.288095e+03
0.4655OOe+03 0.288095e+03
0.4665OOe+03 0.284762e+03
0.4675OOe+03 0.281429e+03
- 159 -

BR13 Response to I 20 Discharge


Time Specific Pressure
(h) (k~a-s~rn~)
0.4685OOe+03 0.281429e+03
0.4695OOe+03 0.274762e+03
0.4705OOe+03 0.274762+03
0.47 15OOe+03 0.274762e+03
0.4725OOe+03 0.271429e+03
0.4735OOe+03 0.268095e+03
0.4745OOe+03 0.268095e+03
0.4755OOe+03 0.268095e+03
0.4765OOe+03 0.268095e+03
0.4775OOe+03 0.274762+03
0.4785OOe+03 0.27476&+03
0.4795OOe+03 0.274762e+03
0.4805OOe+03 0.274762e+03
0.48 15OOe+03 0.274762e+03
0.4825OOe+03 0.274762+03
0.4835OOe+03 0.268095e+03
0.4845OOe+03 0.268095e+03
0.4855OOe+03 0.268095e+03
0.4865OOe+03 0.261429e+03
0.4875OOe+03 0.261429e+03
0.4885OOe+03 0.254762e+03
0.4895OOe+03 0.254762e+03
0.4905OOe+03 0.254762e+03
0.49 15OOe+03 0.254762e+03
0.4925OOe+03 0.254762e+03
0.4935OOe+03 0.254762e+03
0.4945OOe+03 0.254762e+03
0.4955OOe+03 0.254762e+03
0.4965OOe+03 0.248095e+03
0.4975OOe+03 0.25 1429e+03
0.4985OOe+03 0.248095e+03
0.4995OOe+03 0.241429e+03
0.5005OOe+03 0.248095e+03
0.5015OOe+03 0.248095e+03
0.5025OOe+03 0.254762e+03
0.5035OOe+03 0.254762e+03
0.5045OOe+03 0.2547 62e+03
0.5055OOe+03 0.251429e+03
0.5065OOe+03 0.254762e+03
0.5075OOe+03 0.25476%+03
0.5085OOe+03 0.248095e+03
0.5095OOe+03 0.251429e+03
0.5105OOe+03 0.248095e+03
0.5115OOe+03 0.248095e+03
0.5125OOe+03 0.248095e+03
0.5135OOe+03 0.248095e+03
0.5145OOe+03 0.241429e+03
0.5155OOe+03 0.241429e+03
0.5165OOe+03 0.24 1429e+03
0.5 175OOe+03 0.234762e+03
0.5185OOe+03 0.241429e+03
-160-

BR13 Response to BR20 Discharge


Specific Pressure
(k~a-s~m’)
0.241429e+03
0.241429e+03
0.241429e+03
0.23476%+03
0.23476%+03
0.23476&+03
0.23476%+03
0.241429e+03
0.241429e+03
0.241429e+03
0.241429e+03
0.241429e+03
0.244762e+03
0.241429e+03
0.244762e+03
0.241429e+03
0.241429e+03
0.241429e+03
0.234762e+03
0.234762e+03
0.234762e+03
0.241429e+03
0.234762e+03
0.234762e+03
0.23476%+03
0.23476%+03
0.221429e+03
0.228095e+03
O.l85OOoe+03
0.183095e+03
0.18119Oe+03
0.18119Oe+03
0.17738le+03
0.173452e+03
0.17738le+03
0.177381e+03
0.17345%+03
0.173452e+03
0.17345%+03
0.173452e+03
0.17345%+03
O.l73452e+03
0.173452e+03
0.177381e+03
0.17345%+03
0.173452e+03
0.173452e+03
0.171548e+03
0.17345%+03
0.17345%+03
0.177381e+03
- 161 -
BR13 Response to BR20 Discharge-
Time specific pressure
(h) (~a-slrn~)
0.9115OOe+03 0.177381e+03
0.9125OOe+03 0.18119Oe+03
0.9135OOe+03 O.l83095e+03
0.9145OOe+03 0.18119Oe+03
0.9155OOe+03 0.185OOOe+03
0.9165OOe+03 0.185OOOe+03
0.9175OOe+03 0.185OOOe+03
0.9185OOe+03 0.18119Oe+03
0.9195OOe+03 0.179286e+03
0.9205OOe+03 0.177381e+03
0.9215OOe+03 0.17738le+03
0.9225OOe+03 0.17738le+03
0.9235OOe+03 0.175357e+03
0.9245OOe+03 0.173452e+03
0.9255OOe+03 0.173452e+03
0.9265OOe+03 0.17738le+03
0.9275OOe+03 0.17738le+03
0.9285OOe+03 0.173452e+03
0.9295OOe+03 0.181 19Oe+03
0.9305OOe+03 0.183095e+03
0.9315OOe+03 0.185OOOe+03
0.9325OOe+03 0.185OOOe+03
0.9335OOe+03 0.192738e+03
0.9345OOe+03 0.192738e+03
0.9355OOe+03 0.190833e+03
0.9365OOe+03 0.18881Oe+03
0.9375OOe+03 0.18881Oe+03
0.9385OOe+03 0.18881Oe+03
0.9395OOe+03 O.l8881Oe+O3
0.9405OOe+03 O.l85OOOe+O3
0.9415OOe+03 0.185OOOe+03
0.9425OOe+03 0.185OOOe+03
0.9435OOe+03 0.185OOOe+03
0.9445OOe+03 0.185oOOe+03
0.9455OOe+03 0.185OOOe+03
0.9465OOe+03 O.t85OOOe+03
0.9475OOe+03 0.185OOOe+03
0.9485OOe+03 0.185oOOe+03
0.9495OOe+03 0.173452e+03
0.9505OOe+03 0.18119Oe+03
0.9515OOe+03 0.181 19Oe+03
0.9525OOe+03 O.l8119Oe+O3
0.9535OOe+03 0.177381e+03
0.9545OOe+03 0.18119Oe+03
0.9555OOe+03 0.181 19Oe+03
0.9565&+03 0.18119Oe+03
0.9575OOe+03 0.185ooOe+03
0.9585OOe+03 0.18881Oe+O3
0.9595OOe+03 0.18881Oe+03
0.9605OOe+03 0.18881Oe+03
0.9615OOe+03 0.18881Oe+03
- 162 -
BR13 Response to BR20 Discharge-
Time Specific Pressure
(4 (~a-s~rn~)
0.9625OOe+03 0.185oOOe+03
0.9635OOe+03 0.181 19Oe+03
0.9645OOe+03 0.18119Oe+03
0.9655OOe+03 0.17738 1e+03
0.9665OOe+03 0.177 38 1e+03
0.9675OOe+03 0.17738 le+03
0.9685OOe+03 0.181 19Oe+03
0.9695OOe+03 0.177381e+03
0.9705OOe+03 0.17345&+03
0.97 15OOe+03 0.169643e+03
0.9725OOe+03 0.169643e+03
0.9735OOe+03 O.l65833e+03
0.9745OOe+03 0.165833e+03
0.9755OOe+03 0.163929e+03
0.9765OOe+03 0.165833e+03
0.9775OOe+03 0.162024e+03
0.9785OOe+03 0.158095e+03
0.9795OOe+03 0.135OOOe+03
0.9805OOe+03 0.165833e+03
0.98 15OOe+03 0.173452e+03
0.9825OOe+03 O.l73452e+03
0.9835OOe+03 0.17738 le+03
0.9845OOe+03 0.17738 1e+03
0.9855OOe+03 0.17738 le+03
0.9865OOe+03 0.173452e+03
0.9875OOe+03 0.171548e+03
0.9885OOe+03 0.169643e+03
0.9895OOe+03 0.165833e+03
0.9905OOe+03 O.l65833e+03
0.99 15OOe+03 O.l65833e+03
0.9925OOe+03 0.158095e+03
0.9935OOe+03 0.165833e+03
0.9945OOe+03 O.l65833e+03
0.9955OOe+03 0.162024e+03
0.9965OOe+03 0.162024e+03
0.9975OOe+03 0.158095e+03
0.9985OOe+03 0.158095e+03
0.9995OOe+03 0.162024e+03
0.10005Oe+04 0.158095e+03
O.l0015Oe+O4 0.158095e+03
0.10025Oe+O4 0.154286e+03
0.10035Oe+O4 0.158095e+03
0.10045Oe+04 0.158095e+03
0.10055Oe+04 O.l69643e+0?
O.l0065Oe+O4 O.l69643e+O?
0.10075Oe+O4 O.l73452e+O?
0.10085Oe+04 0.17738 le+O?
0.10095Oe+04 0.175357e+OT
0.10 105Oe+O4 0.173452e+O?
0.101 15Oe+O4 0.169643e+O:
0.10125Oe+O4 0.165833e+0:
- 163 -

-
BR13 Response to BR20 Discharge
Time specific pressure
(h) (k~a-sirn~)
0.10135Oe+O4 0.162024e+03
0.10145Oe+O4 0.162024e+03
0.10155Oe+04 0.158095e+03
0.10165Oe+04 0.160000e+03
0.10175Oe+O4 0.158095e+03
0.10185Oe+O4 0.162024e+03
0.10195Oe+04 0.158095e+03
0.10205Oe+04 0.162024e+03
0.10215Oe+04 0.162024e+03
0.10225Oe+04 O.l65833e+03
0.10235Oe+04 0.165833e+03
0.10245Oe+04 0.165833e+03
0.10255Oe+04 0.165833e+03
0.10265Oe+O4 O.l65833e+03
0.10275Oe+04 0.154286e+O3
0.10285Oe+04 0.169643e+03
0.10295Oe+O4 0.154286e+03
0.10305Oe+04 0.18119Oe+03
0.103 15Oe+04 0.181 19Oe+03
0.10325Oe+04 0.17738 le+03
O.l0335Oe+04 0.177 381e+03
0.10345Oe+04 0.177381e+03
0.10785Oe+04 0.150476e+03
0.10795Oe+04 0.154286e+03
0.10805Oe+04 0.158095e+03
0.10815Oe+04 0.158095e+03
0.10825Oe+04 0.162024e+03
0.10835Oe+04 0.158095e+03
0.10845Oe+04 0.158095e+03
0.10855Oe+04 0.154286e+03
0.10865Oe+04 0.15428&+03
0.10875Oe+04 0.154286e+03
0.10885Oe+04 0.158095e+03
0.10895Oe+04 O.l54286e+O3
0.10905Oe+04 O.W286e+O3
0.10915Oe+04 0.15238 le+03
0.10925Oe+04 0.150476e+03
0.10935Oe+04 0.154286e+03
0.10945Oe+04 O.l54286e+O3
0.10955Oe+O4 0.154286e+03
0.11015Oe+04 0.150476e+03
O.l1025Oe+O4 0.150476e+03
0.1 1035Oe+W 0.158095e+03
O.l1045Oe+04 0.162024e+03
0.11055Oe+04 0.160000e+03
0.11065Oe+04 0.158095e+03
0.11075Oe+04 0.158095e+03
0.11085Oe+04 0.158095e+03
0.11095Oe+04 0.154286e+03
0.11 105Oe+O4 0.154286e+03
0.1 1115Oe+O4 0.150476e+03
- 164 -

BR13 Reswnse to BR20 Discharge


1

Time Specific Pressure


(h) (~a-s~m')
0.11 125Oe+O4 0.148452e+03
0.11 135Oe+O4 0.150476e+03
0.11145Oe+O4 0.150476e+03
0.11155Oe+O4 0.150476e+03
0.11165Oe+O4 0.150476e+03
0.11175Oe+O4 0.150476e+03
0.11185Oe+O4 0.15238le+03
0.1119sOe+O4 0.1M286e+03
0.11205Oe+O4 0.1928&+03
0.11215Oe+O4 0.1928&+03
O.l1225Oe+O4 0.154286e+03
0.11235Oe+O4 0.150476e+03
0.11245Oe+O4 0.1928&+03
O.l1255Oe+O4 0.158095e+03
O.l1265Oe+O4 O.l6S833e+03
O.l1275Oe+O4 0.165833e+03
0.11285Oe+O4 O.l65833e+O3
0.11295Oe+O4 O.l65833e+03
0.11305Oe+O4 0.165833e+03
O.l1315Oe+O4 0.162024e+03
0.11325Oe+04 0.158095e+03
0.11335Oe+O4 0.154286e+03
0.11345Oe+O4 0.15O476e+03
0.113S5Oe+O4 0.150476e+03
O.l1365Oe+OQ O.l58095e+03
0.11375Oe+O4 0.158095e+03
0.1 1385Oe+O4 0.15809Se+03
0.11395Oe+04 0.15428&+03
O.l1405Oe+O4 0.154286e+03
0.11415Oe+O4 0.19286e+03
0.11425Oe+O4 0.154286e+03
0.11435Oe+O4 0.150476e+03
0.11445Oe+O4 0.15O476e+03
0.11455Oe+04 0.15O47&+03
O.l14650e+W 0.150476e+03
O.l1475Oe+O4 0.15428&+03
0.11485Oe+O4 O.l58095e+03
0.11495Oe+O4 0.158095e+03
O.l1505Oe+OQ 0.162024e+03
0.11515Oe+O4 0.162024e+03
0.11525Oe+O4 0.162024e+03
0.11535Oe+O4 0.162024e+O?
O.l1545Oe+O4 0.158095e+03
O.l15S5Oe+O4 0.158095e+0?
O.l1565Oe+O4 0.154286e+03
0.11575Oe+O4 0.152381e+O?
0.11585Oe+O4 0.150476e+03
O.l1595Oe+O4 0.148452e+03
0.11605Oe+O4 0.146548e+03
0.11615Oe+04 0.146548e+03
0.11625Oe+O4 0.142738e+03
- 165 -
BR13 Response to BR20 Discharee
Time
(h)
0.11635Oe+O4
0.11645Oe+O4
0.1 1655Oe+O4
0.11665Oe+O4
0.11675Oe+O4
0.1 1685Oe+O4
0.11695Oe+O4
0.11705Oe+O4
0.11715Oe+O4
0.11725Oe+O4
0.11735Oe+O4
0.11745Oe+O4
0.11755Oe+O4
0.11765Oe+O4
0.11775Oe+O4
O.l1785Oe+O4
0.11795Oe+O4
O.l1805Oe+W
0.118 15Oe+O4
0.11825Oe+O4
0.1 1835Oe+O4
0.11845Oe+O4
0.11855Oe+O4
O.l1865Oe+O4
0.1 1875Oe+O4
O.l1885Oe+W
0.11895Oe+O4
0.11905Oe+O4
0.119 15Oe+04
0.1 1925Oe+O4
O.l1935Oe+O4
0.11945Oe+O4
0.11955Oe+O4
0.1 1965Oe+M
0.11975Oe+O4
0.11985Oe+O4
0.11995Oe+O4
0.12OO5Oe+O4
0.12015Oe+O4
0.12025Oe+O4
0.12035Oe+O4
0.12045Oe+O4
0.12055Oe+04
0.120650e+04
0.12075Oe+04
0.12085Oe+04
0.12095Oe+04
0.12105Oe+O4
0.121 15Oe+04
0.12125Oe+04
0.12135Oe+04
- 166-

I BR13 Reswnse to BR20 Discharge


Time specific pressure
(h) (k~a-slrn~)
0.12145Oe+04 0.146548e+03
0.12155Oe+O4 0.146548e+03
0.12165Oe+O4 0.146548e+03
0.12175Oe+O4 0.146548e+03
- 167 -
TABLE B11
TEST C2: BR23 RESPONSE TO BR20 DISCHARGE

-
BR23 Response to BR20 Discharge
Specific Pressure
(k~a-slrn~)
0.272618e+02
0.34166&+02
0.546428e+02
0.614286e+02
0.819046e+02
0.819046e+02
0.955952e+02
0.1025oOe+03
0.116071e+03
0.116071e+03
0.12297&+03
0.136667e+03
0.293809e+03
0.328095e+03
0.328095e+03
0.334881e+03
0.334881e+03
0.33488le+03
0.334881e+03
0.334881e+03
0.341667e+03
0.34857 le+03
0.341667e+03
0.34857 le+03
0.34857 le+03
0.355357e+03
0.355357e+03
0.362262e+03
0.362262e+03
0.36226&+03
0.3657 14e+03
0.362262e+03
0.369048e+03
0.369048e+03
0.369048e+03
0.375952e+03
0.375952e+03
0.37595&+03
0.375952e+03
0.382738e+03
0.382738e+03
0.382738e+03
0.38619 le+03
0.3 82738e+03
0.389643e+03
-168-

BR23 Response to BR20 Discharge-


Time specific Pressure
(4 (~a-slrn~)
0.1015OOe+03 0.389643e+03
0.1025OOe+03 0.389643e+03
0.1035OOe+03 0.389643e+03
0.1045OOe+03 0.389643e+03
0.1055OOe+03 0.396429e+03
0.1065OOe+03 0.396429e+03
0.1075OOe+03 0.403333e+03
0.1085OOe+03 0.403333e+03
0.1095OOe+03 0.403333e+03
0.1105OOe+03 0.403333e+03
0.11 15OOe+03 0.4101 19e+03
0.1125OOe+03 0.403333e+03
0.1135OOe+03 0.410119e+03
0.1145OOe+03 0.410119e+03
0.1155OOe+03 0.4101 19e+03
0.1165OOe+03 0.410119e+03
0.1175OOe+03 0.410119e+03
0.1 185oOe+03 0.416905e+03
0.1195OOe+03 0.4 10119e+03
0.1205OOe+03 0.4 16905e+03
0.1215OOe+03 0.410119e+03
0.1225OOe+03 0.4 16905e+03
0.1235OOe+03 0.4 16905e+03
0.1245OOe+O3 0.420357e+03
0.1255OOe+03 0.423809e+03
0.1265OOe+03 0.423809e+03
0.127 5OOe+03 0.423809e+03
0.1285OOe+03 0.423809e+03
0.1295OOe+03 0.423809e+03
0.1305OOe+03 0.430595e+03
0.1315OOe+03 0.423809e+03
0.1325OOe+03 0.430595e+03
0.1335OOe+03 0.430595e+03
0.1345OOe+03 0.430595e+03
0.1355OOe+03 0.4375OOe+03
0.136500e+03 0.437500e+03
0.1375OOe+03 0.4375OOe+03
0.1385OOe+03 0.4375OOe+03
0.1395OOe+03 0.4375OOe+03
0.1405OOe+03 0.440952e+03
0.1415OOe+03 0.4375OOe+03
0.1425OOe+03 0.4375OOe+03
0.1435OOe+03 0.444286e+03
0.1445OOe+03 0.4375OOe+03
0.1455OOe+03 0.444286e+03
0.1465OOe+03 0.444286e+03
0.1475OOe+03 0.444286e+03
0.1485OOe+03 0.444286e+03
0.1495OOe+03 0.44428&+03
0.1505OOe+03 0.451 19Oe+03
0.15 15OOe+03 0.451 19Oe+03
-169-

-
BR23 Resmnse to BR20 Discharge
Specific Pressure
(k~a-s~rn’)
0.451 19Oe+03
0.451 19Oe+03
0.45797&+03
0.45119Oe+03
0.45797&+03
0.45797&+03
0.45797&+03
0.45797&+03
0.461428e+03
0.464881e+03
0.499048e+03
0.499048e+03
0.499048e+03
0.499048e+03
0.499048e+03
0.505952e+03
0.505952e+03
0.505952e+03
0.505952e+03
0.50595&+03
0.505952e+03
0.505952e+03
0.505952e+03
0.505952e+03
0.505952e+03
0.512738e+03
0.509405e+03
0.505952e+03
0.512738e+03
0.512738e+03
0.512738e+03
0.5 12738e+03
0.509405e+03
0.512738e+03
0.512738e+03
0.5 12738e+03
0.5 12738e+03
0.5 12738e+03
0.5 1619Oe+03
0.512738e+03
0.5 12738e+03
0.5 12738e+O:!
0.5 19643e+03
0.519643e+O?
0.519643e+03
0.512738e+O?
0.519643e+O?
0.519643e+O?
0.519643e+0?
0.519643e+O?
0.519643e+0?
- 170 -
BR23 ResDonse to BR20 Discharge-
Time specific pressure
(h) (k~a-strn~)
0.2355OOe+03 0.5 19643e+03
0.2365OOe+03 0.5 19643e+03
0.2375OOe+03 0.5 19643e+03
0.2385OOe+03 0.522976e+03
0.2395OOe+03 0.519643e+03
0.2405OOe+03 0.519643e+03
0.2415OOe+03 0.519643e+03
0.2425OOe+03 0.519643e+03
0.2435OOe+03 0.522976e+03
0.2A45OOe+03 0.526429e+03
0.2455OOe+03 0.526429e+03
0.2465OOe+03 0.519643e+03
0.2475OOe+03 0.526429e+03
0.2485OOe+03 0.526429e+03
0.2495OOe+03 0.526429e+03
0.2505OOe+03 0.519643e+03
0.25 15OOe+03 0.526429e+03
0.2525OOe+03 0.533333e+03
0.2535OOe+03 0.526429e+03
0.2545OOe+03 0.526429e+03
O,2555OOe+03 0.533333e+03
0.2565OOe+03 0.533333e+03
0.2585OOe+03 0.536667e+03
0.2595OOe+03 0.533333e+03
0.2605OOe+03 0.533333e+03
0.2615OOe+03 0.533333e+03
0.2625&+03 0.533333e+03
0.2635OOe+03 0.533333e+03
0.264500e+03 0.533333e+03
0.2655OOe+03 0.5401 19e+03
0.2665OOe+03 0.5401 19e+03
0.2675OOe+03 0.533333e+03
0.2685OOe+03 0.5401 19e+03
0.2695OOe+03 0.5401 19e+03
0.2705OOe+03 0.547024e+03
0.27 15OOe+03 0.5401 19e+03
0.2725OOe+03 0.54357 le+03
0.2735OOe+03 0.547024e+O?
0.2745OOe+03 0.547024e+O?
0.2755OOe+03 0.547024e+O?
0.2765OOe+03 0.547024e+O?
0.2775OOe+03 0.547024e+O:
0.2785OOe+03 0.547024e+O?
0.2795OOe+03 0.553810e+Of
0.2805OOe+03 0.547024e+O:
0.28 15&+03 0.547024e+O:
0.2825OOe+03 0.5538 10e+O:
0.2835OOe+03 0.553810e+O:
0.2845OOe+03 0.553810e+K
0.2855OOe+03 0.547024e+O:
0.2865OOe+03 0.547024e+0:
- 171 -

Time Specific Pressure


(h)
0.287 5OOe+03
0.2885OOe+03 0.55381Oe+03
0.2895OOe+03 0.55381Oe+03
0.2905OOe+03 0.55381Oe+03
0.2915OOe+03 0.55381Oe+03
0.2925OOe+03 0.55381Oe+03
0.2935OOe+03 0.55381Oe+03
0.2945OOe+03 0.55381Oe+03
0.2955&+03 0.5607 14e+03
0.2965OOei-03 0.557262e+03
0.2975OOe+03 0.560714e+03
0.2985&+03 0.5607 14e+03
0.2995OOe+03 0.55381Oe+03
0.3005OOe+03 0.560714e+03
0.3015OOe+03 0.55381Oe+03
0.3025OOe+03 0.560714e+03
0.3035OOe+03 0.560714e+03
0.3045OOe+03 0.5607 14e+03
0.3055OOe+03 0.5607 14e+03
0.3065OOe+03 0.5607 14e+03
0.3075OOe+03 0.5607 14e+03
0.3085OOe+03 0.564047e+03
0.3095OOe+03 0.5607 14e+03
0.3 105OOe+03 0.5675OOe+03
0.3115OOe+03 0.5675OOe+03
0.3125OOe+03 0.5 675OOe+03
0.3135OOe+03 0.5675OOe+03
0.3145OOe+03 0.5675OOe+03
0.3 155OOe+03 0.5675OOe+03
0.3165OOe+03 0.5675OOe+03
0.3175OOe+03 0.560714e+03
0.3185OOe+03 0.5675OOe+03
0.3195OOe+03 0.5675OOe+03
0.3205OOe+03 0.5675OOe+03
0.3215oOe+03 0.5675OOe+03
0.3225OOe+03 0.5675OOe+03
0.3235OOe+03 0.5675OOe+03
0.3245OOe+03 0.5675OOe+03
0.325 5OOe+03 0.5675OOe+03
0.3265OOe+03 0.574405e+03
0.327 5.OOe+03 0.5675OOe+03
0.3285OOe+03 0.567 5OOe+03
0.3295OOe+03 0.5675OOe+03
0.3305OOe+03 0.5675OOe+03
0.3315OOe+03 0.5675OOe+03
0.3325OOe+03 0.574405e+03
0.3335OOe+03 0.574405e+03
0.3345OOe+03 0.5675OOe+03
0.3355OOe+03 0.574405e+03
0.3365OOe+03 0.574405e+03
0.3375OOe+03 0.574405e+03
- 172 -
BR23 Response to BR20 Discharge-

specific pressure
(k~a-strn~)
0.574405e+03
0.5675OOe+03
0.560714e+03
0.547024e+03
0.533333e+03
0.519643e+03
0.505952e+03
0.492262e+03
0.485357e+03
0.47857 le+03
0.471667e+03
0.46488le+03
0.457976e+03
0.451 19Oe+03
0.44428&+03
0.44428&+03
0.4375OOe+03
0.430595e+03
0.430595e+03
0.423809e+03
0.423809e+03
0.4 16905e+03
0.4101 19e+03
0.410119e+03
0.403333e+03
0.403333e+03
0.396429e+03
0.396429e+03
0.389643e+03
0.389643e+03
0.382738e+03
0.382738e+03
0.375952e+03
0.369048e+03
0.369048e+03
0.369048e+03
0.362262e+03
0.362262e+03
0.362262e+03
0.355357e+03
0.355357e+03
0.34857 le+03
0.287024e+03
0.287024e+03
0.287024e+03
0.287024e+03
0.283571e+03
0.2801 19e+03
0.2801 19e+03
0.273333e+03
0.273333e+03
- 173 -
BW3 Response to BR20 Discharge-

Specific Pressure
(~a-s~rn’)
0.273333e+03
0.273333e+03
0.266548e+03
0.266548e+03
0.266548e+03
0.266548e+03
0.266548e+03
0.266548e+03
0.266548e+03
0.259643e+03
0.252857e+03
0.259643e+03
0.252857e+03
0.252857e+03
0.252857e+03
0.252857e+03
0.252857e+03
0.24595&+03
0.245952e+03
0.249405e+03
0.245952e+03
0.245952e+03
0.245952e+03
0.245952e+03
0.24595&+03
0.245952e+03
0.239167e+03
0.239167e+03
0.242619e+03
0.239 167e+03
0.239 167e+03
0.239167e+03
0.239 167e+03
0.239 167e+03
0.239167e+03
0.239 167e+03
0.232262e+03
0.2357 14e+03
0.232262e+03
0.232262e+03
0.232262e+03
0.232262e+03
0.232262e+03
0.232262e+03
0.232262e+03
0.232262e+03
0.232262e+03
0.23226&+03
0.232262e+03
0.232262e+03
0.22547&+03
- 174 -

BR23 Response to BR20 Discharge


-
Time specific pressure
(h) (k~a-sim~)
0.4685OOe+03 0.22547&+03
0.4695OOe+03 0.22547&+03
0.4705OOe+03 0.21869Oe+03
0.47 15OOe+03 0.21869&+03
0.4725OOe+03 0.21869Oe+03
0.4735OOe+03 0.2 1869Oe+03
0.474 5OOe+03 0.21869Oe+03
0.4755OOe+03 0.21178&+03
0.4765OOe+03 0.21178&+03
0.4775OOe+03 0.21178&+03
0.4785OOe+03 0.21 178&+03
0.4795OOe+03 0.21178&+03
0.4805OOe+03 0.21178&+03
0.4815OOe+03 0.208452e+03
0.4825OOe+03 0.205OOOe+03
0.4835OOe+03 0.205OOOe+03
0.4845OOe+03 0.205OOOe+03
0.4855OOe+03 0.205OOOe+03
0.4865OOe+03 0.205OOOe+03
0.4875OOe+03 0.205OOOe+03
0.48 85OOe+03 0.198095e+O3
0.4895OOe+03 0.198095e+03
0.4905OOe+03 0.198095e+03
0.4915OOe+03 0.198095e+03
0.4925OOe+03 0.198095e+03
0.4935OOe+03 0.198095e+03
0.4945OOe+03 0.198095e+03
0.4955OOe+03 0.198O95e+O3
0.4965OOe+03 0.198095e+03
0.4975OOe+03 0.191309e+03
0.4985OOe+03 0.191309e+03
0.4995OOe+03 0.191309e+03
0.5005OOe+03 0.191309e+03
0.50 15OOe+03 0.191309e+03
0.5025OOe+03 0.191309e+03
0.5035OOe+03 0.191309e+03
0.5045OOe+03 0.184405e+03
0.5055OOe+03 0.191309e+03
0.5065OOe+03 0.191309e+03
0.507 5OOe+03 0.1913We+03
0.5085OOe+03 0.19 1309e+03
0.5095OOe+03 0.191309e+03
0.5105OOe+03 0.191309e+03
0.51 15OOe+03 0.191309e+03
0.5125OOe+03 0.191309e+03
0.5 135OOe+03 0.184405e+03
0.5145OOe+03 0.191309e+03
0.5155OOe+03 0.184405e+03
0.5 165OOe+03 O.l84405e+03
0.5 175OOe+03 0.184405e+03
0.5185OOe+03 0.184405e+03
- 175 -
BR23 Response to 120 Discharge
Time Specific Pressure
(h) ~a-s~rn’)
0.5 195OOe+03 0.184405e+03
0.5205OOe+03 0.184405e+03
0.52 15OOe+03 0.184405e+03
0.5225OOe+03 0.184405e+03
0.5235OOe+03 0.184405e+03
0.5245OOe+03 0.184405e+03
0.5255OOe+03 0.184405e+03
0.5265OOe+03 0.184405e+03
0.5275OOe+03 0.184405e+03
0.5285OOe+03 0.177619e+03
0.5295OOe+03 0.1776 19e+03
0.5305OOe+03 0.184405e+03
0.5315OOe+03 0.177619e+03
0.5325OOe+03 0.1776 19e+03
0.5335OOe+03 0.177619e+03
0.5345OOe+03 0.177619e+03
0.5355OOe+03 0.1776 19e+03
0.5365OOe+03 0.177619e+03
0.5375OOe+03 0.177619e+03
0.5385OOe+03 0.170833e+03
0.5395OOe+03 0.177619e+03
0.5405OOe+03 0.177619e+03
0.54 15OOe+03 0.170833e+03
0.5425OOe+03 0.170833e+03
0.5435OOe+03 0.177 619e+03
0.5445OOe+03 0.177619e+03
0.57 85OOe+03 0.150357e+03
0.5795OOe+03 0.150357e+03
- 176 -
TABLE B12
TEST C3: BR34 RESPONSE TO BR20 DISCHARGE

Time Specific Pressure


(h) (k~a-strn’)
0.25m+01 0.214277e+01
0.450000e+01 0.42859Oe+Ol
0.550000e+01 0.42859Oe+Ol
0.650000e+01 0.214277e+01
0.750000e+01 0.42859Oe+01
0.850000e+01 0.42859Oe+01
0.950000e+01 0.642867e+01
0.105OOOe+O2 0.857 144e+01
0.1 15OOOe+O2 0.857 144e+01
0.125OOOe+M 0.857 144e+01
0.545OOOe+02 0.100357e+03
0.695OOOe+02 0.136667e+03
0.705OOOe+02 0.136667e+03
0.715000e+M 0.140952e+03
0.725OOOe+02 0.140952e+03
0.735OOOe+02 0.145238e+03
0.745OOOe+02 0.145238e+03
0.755OOOe+02 0.145238e+03
0.765OOOe+02 0.149524e+03
0.775OOOe+02 0.149524e+03
0.785OOOe+02 0.15381Oe+03
0.795OOOe+02 0.15381Oe+03
0.805OOOe+O2 0.158095e+03
0.8 15OOOe+02 0.158095e+03
0.825OOOe+02 0.16238 le+03
0.835OOOe+02 0.166667e+03
0.845OOOe+02 0.170833e+03
0.855OOOe+02 0.170833e+03
0.865OOOe+02 0.175 119e+03
0.875OOOe+02 0.1751 19e+0?
0.885OOOe+02 0.179405e+O?
0.895OOOe+02 0.179405e+Or
0.905OOOe+02 O.l8369Oe+O:
0.915OOOe+02 O.l8369Oe+O:
0.925OOOe+02 0.187976e+E
0.935OOOe+02 0.187976e+O:
0.945OOOe+02 0.192262e+O:
0.955OOOe+02 0.192262e+0:
0.965OOOe+02 0.192262e+0:
0.975OOOe+02 0.196548e+0:
0.985OOOe+02 0.200833e+U
0.995OOOe+02 0.200833e+O:
0.1005OOe+03 0.205 119e+0:
0.10 15OOe+03 0.205 119e+0:
0.1025OOe+03 0.205 119e+0:
- in -

-
BR34 Reswnse to B E 0 Discharge
Time specificpressure
(h) (k~a-sirn~)
0.1035OOe+03 0.2051 19e+03
0.1045OOe+03 0.209405e+03
0.1055OOe+03 0.21369le+03
0.1065OOe+03 0.213691e+03
0.1075OOe+03 0.217976e+03
0.1085OOe+03 0.217976e+03
0.1O95OOe+03 0.217976e+03
O.l105OOe+O3 0.222262e+03
0.1 1 15OOe+03 0.222262e+03
0.1125OOe+03 0.226428e+03
0.1135OOe+03 0.230714e+03
O.l145OOe+O3 0.230714e+03
0.1155OOe+03 0.232857e+03
0.1165OOe+03 0.235OOOe+03
O.l175OOe+O3 0.237143e+03
O.l185OOe+O3 0.23928&+03
0.1195OOe+03 0.239286e+03
0.1205OOe+03 0.243571e+03
0.1215OOe+03 0.243571e+03
0.1225OOe+03 0.243571e+03
0.1235OOe+03 0.247857e+03
0.1245OOe+03 0.247857e+03
0.1255OOe+03 0.247857e+03
0.1265OOe+03 0.252143e+03
0.1275OOe+03 0.252143e+03
0.1285OOe+03 0.252143e+03
0.1295OOe+03 0.256429e+03
0.1305OOe+03 0.256429e+03
0.1315OOe+03 0.260714e+03
0.1325OOe+03 0.260714e+03
0.1335OOe+03 0.260714e+03
0.1345OOe+03 0.265OOOe+03
0.1355OOe+03 0.265OOOe+03
0.1365OOe+03 0.265OOOe+03
0.1375OOe+03 0.267143e+03
O.l385OOe+O3 0.269286e+03
0.1395OOe+03 0.269286e+03
0.1405OOe+03 0.269286e+03
0.1415OOe+03 0.273572e+03
0.1425OOe+03 0.273572e+03
0.1435OOe+03 0.277857e+03
0.1445OOe+03 0.277857e+03
0.1455OOe+03 0.277857e+03
O.l465OOe+O3 0.286309e+03
0.1475OOe+03 0.286309e+03
0.1485OOe+03 0.286309e+03
0.1495OOe+03 0.288453e+03
0.1505OOe+03 0.290595e+03
0.1515OOe+03 0.290595e+03
0.1525OOe+03 0.294881e+03
O.l535OOe+O3 0.290595e+03
- 178 -
-
BR34 Resmnse to BR20 Discharge
specific Pressure
(k~a-sirn~)
0.29488le+03
0.29488 le+03
0.299167e+03
0.299167e+03
0.301 31Oe+03
0.30131Oe+03
0.303452e+03
0.307738e+03
0.363333e+03
0.359 167e+03
0.363333e+03
0.363333e+03
0.363333e+03
0.363333e+03
0.363333e+03
0.367619e+03
0.369762e+03
0.367 619e+03
0.371905e+03
0.37 1905e+03
0.371905e+03
0.376191e+03
0.37619 le+03
0.380476e+03
0.380476e+03
0.380476e+03
0.38476%+03
0.384762e+03
0.384762e+03
0.384762e+03
0.38476%+03
0.389048e+03
0.389048e+03
0.389048e+03
0.391191e+03
0.391 191e+03
0.393333e+03
0.393333e+03
0.393333e+03
0.395476e+03
0.397 619e+03
0.397619e+03
0.39761 9e+03
0.397619e+03
0.397 619e+03
0.401905e+03
0.399762e+03
0.40619Oe+03
0.406 19Oe+03
0.408 333e+03
0.410476e+03
- 179 -
BR34 Response to BR20 Discharge -
Time specific pressure
(h) (Wa-sld)
0.2375OOe+03 0.410476e+03
0.2385OOe+03 0.41047&+03
0.2395OOe+03 0.410476e+03
0.2405OOe+03 0.4 14762e+03
0.24 15OOe+03 0.414762e+03
0.2425OOe+03 0.414762e+03
0.2435OOe+03 0.414762e+03
0.2445OOe+03 0.416905e+03
0.2455OOe+03 0.419048e+03
0.2465OOe+03 0.4 19048e+03
0.2475OOe+03 0.421 191e+03
0.2485OOe+03 0.419048e+03
0.2495OOe+03 0.419048e+03
0.2505OOe+03 0.419048e+03
0.2515OOe+03 0.423333e+03
0.2525OOe+03 0.423333e+03
0.2535OOe+03 0.427619e+03
0.2545OOe+03 0.427619e+03
0.2555OOe+03 0.429762e+03
0.2565OOe+03 0.431905e+03
0.2585OOe+03 0.431905e+03
0.2595OOe+03 0.431905e+03
0.2605OOe+03 0.43619 le+03
0.26 15OOe+03 0.4 3619le+03
0.2625OOe+03 0.43619 le+03
0.2635OOe+03 0.43619 le+03
0.2645OOe+03 0.440476e+03
0.2655OOe+03 0.43619 le+03
0.2665OOe+03 0.44047&+03
0.2675OOe+03 0.444762e+03
0.2685OOe+03 0.44476%+03
0.2695OOe+03 0.444762e+03
0.2705OOe+03 0.444762e+03
0.27 15OOe+03 0.444762e+03
0.2725OOe+03 0.444762e+03
0.2735OOe+03 0.449048e+03
0.2745OOe+03 0.449048e+03
0.2755OOe+03 0.449048e+03
0.2765OOe+03 0.453333e+03
0.2775OOe+03 0.453333e+03
0.27 85OOe+03 0.453333e+03
0.2795OOe+03 0.457619e+03
0.2805OOe+03 0.453333e+03
0.28 15OOe+03 0.457 619e+03
0.2825OOe+03 0.461905e+03
0.2835OOe+03 0.461905e+03
0.2845OOe+03 0.464048e+03
0.2855OOe+03 0.466 191e+O3
0.2865OOe+03 0.466191e+03
0.2875OOe+03 0.470476e+O?
0.2885OOe+03 0.46619 le+O?
- 180 -
BR34 Response to BR20 Discharge-
Time specific Pressure
(h) (wa-s/m3)
0.2895OOe+03 0.470476e+03
0.2905OOe+03 0.470476e+03
0.29 15OOe+03 0.470476e+03
0.2925OOe+03 0.470476e+03
0.2935OOe+03 0.474762e+03
0.2945OOe+03 0.470476e+03
0,2!955OOe+O3 0.47476%+03
0.2965OOe+03 0.474762e+03
0.297 5OOe+03 0.47476%+03
0.2985OOe+03 0.479048e+03
0.2995OOe+03 0.479048e+03
0.3005OOe+03 0.479048e+03
0.3015OOe+03 0.479048e+03
0.3025OOe+03 0.483214e+03
0.3035OOe+03 0.479048e+03
0.3045OOe+03 0.483214e+03
0.3055OOe+03 0.485357e+03
0.3065OOe+03 0.4832 14e+03
0.3075OOe+03 0.4875OOe+03
0.3085OOe+03 0.4875OOe+03
0.3095OOe+03 0.491786e+03
0.3 105OOe+03 0.49 1786e+03
0.3115OOe+03 0.49 1786e+03
0.3 125OOe+03 0.491786e+03
0.3135OOe+03 0.491786e+03
0.3 145OOe+03 0.49 1786e+03
0.3 155OOe+03 0.49 1786e+03
0.3 165OOe+03 0.496072e+03
0.3175OOe+03 0.496072e+03
0.3185OOe+03 0.496072e+03
0.3195OOe+03 0.496072e+03
0.3205OOe+03 0.500357e+03
0.32 15OOe+03 0.500357e+03
0.3225OOe+03 0.504643e+03
0.3235OOe+03 0.500357e+03
0.3245OOe+03 0.504643e+03
0.3255OOe+03 0.504643e+03
0.3265OOe+03 0.504643e+03
0.3275OOe+03 0.500357e+03
0.3285OOe+03 0.504643e+03
0.3295OOe+03 0.506786e+03
0.3305OOe+03 0.508929e+03
0.33 15OOe+03 0.508929e+03
0.3325OOe+03 0.508929e+03
0.33 35OOe+03 0.508929e+03
0.3345OOe+03 0.5 13214e+0?
0.3355OOe+03 0.508929e+03
0.3365OOe+03 0.5 13214e+O?
0.3375OOe+03 0.5 13214e+03
0.3385OOe+03 0.5 13214e+0?
0.3395OOe+03 0.5 13214e+03
- 181 -
BR34 ResDonse to BR20 Discharge
1

Time Specific Pressure


(h) (k~a-strn~)
0.3405OOe+03 0.5 13214e+03
0.3415OOe+03 0.513214e+03
0.3425OOe+03 0.5 13214e+O3
0.3435OOe+03 0.5175OOe+03
0.3445OOe+03 0.508929e+03
0.3455OOe+03 0.5 13214e+03
0.3465OOe+03 0.5 13214e+03
0.3475OOe+03 0.5 15357e+03
0.3485OOe+03 0.5 13214e+03
0.3495OOe+03 0.5 13214e+03
0.3505OOe+03 0.513214e+03
0.3515OOe+03 0.508929e+03
0.3525OOe+03 0.508929e+03
0.3535OOe+03 0.508929e+03
0.3545OOe+03 0.504643e+03
0.3555OOe+03 0.504643e+03
0.3565OOe+03 0.504643e+03
0.3575OOe+03 0.504643e+03
0.3585OOe+03 0.500357e+03
0.3595OOe+03 0.500357e+03
0.3605OOe+03 0.500357e+03
0.36 15OOe+03 0.496072e+03
0.3625OOe+03 0.493928e+03
0.3635OOe+03 0.491786e+03
0.3645OOe+03 0.491786e+03
0.3655OOe+03 0.489643e+03
0.3665OOe+03 0.4875OOe+03
0.3675OOe+03 0.4875OOe+03
0.3685OOe+03 0.485357e+03
0.3695OOe+03 0.483214e+03
03705OOe+O3 0.483214e+03
0.37 15OOe+03 0.481191e+03
0.3725OOe+03 0.479048e+03
0.3735OOe+03 0.474762e+03
0.3745OOe+03 0.474762e+03
0.3755OOe+03 0.474762e+03
0.3765OOe+03 0.470476e+03
0.3775OOe+03 0.470476e+03
0.3785OOe+03 0.466 191e+03
0.3795OOe+03 0.466191e+03
0.4075OOe+03 0.419048e+03
0.4085OOe+03 0.419048e+03
0.4095OOe+03 0.4 1476%+03
0.4105OOe+03 0.414762e+03
0.4 115OOe+03 0.4 10476e+03
0.4 125OOe+03 0.410476e+O?
0.4135OOe+03 0.410476e+03
0.4145OOe+03 0.410476e+O?
0.4155OOe+03 0.40619Oe+03
0.4165OOe+03 0.406190e+O?
0.4175OOe+03 0.401905e+0?
- 182 -
BR34 Response to BR20 Discharge
-
Time Specific Pressure
(h) (k~a-strn~)
0.4185OOe+03 0.401905e+03
0.4195OOe+03 0.401905e+03
0.4205OOe+03 0.40 1905e+03
0.42 15OOe+03 0.3976 19e+03
0.4225OOe+03 0.3976 19e+03
0.4235OOe+03 0.3976 19e+03
0.4245OOe+03 0.393333e+03
0.4255OOe+03 0.393333e+03
0.4265&+03 0.389048e+03
0.4275OOe+03 0.389048e+03
0.4285OOe+03 0.389048e+03
0.4295OOe+03 0.389048e+03
0.4305OOe+03 0.384762e+03
0.43 15OOe+03 0.384762e+03
0.4325OOe+03 0.382619e+03
0.4335OOe+03 0.380476e+03
0.4345OOe+03 0.380476e+03
0.4355OOe+03 0.380476e+03
0.4365OOe+03 0.37619 le+03
0.4375OOe+03 0.376191e+03
0.4385OOe+03 0.374048e+03
0.4395OOe+03 0.371905e+03
0.4405OOe+03 0.371905e+03
0.44 15OOe+03 0.36976%+03
0.4425OOe+03 0.3 67619e+03
0.4435OOe+03 0.367619e+03
0.4445OOe+03 0.367619e+03
0.4455OOe+03 0.363333e+03
0.44650&+03 0.3 63333e+03
0.4475OOe+03 0.363333e+03
0.4485OOe+03 0.359 167e+03
0.449500e+03 0.359 167e+03
0.4505OOe+03 0.357024e+03
0.45 15OOe+03 0.354881e+03
0.4525OOe+03 0.354881e+03
0.4535OOe+03 0.354881e+03
0.4545OOe+03 0.350595e+03
0.4555&+03 0.350595e+03
0.4565OOe+03 0.348452e+03
0.4575&+03 0.346309e+03
0.4585OOe+03 0.346309e+03
0.4605OOe+03 0.342024e+03
0.46 15OOe+03 0.342024e+03
0.4625OOe+03 0.342024e+0!
0.4635OOe+03 0.342024e+O?
0.4645OOe+03 0.337738e+O?
0.4655OOe+03 0.333453e+O?
0.4665OOe+03 0.333453e+O?
0.4675OOe+03 0.3377 38e+0?
0.4685OOe+03 0.333453e+O:
0.4695OOe+03 0.337738e+O?
- 183 -

-
Time Specific Pressure
(~a-s~rn~)
0.4705OOe+03 0.333453e+03
0.47 15OOe+03 0.333453e+03
0.4725OOe+03 0.329167e+03
0.4735OOe+03 0.33 13We+O3
0.4745OOe+03 0.329167e+03
0.4765OOe+03 0.320595e+03
0.4795OOe+03 0.320595e+03
0.4805OOe+03 0.320595e+03
0.48 15OOe+03 0.320595e+03
0.4825OOe+03 0.320595e+03
0.4835OOe+03 0.31631Oe+03
0.4845OOe+03 0.31631Oe+03
0.4855OOe+03 0.31631Oe+03
0.4865OOe+03 0.3 14167e+03
0.4875OOe+03 0.312024e+03
0.4885OOe+03 0.312024e+03
0.4895OOe+03 0.307738e+03
0.4905OOe+03 0.307738e+03
0.49 15OOe+03 0.307738e+03
0.4925OOe+03 0.307738e+03
0.4935OOe+03 0.307738e+03
0.4945OOe+03 0.307738e+03
0.4955OOe+03 0.307738e+03
0.4965OOe+03 0.303452e+03
0.4975OOe+03 0.303452e+03
0.4985OOe+03 0.303452e+03
0.4995OOe+03 0.303452e+03
0.5005OOe+03 0.303452e+03
0.5015OOe+03 0.303452e+03
0.5025OOe+03 0.299167e+03
0.5035OOe+03 0.299 167e+03
0.5045OOe+03 0.297024e+03
0.5055OOe+03 0.294881e+03
0.5065OOe+03 0.29488le+03
0.5075OOe+03 0.294881e+03
0.5085OOe+03 0.294881e+03
0.5095OOe+03 0.29488 le+03
0.5 105OOe+03 0.290595e+03
0.5115OOe+03 0.290595e+03
0.5 125OOe+03 0.290595e+O?
0.5 135OOe+03 0.288453e+O?
0.5145OOe+03 0.286309e+O?
0.5 155OOe+03 0.290595e+O?
0.5165OOe+03 0.290595e+0?
0.5175OOe+03 0.286309e+O?
0.5185OOe+03 0.288453e+O:
0.5195OOe+03 0.286309e+O?
0.5205OOe+03 0.286309e+O:
0.5215OOe+03 0.286309e+O:
0.5225OOe+03 0.288453e+O:
0.5235OOe+03 0.286309e+O:
-184-

BR34 ResDonse to BR20 Discharee

You might also like