Craig Geffen y Morse Teoria
Craig Geffen y Morse Teoria
fp = 7.758A/?(|) (7.115)
donde:
Vn = volumen poroso,
Bbl Interferencia entre
bancos de petróleo
A = área del yaci
miento, acres
Etapa 2
h = espesor neto pro . .
N =V pS0/B0 (7.116)
3. Cálculo de la razón de
movilidad, M, previa a la
ruptura del agua usando
los datos de flujo fraccio Producción de agua
nal y la siguiente ecua
ción: i Banco de agua |¡¡) Reglón de gas O Banco de petróleo
2
=}
£
i
<
<
UJ
ce
<
o
o
ce
ce
<
m
UJ
o
<
o
0.1 1 10
RAZÓN DE MOVILIDAD
Figura 7.17. Correlación de Craig, Geffen y Morse para determinar la eficiencia de barrido areal a
la ruptura6.
O
CO
z
UJ
z
o
<
CJ
¡ü
RAZÓN DE MOVILIDAD
(7.121)
donde: W¡f es el agua inyectada al momento del llene, en Bbl; y Sgl, la satura
ción de gas al inicio de la invasión, en fracción.
8. Se calcula el agua inyectada acumulada al momento en que ocurre la ruptura
de agua:
Wm = V o*EAiASwp,n-Swc) (7.122)
0,00708fthAp
q,=
. r V -o . re
— ln — +— ln —
r n /'
*v *r n i r
•'■m *i*
(7.123)
donde:
5,615W,
r„2 = (7.125)
T5,615W;. 1^
Por lo tamo, r, (7.126)
Luego:
2 2
r ¿ =r; (7.128)
(<7, )„+(<?,)„-,
[(* 7 1) p rom ] R (7.130)
(w a - w v ,
(7.131)
\ y 1 t J p ro m \ n
(7.132)
i=1
Tabla 7.9
Resumen de cálculos para la etapa 1
rei2
9 re¡
( q , l _____________________ t,
ció dejado por el gas (segunda etapa) dentro del arreglo de 5 pozos, la forma del beinco
de petróleo cambiará continuamente. Por esta razón, no es posible expresar matemáti
camente el comportamiento de la predicción. Como la duración de esta etapa es más
corta que la de las otras, se debe calcular la tasa de inyección de agua al final de la eta
pa 1 y al comienzo de la etapa 3, suponiendo que q, cambia linealmente entre estos
dos valores. Por lo tanto, el intervalo de tiempo entre la interferencia y el llene, vendrá
dado por:
wif- w u
At = (7.133)
0,5(<7„ +q,f)
Los valores de Wif yWn se conocen por los cálculos iniciales. La tasa de inyección
de agua a la interferencia, qu.corresponde a la tasa de inyección al final de la etapa 1.
La tasa de inyección al llene, q u, así como las tasas de inyección desde el llene hasta la
ruptura, se calculan mediante:
Q, =y*Q„ (7.134)
0,003541* h * (k J SwlrAp
(7.135)
ln — -0,619+0,5sp+0,5s,.
donde:
Qbase ~ tasa básica de inyección de agua, en BPD
d = distancia de la diagonal que une un inyector con el productor adyacen
te, pies
sp = factor de daño en el pozo productor, adimensional
s, = factor de daño en el pozo inyector, adimensional
(ka) Sw¡r = permeabilidad efectiva al petróleo a la saturación de agua irreducible; y
Ap = diferencia de presión de fondo del pozo entre el inyector y el productor,
después del llene en lpc.
La razón de conductancia, y, es un factor determinado experimentalmente por
Caudle y Witte20 que permite estimar el valor corregido de la tasa de inyección me
diante la ecuación 7.134. Esta razón de conductancia se presenta gráficamente en la
Figura 7.20 en función de la razón de movilidad, M, y de la eficiencia de barrido areal,
Ea■
248 Magdalena Paris de Ferrer
Figura 7.20. Razón de conductancia para un arreglo de 5 pozos (según Caudle y Witte20).
(7.137)
7.5- EtaDa
r
4: ComDortamiento desD
*
ués d e la ruDtura
a
del amia
o
ANpu - A ($ w Z E S wc ) * V p (7.139)
250 Magdalena París de Ferrer
T a b la 7.10
R e s u m e n d e c á lc u lo s p a r a la e ta p a 3
w,. e a Y Q,
(Ecuación 7.136) (Figura 7.20) (Ecuación 7.134)
<w x U a \, íq X
f— N
II
c
C
2 ^
a =<*1 W , - lV „ = (l)- (w a yv = —
B0 p B0
C^P),
W
v v ibl
Simplificando:
donde:
A£\
X= ^ (7.142)
Estos cálculos pueden expresarse con base en un barril de inyección total (o pro
ducción, ya que las tasas de inyección y producción se suponen iguales en condiciones
de yacimiento) fijando AW¡ =1. Luego:
^■VptSwZE ~Swc ) . .
--------
v v ibt
C7I43)
W - ( £ . . * . ( 7 . 1 4 4 )
El petróleo producido del área no barrida, A/Vpu, durante el tiempo en que se in
yectan AW¡ barriles de agua, puede estimarse usando la ecuación 7.145, ya que el valor
de X es conocido.
Nótese, sin embargo, que X depende, además, del aumento de la eficiencia de
barrido areal, AEA, lo cual ocurre como resultado de inyectar AW, barriles de agua.
Craig, Geffen y Morse6 encontraron experimentalmente que EA aumenta linealmente
W,
con el logaritmo de ttj—. Esta relación se muestra gráficamente en la Figura 7.21 y pue-
’W
v v ibt
de expresarse por:
'W , '
Ea =0,2749 ln +EAbl (7.146)
W bU
252 Magdalena París de Ferrer
w¡/w¡bt
Figura 7.21. Eficiencia de barrido areal en función de los volúmenes de agua inyectada (según
Craig1)-
f W, V1
X=0,2749 (7.148)
\
W tm
L /.
Luego:
0,2749(Su;„£. - S ) W,
AN (7.149)
w y
\vvitn
A N p s^ O -A / V pJ (7.150)
dL 1
(7.152)
dSw ' Sw2
w ¡h
= s wpbt (7.153)
EM>y P
W ^ G V ^ + C A IV ), después de (7.154)
la ruptura
Q i Q i b t ~^~Qdespuésde (7.155)
la ruptura
?
' después de - iJ * y * r - (7.156)
la ruptura w ih, v p C tA
'w.
u' vw ibt y
'd espu és de
(7.157)
¡a ruptura vP ea
De la ecuación 7.153:
W
y v ib
-QibtE a (7.158)
Vn
por lo tanto:
¿después de : Q i b l E Abt I
(7.159)
la ruptura
w, ,( w , '
'V„„ d
t
(7.160)
y¿ ib l 1 C /4
w,
VV, ,
Wr d w
Qi_
Q lb t
| \vyibt y
(7.161)
EÁ+0,2749 ln
VV,
Una solución de la ecuación 7.161 como una función de EAbl y ttt- se presenta en
y y ibt
T a b la 7.11
w,
E/tos,%
W:hl
50 51 52 53 54 55 56 57 58 59
1,0 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000
1,2 1,190 1,191 1,191 1,191 1,191 1,191 1,191 1,191 1,192 1,192
1,4 1,365 1,366 1,366 1,367 1,368 1,368 1,369 1,369 1,370 1,370
1,6 1,529 1,530 1,531 1,532 1,533 1,535 1,536 1,536 1,537 1,538
1,8 1,684 1,686 1,688 1,689 1,691 1,693 1,694 1,696 1,697 1,699
2,0 1,832 1,834 1,837 1,839 1,842 1,844 1,846 1,849 1,851 1,853
2,2 1,974 1,977 1,981 1,984 1,987 1,990 1,993 1,996 1,999 2,001
2,4 2,111 2,115 2,119 2,124 2,127 2,131 2,135 2,139 2,142 2,146
2,6 2,244 2,249 2,254 2,259 2,264 2,268 2,273 2,277 2,282 2,286
2,8 2,373 2,379 2,385 2,391 2,397 2,402 2,407 2,413 2,418 2,422
3,0 2,500 2,507 2,513 2,520 2,526 2,533 2,539 2,545 2,551 2,556
3,2 2,623 2,631 2,639 2,646 2,653 2,660 2,667 2,674 2,681 2,687
3,4 2,744 2,752 2,761 2,770 2,778 2,786 2,793 2,801 2,808 2,816
3,6 2,862 2,872 2,881 2,891 2,900 2,909 2,917 2,926 2,934 2,942
3,8 2,978 2,989 3,000 3,010 3,020 3,030 3,039 3,048 3,057 3,066
4,0 3,093 3,105 3,116 3,127 3,138 3,149 3,159 3,169 3,179 3,189
4,2 3,205 3,218 3,231 3,243 3,254 3,266 3,277 3,288 3,299 3,309
4,4 3,316 3,330 3,343 3,357 3,369 3,382 3,394 3,406 3,417 3,428
4,6 3,426 3,441 3,455 3,469 3,483 3,496 3,509 3,521 3,534 3,546
4,8 3,534 3,550 3,565 3,580 3,594 3,609 3,622 3,636 3,649
5,0 3,641 3,657 3,674 3,689 3,705 3,720 3,735
5,2 3,746 3,764 3,781 3,798 3,814 3,830
5,4 3,851 3,869 3,887 3,905 3,922
6,2 4,257
Valores ^ w —¡
de— ¡ para los cuales E A = 100%
W ibt
6,164 5,944 5,732 5,527 5,330 5,139 4,956 4,779 4,608 4,443
256 Magdalena París de Ferrer
wt E
u A bo %
W» 60 61 62 63 64 65 66 67 68 69
1,0 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000
1,2 1,192 1,192 1,192 1,192 1,192 1,192 1,193 1,193 1,193 1,193
1,4 1,371 1,371 1,371 1,372 1,372 1,373 1,373 1,373 1,374 1,374
1,6 1,539 1,540 1,541 1,542 1,543 1,543 1,544 1,545 1,546 1,546
1,8 1,700 1,702 1,703 1,704 1,706 1,707 1,708 1,709 1,710 1,711
2,0 1,855 1,857 1,859 1,861 1,862 1,864 1,866 1,868 1,869 1,871
2,2 2,004 2,007 2,009 2,012 2,014 2,016 2,019 2,021 2,023 2,025
2,4 2,149 2,152 2,155 2,158 2,161 2,164 2,167 2,170 2,173 2,175
2,6 2,290 2,294 2,298 2,301 2,305 2,308 2,312 2,315 2,319 2,322
2,8 2,427 2,432 2,436 2,441 2,445 2,449 2,453 2,457 2,461 2,465
3,0 2,562 2,567 2,572 2,577 2,582 2,587 2,592 2,597 2,601 2,606
3,2 2,693 2,700 2,705 2,711 2,717 2,723 2,728 2,733 2,738 2,744
3,4 2,823 2,830 2,836 2,843 2,849 2,855 2,862 2,867 2,873
3,6 2,950 2,957 2,965 2,972 2,979 2,986 2,993
3,8 3,075 3,083 3,091 3,099 3,107
4,0 3,198 3,207 3,216 3,225
4,2 3,319 3,329
A
*1i' 3,439
W
Valores de — - para los cuales E A = 100%
w ih,
4,285 4,132 3,984 3,842 3,704 3,572 3,444 3,321 3,203 3,088
w, F %
L‘ A bo
W,b, 70 71 72 73 74 75 76 77 78 79
1,0 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000
1,2 1,193 1,193 1,193 1,193 1,193 1,193 1,193 1,194 1,194 1,194
1,4 1,374 1,375 1,375 1,375 1,376 1,376 1,376 1,377 1,377 1,377
1,6 1,547 1,548 1,548 1,549 1,550 1,550 1,551 1,551 1,552 1,552
1,8 1,713 1,714 1,715 1,716 1,717 1,718 1,719 1,720 1,720 1,721
2,0 1,872 1,874 1,875 1,877 1,878 1,880 1,881 1,882 1,884 1,885
2,2 2,027 2,029 2,031 2,033 2,035 2,037 2,039 2,040 2,042 2,044
Inyección de agua y gas en yacimientos petrolíferos 257
w, E A b lt %
W¡bl 70 71 72 73 74 75 76 77 78 79
2,4 2,178 2,180 2,183 2,185 2,188 2,190 2,192 2,195 2,197
2,6 2,325 2,328 2,331 2,334 2,337 2,340
2,8 2,469 2,473 2,476 2,480
3,0 2,610 2,614
W
Valores de — - para los cuales E A = 100%
W ibt
2,978 2,872 2,769 2,670 2,575 2,483 2,394 2,309 2,226 2,147
wi E,1b t i %
/u
W* 80 81 82 83 84 85 86 87 88 89
1,0 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000
1,2 1,194 1,194 1,194 1,194 1,194 1,194 1,194 1,194 1,194 1,194
1,4 1,377 1,378 1,378 1,378 1,378 1,379 1,379 1,379 1,379 1,379
1,6 1,553 1,553 1,554 1,555 1,555 1,555 1,556 1,556 1,557 1,557
1,8 1,722 1,723 1,724 1,725 1,725 1,726 1,727 1,728
2,0 1,886 1,887 1,888 1,890
o o n
W
Valores de — - para los cuales E , = 100%
W ibt
2,070 1,996 1,925 1,856 1,790 1,726 1,664 1,605 1,547 1,492
W, E fibtt %
/ü
w lhl 90 91 92 93 94 95 96 97 98 99
1,0 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000
1,2 1,194 1,195 1,195 1,195 1,195 1,195 1,195 1,195 1,195 1,195
1,4 1,380 1,380 1,380 1,380 1,381
1,6 1,558
Q,
Una vez determinado ~r—de las tablas, se puede calcular Q, el cual al usarse con
la curva de flujo fraccional permite definir fw2. Finalmente se calculan fo2 y A/Vps usan
do las ecuaciones 7.151 y 7.150, respectivamente.
El incremento del agua producida con base en un barril, AWps, se obtiene de:
1-A/V -A N pu
(7163)
RAPcn=RAP I 2- (7.164)
q ^ p s +A/vpu)
Qo = -------- g --------- (7.165)
<7,(l-A/Vps-AyVpu)
Q w = --------- R---------- (7.166)
Kp[
oí
£ ,n (-S ; p- 5 im
lu u
^. )■ -' 5« sí]
i i ..
(7.168)
Q,
Jk (7 .170)
'■” r o J Swc ^ ir
W , Qi
w i ea
w , bl (Ecuación 7.146) Qibl
(Tabla 7.10)
Q,
dfw $w 2 fo 2
w r v ib t U J u,
1,0/(5)
SLP X AN pu AN ps
RAP r a p cn N p
M
(Ecuación 7.163) (Ecuación 7.164) (Ecuación 7.167) (Ecuación 7.170)
Inyección de agua y gas en yacimientos petrolíferos 261
Y Q,
(Figura 7.20) (Ecuación 7.134)
<■*«»
A W, Aj _ (20)
II
M
(19)
Qo Qu Wp
(Ecuación 7.165) (Ecuación 7.166) (Ecuación 7.169)
3. Los datos de permeabilidades relativas son los mismos para todos los estratos.
4. Las tasa de inyección y producción asociadas con cada estrato son propocio-
nales ak*h.
Si se ha realizado la predicción del estrato 1 usando los cálculos anteriores, el
tiempo requerido para inyectar el mismo número de volúmenes porosos de agua hasta
el estrato n, tiempon, tn, como fue inyectado para el estrato 1 durante el tiempo /,, será:
(7.171)
N -N (7.172)
W>),
W =W (7.173)
p" P' W ,
262 Magdalena Paris de Ferrer
Win=W n T ^ (7-174)
m i
q^ q« J m [ (7175)
f 7 1 7 C *
Í 7 I 7 6 )
f 7 1 7 7 ^
Procedimiento:
1. Predecir el comportamiento del estrato 1 usando el método previamente des
crito.
2. Representar los datos de N p, W¡, Wp, q, , qa y qwen función de tiempo para el
estrato 1.
3 ^
‘ Obtener los valores de —, tyh y kh para todos los estratos.
Tabla 7.13
Resumen de las predicciones realizadas en el estrato 1
5. Considérense luego los estratos remanentes del yacimiento, por ejemplo el es
trato n. Correspondiendo al tiempo í, , seleccionar en el paso 4 para el estrato 1
y usar la ecuación 7.171 para calcular los tiempos t„ que se necesitan para el
estrato n, con el mismo volumen poroso de agua inyectada en los dos estratos.
Estos cálculos se ilustran en la Tabla 7.14.
Inyección de agua y gas en yacimientos petrolíferos 263
T a b la 7.14
Resumen de las predicciones realizadas en el estrato n
Estrato 1 Estrato n
co.
ca (U
(Os (O s
P ro b le m a s79-34'38
donde: S
Parsons, si se considera que el yacimiento está formado por 8 estratos con las
siguientes propiedades:
Otros datos:
Área productiva, acres-pie 100
Saturación de agua connata, % 24
Saturación de gas inicial, % 17
Porosidad promedio, % 19
Espesor neto de arena, pies 10,0
km en el banco de petróleo, fracción 0,80
km detrás del banco de petróleo, fracción 0,20
Factor volumétrico inicial del petróleo, BY/BN 1,215
Factor volumétrico del petróleo, BY/BN 1,073
Temperatura del yacimiento, °F 121
Saturación de petróleo residual, % 22,5
Recuperación primaria, BN/acre-pie 111,5
Viscosidad del agua, cp 0,82
Viscosidad del petróleo, cp 4,34
Tasa de inyección, BPD 100
Ea,% 85
Factor del agua en la formación, BY/BN 1,0
Realizar la predicción aplicando el Método de Stiles.
¿S2
J
km= o m - s wDy
donde S,„n
WD
=1 C _ O
1 °o r °w i
¡ :
4. En el sistema mostrado en la Figu- ¡ j
ra 7.23, la distancia entre los po- O j O j O
zos inyectores es de 1.000 pies; el i ¡
radio de drenaje es de 10 pies y el .............. 4 "............. Y ...............
radio de ios pozos 0,25 pies. La ¡ ¡
permeabilidad efectiva del petró- O ¡ O ¡ o
leo a la saturación de agua irredu
cible es 100 md; la saturación de
agua irreducible es 10% y la satu
...... t ...........t .......
: i
ración promedio de agua detrás o : o o
del frente de invasión es 56,3%. i i
Se conoce, además, la siguiente Figura 7.23. Localización de los pozos en el
información: yacimiento (según Willhite*).
Profundidad, Permeabilidad,
pies md
2.050-2.051 35
2.051-2.052 51
2.052-2.053 27
2.053-2054 116
2.054-2.055 60
2.055-2.056 237
2.056-2.057 519
2.057-2.058 98
2.058-2.059 281
2.059-2.060 164
Estíme qué porcentaje de la capacidad de flujo total contienen los 3 pies más
permeables de la formación.
km=0,621 ( 1 - S ^ )1’638
Inyección de agua y gas en yacimientos petrolíferos 269
K Espesor,
md m
187 1
63 1
310 1
k
K rw
=0 78£372
v >Í O é * w D
w» w, /wibl ea 0, $w2
Bbl
1,00 0,336 0,664
68.658 1,16 0,80 0,387 0,670
77.906 1,31 0,83 0,435 0,675
87.155 1,47 0,87 0,482 0,679
96.403 1,62 0,89 0,527 0,683
105.652 1,78 0,92 0,571 0,686
114.900
124.149 2,09 0,96 0,656 0,691
Inyección de agua y gas en yacimientos petrolíferos 271
fo2 wp N P, Q, RAP
VP VP Bbl
Estime:
a) La razón de movilidad y el agua inyectada a la ruptura
b) La tasa de inyección para un valor de agua inyectada de 114.900 Bbl
c) La tasa de producción de petróleo cuando se considera como límite econó
mico una RAP = 4
d) Complete la tabla usando el método más apropiado.
11. Se está llevando a cabo una inyección de agua en un yacimiento formado por
un arreglo de 5 pozos. Se conocen además los siguientes datos:
Área del arreglo, km2 12,5
Factor volumétrico del agua, m3/m3a CN 1
Factor volumétrico del petróleo, m3/m3a CN 1,25
Espesor de la formación, m 10
Permeabilidad relativa al petróleo a S w¡, fracción 0,63
Permeabilidad relativa al agua a Swp, fracción 0,21
Porosidad, % 20
Viscosidad del petróleo, cp 5
Viscosidad del agua, cp 0,8
Saturación de petróleo residual, fracción 0,3
Saturación de petróleo inicial, fracción 0,7
Tasa de inyección, m3/día 500
Calcule:
a) Petróleo recuperado a la ruptura
Magdalena París de Ferrer
12. Se está llevando a cabo una inyección de agua en un yacimiento formado por
un arreglo de 9 pozos en 10 acres. Los pozos serán operados de tal forma que
los pozos laterales y de los vértices tengan la misma tasa. La porosidad es 0,20 y
la razón de movilidad es 1. El yacimiento tiene una saturación de petróleo ini
cial de 0,70 y la saturación residual es de 0,26.
Considerando un espesor de 1 pie y el Ba =1, complete la siguiente tabla utili
zando las Figuras 7.25, 7.26 y 7.27.
K Ea N p, N„, Q„ fw RAP
___________________________ Bbl__________VP__________ VP___________________________
0,53
0,60
0,70
0,80
0,90
1,00
1,20
1,40________________________________________________________________________________
RAP = 1
fsw = 0,95
Algunas fórmulas útiles
w,
Vri~Vp{\-Swir- S or)
f _ qJcu +2<7A,
Qc +2<7s
RAP = 1
RAP= 1
fsw = 0,95
RAZÓN DE MOVILIDAD
Figura 7.27. Eficiencia de barrido en función de M en un arreglo de 9 pozos para diferen
tes fsw (según Willhite38).
274 Magdalena París de Ferrer
Referencias bibliográficas
1. Craig, F.F., Jr.: The Reservoir Engineering Aspects o f Waterñooding, Monograph Series, SPE,
Richardson, TX (1971) 3.
2. Buckley, S.E. y Leverett, M.C.: Mechanisms of Fluid Displacement in Sands, Trans., AIME
(1942) 146, 107-116.
3. Dykstra, H. y Parsons, R.L.: The Prediction of Oil Recovery by Water Flood, Secondary Reco
very o fO il in The United States, API, N e w York City (1950) 160-174.
4. Johnson, C.E. Jr.: Prediction of Oil Recovery by Waterflood - A Simplified Graphical Treat-
ment of The Dykstra-Parsons Method, Trans., AIME (1956) 207,345-346.
6. Craig, F.; Geffen, T. y Morse, R.A.: Oil Recovery Performance of Pattem Gas or W ater Injec
tion Operations from Model Tests, Trans., AIME (1955) 204,7-15.
10. Scientific Software Intercomp (SSI): Notebook oflm p roved Oil Recovery by Waterñooding
and Gas Injection, Denver (1981).
11. Roberts, T.G.: A Permeability Block Method of Calculating a Water Drive Recovery Factor,
Pet. Eng. (1959) 31, B45-48.
12. Higgins, R.V. y Leighton, A.J.: Computer Prediction of W ater Drive of Oil and Gas Mixtures
through Irregularly Bounded Porous Media-Three Phase Flow. JPT (Sept. 1962) 1048-1054.
13. Rapoport, L.A., Carpenter, C.W . y Leas, W.J.: Laboratory Studies of Five-Spot Waterflood
Performance, Trans., AIME (1958) 213,113-120.
14. Yuster, S.T. y Calhoun, J.C., Jr.: Behavior of W ater Injection Wells, Oil Weekly (Dic. 18 y 25,
1944) 44-47.
15. Suder, F.E. y Calhoun, J.C., Jr.: Waterflood Calculations, Drill and Prod. Prac., API (1949)
260-270.
16. Prats, M., Matthews, C.S., Jewett, R.L. y Baker, J.D.: Prediction of Injection Rate and Produc
tion History for Multifluid Five-Spot Floods, Trans., AIME (1959) 216, 98-105.
17. Felsenthal, M. y Yuster, S.T.: A Study of the Effect of Viscosity in Oil Recovery by Waterfloo-
ding, trabajo N° 163-G presentado en la reunión del SPE W est Coast en Los Angeles (Oct.
25-26,1951).
18. Muskat, M.: Flow o f Hom ogeneous Fluids Through Porous Systems, J.W. Edwards, Inc., Ann
Arbor, Mich. (1946).
19. Hurst, W .: Determination of Performance Curves in Five Spot Waterflood, Pet. Eng. (1953)
25, B40-46.