Leaver, Jonathan D. (1986
Leaver, Jonathan D. (1986
Leaver, Jonathan D. (1986
March 1986
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-
...
-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
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
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
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
Page
ABSTRACT ..................................................................................................................... ii
.
1 INTRODUCTION ............................................................................................. ..........
+ 1
.
3 BAROMETRIC AND EARTH TIDE EFFECTS ..................................................... 13
3.1 INTRODUCTION................................................................................................... 13
4. COMPRESSIBILITY .................................................................................................. 20
5. EARLY TIME DATA MATCHING (Sageev et al . 1986) ....................................... 22
.
7 ANALYSIS OF INTERFERENCE TESTS DATA IN THE OHAAKI
-iv-
7.4.2 Test C2: BR23 Response to BR20 Discharge .......................................... 103
-V-
LIST OF FIGURES
Page
2.3 Log-log Drawdown and Buildup Response for Reservoir with Linear
I
2.5 Semi-log Drawdown Response for Reservoir with Linear Boundary
2.6 Semi-log Buildup Response for Reservoir with Linear No-Flow Boundary
3.2 Barometric and Earth Tide Fluctuations During the "C"Tests .............................. 15
Sine Waves On a Linear Barometric Trend (after Sageev et al. 1986) ..,........... 16
-vi-
3.6 Test C3:Log-log Matches of Data to Stallman Type Curves ................................ 19
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.3 Log-log Match of the Line Source To the Curve for r2/rl = 1.1
5.8 Hypothetical Cases for Missing Pressure Data for r2/rl = 2.0
5.9 Four Possible Log-log matches of Test C3 data to the Stallman Type Curves
-vii-
6.4 West Bank Well Pressure Profiles (after Bixley 1982) ......................................... 36
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
7.1 Source and Observation Well Location for Interference Tests ............................... 42
7.7 Test B1: Incorrect Semi-log Match of Buildup Data to Fox Type Curve
7.8 Test B1: Correct Semi-log Match of Buildup Data to Fox Type Curve
7.9 Test B1: Inference Ellipse Location (after Leaver et a1. 1985) ............................ 52
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.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.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.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.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.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.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.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
92
-xii-
B4 Test B4: Specific Pressure vs Time Data ............................................................. 142
B5 Test B6: Specific Pressure vs Time Da............................................................... 144
1
I,
-1
1. INTRODUCTION
Interference testing in geothermal fields is one of the principal tools used to determine bulk
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
Theis [1935] published the line source solution based on heat transfer analogy and used time1
(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
"The eflect of boundaries can be considered by more elaborate analyses lonce they
-2-
are located."
Jacob [1941] circumspectly used the term "interference testing" to describe the1 change in pres1
on of the pa0 of one of two wells when both are $owing freely, with He pa0 of ~
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
could be described by the Diffusivity equation. The use of Laplace transforms was shown
work on interference between oil fields showed how to determine pressure m a t i o n s due
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
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
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
The solution of the pressure response at any point in an infinite reservoir due 'to the dischxg
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).
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
where:
-7-
Substituting the exponential integral gives: Figure 2.2 is a portion of the figure presented by
Stallman [19521.
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
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:
+ E i [ z -4
A dimensionless pressure-time log-log graph is shown in Figure 2.3.
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
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).
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-
(2.10)
(2.11)
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).
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
- 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
Barometric variations can produce reservoir responses of more than 8 kPa while earth tide
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:
(3.4)
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
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)
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
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.
-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
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
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
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
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
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
Perrine [1956] used the following expression for two phase total system compreskibility:
Ct = S d W+ sgcg+ Cf (4.1)
R m e y [1964] used the correlations of Standing [1952] to obtain approximatiod to the partial
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
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
Fig. 4.1 Variation in water compressibility with temperature and pressure (after Home la80).
1.5
1.o
0.5
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
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
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
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
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
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
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
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).
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-
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)
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
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
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
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-
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).
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 -
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
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
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-
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).
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
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-
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.
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
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
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
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
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
Two types of instruments were used to measure pressures in the interference test$. These were:
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
- 42 -
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 -
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
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
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
10
s
0
$ 1
h
0
-:
.- 0.1
C
e
.-
0
E
I I I
0.01 Y
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:
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-
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
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
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:
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
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).
Type curve matching the data requires careful interpretation. A simple best fit /of the data on
-
e(3.4-o.o)
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
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
The corresponding inference ellipse to which the postulated linear boundab is tangent i
- 52 -
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
variations in flow rate are detailed in Table 7.3. The flow rate variation for th& test period u
,
For a minimum pressure significance level of 3 kPa (refer $3.4) and an average( early time f l o i
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
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
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
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
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
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
For a minimum pressure significance level of 3 f f u (refer §3.4) and an average etrly time flow
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
Discharge Temp. -
Discharge Enthalpy 1180 W/kg
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.
Matching data between times of 8 and 27 hours yields a match point of:
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
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:
The corresponding inference ellipse to which the postulated linear no-flow boundlary is tangent
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.
- 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
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
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.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
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
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:
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.
(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
The data appear to show oscillations about the theoretical curves which may be the result of
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
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
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
Discharge Enthalpy
Drawdown Period 6h
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
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
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
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-
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
f::0.1
h
*
-3
.-5
p.01
.-
0
0.001
0.01 0.1 I 10 100
Dinensionleso fine
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.
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
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
1 10
Horner Time
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.
For a minimum pressure significance level of 3 &Pa (refer §3.4)and an average early time flow
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
Discharge Temp. -
Discharge Enthalpy 1270 W/kg
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
The data show minor oscillations about the theoretical curves which may be the result of earth
A semi-log graph of the data is shown in Figure 7.28. The lower curve is
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:
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)
variations in flow rate are detailed in Table 7.10. A 32% LOP 1.A flow rate occur&( in the final
For a minimum pressure significance level of 3 kPa (refer $3.4)and an average ebrly time flow
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
Discharge Temp.
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.
The data show an excellent fit to the Theis line source solution but this match is non-unique
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-
(7.21)
For a minimum pressure significance level of 3 kPa (refer §3.4)and an average early time flow
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
Discharge Temp.
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
ps = 10 kPa-slm3 p~ = 0.048
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
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
Semi-log plots for the r2/tl = 1.5 match and the line source match are shown in Figures 7.36
and
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 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
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
For a minimum pressure significance level of 3 &Pa (refer $3.4) and an average early time flow
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
Permeable Depth -
1015 1055 915 * 30 m
Discharge Temp.
Wu-slm3 and chose after 100 hours of the test where there is a slight indication of a no-flow
boundary.
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
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
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
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:
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
T i r (kl
-1. mtch point0
p, = 0 p; 0 9.40
I , . . . , .,.I , , , , , , . . I , , , , ,
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
- 96 -
ower
The data for this test are presented in Table B10 and plotted in Cartesian form in Figure 7.45.
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
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
Discharge Enthalpy
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
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
Drawdown data again show oscillations about the theoretical curves which may ble the result of
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
The corresponding inference ellipse to which the first postulated linear no-flow boundary is
E
a
- 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.
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
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
- 104-
TABLE 7.15
TEST C2: BR23 RESPONSE TO BR20 DISCHARGE
TEST SPECIFICATIONS
Permeable Depth -
1015 1055 m 815 m, 945 m, 1045 rn
Discharge Enthalpy
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
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
Drawdown data again show oscillations about the theoretical curves which may be the result of
-
The Homer plot for r2/rl 4.0 (Figure 7.53) shows a fair match of the early time buildup data
- 106 -
0.1
1 10 100 lo00
tiornor Time
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
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:
The corresponding inference ellipse to which the postulated linear no-flow boundary is tangent
is presented in Figure 7.54.
I / I
- 108 -
The data for this test are presented in Table B12 and plotted in Cartesian form in Figure 7.55.
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
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
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
Discharge Enthalpy -
Recording Meter ~
quartz crystal gauge
1.5.
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
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
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
(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-
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
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
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.
Horner plots for r2/r1 = 1.5 a m 13r the Theis line source solution are presented ixb Figures 7.60
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
8. RESULTS
Transmissivities and storativities are calculated using Equations 2.13 and 2.14. Results are
TABLE 8.1
(m) (tdkp~io~)
3
B2 23 13 279 2.50
B4 23 13 279
B6 23 20 411
B8 34 23 795 40.00
B9 34 19 727 5.20
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.
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
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
which showed both pressure support and no-flow boundary trends in Tests B1, 2,B9 and C3
resistivity
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-
9. DISCUSSION
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
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
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-
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).
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
- The principal results obtained were that the average test transmissivity was 80 d-m and tha
two high storativity values from Tests B7 and B8 are excluded. For a porosity of 20% an4
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).
(i) Atmospheric earth tide and barometric pressure measurements measured to 1 Pa which al.
(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~
~
(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-
- 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
permeability (durcy)
P pressure ( P u )
PD dimensionless pressure
storage coefficient
vapor saturation
water saturation
time (s)
dimensionless time
porosity
formation density ( k g / d )
- 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
Homogeneous Liquids in Fissured Rocks (Strata)," J. Appl. Math. Mech., 0.S.S.R. (1960),
pp 1286-1303.
Workshop (1980).
27 11-27 18.
I
Braithwaite R L : "Perrological Report: Broadlands Drillhole Br34," Gedlogical Survey4
(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
1
Reservoirs With Double Porosity Behaviour-Theory and Field Examples,', SPE 57th An-1I
i
I
(1977).
I
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 : "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
,
Grant M A : "Broadlands Reservoir Structure," AppZied Mathematics Division1 D.S.I.R., Well
Grant M A, Iles D W : "A Simple Model of the Connection Between Easteb and Wester;
I
lands Geothermal Field," Proc., 2nd UN. Cogerence on the DeveZopdpr and Use of
Ohaaki Geothermal System, New Zealand: Their Chemistry, Distributioq and Corrosiv$
Nature," Trans., Geothermal Resources Council, Vol. 9, Part 2 (Aug. 19851, pp 245-250. '
Horner D R : "Pressure Build-Up in Wells," Proc. Third World Pet. Congress1 Section 11, E
dI
Repeat Gravity Measurements," Proc., 6th N.Z. Geothermal Workshop 1984), pp 195
200.
Jacob C E : "Correlation of Ground-water Levels and Precipitation on Long Island New York,
Reservoirs With and Without Fluid Banks," J. Pet. Tech. (Sept. 1972), pp 1t147-1156.
Fractured Reservoirs With Uniform Fracture Distribution," SOC. Pet. Eng. 9. (Dec. 1969),
pp 463-472.
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.: "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.
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
,
~
p 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).
~
Standing M B : Volumerric and Phase Behavior of Oil Field Hydrocarbon Systems, Reinhold1
,
Pub. Co., New York (1952). I
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.
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-
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.
1963), pp 245-255.
I
Young A : “Tidal Phenomena at Inland Boreholes Near Cradock,” Trans.,Ropal SOC.South
I
- 135 -
APPENDIX A
(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
TABLE B1
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 -
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
- 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
- 157 -
-
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 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 -
-
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 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 -
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 -
-
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
(~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-