1 IPR Oil y Gas PI

Comportamiento de las Formaciones Productoras
Bibliografía
Producción y mantenimiento de pozos petroleros, T.E.W. Nind
The Technology of artificial Lift Methods, Kermit E. Brown
Petroleum Engineers Handbook, Cap 32
Tarek Ahmed, cap 7 y 8
Fundamentals of Reservoir Engineering, L.P.Dake
1
Introducción
Un pozo productor de petróleo es solo una parte de un

sistema complejo el cual comprende:
El Yacimiento
El pozo
Las instalaciones de superficie.
Técnicas de Producción
Pozos Fluyentes
Bombeo Neumático
Bombeo Mecánico
Bombeo Electrocentrífugo
Bombeo por Cavidades Progresivas
Bombeo Hidráulico
2
Introducción
Cada Técnica tiene un amplio rango de aplicación
Diagrama de Babson
3
Introducción
Otros parámetros a considerar son:
1. Profundidad del pozo
2. RGL actual y futura
3. Problemas (arena, parafinas, gas, tipo de fluido…)
4. Desviación del pozo
5. Diámetro de la tubería de revestimiento
6. Comportamiento del flujo de entrada (IPR, IP, AOF)
7. Presión de formación y declinación
8. Planes de recuperación (secundaria, terciaria,
mantenimiento de presión)
9. Estimación de acumulada, tiempo de vida
10. Disponibilidad de gas a alta presión
11. Viscosidad del crudo
12. WOR actual y futura
13. Criterios económicos y política de la empresa
14. Otros
4
Índice de Productividad
Índice de productividad, IP - J
Una simple aproximación para describir el comportamiento del
flujo de entrada de un pozo petrolero es el uso del concepto de índice
de productividad (IP).
Definiendo el abatimiento como: ∆ P = p r − p wf
bl
Se define el IP como: IP = J =
q
[= ] día
pr − p wf psi
q= ql=qo+qw
Se define el IP específico como:
bl
IP = J =
q
[= ] día
h ⋅( pr − p psi ⋅ ft
s s
wf )
5
Índice de productividad, IP – J
Desarrollado asumiendo la siguiente simplificación:
• Flujo radial alrededor del pozo

• Flujo de una sola fase ligeramente compresible
• Distribución de la permeabilidad homogénea
• La formación esta saturada completamente con el fluido
Para estas condiciones, la ecuación de Darcy queda

resumida a:
0 . 007082 ⋅ k ⋅ h ⋅ ( p r − p wf )
q=
µ ⋅ β o ⋅  Ln  re r  − 3 / 4 + S 
  w  
6
0 . 007082 ⋅ k ⋅ h ⋅ ( p r − p wf )
q=
µ ⋅ β o ⋅  Ln  re r  − 3 / 4 + S 
  w  
La mayoría de los parámetros a la derecha son constantes (en un
tiempo determinado, para un intervalo de presión acotado), lo que permite
juntarlos en un solo coeficiente llamado Índice de Productividad (IP)
q = IP ⋅ ( p r − p wf ) q = J ⋅ ( p r − p wf )
Esta ecuación nos dice que el flujo de líquido que entra en el
pozo es directamente proporcional al abatimiento de presión .
(El máximo caudal es el Potencial Absoluto a Pozo abierto (Absolute Open

Flow Potential –AOFP)
7
IPR (Inflow Performance Relationship)
Este se define como la relación funcional entre el
caudal de producción y la presión dinámica de fondo.
Gilbert, en 1954 fue el primero en proponer esta relación.
El IPR se define en el rango de presión desde la
presión promedio del Reservorio y la presión atmosférica.
El caudal correspondiente a la presión atmosférica,
como presión dinámica de fondo, se define como
Potencial del Pozo (Qmax, AOF),
El caudal a la presión promedio del reservorio en
fondo, es siempre cero.
Una relación de performance de influjo típica (IPR)
es:
8
IPR (Vogel)
3000
2700
2400
2100
1800
Pwf (psi)
1500
1200
900
600
300
0
0 300 600 900 1200 1500 1800
q (bl/día)
9
Factores que influyen en la forma del IPR

En zonas de Permeabilidad constante
• Una Fase
• Dos Fases
• Gas
10

Formación Estratificada
11
Flujo Turbulento
12
Factores que influyen en la
forma del IPR
Tipo de Reservorio.
13
Empleo del IPR

Muskat propone que el IP es excelente para determinar problemas
como:
1. Comparación antes y después de un tratamiento.
2. GOR Estable y J declinando indica taponamiento en el pozo.
3. Una declinación del IP es normal durante la depleción del reservorio
y debe ir acompañada del crecimiento del GOR o de la relación agua
petróleo. Si no pasa esto debe considerarse que el pozo se esta
tapando.
4. Un incremento rápido en la producción de agua debe ir acompañado
de una declinación del IP. Si este se mantiene entonces el agua
viene de otro estrato
5. GOR incrementándose marcadamente sin declinación del IP, indica
la entrada de un gas extraño
14
Predicción del IPR

Se presentan dos problemas diferentes
1. Determinar la forma del IPR actual
2. La estimación de la IPR futura, cuando la

producción continua
15
Predicción del IPR de Pozos de
IPR de Vogel
petróleo
En el caso del flujo bifásico, donde la presión promedio del

reservorio (pr) está por debajo de la presión del punto de burbuja,
es recomendado el IPR de Vogel. Vogel desarrolló una ecuación
empírica:
2
q  p wf   p wf 
= 1 − 0 , 2 ⋅   − 0 , 8 ⋅  
q′  pr   pr 
Para cualquier caso, se necesita realizar un ensayo de
recupero de presión (build up) para conocer pr y también se
debe realizar un ensayo de flujo por orificio calibrado, para
conocer el caudal entregado por la formación a una determinada
pwf.
16
Predicción del IPR de Pozos de
IPR de Vogel
petróleo
IPR (Vogel)
3000
2700
2400
2100 2
q  p wf   p wf 
1800
= 1 − 0 , 2 ⋅   − 0 , 8 ⋅  
Pwf (psi)
1500
q′  pr   pr 
1200
900
600
300
0
0 300 600 900 1200 1500 1800
q (bl/día)
17
Predicción del IPR
IPR de Vogel
Aplicación
Reservorios con empuje por gas disuelto (Aplicable a otros con
reservas)
Presión promedio debajo del PB
Dos fases fluyendo en el reservorio (gas y petróleo) (aplicable a tres
fases con reservas)
Amplio rango de PVT de petróleos.
Reservorios con distinto KR
No aplica
o Reservorios donde el daño da una IPR recta
o Reservorios con crudos muy viscosos
o Reservorios con Presión promedio sobre el PB
18
Predicción del IPR
IPR de Standing
Standing extiende el trabajo de Vogel considerando el
daño o mejora del pozo
Vogel asume eficiencia al flujo (FE) igual a 1, Standing
propone otras eficiencia ∆ p ideal p r − p ′wf
Define FE FE = =
∆p p − p
real r wf
q 0 . 007082 ⋅ k ⋅ h
IP = =
IP sin daño ′ )
( p r − p wf µ ⋅ β o ⋅ Ln  re r 
 w 
q 0 . 007082 ⋅ k ⋅ h
IP con daño IP =
( p r − p wf )
=
µ ⋅ β o ⋅  Ln  re r  + S 
  w  
 re 
Ln  

 rw 
FE =
  r  
 Ln  e  + S 

  rw  
19
Predicción del IPR
IPR de Standing
20
Predicción del IPR
IPR de Standing
Ej 7 pag 20
Ej 8 pag 20
Ejemplo
Q = 70 bl/d; Pp = 2400 psi; Pwf = 1800 psi; FE = 0.7
a) Hallar qmax 0.7
b) Hallar qmax 1.3
c) Q cuando Pwf = 1200 psi
21
Predicción del IPR Futuro
IPR de Standing
2
q  p wf   p wf  q  p wf   p wf 
= 1 − 0 , 2 ⋅   − 0 ,8 ⋅   =  1 −  ⋅  1 + 0,8 ⋅ 
q′  pr   pr  q′  pr   pr 
q q′  p wf  1,8 ⋅ q ′
IP = = 
 1 + 0,8 ⋅  Si pwf tiende a ps J =
*
p r − p wf pr  pr  pr
De la ecuación de flujo radial

1,8 ⋅ IP
J* = 0 . 007082 ⋅ k o ⋅ h
 p wf  J *
=
 1 + 0 ,8 ⋅ 
pr µ o ⋅ β o ⋅  Ln  re r  − 3 / 4 + S 
    w  
Donde Βo, µo y ko se toman a la
presión promedio del reservorio y
considerando una saturación media.
Sobre todas las cosas cuida tu mente,
Los demás parámetros son constantes
porque ella es fuente de vida 22
IPR de Standing
 ko 
 
J *
 µo ⋅β o  1,8 ⋅ q ′f
=
f f
*
= ⇒ J * f ⋅ pr f
J *
p  ko  J f
q ′f = ⇒
  pr f
1 .8
 µo ⋅β o  p
  p wf   p wf 
2

qo f = q ′f ⋅ 1 − 0 , 2 ⋅   − 0 ,8 ⋅   
  pr   pr  
  f   f  
23
IPR de Standing
24
IPR de Standing
IPR futuras (Standing)

2500
IPR a 1900 psi
2000 IPR a 1700 psi
IPR a 1500 psi
1500
Pwf (psi)
IPR a 1300 psi
1000 IPR a 1100 psi
500
0
0 200 400 600 800 1000
q (bl/día)
25
Predicción del IPR
En el caso del flujo bifásico, donde la presión promedio del
reservorio (pr) está por encima de la presión del punto de burbuja,
Para cualquier caso, se necesita realizar un ensayo de

recupero de presión (build up) para conocer pr y también se
debe realizar un ensayo de flujo por orificio calibrado, para
conocer el caudal entregado por la formación a una determinada
pwf.
IPR actual
4000
3500
3000
2500
Pwf (psi)
2000
1500
1000
500
0
0 200 400 600 800 1000 1200 1400
q (bl/día)
26
Predicción del IPR
Conocida pr y q para pwf > pb
q
J =
p r − p wf qb = J ⋅ (pr − pb )
∗ 1 , 8 ⋅ q v′ J ∗ ⋅ pb
J = J = q v′ =
pb 1,8
Si pwf > Pb q = J ⋅ (pr − pb )
  p wf   wf  
p
2
Si pwf < Pb
q = q b + q v′ ⋅ 1 − 0 , 2 ⋅   − 0 , 8 ⋅   
  pb   p b  
27
Predicción del IPR
Conocida pr y q para pwf > pb
IPR actual
4000
3500
3000
2500
Pwf (psi)
2000
1500
1000
500
0
0 200 400 600 800 1000 1200 1400
q (bl/día)
28
Predicción del IPR
Conocida pr y q para pwf < pb
q = qb + qv
  p wf   p wf 
2

q = q b
′
+ q v ⋅ 1 − 0 , 2 ⋅ 
  − 0 , 8 ⋅   
  p b   p b  
1, 8 ⋅ q v′ qb
q = q b + q v′ ⋅ a J ∗
= J = J = J ∗
pb pr − pb
pr − pb
q b = 1,8 ⋅ ⋅ q v′
pb
29
Predicción del IPR
Conocida pp y q para pwf < pb
pr − pb
q b = 1,8 ⋅ ⋅ q v′ q = q b + q v′ ⋅ a
pb
pr − pb
q = 1,8 ⋅ ⋅ q v′ + q v′ ⋅ a
q pb
⇒ q v′ =
p − pb
1 ,8 ⋅ r + a
pb
Si pwf < Pb   p wf   p wf 
2

q = q b + q v′ ⋅  1 − 0 , 2 ⋅   − 0 , 8 ⋅   
  pb   pb  
30
Predicción del IPR
Conocida pp y q para pwf < pb
IPR actual
3500
3000
2500
2000
Pwf (psi)
1500
1000
500
0
0 200 400 600 800 1000 1200
q (bl/día)
31
Predicción del IPR
IPR Log-Log
Método Propuesto por Gilbert
Modifica la ecuación
q = IP ⋅ ( p r − p wf ) por q = J ⋅ ( p r − p wf ) n
Donde n, no necesariamente es 1
Se grafica el caudal en función del abatimiento en escala log-log,

como línea recta, y la pendiente se desvía suavemente de 1
Para emplear este método son necesarios al menos tres ensayos
de pozo a distinto q.
32
Predicción del IPR
IPR Log-Log
y = 13.22533x0.61893
q vs abatimiento
100000
10000
1000
q (bl)
100
10
0
1 10 100 1000 10000
Abatimiento (psi)
Q = 0 Bl Q = 2 0 0 0 0 Bl Q = 4 0 0 0 0 Bl Q = 6 0 0 0 0 Bl Q = 8 0 0 0 0 Bl Q = 10 0 0 0 0 Bl Q = 12 0 0 0 0 Bl
33
IPR Log-Log
q = J ⋅ ( pr − pwf ) n
1. Graficar los datos de ensayo de pozo y calcular J y n

2. Graficar IP vs acumulada de la historia del pozo. (Los IP
deben ser medidos a abatimientos semejantes)
3. Graficar las Pr vs acumulada de la historia del pozo.
IP a distintas acumuladas -0.0000114x
PS vs Q y = 2794.3534482e
-0.0000499x
y = 10.1562233e
100 10000
10
IP (bl/día/psi)
1 Ps (psi) 1000
0.1
100
0.01
0 20000 40000 60000 80000 100000 120000 140000
0 20000 40000 60000 80000 100000 120000
Q (bl)
Q (bl)
34
IPR Log-Log
q = J ⋅ ( pr − pwf ) n
4. Para la Ps futura deseada estimar el IP futuro

5. Con el IP futuro al abatimiento dado, calcular el qf
6. Con qf calcular Jf
qf
IP f = q f = IP f ⋅ ( p r f − p wf )
( p r f − p wf )
qf
Jf =
( pr f − pwf ) n
35
IPR Log-Log
q = J ⋅ ( pr − pwf ) n
7. Con Jf realizar la IPR futura.
IPRs a distintas acumuladas

1600
Q = 60000
1400
Q = 80000
1200 Q = 120000
1000 Q = 130000
Pwf (psi0
800
600
400
200
0
0 100 200 300 400 500
q (bl/día)
36
Predicción del IPR
IPR de Fetkovich
0 . 007082 ⋅ k ⋅ h pe
f ( p ) ⋅ dp
Parte de la ecuación
qo =
Ln  e 
r ∫ p wf
 rw 
0 . 007082 ⋅ k ⋅ h  pb k ro ( p , s ) pe k ro ( p ) 
qo =
Ln  e 
r

 ∫ p wf µ o ⋅β o
⋅ dp +
∫ pb µ o ⋅β o
⋅ dp 

 rw 
Para flujo de una sola fase sobre la presión de saturación
qo =
0 . 007082 ⋅ k ⋅ h
⋅
(
p e − p wf )
  re   µ o ⋅β o µo y βo medidos a
Ln 
  rw   + S  Pprom
 
37
Predicción del IPR
IPR de Fetkovich
Para flujo de dos fase bajo la presión de saturación
0 . 007082 ⋅ k ⋅ h  pb k ro ( p , s ) 
qo =

Ln 
re 



∫ p wf µ o⋅β o
⋅ dp 

 rw
Fetkovich supone que la función de la presión es lineal y pasa por el origen,

entonces hace
 k ro   k ro  p
  =  ⋅
µ β   
 o ⋅ o  p  µ o ⋅β o  pi pi
Luego
0 . 007082 ⋅ k ⋅ h  pi  k ro  p 
qo =
Ln  e 
r

 ∫ p wf
 µ β  ⋅ p ⋅ dp 
 
 o⋅ o i i 
 rw 
38
Predicción del IPR
IPR de Fetkovich
qo =
0 . 007082 ⋅ k ⋅ h  k (p, s )
⋅  ro  ⋅

1
⋅ (
p 2
− p wf 2
)
 µ o ⋅β o  i 2 ⋅ pi
i
Ln  e 
r
 rw 
Luego
(
′ ⋅ p i 2 − p wf
q o = J oi 2
)
39
Predicción del IPR
IPR de Fetkovich
IPR a 1900 psi (Fetcovich)

2500
2000
1500
Pwf (psi)
1000
500
0
0 100 200 300 400 500 600 700
q (bl/día)
40
IPR de Fetkovich
IPR Actual q o = J oi (
′ ⋅ p i 2 − p wf 2
)
IPR Futura (
′ ⋅ pf
q o = J of 2
− p wf 2
)
pf
Donde ′ = J oi
J of ′ ⋅
pi
41
IPR de Gas
Ec de Darcy para gas
k ⋅ h ⋅ (Ψ − Ψ )
=
r wf
q g
  r  3 
1422 ⋅ T ⋅  ln  e  − + S 
  rw  4 
qg caudal de gas (mft3/d)
k permeabilidad (mD)
2 P
Ψ Pseudo presión promedio (psi2/cp)
P
Ψ = ∫ 0 µ g ⋅ Z
⋅ dp
T temperatura (R)
S factor de daño
H espesor (ft)
re radio de drene (ft)
rw radio del pozo (ft)
Definiendo J
k ⋅h
qg J =
J =   re  3 
Ψ r − Ψ wf 1422 ⋅ T ⋅  ln   − + S
  rw  4 
42
IPR de Gas
q g = J ⋅ (Ψ r − Ψ wf )
qg
J =
Ψ r − Ψ wf
q g (max) = J ⋅Ψr
J indice de productividad (mft3/d/psi2/cp)
1
Ψ wf = Ψ r − ⋅qg
J
43
IPR de Gas
La ecuación de Darcy puede escribirse
k ⋅ h Pr 2 P
q g =
  r  3 
⋅ ∫ P wf µ g ⋅ Z
⋅ dp
1422 ⋅ T ⋅  ln  e  − + S 
  rw  4 
Z ⋅T
P/µg.Z es proporcional a 1/µg.Bg B g = 0 , 00504 ⋅
P
Bg Factor de volumen de gas (bl/ft3)
Z factor de compresibilidad
T temperatura (R)
Reemplazando queda
7 , 08 E − 6 ⋅ k ⋅ h Pr 1
q g =
 r  3
⋅ ∫ P wf µ ⋅ B
⋅ dp
ln  e  − + S g g
 rw  4
44
IPR de Gas
1 2 P
Graficando o En función de la presión
µ g ⋅ B g µ g ⋅ Z
45
IPR de Gas
Región III, Presión Alta
7 , 08 E − 6 ⋅k ⋅h Pr 1
q g =
 re  3
⋅ ∫ P wf µ ⋅B
⋅ dp
ln   − + S g g
 rw  4
Para Pwf > 3000 psi la función µg.Bg es cte

Por lo tanto
7 , 08 E − 6 ⋅ k ⋅ h ⋅ (P r − P wf )
q =
 
g
 r  3
µ g ⋅ B g ⋅  ln  e  − + S 
  rw  4 
46
IPR de Gas
Región II, Presión Intermedia
7 , 08 E − 6 ⋅k ⋅h Pr 1
q g =
 re  3
⋅ ∫ P wf µ ⋅B
⋅ dp
ln   − + S g g
 rw  4
Para 2000< Pwf < 3000 psi la función µg.Bg no

es cte. Por lo tanto debe emplearse
k ⋅ h ⋅ (Ψ − Ψ )
=
r wf
q g
  r  3 
1422 ⋅ T ⋅  ln  e  − + S 
  rw  4 
47
IPR de Gas
Región I, Presión Baja
7 , 08 E − 6 ⋅k ⋅h Pr 1
q g =
 re  3
⋅ ∫ P wf µ ⋅B
⋅ dp
ln   − + S g g
 rw  4
Para Pwf < 2000 psi la función µg.Bg tiene un

comportamiento lineal
Por lo tanto debe emplearse
k ⋅ h ⋅ ( Pr − P wf
2 2
)
q g =
  r  
⋅ (µ ) 3
1422 ⋅T g ⋅ Z prom
⋅  ln  e  − + S 
  rw  4 
µg.Z tiene que ser evaluado a la presión promedio, se recomienda calcular esta como
+ P wf
2 2
Pr
P prom =
2
48
IPR de Gas
Modificaciones por Turbulencia
I) Aproximación por Cuadrado de la Presión
k ⋅ h ⋅ ( Pr − P wf
2 2
)
q g =
  r  
⋅ T ⋅ (µ ) 3
1422 g ⋅ Z prom
⋅  ln  e  − + S + D .q g 
  rw  4 
D es el factor de Turbulencia
F . K .h  β .T . γ g 
D = F = 3 , 161 E − 12 .  
1422 .T  µ .h 2 .r 
 g w 
49
IPR de Gas
I) Aproximación por Cuadrado de la Presión
F . K .h  β .T . γ g 
D = F = 3 , 161 E − 12 .  
 µ .h 2 .r 
1422 .T  g w 
D es el factor de Turbulencia β = 1 , 88 E − 10 . K − 1 , 47
.φ − 0 , 53
F coeficiente de flujo No-Darcy
β parámetro de turbulencia
50
IPR de Gas
II) Aproximación por Presión
7 , 08 E − 6 . k ⋅ h ⋅ ( P r − P wf )
q g =
  r  
(µ g ⋅ B g ) prom
⋅ T .  ln  e  −
3
+ S + D .q g 
  rw  4 
51
IPR de Gas
III) Aproximación por Pseudo-Presión
k ⋅ h ⋅ (Ψ − Ψ )
q g = r wf
  r  3 
1422 ⋅ T ⋅  ln  e  − + S + D .q g 
  rw  4 
Las tres aproximaciones son cuadráticas en qg y aceptan dos tratamientos

distintos
1)Tratamiento simplificado
2)Tratmiento Laminar-Inercial-Turbulento (Tratamiento LIT)
52
Predicción del IPR de Gas
1) Tratamiento simplificado
qg = C ⋅ ( Pr − Pwf )
2 2 n
Aplicando Logaritmo Log(qg ) = Log(C ) + n ⋅ Log Pr − Pwf ( 2 2

)
Isocronal
10000
(Ps2 - Pwf2)/1000
1000
100
10
10 100 1000
q (bl/día)
53
qg = C ⋅ ( Pr − Pwf )
2 2 n
El exponente n intenta incluir los efectos de turbulencia

Varía entre 0,5 y 1 (1 completamente laminar, 0,5 turbulento)
El coeficiente C incluye las propiedades del Reservorio, del fluido

y la geometría del flujo
El exponente n y el coeficiente C pueden determinarse por un ensayo

Isocronal
Isocronal modificado
Back-pressure test (Ensayo de flujo estabilizado)
54
qg = C ⋅ ( Pr − Pwf )
2 2 n
IPR de Gas
2500
2000
Pwf (psi)
1500
1000
500
0
0 2000 4000 6000 8000 10000
qg (mcft/d)
55
2) Tratamiento Laminar-Inercial-Turbulento
Puede ser tratado de tres formas distintas

2a) Forma de Presión cuadrática
Pr − Pwf = a.qg + b.qg

2 2 2
1422.T .µ g .Z   re  3 
a= .ln  − + S 
k.h   rw  4 
1422.T .µ g .Z El término a.qg da cuenta del flujo laminar
b= .D
k.h El término b.qg2 da cuenta de los efectos inerciales y
de turbulencia
56
Determinación de a y b
Pr − Pwf
2 2
= a + b.qg qg =
2
(
− a + a + 4.b. Pr − Pwf
2 2
)
qg 2.b
Esta ecuación es válida para
Una sola fase fluyendo

Reservorio Homogéneo
Permeabilidad independiente de la presión
µg.Z constante
Es aplicable para presiones menores a 2000 psi
57
qg =
− a + a 2
(
+ 4 .b . P r
2
− P wf
2
)
2 .b
IPR de Gas
2500
2000
Pwf (psi)
1500
1000
500
0
0 2000 4000 6000 8000 10000
qg (mcft/d)
58
2b) Forma de Presión
Pr − Pwf = a.qg + b.qg

2
141,2E − 3.µ g .Bg   re  3 

a= .ln  − + S 
k.h   rw  4 
141,2E − 3µ g .Bg
b= .D El término a.qg da cuenta del flujo laminar
k.h
El término b.qg2 da cuenta de los efectos inerciales y
de turbulencia
59
2b) Forma de Presión

Pr − Pwf
= a + b.qg − a + a 2 + 4.b.(Pr − Pwf )
qg qg =
2.b

µg.Z constante
Es aplicable para presiones mayores a 3000 psi
60
2b) Forma de Presión − a + a 2 + 4 .b . (Pr − P wf )

qg =
2 .b
IPR de Gas
2500
2000
Pwf (psi)
1500
1000
500
0
0 2000 4000 6000 8000 10000 12000
qg (mcft/d)
61
2c) Forma de Pseudo-Presión
Ψr − Ψwf = a.qg + b.qg

2
1422   re  3 
a= .ln  − + S 
k.h   rw  4 
1422
b= .D El término a.qg da cuenta del flujo laminar
k.h
El término b.qg2 da cuenta de los efectos inerciales y
de turbulencia
62
2c) Forma de Pseudo-Presión
Ψr − Ψwf
= a + b.qg − a + a 2 + 4.b.(Ψr − Ψwf )
qg qg =
2.b

Es aplicable para cualquier rango de presiones
63
2c) Forma de Pseudo-Presión − a + a 2 + 4 .b . (Ψ r − Ψ wf )
qg =
2 .b
IPR de Gas
2500
2000
Pwf (psi)
1500
1000
500
0
0 2000 4000 6000 8000 10000 12000 14000
qg (mcft/d)
64
Comparación
IPR de Gas
2500
2000 qg 1
qg 2a
Pwf (psi)
1500
qg 2b
1000 qg 2c
500
0
0 2000 4000 6000 8000 10000 12000 14000
qg (mcft/d)
65
Back-Pressure test
1. Cerrar el pozo suficiente tiempo para que la presión de

formación se ecualice a la presión promedio volumétrica.
2. Tener el pozo en producción suficiente tiempo para que la Pwf

alcance el estado pseudo-estacionario
3. Repetir el paso dos para tres o cuatro caudales distintos
1200.Φ.S g .µ g .re
2
ts =
k.Pr
66
IPR Actual q g = C a ⋅ Pr − Pwf ( 2 2 n

)
IPR Futura q g = C f ⋅ Prf − Pwf ( 2 2 n

)
µ .Z µ .B . P rf
= C ⋅ = C ⋅
f f f gf
Donde C C
µ µ
f a f a
a .Z a a .B ga . P ra
67
IPR Actual q g = C a (
⋅ Pr
2
− P wf
2
)
n
IPR Futura q g = C f (
⋅ P rf
2
− P wf
2
)
n
IPR de Gas
2500
2000
Pwf (psi)
1500
1000
500
0
0 2000 4000 6000 8000 10000
qg (mcft/d)
68
2)Tratamiento LIT
2a) Cuadrado de la Presión
Pr − Pwf = a .q g + b .q g
2 2 2
IPR Actual
Pr − Pwf = a f .q g + b f . q g
2 2 2
IPR Futura
µ .Z µ .B . P rf
= a ⋅
f gf
= a ⋅ a
f f
Donde a f a
µ
µ
f a
µ .Z µ .B . P rf
= ba ⋅
f gf
= ba ⋅ b
f f
b
µ
f
µ
f
69
2)Tratamiento LIT
2a) Cuadrado de la Presión
Pr − Pwf = a f .qg + b f .qg
2 2 2
IPR Actual Pr − Pwf = a .q g + b.q g
2 2 2
IPR Futura
IPR de Gas
2500
2000
Pwf (psi)
1500
1000
500
0
0 2000 4000 6000 8000 10000
qg (mcft/d)
70
2) Tratamiento LIT
2b) Presión
Pra − Pwf = a .q g + b .q g
2
IPR Actual
Ψ rf − Pwf = a .q g + b .q g
IPR Futura 2
µ .Z µ .B . P rf
= a ⋅
f gf
= a ⋅ a
f f
Donde a f a
µ
µ
f a
µ .Z µ .B . P rf
= ba ⋅
f gf
= ba ⋅ b
f f
b
µ
f
µ
f
71
2) Tratamiento LIT
2b) Presión
Pra − Pwf = a .q g + b .q g IPR Futura Ψ rf − Pwf = a .q g + b .q g
2 2
IPR Actual
IPR de Gas
2500
2000
1500
Pwf (psi)
1000
500
0
0 2000 4000 6000 8000 10000 12000
qg (mcft/d)
72
Tratamiento LIT
c) Pseudo-Presión
Ψ ra − Ψ wf = a .q g + b .q g
2
IPR Actual
Los coeficientes a y b en este método son escencialmente

independientes de la presión del reservorio.,
Ψ rf − Pwf = a .q g + b .q g
2
IPR Futura
73
Tratamiento LIT
c) Pseudo-Presión
IPR Actual Ψra −Ψwf =a.qg +b.qg

2
IPR Futura Ψrf − Pwf = a.q g + b.q g
2
IPR de Gas
2500
2000
Pwf (psi)
1500
1000
500
0
0 2000 4000 6000 8000 10000 12000 14000
qg (mcft/d)
74
Comparación
IPR de Gas
2500
2000
qg 1
qg 2a
Pwf (psi)
1500
qg 2b
1000 qg 2c
500
0
0 2000 4000 6000 8000 10000 12000 14000
qg (mcft/d)
75
Comparación
IPR de Gas
2500
2000
qg 1
qg 2a
Pwf (psi)
1500
qg 2b
1000 qg 2c
500
0
0 2000 4000 6000 8000 10000 12000 14000
qg (mcft/d)
76
Preguntas
1. Defina IP, que condiciones deben cumplirse para su definición?
2. ¿Qué parámetros deben considerarse constantes en la ec de flujo radial
para que q=IP*∆∆P?
3. Defina IPR, ¿Qué es el potencial de la formación?
4. ¿Qué factores influyen en la forma del IPR? ¿Cómo influyen?
5. ¿Qué es una IPR compuesta?
6. ¿Ejemplifique y explique algunas aplicaciones del IPR?
7. ¿Bajo que condiciones se puede aplicar la ecuación de Vogel para
calcular una IPR?
8. ¿Cómo determina una IPR de Vogel?
9. ¿Qué es eficiencia al Flujo?
10. ¿Cómo define Standing Una IPR para una eficiencia distinta de 1?
11. ¿Qué consideraciones hace Standing para calcular una IPR futura? ¿Por
qué correlaciona los J* y no otro J?
12. ¿Cuáles son las consideraciones para calcular una IPR por el método
propuesto por Gilbert?
13. Explicar los pasos necesarios para determinar una IPR Futura, empleando
el método de Gilbert.
14. ¿En qué se fundamenta la determinación de un IPR de Fetkovich?
15. ¿Cómo determina una IPR futura por el método de Fetkovich?
16. ¿Cómo determina una IPR para un pozo de gas?
17. ¿Qué determina el exponente n en el método simplificado?
18. ¿Cómo emplea un ensayo isocronal?
19. ¿Cómo puede determinar el grado de turbulencia en el flujo?
77

1 IPR Oil y Gas PI

Cargado por

Copyright:

Formatos disponibles

1 IPR Oil y Gas PI

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

1 IPR Oil y Gas PI

Cargado por

Copyright:

Formatos disponibles

Comportamiento de las Formaciones Productoras

Producción y mantenimiento de pozos petroleros, T.E.W. Nind

The Technology of artificial Lift Methods, Kermit E. Brown

Petroleum Engineers Handbook, Cap 32

Tarek Ahmed, cap 7 y 8

Fundamentals of Reservoir Engineering, L.P.Dake

Un pozo productor de petróleo es solo una parte de un

Se define el IP específico como:

Desarrollado asumiendo la siguiente simplificación:

• Flujo radial alrededor del pozo

Para estas condiciones, la ecuación de Darcy queda

(El máximo caudal es el Potencial Absoluto a Pozo abierto (Absolute Open

Factores que influyen en la forma del IPR

Factores que influyen en la forma del IPR

Factores que influyen en la forma del IPR

Empleo del IPR

Predicción del IPR

1. Determinar la forma del IPR actual

2. La estimación de la IPR futura, cuando la

En el caso del flujo bifásico, donde la presión promedio del

o Reservorios con crudos muy viscosos

o Reservorios con Presión promedio sobre el PB

De la ecuación de flujo radial

IPR futuras (Standing)

IPR a 1300 psi

1000 IPR a 1100 psi

Para cualquier caso, se necesita realizar un ensayo de

Si pwf > Pb q = J ⋅ (pr − pb )

Si pwf > Pb q = J ⋅ (pr − pb )

Se grafica el caudal en función del abatimiento en escala log-log,

1. Graficar los datos de ensayo de pozo y calcular J y n

4. Para la Ps futura deseada estimar el IP futuro

7. Con Jf realizar la IPR futura.

IPRs a distintas acumuladas

Fetkovich supone que la función de la presión es lineal y pasa por el origen,

IPR a 1900 psi (Fetcovich)

J indice de productividad (mft3/d/psi2/cp)

Para Pwf > 3000 psi la función µg.Bg es cte

Para 2000< Pwf < 3000 psi la función µg.Bg no

Para Pwf < 2000 psi la función µg.Bg tiene un

I) Aproximación por Cuadrado de la Presión

I) Aproximación por Cuadrado de la Presión

F coeficiente de flujo No-Darcy

II) Aproximación por Presión

III) Aproximación por Pseudo-Presión

Las tres aproximaciones son cuadráticas en qg y aceptan dos tratamientos

2)Tratmiento Laminar-Inercial-Turbulento (Tratamiento LIT)

Aplicando Logaritmo Log(qg ) = Log(C ) + n ⋅ Log Pr − Pwf ( 2 2

El exponente n intenta incluir los efectos de turbulencia

El coeficiente C incluye las propiedades del Reservorio, del fluido

El exponente n y el coeficiente C pueden determinarse por un ensayo

Puede ser tratado de tres formas distintas

Pr − Pwf = a.qg + b.qg

Una sola fase fluyendo

Es aplicable para presiones menores a 2000 psi

2b) Forma de Presión

Pr − Pwf = a.qg + b.qg

141,2E − 3.µ g .Bg   re  3 

2b) Forma de Presión