Mecanica de Yacimientos
Mecanica de Yacimientos
Apuntes
Mecnica de
Yacimientos
Realizado por:
Irvin Garduza Antonio
INDICE
INTRODUCCIN A LA INGENIERA DE YACIMIENTOS............................................5
DEFINICIN DE RESERVAS................................................................................... 6
CLASIFICACIN DE YACIMIENTOS......................................................................12
PROPIEDADES DE LAS ROCAS, GAS, ACEITE Y AGUA........................................14
ECUACIN GENERAL DE BALANCE DE MATERIA................................................59
YACIMIENTO CON ENTRADA DE AGUA O INYECCIN DE AGUA..........................71
YACIMIENTO DE GAS......................................................................................... 83
YACIMIENTOS DE ACEITE BAJO SATURADO......................................................113
LEY DE DARCY................................................................................................. 123
BIBLIOGRAFIA.................................................................................................. 135
Yacimiento
Acumulacin de hidrocarburos tcnica y econmicamente explotable.
Yacimiento de gas
Los yacimientos de gas natural se definen por la ubicacin de la presin y la
temperatura inicial en un diagrama de presin, temperatura, se subdividen en
yacimientos de gas seco, gas hmedo y gas y condensado.
- Condensacin retrgrada
- Gas seco
Yacimientos de aceite
- Con capa de gas
- Con gas disuelto
DEFINICIN DE RESERVAS
CLASIFICACION DE RESERVAS
Estimacin Alta
Estimacin Central
Estimacin Alta
Reservas
Estimacin Baja
1P
Probada
2P
Probada
+
Probable
3P
Probada
+
Probable
+
Produccin
Estimacin
Central
Recursos Contingentes
Estimacin Baja
No Recuperable
Recurso Prospectivo
No Recuperable
Incertidumbre
Posible
Recursos contingentes
Son aquellas cantidades de hidrocarburos que son estimadas, a una fecha dada,
las cuales potencialmente son recuperables de acumulaciones conocidas pero que
bajo las condiciones econmicas de evaluacin a esa fecha, no son consideradas
comercialmente explotables.
Los recursos contingentes pueden incluir, por ejemplo, acumulaciones donde no
exista un mercado para comercializar lo producido, o donde la recuperacin
deseada de hidrocarburos depende del desarrollo de nuevas tecnologas, o bien,
donde la evaluacin de la acumulacin no se ha concluido.
Reservas
Originales
(Recurso
Econmico)
Reservas
Probadas
Originales
Produccin
Acumulada
Reservas no
Probadas
Reservas
Probadas
Desarroladas
Reservas
Probables
Reservas
Posibles
No Recuperable
Reserva
Se definen como aquellas cantidades de hidrocarburos que se prev sern
recuperadas comercialmente de acumulaciones conocidas a una fecha dada.
Clasificacin de las reservas
Todas las reservas estimadas involucran algn grado de incertidumbre. La
incertidumbre depende principalmente de la cantidad y calidad de la informacin
7
las que sustentan los proyectos de inversin, y por ello la importancia de adoptar
definiciones emitidas por la SEC.
Reservas probadas desarrolladas
Son aquellas reservas que se espera sean recuperadas de pozos existentes,
incluyendo las reservas detrs de la tubera de revestimiento, que pueden ser
8
Reservas Probables o 2P
Son aquellas reservas no probadas en donde el anlisis de la informacin
geolgica y de ingeniera del yacimiento sugiere que son ms factibles de ser
comercialmente recuperables, que de lo contrario. Si se emplean mtodos
probabilistas para su evaluacin, existir una probabilidad de al menos 50 por
9
CLASIFICACIN DE YACIMIENTOS
10
11
12
Porosidad :
El porcentaje de volumen de poros o espacio poroso, o el volumen de roca
que puede contener fluidos. La porosidad puede ser un relicto de la depositacin
(porosidad primaria, tal como el espacio existente entre los granos que no fueron
completamente compactados) o puede desarrollarse a travs de la alteracin de
las rocas (porosidad secundaria, tal como sucede cuando los granos de feldespato
o los fsiles se disuelven preferentemente a partir de las areniscas). La porosidad
puede generarse a travs del desarrollo de fracturas, en cuyo caso se denomina
porosidad de fractura. La porosidad efectiva es el volumen de poros
interconectados, presentes en una roca, que contribuye al flujo de fluidos en un
yacimiento. Excluye los poros aislados. La porosidad total es el espacio poroso
total presente en la roca, sin importar si contribuye o no al flujo de fluidos. Por
consiguiente, la porosidad efectiva normalmente es menor que la porosidad total.
Los yacimientos de gas de lutita tienden a exhibir una porosidad relativamente
alta, pero la alineacin de los granos laminares, tales como las arcillas, hace que
su permeabilidad sea muy baja.
Porosidad=
Ejemplos:
r=1.16 cm
13
numero de esferas : 27
4
4
volumen de laesfera= r 3= (1.16)3=6.538 cm3
3
3
volumen de la esfera = 6.538
cm 3
(6.96)3=337.153 cm 3
337.153176.533 160.62
=
=0.457 ( 100 )=47.64
337.153
337.153
r= 0.25cm
Nmero de esferas : 216
4 3 4
3
3
volumen de laesfera= r = (0.25) =0.065 cm
3
3
volumen de la esfera = 0.065
cm 3
(3)3=27 cm 3
2714.137 12.863
=
=0.4764 ( 100 )=47.64
27
27
14
Volumen de la esfera=( 8 )
4
3
3
R =33.51 R
3
R3
Vporoso 21.92 R3
=
=0.3954 ( 100 ) =39.54
Vtotal 55.426 R3
Volumen de la esfera=( 8 )
4
3
3
R =33.51 R
3
15
45.25 R
=
-33.51 R
= 11.74 R
Vporoso 11.74 R3
=
=0.2594 ( 100 )=25.94
3
Vtotal
45.25 R
3
( 8) 4 R
3 2
Vporoso VtotalVesferas
=
=
=1
=0.4764 ( 100 )=47.64
Vtotal
Vtotal
R 3
( 4)
2
( )
( )
DIAGRAMA DE FASES
16
COMPRESIBILIDAD DE LA ROCA
C=
ISOTRMICA
1 dV
V dP
Donde:
C= compresibilidad de la roca
V= volumen
P= presin
Cf =
97.3200 10 6
1.42854
(1+55.8721)
Cf =
0.853531
(1+2.47664 106 )0.92990
SATURACIN DE FLUIDOS
18
Saturacion de aceite So
Saturacion de gas Sg
So+Sg+Sw= 1 100%
Saturacion de agua Sw
Sw=
So=
Aceit
e
Agua
Sw= 50%
So=
Sw= 50%
50%
So= 50%
Si la roca es mojada por agua es por la afinidad ionica con la que el gua se pega
con el aceite, en cambio con el aceite se pegara con el agua.
Medicin: En base al laboratorio, en un ncleo, en base a registros electricos.
19
V=
Ft
R=10.73
n=lb-mol
T=R
Ejemplo:
Ft
Ft
de etano y se
V= 500
Ft 3
P= 100 psia
T= 100 F + 460= 560 R
Mw= 30.07 mol
lb
lbmol
PV=nRT
PV
=n
RT
(100)(500 Ft ) 50000
=
=8.32 moles
(10.73)(560) 6008.8
20
lb
lbmol )= 250.21
lb
moleculas
SCF=?
Ft 3
( 8.32mol lbmol ) 10.73 psia
(520 R)
lbmolR
V= nRT =
=3158Cf
Psi
14.7 Psi
Scf
=3157 scf
lbmol
GRAVEDAD ESPECFICA
m
Densidad== v
P=0.052 h
Definiendo que la densidad =
m
v
Densidad
21
m nMw
=
v
v
PV
Mw
RT
PMw
=
=
V
RT
PV =nRT
aire=
n=
PV
RT
P 28.97
RT
PMw
g
RT
Mw
=
=
g= aire P 28.97 28.97
RT
z=
z
FactorVactual=Videal
de compresibilidad,
factor de desviacin z
Videal=
Vactual
z
PV =nRT
=nRT
( Vactual
z )
22
PVactual=znRT
Ejemplo:
Un gas tiene un volumen de 364.6 pies3 a 213 F y 3250 psia, a 14.8 psia y 82 F
tiene un volumen de 70860 pies3. Calcular el factor de desviacin.
A condiciones estndar el gas se comporta como ideal.
213F+460= 673R
82F+460= 542R
n=
PVactual=znRT
n=
PVactual
zRT
n=180.32lbmol
z=
PVact
nRT
23
z=
(3250)(364.6)
=0.9099
(180.32)(10.73)(673)
24
Ejemplo:
Calcular el factor de desviacin z del Campo Bell, el cual tiene las siguientes
caractersticas:
Gravedad especifica= 0.665
Temperatura de yacimiento= 213F=673R
Presin de yacimiento= 3250 psia
Ppc=668.09
Tpc= 169.2+349.5 (0.665)-74.0 (0.665)2
25
Tpc=368.89
Ppr=
3250 psia
=4.8
668.09 psia
Tpr=
673 R
=1.8
368.89 R
26
Drachuk Abou-Kassem
Donde:
r=0.27
Ppr
( z Tpr)
Mtodo de la Secante:
28
x n+1=x nf n
x nx n1
f nf n1
x=z
Mtodo de Newton-Raphson:
x n+1=x n
fn
F'n
f(z)
29
11
r )
Component
e
Fracc.
Mol
Metano
Etano
Propano
Butano
Pentano
CO2
N2
0.8612
0.0591
0.0358
0.0172
0.005
0.001
0.0207
1.000
3
Mw
(lb/lbmol)
16.04
30.07
44.09
58.12
72.15
44.01
28.02
Pc
(psia)
Tc
(R)
673
708
617
550
490
1070
492
343
550
666
766
846
548
227
(2)*(3)
Mwi
(lb/lbmol)
13.81
1.78
1.58
1.00
0.36
0.04
0.58
19.15
(2)*(4)
(2)*(5)
Pci
(R)
Tci
(R)
579.59
41.84
22.09
9.46
2.45
1.07
10.18
666.68
295.39
32.51
23.84
13.18
4.23
0.55
4.70
374.40
Gravedad Especfica:
g=
Mw 19.15
=
=0.661
28.97 28.97
3250
=4.87
666.68
Ty= 673R
Tpr=
673
=1.80
374.40
30
Tpc= Tpc-
P' pc=
Ppc T pc
( Tpc+ B ( 1B ) )
0.97%
154< P (psia)< 7026
40< T (F) < 300
CO2 < 54.4% mol
H2S < 73.8% mol
Bg=
Vol . gas @c . y . 3
bbls
ft /SCF
Vol. gas @ c . s .
SCF
31
Deduccin de la ecuacin:
PV=znRT
V=
znRT
P
[email protected].=14.71 psia
[email protected].=60F=520R
znRT
P
P @c . s . ZT ( 14.71 ) zT
zT
Bg=
=
=
=0.02829
nRT @ c . s . T @ c . s . P
520 P
P
P@ c . s .
znRT
P
P @c . s . ZT ( 14.71 ) zT
zT
Bg=
=
=
=0.00504
nRT @ c . s . T @ c . s . P
520 P
P
P@ c . s .
32
Ejemplo:
Calcula el factor de volumen del gas dadas las siguientes condiciones:
Bg=?
Z=0.91
T= 673R
P=3250 psia
Bg=0.02829
zT
P
Bg=0.02829
( o .91 )( 673 R )
3250 psia
Bg=0.02829
612.43 R
3250 psia
Bg=0.02829 ( 0.18844 )
Bg=5.3309676 103
ft 3
SCF
Bg=0.00504
( 0.91 )( 673 R )
3250 psia
Bg=0.00504
612.43
3250
Bg=0.00504 ( 0.18844 )
Bg=9.497376 104
bls
SCF
33
m
v
PV=znRT
V=
znRT
P
Mw= g 28.97
g=
g=
g 28.97
znRT
P
28.97 g P
zRT
lb
= ft 3
Ejemplo:
Calcule la densidad del gas con los siguientes datos:
g = 0.665
P= 3250 psia
T= 673R
34
R= 10.73
Z= 0.91
g=
28.97(0.665)(3250)
62611.4125
lb
=
=9.52 3
( o .91)(673 R)( 10.73) 6571.3739
ft
z nRT
P
V =C
Ct
z
P
[ ]
dV d
z
=
C
dP dP
P
] [ ]
dz
dP
dz
Pz
Pz
dV
dP
dP
dP
=C
= nRT
dP
P2
P2
dV nRT P dz znRT
=
2
dP
P2 dP
P
35
dV nRT dz znRT 1
=
dP
P dP
P P
dV znRT 1 dz znRT 1
=
dP
P z dP
P P
dV
1 dz
1
1 dz 1
=v
v =v
dP
z dP
P
2 dP P
1 dV 1 dz 1
=
=
V dP z dP P
De la deferencia de Compresibilidad:
Cg=
1 dV
V dP
Cg=
1 dz 1
z dP P
1 1 dz
Cg=
( Psia1 )
P z dP
36
P(psia)
1000
Cg=
dz
dP
0.83
127 106
1
1
(127 106 )
1000 0.83
P(psia)
dz
dP
2500
0.73
4000
0.90
110 1 0
37
Cg=
1
1
( 0)
2500 0.73
Cg=
1
1
( 110 106 )
400 0.9
Viscosidad
CO2= 5% mol
1= 0.0117 cp.
0.0003 cp.
T= 0.0117+0.0003= 0.0120 cp.
Tpr=
672
=1.6
420
Ppr=
2680
=4
670
=( 104 ) k e ( x )
Donde:
3
=1.4935 10
k=
PMw
zT
x=3.5+
986
+0.01 Mw
T
P; psia
T; R
M; lb-lb mol
39
y=240.2 x
P=2680 psia
T= 212F
g= 0.9
Ppc= 756.8-131.0(0.9)-3.6(0.9)2= 635.984 psia
Tpc= 169.2+349.5(0.9)-74.0(0.9)2= 423.81R
Ppr=
2680
=4.2139 psia
635.984
Tpr=
672
=1.58 R
423.81
Z= 0.82
Mw= 0.9 x 28.97= 26.073
=1.4935 103
=0.189386
( 2680 )( 26.073 )
69875.64
=1.4935 103
551.04
( 0.82 )( 672 )
g
cc
1.5
9.4 +0.02 ( 26.073 ) ] ( 672 )
[
172834.125
k=
=
1376.387
[ 209+19 ( 26.073 ) + ( 672 ) ]
40
x=3.5+
= 125.57
986
+0.01 ( 26.073 )
672
X= 5.228
Y= 2.4 - 0.2 (5.228)
Y= 1.3544
4
[( 5.228) (0.189386 )
g =( 10 ) ( 125.57 ) e
1.3544
g =0.02174279229 cp .
g =0.0217+0.0192=0.0025 cp .
Rel . Dif .=
0.0025
=0.115=11.5
0.0217
Separacin Flash
42
Rso= g
P
y
18 ( 10 )
1.204
Donde:
y g =0.00091T 0.0125 o
T= F
o= F
P= psia
Rso=
SCF
STB
0.0125 o
P 10
Rso= g
18 100.00091T
1
0.83
Aceite: 18API
P= 2500 psia
T= 212F
43
g=0.9
yg= 0.00091 (212F) 0.0125 (18)
yg= - 0.03208
Rso=0.9
2500 psia
0.03208
18 ( 10 )
Rso=373.79
1.204
SCF
STB
1
0.83
bbl
( STB
)
44
45
Bo=0.972+0.000147 F
F=Rs g
o
0.5
( )
o=
1.75
+1.25 T
141.5
131.5+ API
Bo=Bobe Co ( Pb P)
T= F
46
SCF
STB
Bg=0.00504
zT bbls
P SCF
Bg=0.00504
0.89 ( 620 )
bbs
=0.002317
1200
SCF
( )
Conociendo z=0.89
Bg=0.02829
( 0.89 )( 620 )
ft 3
=0.013
1200
SCF
47
Vg=230 0.013=2.99 ft 3
2.99 ft
=0.533 bls.
5.615
1 dv
v dP
Bo
Pb
Pi
Co=
48
Ecuacin de Egbogah
Beggs y Robinson
o= A odB
Donde:
A=10.715 ( Rso+100 )0.515
B=5.44 ( Rso+150 )0.338
o= ob
P
Pb
( )
P)
49
Mc Cain
Bw =( 1+ V WT ) ( 1+ V )
bbl
( STB
)
Donde:
2
50
Correlacion de Mc Cain
Rsw
=10(0.0840655 ST
R swp
0.285854
Donde:
S= Salinidad % en peso
Rswp= relacion gas disuelto-agua pura (scf/STB)
Rswp= A + BP + CP2
A=8.158396.12265 102 T +1.91663 104 T 22.1654 107 T 3
51
Compresibilidad isotrmica Cw
Cw=
Donde:
CNaCl= Salinidad g NaCl/litro
T= F
Cw=
1 Bw
Bg Rswp
+
Bw P T Bw P
R swp
=B+2 Cp
P T
Correlacion de Mc Cain:
52
w 1= AT
Donde:
A=109.574840564 S+ 0.313314 S 2 +8.72213 103 S 3
B=1.12166+2.63951 102 S6.79461 104 S25.47119 105 S 3+ 1.55586 106 S 4
T= F
S= Salinidad % peso slido
w
=0.999 4+ 4.0295 105 P+3.1062 109 P2
w1
Ejercicios:
53
V=
V=
nRT
P
b) 14.7 psia
32F
V=
nRT
P
V=
c) 14.7 psia
10 oz@ 80F
Cambiar oz a unidades de presin para eso le sumamos 0.65+ 14.7= 15.35 psia
nRT
V= P
(1)( 10.73)(80+ 460)
=377.47 ft 3
V=
15.35
d) 15.025 psia
60 F
nRT
V= P
V=
3
2) Un tanque tiene 500 ft , 10 lb de metano y 20 lb de etano a 90F
m
n= M
m= peso (lbs)
lb
M= peso mol ( lbmol )
n= lb-mol
10
20
+
=1.29 lbmol
Nt=nm + ne= 16.04 30.07
P=
55
metano=
0.623
=0.484
1.29
0.6651
=0.516
etano= 1.29
lb
lbmol
56
que
contien aire a presin atmosfrica 14.7 psia y 75F. cual ser la presin final
del tanque sellado cuando todo el hielo seco se ha enfriado a 45F.
PV=nRT
lb
co2= 44. 01 lbmol
(14.7)(50)
PV
=
n= RT (10.73)(75+ 460) =0.128
10
lb
=0.227
nco2= 44.01
lbmol
57
nt= naire+nco2
nt= 0.128 + 0.227=0.355 lb-mol
PV=nRT
nRT (O .355 ) ( 10.73 )( 45+460 )
=
=38.47 psi
P= V
50 ft 3
lb
lbmol
Mg=yiM1
18.83=ym(16.04)+ye(30.07)
Ym+ye=1
58
18.83=(1-ye)(16.04)+ye(30.07)
18.83=16.04-16.04ye+30.07ye
18.83-16.04=-16.04ye+30.07ye
2.79=14.03ye
2.79
= ye
14.03
Ye=0.198 lb-mol
Ym=1-ye
Ym=1-0.198
Ym=0.802 lb-mol
Peso en volumen
Ye=19.8%
Ym=80.2%
Component
e
Metano
Etano
Yi
0.802
0.198
Masa Total=
M1
Masa Mi
Peso Wi
Porcentaje en
peso
16.04
30.07
12.864
5.953
18.817 lb
0.6836
0.3164
1
68.36%
31.64%
Conclusin:
El etano es ms pesado y el metano es ms ligero
59
3
6) Un tanque contiene 50 ft
con otro tanque que contiene gas a 25psia y 50f. cuando la valvula de en
medio de los dos es abierta, la presion se estabiliza a 35 psia a 50F. Cul
es el volumen en el segundo tanque?
N1
N2
3
V1=50 ft
V2=?
Pf= 35psia
P1= 50 psia
P2= 25 psia
tf= 50F
T1= 50 F
T1= 50 F
60
PV=nRT
P 1 V 1 P 2 V 2 Pf (VI +V 2)
+
=
RT 1
RT 2
RTF
P1V1+P2V2=Pf(VI+V2)
(50)(50)+25V2=35(50+V2)
2500+25V2=1750+35V2
2500-1750=35V2- 25V2
750=10V2
3
V2= 75 ft
dolares
MCF
61
(14.4)(1000) 15.025(V 2)
=
80+460
60+460
80 15.025 V 2
=
3
520
13866.66=15.025V2
13866.66
=V 2
15.025
3
V2=922.91 ft
Tanque 1
Tanque 2
V1
V2
Aceite
Agua
Fuente de Agua
Tanque 1
P1=Pi
T 1 =T
V w 1=WB wi=WB twi
V g 1=GB gi
W; volumen de agua inicial en STB
V o 1=NB oi =NB ti
N; volumen de aceite inicial en STB
Tanque 2
P2=0
T 2 =T
V w 2=0
V o 2=0
V g 2=0
63
A ( X)
V =
Aceite
Aceite
Agua
Agua
Fuente de agua
V total=V t =V 1V +V 2
Nuevas condiciones de presin
La misma temperatura
GB g
WBtw +We
64
Expansin
del aceite
inicial y gas
asociado
Np.-
Expansin
del gas
inicial
Expansi
n del
Entrad
a de
Variaci
n del
volume
n
N ( Bt Bt i ) +G ( BgBgi ) +W ( Bt w Bt wi ) +We+ V =V 2
N ( Bt Bt i ) +G ( BgBgi ) +W ( Bt w Bt wi ) +We+ V =Np Bo+ Bw+ ( GpNp Rs Rsw ) Bg
Mecanismos:
1. Expansin del aceite
2. Expansin del gas
3. Disminucin del volumen poroso en la zona aceite-gas y agua congnita
por compresibilidad
4. Empuje por entrada de agua
66
NOMENCLATURA
rb/scf
rb/scf
rb/stb
rb/stb
Bt.-
rb/stb
Bti.-
rb/stb
rb/stb
rb/stb
Bw.-
rb/stb
rb/stb
Cf.-
psi -1
psi -1
G.-
scf
GI.-
scf
Gp.-
scf
m.Relacin del Volumen del casquete de Gas Inicial con el Volumen de Aceite
Inicial
adimensional
N.-
stb
stb
67
P.-
Presin de Yacimiento
psia
Pi.-
psia
scf/ stb o
scf/stb o
scf/stb
scf/stb w
scf/stbw
adimensional
adimensional
adimensional
T.-
Temperatura
rb
rb
rb
rb
rb
rb
W.-
rb
WI.-
Stb w
Stb w
68
1 dVp
Vp dP
= Cf dP
dVp
Vp Vp
P
i
ln Vp
Vpi
Cf Vp
Pi
ln VplnVpi = Cf Vp
Pi
ln
Vp
= Cf dP
Vpi P
( )
Cf dP
Vp
=e
Vpi
Pi
[ ]
Vp=Vp e
P
Cf dP
Pi
[ ]]
Vp=Vp Vp=Vp [ 1e
P
Cf dP
Pi
69
1ra. Suposicin:
Swi es la misma en la zona de aceite y zona de gas.
Voi=Vpoi Soi Vpoi=
Voi
So i
Vgi=Vpgi Sg i Vpgi=
Vgi
Sgi
Dado que:
Vpi=Vpoi +Vpgi Adem s Soi =1Sw i Sgi=1Sw i
Tenemos:
Vpi=
Voi
Vgi
N Bt i G B gi
+
=
+
1Swi 1Sw i 1Sw i 1Swi
Vpi=
N Bt i +G Bg i
1Sw i
Vp=
N Bt i+ G Bgi
1e
1Sw i
Cf dP
Pi
]]
70
Cf dP
N Bt i+G Bgi
Vp=VpiVp=
1e
1Sw i
Pi
]]
N Bt i +G Bg i
1e
1Sw i
Cf dP
Pi
1.- Expansin
2.- Entrada de Agua
3.- Variacin del Volumen Poroso
4.- Drenaje del Yacimiento
Pi
(Cf ( PPi ))
1+ Cf ( PPi )
e x =1+ x +
x2 x3
+ +
2! 3 !
N ( Bt Bt i ) +G ( BgBgi ) +W ( Bt w Bt wi ) +We+
N Bt i +G Bg i
1e
1Sw i
Cf dP
Pi
71
1 dVw 1 dBw
=
( No hay cambio de fase )
Vw dP Bw dP
Entonces:
Bw
dBw
dBw
=Cw dP
= Cw dP
Bw
Bw Bw
P
Bw
Bw
ln
= Cw dP=Cw dP
Bw i
Bw
P
Bw=Bwi e
( Cw ( Pi P ))
Cw es muy pequeo
Bw Bw i [ 1+Cw ( PiP ) ]
Bw Bw i Bw i Cw ( PiP )
72
W ( Bt w Bt wi ) W ( BwBw i ) W Bw i Cw ( Pi P )
Notando que:
W Bw i =Vpi Swi
Adems, considerando que Sw es la misma en la zona de gas y de aceite:
W =Vpi
W=
Swi
Bw i
N Bt i+ G Bgi Sw i
1Sw i Bw i
W ( Bt w Bt wi ) W Bw i Cw ( PiP )
N Bt i +G Bgi Sw i
N Bt i +G B
Bw i Cw ( PiP ) N ( Bt Bt i ) +G ( BgBgi ) +
1Sw i Bw i
1Swi
Agrupando:
N ( Bt Bt i ) +G ( BgBgi ) +
m=
N Bt i +G Bgi
( Swi Cw +Cf ) P+We=Np Bo+ [ GpNp Rs ] Bg+ Bw
1Swi
Bt
Vol. Gas Libre Inicial GBg i
=
G=m N i
Vol . Aceite Inicial
NBt i
Bg i
Rp=
Gp
Gp=Rp Np
Np
73
N ( Bt Bt i ) +m N
Bt i
+
Bg i
N Bt i + m N
Bt i
Bgi
Bgi
1Sw i
Bt i
Swi Cw +Cf
( BgBgi ) + ( 1+ m ) N Bt i
P+ We=Np ( Bt + ( RpRsi ) Bg ) +WpBw
Bg i
1Swi
N ( Bt Bt i ) +m N Bt i
Sw i Cw +Cf
BgBgi
+ ( 1+m ) N Bt i
P+We+W I Bw I +GI BgI =Np ( Bt + ( RpRsi ) Bg ) +
Bgi
1Sw i
Sw i Cw +Cf
m NBti (
( 1+m ) N Bt i
P
)
BgBgi
+G
Bg
I
I
1Sw i
N ( Bt Bt i )
WeWpBw+W I Bw I
Bgi
+
+
+
=1
Fo
Fo
Fo
Fo
74
N ( Bt Bt i )
= Oil Drive Index o ndice de Empuje de Aceite
Fo
(ODI), porcentaje de produccin debido a la expansin del aceite y gas en
solucin.
m NBti
( bgbgi ) +GiBgi
Bgi
=
Fo
( 1+m ) N Bt i
Sw i Cw+Cf
P
1Swi
Fo
del
aceite
drenado
por
la
Fo= Np ( Bt + ( RpRsi ) Bg )
Eo=( Bt Bt i )
Swi Cw+Cf
(PiP) Incremento en volumen en la zona congnita de
1Swi
Se reduce a:
75
NEo+
m NBti
Eg+ ( 1+m ) N Bt i EGW +We+W I B WI +G I Bg I =Fo+WpB
Bgi
Bt i
Eg+ ( 1+m ) Bt i EGW
Bgi
Yo=FoWe+ Bw W I BWI G I B gI
Xo=Eo+ m
Bt i
Eg+ ( 1+ m) Bt i E GW
Bgi
Eo
Fo
m NBti
Eg+G I Bg I
Bgi
GDI =
Fo
WGI =
CDI=
( 1+ m ) N Bt i EGw
Fo
Bt i
Eg+ ( 1+m ) Bt i EGW
Bgi
]
76
FoYo=WeWpBw+W I B WI + GI Bg I
Fo=( FoYo ) + N Xo
Fo FoYo
=
+ N Recta con pendiente unitaria.
Xo
Xo
Yo = Yo + We ; Yo We= Yo
Yo= NXo
Por lo tanto:
Yo We = NXo
Yo = We + NXo
Dividiendo entre Xo:
Y o We
= + N Ecuacin de una recta con ordenada al origen
Xo Xo
We= k f (p,t)
Dnde:
77
Yo= N
NBti
Eg)=NEo+GEg
( Eo+ mBgi
Yo
EG
=G
+ N Recta con pendiente G y ordenada al origen N.
Eo
Eo
Yacimiento de Gas
G=
43560 Vb ( 1Sw )
Bg
Ejemplo:
El campo Bell, es un campo que tiene un rea de 1500 acres y el espesor de 40 ft.
Cul es el volumen total en acre-pie?
1500
ac.
40 ft
79
A= 1500 ac.
h= 40 ft
Vb= 1500 x 40= 60 000 ac-ft
= 22%
Swi = 23%
Bgi @ 3250 psia= 0.00533
G=
Pirmide:
Vb=
h
( A + A n+1 + A n A n+1 )( acft )
3 n
An+1
An
80
Trapezoide:
h
Vb= ( A n + An +1 )
2
Para una serie de trapezoides:
h
Vb= ( A 0 +2 A 1+ 2 A2 ++ 2 A n1 + A n ) +t avg A n
2
t avg = Espesor Promedio de las Isopacas; ft
Cuando la relacin de reas de dos isopacas sucesivas es menor a 0.5, se
aplica la frmula de pirmide.
Ejemplo:
Calcula el volumen neto de un yacimiento ideal del siguiente mapa de isopacas.
81
1 pg= 1000 ft
1 pg2= 22.96 acres
rea
Productiv map
a
a
A0
5''
A1
456''
A2
4.1''
A3
3.57''
A4
2.92''
A5
2.02''
A6
0''
1000 ( 1000 )
43560
rea del
Relaci
rea
Planmetro
n de
Acres
pg2
reas
19.64
450
16.34
375
0.83
13.19
303
0.8
10.05
231
0.76
6.69
154
0.67
3.22
74
0.48
0
0
0
Interval
oH
Ec. Que
Aplica
Vb acft
5
5
5
5
5
4
Trapezoide
Trapezoide
Trapezoide
Trapezoide
Pirmide
Pirmide
2063
1695
1335
963
558
99
6713
82
5
Vb= ( 450+375 ) =2063 acft
2
5
Vb= ( 154+74+ 154 74 ) =558 acft
3
4
Vb= ( 74 )=99 acft
3
Pi
0
P i Ai
Presin Promedio Areal=
Ai
0
Pi A i hi
Presin Promedio Volumen=
A i hi
0
Ejemplo:
Calculo de la Presin Promedio del Yacimiento.
Pozo
No.
1
2
Presin
Psia
2750
2680
Drenaje rea
acres
160
125
PxA
440,000
335,000
Est. Sd.
Thick
20
25
PxAxh
Axh
8,800,000
8,375,000
3200
3125
83
3
4
2840
2700
10970
190
539,600
26
14,029,600
145
620
391,500
1,706,100
31
12,136,500
43,341,100
4940
4495
15760
10,970
=2743 Psia
4
1,706,100
=2752 Psia
620
43,341,100
=2750 Psia
15,760
rea
A
B
C
Acres
25.5
15.1
30.5
Presin (Psi)
2750
2750
2850
h (ft)
25
15
25
Axh
637.5
226.5
1262.5
PxAxh
1,753,125
622,875
3,598,125
84
30.2
D
Presin Media=
2850
15
453
1,291,050
2579.5
7,265,175
7,265,175
=2816.5 Psia
2579.5
G=
Ga=
RecobroUnitario=43560 ( 1Sw )
Factor de Recobro=
100 ( GGa )
=
G
1
1
Bgi Bg a
100
1
1
Bgi Bg a
1
Bg i
=100 1
Bgi
Bga
85
Ejemplo:
Calcular el Factor de Recuperacin del Campo Bell cuando la presin baja a 500
psi.
Bgi = 0.00533 ft3/SCF
Bga = 0.03623 ft3/SCF
Factor de Recuperacin=100 1
0.00533
= 85
0.03623
86
Agua Congnita=43569 Sw i ft
Factor de Recuperacin
1Swi Sgr
SCF
=43560
acft
Bgi
Bga
100
Factor de Recobro=
1
scf
Bg i
1Sw i Sgr
Bgi
Bg a
1Sw i
Bgi
Ejemplo:
1. El Campo Bell produce bajo un acufero productivo y se estabiliza a
1500 Psia. El gas de saturacin residual es 24% y el factor de volumen
a 1500 lb es de 0.01122 ft/SCF. Cul ser entonces la recuperacin
unitaria en ft3-ac-ft y luego el factor de recobro en %?
Factor de Recuperacin=43560 ( 0.22 )
10.23
0.24
M SCF
=1180
0.00533 0.01122
acft
87
100
Factor de Recobro=
10.23
0.24
0.00533 0.01122
=85
10.23
0.00533
Fac . Rec .
SCF
=43560 ( 1Swi Sg r )
acft
Fac . Rec .
SCF
SCF
=43560 ( 0.22 ) ( 10.230.24 ) =43560 ( 0.22 ) ( 0.53 )=5079
acft
acft
Fac . Rec . =
Fac . Rec . =
100 ( 1SwiSgr )
( 1Swi )
100 ( 10.230.24 )
=68.83 69
( 10.23 )
2. Calcule la reserva inicial de gas que tiene un rea de 160 acres del
Campo Bell por depresionamiento volumtrico, bajo empuje de agua
parcial y empuje de agua total.
= 0.22
Swi = 23%
Sgr = 34%
Bgi =0.00533 ft3/SCF @ Pi =3250 Psia
88
Volumen Poroso
Vp=43560 A h=43560 ( 160 )( 40 )( 0.22 )
6
Vp=61.33248 10 ft
P= 500 Psia
G 3=
89
F . R .=85
N ( Bt Bt i ) +mN
Bt i
Swi C w +Cf
BgBgi ) + ( 1+m ) N Bt i
P+We=Np [ Bt + ( RpRso i ) Bg ] +WpBw
(
Bgi
1Sw i
YACIMIENTO DE GAS
Np Rp=Gp
m=
G Bgi
N mBt i=G Bgi
N Bt i
N ( Bt Bt i ) +
G Bgi
Swi Cw +Cf
BgBgi ) + [ N Bt i+ G Bgi ]
P+ We=NpBt + ( NpRpNp Rso i ) Bg+WpBw
(
Bgi
1Swi
90
Swi Cw +Cf
P+ We=NpBt ( GpNp Rso i ) Bg+WpBw
1Sw i
N y Np= 0
G ( BgBgi ) +G Bgi
Considerando que
Sw i C w +Cf
P+We=Gp Bg+ Bw 2.10
1Swi
Cw 0Cf 0
Bg=
P @C . S ZT
T @ C .S P
Bg i=
P@ C . S Z i T
T @ C .S Pi
P @C . S ZT P@ C . S Z i T
P
ZT
=Gp @ C . S
T @ C .S P
T @ C .S Pi
T @C . S P
) (
( PZ )G ( P )=Gp ( ZP )
Zi
i
91
Z
Z
Z
Gp
=G i
P
P
Pi
( ) ( ) ( )
Z
Z
( GGp )=G i
P
Pi
( )
( GGp )=
G Zi P
Pi Z
P Pi (GGp )
=
Z
G Zi
P
P Pi
=
Gp+ i
Z G Zi
Zi
y=mx +b
m=
b=
Pi
G Zi
Pi
Zi
Si P/Z 0
0=
Pi
P
Gp+ i
G Zi
Zi
Pi
P Gp
Gp= i
=1
G Zi
Zi G
92
Gp=G
93
G Bgi
V i=G Bgi
Bg
1 =Gp Bg+ Bw We
Bg i
Bg=
P @C . S Z f T
T@C . S Pf
Bg i=
P@ C . S Z i T
T@C . S ZiT
P@ C .S Z f T
T @C . S Pf
P
Z T
V i=
1 = @ C . S f Gp+BwWpWe
P@ C. S Z i T
T @ C .S Pf
T @C . S Z i T
V i=
V i=
P@ C. S Z f T T @ C . S Pi
P@ C . S Z f T
1 =
Gp+ BwWpWe
T @ C .S P+ P@ C . S Z i T
T @ C .S Pf
Z f Pi
P@ C . S Z f T Gp
1 =
Pf Zi
T @ C. S Pf
P@ C . S Gp Pf
Z f Pi
=
Vi
1
T @ C .S
Zf T
Pf Z i
P@ C . S Gp Pf V i Z f Pi Pf V i
=
T @ C .S
Z f T Pf Pi Z f T
P@ C . S Gp Pi V i Pf V i
=
T @ C .S
ZiT
Zf T
94
Presin de abandono
Pt =f (Q , , , H 2 O )
Ejemplo:
Calcule el volumen original de gas y la reserva inicial del yacimiento de gas con
datos de presin-produccin, se trata de un yacimiento volumtrico.
Pi = 3250 Psia
Ty= 213 F= 673R
T@ C.S= 60F= 520 R
Gp= 1.00 x 109 SCF
Pf= 2864 Psia
Z@ 3250 Psia= 0.91
Z@ 2864 Psia= 0.888
Z@ 500 Psia= 0.951
P@ C.S = 15.025 Psia
P@ C . S Gp Pi V i Pf V i
=
T @ C .S
ZiT
Zf T
3250V i
2864 V i
15.025 ( 1 109 )
=
520
0.91 ( 673 ) 0.888 ( 673 )
28894230.77=5.30672 V i4.792331V i
95
28894230.77=0.51143979323V i
V i=
28894230.77
0.51143979323
V i=56.1687 MM ft 3
G=
G=
Vb ( 1Sw i )
Bgi
De la ecuacin:
G ( BgBgi ) +We=Gp Bg+ Bw
We=Gp BgG ( BgBgi ) +Bw
We=935.4 106 ( 0.0057 ) 10181 106 ( 0.00570.005262 ) +1.03 ( 15,200 ) 5.615
pies3
bls
We=960 305 ft 3
Sw f =
Sw f =
Sw f =
( Vb w Swi ) +We Bw
Vbw
13.04 10 ( 0.172 )
Sw f =0.64 o 64
97
Sg r=1Sw f =10.64=36
% Mol
Pc(Psia)
Tc
0.75
0.2
0.05
673.1
708.3
440.1
343.2
549.9
914.2
yPc
yTc
504.8
254.4
141.7
110
22
45.7
Ppc=668. Tpc=413.
5
1
Pi =4200 Psia
T= 180F= 40R
Datos de Produccin:
P/Z (Psia) Gp (MMM SCF)
4600
0
3700
1
2800
2
98
Gp @ 2000 Psia=?
Ppr=
2000
=2.94
668.5
Tpr=
640
=1.55
413.1
z= 0.8
P 2000
=
=2500 Psia
Z
0.8
99
y y 1=
x1=0
y1= 4600
x2= 2
y2= 2800
x 2x 1
( xx 1)
y 2 y 1
y4600=
28004600
( x0 )
20
y4600=
1800
x
2
y4600=900 x
P
=900 Gp+4600
Z
2500=900 Gp+4600
Gp=2.33 MMM SCF
P/Z= 0
Gp= G
0=900 G+ 4600
G=5.11 MM SCF
Recordando V i=G Bgi
Bg i=0.02829
Zi T
Pi
100
V i=G 0.02829
1
T
Pi
Zi
1
( 4600
)( 640 R )( 0.02829 )
V i=20.11 10 6 ft 3
V i= A h ( 1Sw i)
A=
Vi
20.11 106
=
=15.962 106 ft 3
h ( 1Sw i ) ( 0.15 )( 0.12 ) ( 10.3 )
A=
15.962 106 ft 2
=366.4 ac
43560
101
Proceso de Separacin
( Gp )surf
( Gp)Ps
( Gp )SS
( Gp)ST
Pa
Np
Ps1
Ps2
Lquid
o
102
( R )Ps=
R=
( Gp ) Ps
Np
( R )ss=
( Gp )ss
Np
( R )ST =
( Gp )ST
Np
( Gp )surf
=( R )Ps +(R) ss+( R)ST
Np
( g)SS +( R) ST ( g )ST
( g )Ps +( R)SS
R
( R)Ps
( g)surf =
n=
m
M
n(lbmol)
M (lb)
lb
( lbmol
)
R( SCF)
( 28.97 ) ( g ) surf (lbs)
SCF
379.4
lbmol
lbm
ft 3
5.615
3
bl
ft
103
141.5
( n g )= 379.4
28.97
R ( g )surf +350.2 ( o ) STo
379.4
( no )ST =
5954
lb
( o ) ST
R
+350
379.4
( M o ) ST
28.97
R
+350.2 ( o )ST
mwf 379.4 ( g ) surf
M wf =
=
n wf
R
+350 o
379.4
M o ST
( )
M wf
( g )wf = 28.97
( g )wf =
R ( g ) surf + 4586 ( o ) ST
R+132869
o
Mo
( )
ST
104
V=
V = ( ) o=
( )
mo
Mo
( )
=132869
T sc 350.2 ( o )ST
60+ 460
=
( 10.73 )
Psc
14.7
( M o ) ST
ST
o
Mo
nRT
P
( )
ST
SCF
stb
mw T sc 350.2
60+ 460
R
=
( 10.73 )
M w P sc
18
14.7
( )w =7382
scf
stb w
( n g )P
( n g ) P+ ( no ) P
( n g )P = 379.4
mo
( )
( no )P = M
o
ST
=350.2
o
Mo
( )
( ) o
379.4
=
379.4
ST 379.4
105
( )o
( no )P = 379.4
f g=
R
379.4
( )o
R
+
379.4 379.4
R
=
R+ ( )o
R
R+132,869
IGIP=f g G
o
Mo
( )
IOIP=
ST
fgG
R
Ejemplo:
Calcular el IGIP y el IOIP para un yacimiento de gas y condensado y encontrarlo
de forma genrica (acre-pie).
Pi = 2500 Psia
Ty = 215F= 675R
=25%
Swi = 30%
Qo=242 stb
o @ 60F= 48API
Qg sep= 3100 M SCFD
g sep= 0.65
Qg= 120 M SCFD
g ST= 1.20
5954
lb
106
R1=
R 2=
Gp 1 3 100 000
SCF
=
=12 809.91
Np
242
bl
120 000
SCF
=495.86
242
bl
( g ) surf =
( g )wf =
SCF
bl
R ( g ) surf + 4586 ( o ) ST
R+132869
o
Mo
( )
ST
R
13305.77
( n g )P = 379.4 = 379.4 =35.07 lbmol
mo
( )
( no )P = M
o
=350.2
ST
0.788
( 151.43
)=1.8 16 lbmol
f g=
( n g )P
35.07
=
=0.951
35.07+1.816
( n g ) P + ( no ) P
G=
Bg i=0.02829
ZT
P
2
107
P pr=
2740
=4.3
636.6
T pr =
675
=1.596
422.7
De grficas z= 0.825
Bg i=0.02829
G=
3
0.825 ( 675 )
ft @c . y
=5.7496 103
2740
SCF
5.7496 10
=1.326 MM
SCF
acft
IOIP=
SCF
acft
f g G 0.951(1.326 106)
bl
=
=94.77
R
13305.77
acft
Gp r + Gp2 3100+120
=
=3386 M SCFD
fg
0.951
PV =nRT
PV
ZRT
P y V y P sc Z y R T y
=
Zy RT y
R T sc P y
V y = Gp
P sc Z y T y
(14.7)(0.825)(675)
=3386000
=19454 ft 3
T sc P y
(520)(2740)
Ejemplo:
108
Tsup= 70 F
Pws=Pwh e
T avg=
0.01875 g D
Z avg T avg
70+175
=123 F
2
Tpc=377.59=380
Tpr=
583
=1.53
380
Ppc=663.33=670
Ppr=
P avg
670
Pws=Pwh+ 0.25
Pws=3400 e
Pwh D
3400 7500
=3400+0.25
100 100
100 100
0.01875(0.7 )(7500)
0.85 (123+ 460)
Pws
(Psia)
4038
=4147
Pavg
(Psia)
3719
Ppr
Zavg
Pws
5.55
0.85
4147
109
4147
4142
3774
3771
5.63
5.62
0.855
0.855
4142
4142
Problemas:
3.3 La arena M es un pequeo yacimiento con una presin de 3200 Psia y una
temperatura de 220 F.
P (Psia)
3200
2925
2525
2125
a)
Gp (MM
scf)
0
79
221
452
G ( BgBgi ) +
Bg (rcf/scf)
5.2622 x 10-3
5.7004 x 10-3
6.5311 x 10-3
7.736 x 10-3
G (MM
SCF)
1028
1137
1413
P/Z
We (bls)
756
27000
174237
G Bgi
( Cw Sw+Cf ) P+We=Gp Bg+ Bw
1Swi
G ( BgBgi ) =Gp Bg
110
G=
Gp Bg
BgBgi
G=
G=
G=
c)
G= 1018 MM SCF
P/Z
Gp
G ( BgBgi ) +We=Gp Bg
We=Gp BgG ( BgBgi )
We=79 106 ( 5.7004 103 ) 1018 106 ( 5.7004 1035.2622 103 ) =4244
We=22110 6 ( 6.5311 103 )1018 10 6 ( 6.5311 1035.2622 103 ) =151632
We=452 106 ( 7.7360 103 )1018 106 ( 7.7360 1035.2622 103 )=978343.6
111
Bg=0.02829
ZT
P
P 0.02829 T
=
Z
Bg
Z=
Bg P
0.02829 T
3.4 El Campo de gas Sabine tiene un yacimiento que tiene las siguientes
caractersticas:
Pi = 1700 Psia
T= 160 F
Gp= 5.00 BCF
Vgi=( Vp ) HC =G Bgi
Bg=0.02829
G=
ZT
P
Gp Bg
BgBgi
P (Psia)
Bg (rcf/scf)
1700
1550
1400
1100
500
200
0.85
0.86
0.868
0.89
0.944
0.973
8.7635 x 10-3
9.7246 x 10-3
10.8667 x 10-3
14.1809 x 10-3
33.0908 x 10-3
85.2685 x 10-3
Gp (MMM
SCF)
5
9.792
19.327
37.193
45.391
P/Z
(Psia)
2000
1802
1612.9
1236
529.6
205.5
112
G=
5 10 9 ( 9.7246 103 )
Gp Bg
=
=50.59 10 9
BgBgi 9.7246 1038.7635 103
Gp=
G ( BgBgi )
Bg
Z @ 750 =0.92
P 750
=
=815.2
Z 0.92
Cada de presin aproximada por cada mil millones:
36.6 Psia
Valor mnimo de la reserva inicial:
5
+
Pi=1700 Psia
12 Psia
+
P=1550 Psia
Z i=8.949
Z =0.86
Bg i=8.6918 10
Bg=9.8119 103
Gp ( Bg )
4.75 109 ( 9. 8119 103 )
G m n=
=
=41.609 109 SCf
3
3
BgBg i ( 9.8119 10 8.6918 10 )
Gm n=41.6 BCF
3.5 Si de cualquier manera la produccin de los 5000 BCF del ejemplo anterior, 4
bls de agua entraron al yacimiento, la presin cay a 1150 Psia.
G ( BgBgi ) +We=Gp Bg
G=
9
3
6
Gp BgWe 5 10 ( 9.7246 10 )( 4 10 5.615 )
9
=
=27.22 10 =27.22 BCF
3
3
( BgBgi )
(9.7246 10 8.7635 10 )
Error=
50.5927.22
=0.85=85
27.22
114
Gp (BCF)
200
420
Ppr (Psia)
7.53
6.03
3.77
Z
1.02
0.96
0.91
P/Z
4901.96
4166.66
2747.25
T= 237F
g=0.7
Ppc= 663.33 Psia
Tpc= 377.59 R
P
=mGp+ b
Z
4901.96 = m 200 + b
4166.66 = m 420 + b
420 m + b = 4166.66 (-200)
200 m + b = 4901.96 (420)
1225491.2
=5570.41
220
115
1403.75
=3.342
420
Gp=
2747.255570.41
=844.75 BCF
3.342
P
=mGp+ b
Z
P
=m ( 0 ) +5570.41=5570.41
Z
Calcular el IGIP y el IOIP para un yacimiento que tiene gas y lquido, tenemos los
siguientes datos de produccin:
Py= 4350 Psia
Ty= 217F= 677R
HC= (1-Swi)= 17.4%
Cond. Std.= 15.024 Psia y 60F
Qg sep= 842 600 SCFD
Qo ST= 31.1 STBD
Mo cf+= 185 lb/lb-mol
g cf+= 0.8343
116
[email protected]= 0.7675
Bo@ sup= 1.235 bl/ STB
Vol. Molar@ 15.025 y 60F= 371.2 SCF/lb-mol
(1)
Comp
(2)
Gas
(3)
Lquido
CO2
0.012
0.94
0.2024
C2
0.031
0.0484
C3
0.01
0.0312
i-C4
0.002
0.0113
n-C4
0.002
0.0196
n-C5
6E04
3E04
(5)Mol
sep (3)(4)
(7) (3)
(6)
(8) (2)
(59.11)
(9) (3)
(2.107)
(10)
(8)+(9)
(11)
(10)//
(10)
(12)
Pc
(Psia
)
(13)
0.709
0.709
0.0116
0.0266
55.587
0.4265
56.014
0.0086
1.803
0.102
C1
i-C5
(4)
Mo
0.0159
0.017
C6
0.001
0.0384
C7 +
7E04
0.6158
16.0
4
30.0
7
44.0
9
58.1
2
58.1
2
72.1
5
72.1
5
86.1
7
185
3.247
1.455
0.131
7
0.177
1
(15)
(11)
(12)
(14)
Tc
(R)
1070
12.41
548
6.356
0.915
673
615.8
343
313.85
1.905
0.0311
708
22.02
550
17.11
1.376
0.248
0.0077
0.562
0.0657
0.6277
0.0102
617
6.293
666
6.79
0.657
0.294
8
0.0033
0.142
0.0238
0.1658
0.0027
529
1.428
735
1.98
1.139
0.284
0.0056
0.136
0.0413
0.1773
0.0029
550
1.595
766
2.22
0.0052
0.035
0.0335
0.0685
0.0011
484
0.532
830
0.91
0.0056
0.018
0.0358
0.0538
0.0009
490
0.441
846
0.76
0.0142
0.077
0.0809
0.1579
0.0026
440
1.144
914
2.38
0.3902
0.041
1.2975
1.3385
0.0219
300
6.57
1227
26.87
0.4671
59.11
2.107
61.217
1.147
1.227
3.309
113.923
127.478
0.324
8
0.326
4
0.370
6
0.633
6
668.2
379.23
(n0 )sep=
(11)
(14)
117
(n0 )sep=
1
1
lbmol
=
=2.14
bl
(7) 0.46706
(Qg )sep
SCF @ c . s
842600
=
=21 938
QO BO 31.1(1.235)
bl @ c .s
(n g) sep=
R sep 21 938
lbmol
=
=59.11
371.2 371.2
bl
+
7
(350.2)
C1 =
185
bl
=0.6336
lbmol
0.8343(350.2)
C7+
M0
+=
lbmol
C7
bl
16.04
=0.1317
0.348(350.2)
Tpr=
667
=1.785
379.23
Ppr=
4350
=6.509
668.2
Z=0.963
ft 3
acft
118
PV =ZnRT
n=
PV 4350 ( 7579.44 )
lbmol
=
=4713.15
ZRT ( 0.963 ) ( 10.73 )
acft
fg=
( ng ) sep
59.11
=
=0.966
( ng ) sep+ ( no ) sep 59.11+ 2.107
IGIP=fgG
V i=G Bgi
G=
Vi
Bg i
Bg i=
ft
acft
6 SCF
=1.75 10
3
acft
4.33 10
7579.44
G=
SCF
SCF
=4713 371.2=1.75 10 6
lbmol
acft
119
SCF
fg G
acft
STB
IOIP=
=
=62.37
R
bl
acft
SCF
@ c .sep
21938
1.235
Bl @ c .sep
STB
1.69 10
IOIP=
Produccin de Gas
Gp=
Qg sep 842600
=
=872.256 SCFD
fg
0.966
n=
P
Gp
P V
PV
= @ c. s
= y @c . y
ZRT ( 1 ) R T @ c . s Z @ c. y R T y
V y=
P@ c. y GpZ T y
T @ c. y P y
V y=
15.025(872256)(0.963)(677)
4350(520)
V y =3777
ft 3
day
Ejercicios:
120
3.7
y= mx + b
P
=
z
m Gp + b
M= -3.342
b= 5570.36
P
=
z
-3.342 Gp + 5570.36
Gp @2500 Psia =
R@ 2500 =
844.75 BCF
P/Z 0
Gp G
0= -3.342G + 5570.36
G=
5570.36
3.342 =1667 BCF
P
Z
Gp 0
Pi
Zi
Pi
Z i =-3.342(0)+5570.36
Pi
Z i =5570.36
Pi=Zi
5570.36
121
P
600
0
500
0
400
0
250
0
Gp
0
20
0
42
0
84
4
Pi
600
0
606
1
608
8
610
0
610
5
Pr
9.03
Z
1.088
Pi*
6060.55
9.13
1.093
6088
9.17
1.095
6100
9.18
1.096
6105
9.19
1.096
6105
Pi= 6105
Gp (MMM
SCF)
x
0.4
1.2
3.2
4.3
5.5
7.5
8.7
y= mx + b
122
y y 1=
y6553=
( y 2 y 1 )
(xx 1 )
( x 2x 1 )
( 608065538 )
( x0.393)
( 8.7490.393 )
y-6553=-56.6060(x-0.393)
y-6553=-56.6060x+22.2461
y=-56.6060x+22.2461+6553
P
Z =-56.8893 Gp+6569.474
P
Z 0
G=
Gp G
6569.474
56.8893 =115.44 BCF
P
Z =-56.8893 Gp+6569.474
1000=-56.8893 Gp+6569.474
Gp=
6569.4741000
=97.90 BCF
56.8893
123
R=
Gp @ 1000 97.90
=
=0.848=84.8
G
115.44
Gp=97.950=47.9
P
=56.8893 ( 47.9 ) +6569.474
Z
P
=3844.47
Z
3.15
qg sep = 6 MM scf
qo ST = 100 STB
qg ST = 21 M scf
qw ST = 10 bls
P1 = 6000 Psia
Py = 2000 Psia
Ty = 225F
Cantidad de vapor @ 6000 Psia y 225F = 0.86 bls/MM scf
API = 50
124
( o ) =132869
o=
scf
=714
( Mo ) =132869 ( 0.7796
)
145
stb
o
ST
141.5
=0.7796
131.5+ 50
( Mo ) ST =
5954
=145
API 8.811
Np
N
Agua Intersticial
Aceite
7758 Sw
bls
acft
7758 (1Sw)
125
Vol. Aceite@ST
7758 (1Sw)
Boi
Gas
7758 Sg
Ac. Yacimiento
7758 (1SwSg)
Vol. Aceite@ST
7758 (1SwSg)
Bo
Recobro@ST
Ree=1
7758
(1Sw) (1SwSg)
Boi
Bo
1SwSg Boi
(1Sw)
Bo
0.679
1.283
126
300 Psia
73F
69551 scf/stb
( g )wf =
76F
o = 55.9 API
366 scf/stb
R ( g )surf +4586 ( o ) ST
R+132 869
( Mo )
o
ST
( g ) surf =
scf
stb
141.5
141.5
127
5954
5954
lb
Rec=1
(1Sw) (1SwSg)
Boi
Bo
( 1SwSg ) Bo i
( 1Sw ) Bo
Recobro=
7758 Sor
7758 Sor
Boi
Fac . Recobro=
(1SwSor)
(1Sw)
128
Ejemplo:
K= 100 md
Sw= 25%
o=2.0 cp
= 0.20
h= 10 ft
RF=0.114+ 0.272 log ( 1000 ) +0.256 ( 0.25 )0.136 log ( 2.0 )1.538 ( 0.20 )0.00035 ( 10 )=0.642=64.2
Cw Sw i Cf
P+W e=Np [ Bt +( RpRsoi) ] + Bw
1Sw i
Yacimiento Volumtrico
El volumen de aceite permanece constante
N ( Bt Bt i )=Np [ Bt ( RpRsoi ) Bg ] Ec 5.6
Py > Pb
Rp=Rso=Rsoi
N ( Bt Bt i )=Np Bt
129
Factor de Recobro=
Np ( BtBt i )
=
N
Bt
para P > Pb
Ejemplo:
Se desea calcular el volumen original de aceite.
Pb= 3550 Psia
Pi= 4400 Psia
Bti= Boi @ 4400= 1.572 bls/STB
Bt@ 3550= 1.6 bls/STB
FR=
( Bt Bt i) 1.61.572
=
=0.0175=1.75
Bt
1.6
Para
N=
Np Bt
=38.8 MM STB
BtBt i
P y pb
Np [ Bt ( RpRso i ) Bg ]
BtBt i
130
RF=
RP=
Bt Bt i
Np
=
N [ Bt ( RpRsoi ) Bg ]
Gp
N
Gp 0
Rp 0
SCF
STB
B oi =1.572
SCF
STB
bl
STB
Rso@ 2800=900
SCF
STB
Bo @2800 =1.52
SCF
STB
Z @ 2800 =0.87
Bg=0.02829
ZT
1
P 5.615
)
131
Bg=0.02829
Bg= 0.00102
Bt=Bo+Bg (
(0.87)(650) 1
2800
5.615
bl
STB
RsoiRso
RF=
bls
STB
1.7241.572
=0.0383=3.83
1.724+(33001100)(0.00102)
Si reinyecto
2
3
del gas
1
3
SCF
STB
1.7241.572
=0.0881=8.81
1.724+(11001100)(0.00102)
RF=4.98
132
133
Yacimiento
Tubera vertical
q 0=
q 0=
Ko h( PePwf )
141.2 oo ln + s
rw
q0
PwsPwf
Ko h( PePwf )
141.2 oo ln + s
rw
q0
=
( PePwf )
J=
Produccin y productividad
q 0=
Ko h
141.2 0 0 ln
+s
rw
Productividad
q0
PwsPwf
135
Ecuacin de Vogel
q0
q0 max
=10.2
Pwf
Pwf
0.8
Pws
Pwf
( )
GAS
Q0
2
2
= ( Pws Pwf )
Rgimen tranciente: No est en estado estable, hay cambios hasta que los
efectos del fenmeno se manifiestan en la frontera, por ejemplo tirar una
piedra en un lago.
Estado estacionario: Se dice que un sistema fsico est en rgimen
estacionario cuando las caractersticas del mismo no varan con el tiempo,
las propiedades se mantienen constante, es decir se encuentra en equilibrio
termodinmico.
Rgimen tranciente tardo: Existe una variacin con respecto al tiempo,
cuando llega a la frontera.
Pseudo estacionario: Se encuentra en un equilibrio dinmico, la variacin
con respecto no es cero pero es constante.
136
LEY DE DARCY
Es cuando se mueve de un fluido con una gravedad de 1cp a una velocidad de
1
atm
cm
=0.001127
k dp
0.433 cos
ds
= velocidad aparente
atm
.
cm
bbls
2
dayft
k= permeabilidad (md)
= viscosidad del fluido (cp)
P= presin psia
S= distancia a lo largo del flujo ft
= gravedad especifica ( resp. Agua)
= angulo medido en contra de las manecillas del reloj de la vertical a la
direccin positiva de s
137
FLUIDOS INCOMPRESIBLES
C=
1 dv
=0
V dp
C=
1 dv
V dp
C( PR-P)= Ln
e C(PR P) =
PR
vr
Cdp dvv
( PRv )
V
VR
e
(
C
(PRP))
V= VR
X X
X
x
PERO e =1+ X+ 2 ! + 3! + n !
138
Donde la
C ( PRP ) =X
ZnRT
P
1 1 dz
Cg=
P Z dp
Tipos de flujos:
139
DIFUSIVIDAD DE LA FORMACION
n=
k
ct
Ct =C g C g +C o S O +C w C w +C f
140
T pss
horas
2
1200 1200
T pss =
n
k
= 1.5 cp
ct=
15 10 psi
re= 1000 ft
k= 100 md
= 20%
141
GAS
g= 0.015 cp
ct=
400 10 6 psi1
re= 1000 ft
k= 100 md
= 20%
T pss =
142
V=
Q qB
k dp
= =0.001127
A Ac
dx
Q= qB
componente gravitacional=0
qB
k dp
dx= 0.001127
Ac
dx
P2
qB
k
dx=0.001127 dp
Ac 0
P1
qB
k
( L0 ) =0.001127 (P1 P2 )
Ac
143
qB
k
L=0.001127 ( P 1P 2)
Ac
q = 0.001127
kAc( P1 P 2)
BL
P= 100 Psia
K= 250
=
2.5 cp
bbl
STB
B= 1.127
L= 450 FT
A= 45 ft
q = 0.001127
kAc( P1 P 2)
BL
250 ( 45 ) (100)
stb
0.001127
=1
Q=
day
1.127 ( 2.5 ) (450)
144
Q=qb=qR(1+C (PrP) )
V=Vr(1+C(Pr-P))
L
P2
qR
k
dp
dx=0.001127
Ac 0
P 1 1+C ( PrP)
qR=
Calcular qR
P= 100 Psia
K= 250 md
=
2.5 cp
145
L= 450 FT
6
1
C= 65 10 psi
A= 45 ft
P2=P R
P1=P R
qR=
qR=
6
0.001127 ( 250 ) (45) 1+ 65 10 (100)
ln
1
( 2.5 ) ( 450 ) (65 106)
qR=1.123
Ley de Darcy
146
V=
qB
qB
k dp
=
=0.001127
Ac 2 rh
dr
r1
r2
con
P1
P2
qB
k
dr = 0.001127 dp
2 r h
r2
P2
qB
dr
k
=0.001127 dp
2 h r1 r
P1
r
qB
k
ln 2 =0.001127 (P 2P 1)
2h
r1
( )
1
P2P
0.00708 kh
q=
147
1
P2P
0.00708 kh
q=
w
Pe P
0.00708 kh
q=
V=
Vr
V=
(1+C(
Pr
-P))
r2
P2
qR
dr = (1+C (dpP P))
0.00708 k h r 1 r
P1
r
( 1+ C(Pr P1))
C
ln
qR
( ln r 2ln r 1 ) =
0.00708 k h
qR=
0.60708 kh 1+ C( Pr P2)
ln
1+ C( Pr P1)
r2
C ln
r1
( )
149
qBg=
q Pcs T 2
5.615 T cs P
V=
q Pcs T 2
k dp
=0.001127
5.615 T cs P A c
dx
L
P2
q Pcs T 2
dx= Pdp
0.00632811 kT cs A c 0
P1
q=
0.003164 kT cs A c ( P12P22)
Pcs T 2 L
150
BIBLIOGRAFIA
Imgenes sacadas del libro: Applied Petroleum Reservoir Engineering, segunda
edicin, B.C. Craft y M.F. Hawkins
151