Universidad Estatal Península de Santa Elena

Facultad de Ciencias de la Ingeniería

Ingeniería en Petróleo `

Asignatura
Flujo de Fluidos de un Medio Poroso

Profesor:
Ing. Fidel Chuchuca

Tema
Ejercicios Propuestos Capitulo 2

Nombre
Franco Gonzabay Kimberly Dayanna
Llerena Zambrano Ivanna
Mendoza Balderramo Juan Pablo
Tomala Ortega Lidia Verónica
Curso
4/1 Petróleo
 EJERCICIOS PROPUESTOS
Ejercicios 2.1
Verifique la constante 1.127 en la ecuación siguiente

k . A Δp
q=−1.127 , Ec.2.22
μ Δx

1. Considerando las siguientes unidades

d 1
q=bblres /d =
md 1000
cm
k= Darcy pie
=30.48

atm 1
A= pie2 =
Psi 14,7

dp/dx = psi/pie

u = cp q¿

cm2
−1 ( D ) ×30.48 Pie 2( )
Pie 2 1( Psi)
q= ×
( 1.84 βo )( 30.48 pie ) cp 14.7( Pie)

KA ∆ p
q=−1.127 .
u ∆x

2. Encuentre el valor de la nueva constante de ecuación 2.22 si,

A=m 2
dp /dx=atm /m
k . A Δp
q=−1.127
μ Δx
cm
1 ( D ) ×3048 m2 ( ¿ ¿ 2)
m2 atm
q=− ×1 ¿
( 1.84 βo ) ( 30.48 m ) cp m

KA dp
q=−54.347
u dx
Ejercicio 2.2
Un yacimiento en contacto con un acuífero activo produce con una
presión de fondo fluyente ( Pwf ) de 900 psi. Usar los siguientes datos del
yacimiento para calcular la tasa de flujo en STB por día. (STB/día)
Datos
 Porosidad. ϕ = 25%
 Factor Volumétrico de formación del petróleo, βo = 1.25rb/STB
 Espesor neto de formación, h = 30 pies
 Viscosidad del petróleo, μ = 19 cp
 Radio de pozo, r w . =0.49pies
 Radio externo, r e = 1000 pies
 Presión externa del yacimiento, Pe = 1110 psi
 Permeabilidad, k = 130 mD

Asumiendo que el yacimiento es abierto y el régimen de flujo es continuo,

la ecuación de Darcy en unidades de campo, empleando la presión externa
del yacimiento es

1.127 x 10−3 . 2 πkh ( Pe −Pwf )

q= Ec. 2.36
βo μ ( ln r e /r w )

1.127 x 10−3 . 2 π ( 130 mD ) ( 30 pies ) ( 1110 psi−900 psi )

q=
rb 1000 pies
(
1.25
STB ) (
( 19 cp ) ln
0.49 pies )

STB
Respuesta; q=32.04
dia
Ejercicio 2.3
Un yacimiento produce con una presión de fondo fluyente ( Pwf ) de 1088
psi. Usar los siguientes datos del yacimiento para calcular la tasa de
flujo en STB por día. (STB/día)
Datos
 Porosidad. ϕ = 28%
 Factor Volumétrico de formación del petroleo, βo = 1.4rb/STB
 Espesor neto de formación, h = 49pies
 Viscosidad del petroleo, μ = 21 cp
 Radio de pozo, r w . =0.49pies
 Radio externo, r e = 820pies
 Presión promedio del yacimiento, P = 1204 psi
 Permeabilidad, k = 125 mD

Asumiendo que el yacimiento es abierto y el régimen de flujo es continuo,

la ecuación de Darcy en unidades de campo, empleando la presión externa
del yacimiento es

1.127 x 10−3 . 2 πkh ( P−P wf )

q=
1 Ec. 4.97
(
βo μ ln r e /r w −
2 )
1.127 x 10−3 . 2 π ( 125 mD ) ( 49 pies )( 1204 psi−1088 psi )
q=
820 pies 1
(
( 1.4 rb /STB ) ( 21 cp ) ln
0.49 pies 2
−)

STB
Respuesta; q=23.06
dia
Ejercicio 2.4
Un yacimiento circular y pozo en el centro con las siguientes
características, produce con una presión de fondo fluyente ( Pwf ) = 1100
psi. Asumir que la presión del yacimiento se encuentra soportada por
un acuífero (Flujo continuo).
Datos
 Porosidad. ϕ = 25%
 Factor Volumétrico de formación del petroleo, βo = 1.28rb/STB
 Espesor neto de formación, h = 60 pies
 Viscosidad del petroleo, μ = 21 cp
 Radio de pozo, r w . =0.49 pies
 Radio externo, r e = 1000 pies
 Presión inicial del yacimiento, Pi = 1250 psi
 Permeabilidad, k = 125 mD

a) Calcular el caudal de flujo Q en STB/día.

b) Calcular y graficar la presión del yacimiento a diferente distancias r
con respecto al pozo ( P vs r ) , asumiendo que la presión de fondo
fluyente Pwf se mantiene constante. Emplear intervalos de 50 pies.
c) ¿En qué sección del yacimiento ocurre la mayor caída de presión?
Justifique el motivo de su respuesta.

Solución
a) Para un yacimiento abierto con empuje hidráulico, la ecuación de
flujo en unidades de campo es

1.127 x 10−3 . 2 πkh ( Pe −Pwf )

q= Ec. 2.36
βo μ ( ln r e /r w )

1.127 x 10−3 . 2 π ( 125 mD ) ( 60 pies ) ( 1250 psi−1100 psi )

q=
rb 1000 pies
(
1.28
STB ) (
( 21 cp ) ln
0.49 pies )
 STB
Respuesta; q=38.89
dia

b) En flujo continuo, las presiones se mantienen constantes a cualquier

tiempo ( dp /dt ) =0.Recordando que las soluciones de las ecuaciones de
flujo son para tasa de flujo al pozo constante y que en este caso es q
= 38.89 STB/día, entonces,

1.127 x 10−3 . 2 πkh ( Pe −P )

q=
βo μ ( ln r e /r )

qμβo P e−P
=
2 πkh ln ( r e /r )

qμβo r
P=P e−
2 πkh ( )
ln e
r

qμβo r
P=P e− −3
1.127 x 10 .2 πkh ( )
ln e
r

qμβo re
P=P e− −3
7.08 x 10 kh ( )
ln
r

141.2 x qμβo r
P=P e−
kh ( )
ln e
r

La presión P al radio de 50 pies es

( 38.89 STB/dia )( 21 cp )( 1.28 rb / STB ) 1000 pies

P=1250 psi−141.2 x
(125 mD ) ( 60 pies ) ( ln
50 pies )
P=1191.04 psi

La presión P al radio de 100 pies es

( 38.89 STB/dia )( 21 cp )( 1.28 rb / STB ) 1000 pies

P=1250 psi−141.2 x
(125 mD ) ( 60 pies ) ( ln
100 pies )
P=1204.68 psi
 La presión P al radio de 150 pies es
STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 pies
P=1250−141,2
( 125 mD ) ( 60 pies ) ( ln (
150 pies ))

P=1212,66 PSI

La presión P al radio de 200 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
200 ft ))

P=1218.35 PSI

La presión P al radio de 250 pies es

STB rb
( 38,89
dia )
( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
250 ft ))

P=1222.71 PSI

La presión P al radio de 300 pies es

STB rb
( 38,89
dia )
( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
300 ft ))

P=1226,30 PSI

La presión P al radio de 350 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb )
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln ( 1000 ft
350 ft ))

P=1229,33 PSI

La presión P al radio de 400 pies es

 STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( 400 ft ))
ln (
P=1231,96 PSI

La presión P al radio de 450 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
450 ft ))

P=1234,28 PSI

La presión P al radio de 500 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
500 ft ))

P=1236,35 PSI

La presión P al radio de 550 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
550 ft ))

P=1238,23 PSI

La presión P al radio de 600 pies es

STB rb
( 38,89
dia )
( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
600 ft ))

P=1239,94 PSI

La presión P al radio de 650 pies es

 STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( 650 ft ))
ln (
P=1241,52 PSI

La presión P al radio de 700 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
700 ft ))

P=1242,98 PSI

La presión P al radio de 750 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
750 ft ))

p=1244,33 PSI

La presión P al radio de 800 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
800 ft ))

P=1245,60 PSI

La presión P al radio de 850 pies es

STB rb
( 38,89
dia )
( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
850 ft ))

P=1246,80 PSI

La presión P al radio de 900 pies es

 STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( 900 ft ))
ln (
P=1247,92 PSI

La presión P al radio de 950 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
9500 ft ))

P=1248,99 PSI

La presión P al radio de 1000 pies es

STB rb
( 38,89
dia ) ( 21 cp ) (1,28
stb ) 1000 ft
P=1250−141,2
( 125 md ) ( 60 ft ) ( ln (
1000 ft ))

P=1250 PSI

Tabla 2.2 Valores de P vs. r, calculados

r( pies ) P( psi ) r ( pies ) P ( psi )

0.49 1100.01 550 1238.23
50 1191.04 600 1239.95
100 1204.68 650 1241.52
150 1212.66 700 1242.98
200 1218.33 750 1244.34
250 1222.72 800 1245.61
300 1226.30 850 1246.80
350 1229.33 900 1247.92
400 1231.97 950 1248.99
450 1234.28 1000 1250.00
500 1236.36
Gráfica de P vs R
c) ¿En qué sección del yacimiento ocurre la mayor caída de presión?
Justifique el motivo de su respuesta.
Al momento de realizar el cálculo de variación de presión entre dos
puntos consecutivos nos daremos cuenta que en cuanto a mayor sea
el radio menor será la presión, por ende la mayor caída de presión
será en el radio más lejano al pozo. Esto también se ve reflejado al
momento de hallar la velocidad del yacimiento que es el cociente
entre cauda y área
La variación de presión se trabajará con respecto a esta tabla
 R (FT) P (PSI) R (FT) P (PSI)
0,49 1100 550 1238,23
50 1191,04 600 1239,94
100 1204,68 650 1241,52
150 1212,66 700 1242,98
200 1218,35 750 1244,33
250 1222,71 800 1245,60
300 1226,30 850 1246,80
350 1229,33 900 1247,92
ΔP = 400 1231,96 950 1248,99
450 1234,28 1000 1250
500 1236,35
P50−P0,49=1191,04 −1100=91,04 PSI
ΔP = P100−P50=1204,68−1191,04 =13,64 PSI
ΔP = P150 −P100 =1212,66−1204,68=7,98 PSI
ΔP = P200 −P150 =1218,35−1212,66=5,69 PSI
ΔP = P250 −P200 =1222,71−1218,35=4,36 PSI
ΔP = P300 −P250 =1226,30−1222,71=3,59 PSI
ΔP = P350 −P300 =1229,33−1226,30=3,03 PSI
ΔP = P400 −P350 =1231,96−1229,33=2,63 PSI
ΔP = P450 −P400 =1234,28−1231,96=2,32 PSI
ΔP = P500 −P450 =1236,35−1234,28=2,07 PSI
ΔP = P550 −P500 =1238,23−1236,35=1,88 PSI
ΔP = P600 −P550 =1239,94−1238,23=1,71 PSI
ΔP = P650 −P600 =1241,52−1239,94=1,58 PSI
ΔP = P700 −P650 =1242,98−1241,52=1,46 PSI
ΔP = P750 −P700 =1244,33−1242,98=1,35 PSI
ΔP = P800 −P750 =1245,60−1244,33=1,27 PSI
ΔP = P850 −P800 =1246,80−1245,60=1,20 PSI
ΔP = P900 −P850 =1247,92−1246,80=1,12 PSI
ΔP = P950 −P900 =1248,99−1247,92=1,07 PSI
ΔP = P1000 −P950 =1250−1248,99=1,01 PSI
 A continuación se detalla una tabla ubicando los valores de la
variación de presión

R (FT) P (PSI) ΔP(PSI R (FT) P (PSI) ΔP(PSI)

)
0,49 1100 91,04 550 1238,23 1,71
50 1191,04 13,64 600 1239,94 1,58
100 1204,68 7,98 650 1241,52 1,46
150 1212,66 5,69 700 1242,98 1,35
200 1218,35 4,36 750 1244,33 1,27
250 1222,71 3,59 800 1245,60 1,20
300 1226,30 3,03 850 1246,80 1,12
350 1229,33 2,63 900 1247,92 1,07
400 1231,96 2,32 950 1248,99 1,01
450 1234,28 2,07 1000 1250
500 1236,35 1,88

Ejercicio 2.5
Usando la ecuación de Wyllie, calcule los datos de permeabilidad
relativa para una muestra de núcleo cuyos datos de saturación consta
en la tabla 2.5 Asuma que el núcleo es una arenisca bien cementada y
que el agua irreducible satura 30% a la muestra.

Slider, H.C. Wordlwide Practical Petroleum Reservoir Enginieering

Methds, 1983, PennWell Books.
 Swi =0.30
So
 Ec. De Wyllie para k ro y k rg⋮ S = 1−S
¿

wi

Tabla 2.5 Permeabilidades relativas en distintas tipos de arenas en función

de saturación

k ro k rg

Arena no consolidada - bien clasificada ( S¿ )3.0 ( 1−S ¿ )3

Arena no consolidada - pobremente ( S¿ )3.5 ( 1−S ¿ )2 ( 1−S¿1.5 )

clasificada

Arena cementada, caliza oolítica y caliza ( S¿ )4.0 ( 1−S ¿ )2 ( 1−S¿2 )

vugular

Sw −S wi
 Ec. de Wyllie para k ro y k rw⋮ S = 1−S
¿

wi

k ro k rw

Arena no consolidada - bien clasificada ( 1−S ¿ )3 ( S¿ )3.0

Arena no consolidada - pobremente ( 1−S ¿ )2 ( 1−S¿1.5 ) ( S¿ )3.5

clasificada

Arena cementada, caliza oolítica y caliza ( 1−S ¿ )2 ( 1−S¿2.0 ) ( S¿ )4.0

vugular

Tabla 2.6 Calculo de la permeabilidades para el ejercicio 2.5 propuesto.

 Saturación k rw k ro
Water Oil S¿ S¿ 4.0 ( 1−S ¿ )2 ( 1−S¿2 )

100 0 1 1 -

90 10 0.857 0.539 0.0054

80 20 0.714 0.26 0.04

60 40 0.429 0.034 0.266

40 60 0.143 0.0004 0.719

30 70 - 0 1

Swi=0.30
¿ S w −Swi
S=
1−Swi
¿
( 1−S wi ) S + S wi=S w
( 1−0,30 ) 1+0,30=1
1−0,30
S¿ = =1
1−0,30

k rw=S ¿4 =(1)4 =1
k ro=¿

¿ S w −Swi
S=
1−Swi
¿
( 1−S wi ) S + S wi=S w
( 1−0,30 ) 0,857+ 0,30=0,9
0,9−0,30
S¿ = =0,857
1−0,30

k rw=S ¿4 =(0,857) 4=0,539

k ro=¿

¿ S w −Swi
S=
1−Swi
 ¿
( 1−S wi ) S + S wi=S w
( 1−0,30 ) 0,714 +0,30=0,8
0,8−0,30
S¿ = =0,714
1−0,30

k rw=S ¿4 =(0,714 )4=0,26

k ro=¿

S w −Swi
S¿ =
1−Swi
¿
( 1−S wi ) S + S wi=S w
( 1−0,30 ) 0,429+ 0,30=0,6
0,6−0,30
S¿ = =0,429
1−0,30

k rw=S ¿4 =(0,429)4 =0,034

k ro=¿

S w −Swi
S¿ =
1−Swi
¿
( 1−S wi ) S + S wi=S w
( 1−0,30 ) 0,143+ 0,30=0.4
0,4−0,30
S¿ = =0,143
1−0,30

k rw=S ¿4 =(0,143)4 =0,0004

k ro=¿

S w −Swi
S¿ =
1−Swi
¿
( 1−S wi ) S + S wi=S w
( 1−0,30 ) 0+ 0,30=0,3
0,3−0,30
S¿ = =0
1−0,30

k rw=S ¿4 =(0)4 =0
k ro=¿
Ejercicio 2.6
Un pozo drena un yacimiento radial que consiste de dos capas
( camas del ingles beds ) con permeabilidades sin daño de 100 y 75 md.
Cerca del pozo la formación está dañada lo cual afecta a ambas zonas
hasta un radio de 33 pies, resultando en una reducción de la
permeabilidad a 60 y 30 md, respectivamente. Si el radio del pozo es
0.33 pies y no existe movimiento vertical de los fluidos, cuál es la tasa
de flujo continuo (steady state) de gas cuando la saturación de agua es
30% y pe y pw son 1000 y 800 psia, respectivamente. Asuma que no
existe caída presión adicional ocasionada por la turbulencia. Los
espesores de las zonas son 2 y 3 pies, respectivamente. El radio externo
de drenaje es 3300 pies, y la Fig.2.15 es aplicable
( k ro=k rw en la Fig ) No fluye agua .
Slider, H.C. Wordlwide Practical Petroleum Reservoir Enginieering
Methds, 1983, PennWell Books
 Viscosidad del gas = 0.01 cp
 Factor de desviación del gas = 0.97
 Temperatura del yacimiento = 140°F

Nota; Antes de empezar a resolver haga un dibujo de la distribución de

las capas, con sus permeabilidades y espesores.

K . K r=K i

K i= ( 0.37 ) ( 0.100 )

K i=0.037

q i=¿q + q ¿
1 2

q g=¿ kh¿¿¿

re k
S=ln ( −1)
rw ks

3300 ft 0.100
s1=ln (
33 ft 0.060
−1 )
 3300 ft
s1=ln ( 0.667 )
33 ft
s1=( 4.6 ) ( 0.667 )

s1=(3.07)
q g=¿ kh¿¿¿

q g=¿ 0.037d (2 ft)¿ ¿¿

q 0.037 d (2 ft)(360000)
g 1=¿ ¿
(0.01)(1.38)(600) [ ln 100+3.07 ]

q g 1=¿694,09 f t / s ¿ 3

CAPA 1
K= 100 md Pw=800 psia
Rw= 33 ft; 0,33 ft h= 2 ft
Re= 3300 ft μ=¿0.01 cp

Kdañ= 60 md z= 0,97
Sw= 30% t= 140 °F
Pe= 1000 psia
q i=¿q + q ¿
1 2

q 0,703 kh (Pe 2−P wf 2)

g=¿ ¿
re
[
μ zt ln
rw
+s
]
re k
S = Ln r ( k −1)
w s

3300 ft 0.100
s1=ln
33 ft 0.060
−1 ( )
3300 ft
s1=ln ( 0.667 )
33 ft
s1=( 4.6 ) ( 0.667 )

s1=(3.07)
q 0,703 kh (Pe 2−P wf 2)
g=¿ ¿
re
[
μ zt ln
rw
+s
]

q g=¿ 0,703(0.1 d)(2 ft)¿¿ ¿

q (0.703)(0.1)(2 ft)(360000 ft 2)
g 1=¿ ¿
(0.01cp)(0,97)(600° R) [ ln100+3.07 ]

q g 1=¿694,10 f t / s ¿ 3

CAPA 2
K= 75 md Pw=800 ft
Rw= 33 ft h= 3 ft
Re= 3300 ft μ=¿0.01 cp

Rdañ= 30 md z= 0,97
Sw= 30% t= 140 °F
Pe= 1000 ft
q i=¿q + q ¿
1 2

q 0,703 kh (Pe 2−P wf 2)

g=¿ ¿
re
[
μ zt ln
rw
+s
]
re k
S = Ln r ( k −1)
w s

3300 ft 0.060
s1=ln
33 ft 0.030
−1 ( )
3300 ft
s1=ln ( 1)
33 ft
s1=( 4.6 ) ( 1 )

s1=(4.6)

q 0,703 kh (Pe 2−P wf 2)

g=¿ ¿
re
[
μ zt ln
rw
+s
]
q g=¿ 0.703(0.06 d)(3 ft)¿¿ ¿

q 0.703 (0.06 )(3 ft)(360000 ft 2)

g 1=¿ ¿
(0.01cp)(0.97)(600° R) [ ln100+ 4.6 ]

q g 1=¿1683,01 f t / s¿
3