Taller de Productividad de Pozos Capítulo SPE

Ricardo Posadas Mondragn Productividad de Pozos
Captulo Estudiantil SPE, UNAM

Productividad de Pozos

Ing. Ricardo Posadas Mondragn
Activo de Produccin Cantarell

correo: ProductividadUNAM@hotmail.com
Contenido
1) Introduccin a la ingeniera de produccin
2) Comportamiento de afluencia
3) Factor de dao y su relacin con el comportamiento
de afluencia
4) Curvas de declinacin
5) Registros de produccin
6) Anlisis integral del pozo

Contenido
1) Introduccin a la ingeniera de produccin
1) Funciones del ingeniero de produccin
2) Descripcin de los sistemas de produccin


1) Solucin de ecuacin de difusin
2) Yacimientos bajosaturados
3) Yacimientos saturados
4) Yacimientos de gas
2) IPR futuro
3) Pruebas en pozos

* IPR Pozos Horizontales
Contenido
3) Factor de dao y su relacin con el comportamiento de afluencia
1) Factores que provocan dao
2) Obtencin del factor de dao a partir de pruebas de presin
3) Anlisis de los factores de Pseudodao
4) Efecto del dao sobre el comportamiento de Afluencia

4) Curvas de Declinacin
1) Declinacin Exponencial
2) Declinacin Hiperblica
3) Declinacin Armnica

5) Registros de Produccin
1) Registros de presin
2) Registros de temperatura
3) Registros de gasto a condiciones de pozo
4) Combinacin de registros
Contenido
6) Anlisis integral del pozo
1) Importancia de la caracterizacin del fluido
2) Flujo en el yacimiento
3) Flujo en la tubera de produccin
4) Flujo en estranguladores
5) Metodologa de anlisis nodal

Bibliografa:

1) Golan, M. , Whitson, C.H., Well Performance, New Jersey, Prentice
Hall, 1991.
2) Economides, H.,Petroleum Production Systems, Tulsa, Pennwell, 1994.
3) Beggs, H.D., Production optimization using nodal analysis, Tulsa,
OGCI, 1991.
4) Cholet, H., Well production practical handbook, France, IFP,2000.
I.- INTRODUCCIN A LA INGENIERA DE
PRODUCCIN
Analizar e interpretar los datos de produccin de pozos e
instalaciones, para garantizar las condiciones ptimas de
operacin de pozos, garantizando la continuidad operativa del
sistema de produccin.

Condiciones Optimas de Pozo

Mxima produccin posible del pozo, sin generar un efecto
negativo en la recuperacin final del Yacimiento.
Funcin del ingeniero de produccin
D
e
s
c
u
b
r
i
m
i
e
n
t
o

D
e
l
i
m
i
t
a
c
i
n

D
e
s
a
r
r
o
l
l
o

D
e
c
l
i
n
a
c
i
n

R
S
,

R
M

tiempo
Qo
Etapas del yacimiento
Componentes de un sistema de produccin

Yacimiento:
Trampa de hidrocarburos, con condiciones para permitir el flujo de
fluidos a travs del medio poroso hacia el pozo.

Pozo:
Medio de comunicacin entre el yacimiento y la superficie.,
conformado por tubera de produccin y dispositivos de terminacin
para garantizar el flujo de los hidrocarburos hacia la superficie.

Principales componentes del sistema de
produccin
Estrangulador:
Regulador del dimetro de flujo en superficie con la finalidad de
regular presiones y gastos.

Lnea de Flujo o descarga: Tubera para el transporte de los
hidrocarburos del pozo al separador.

Separador:
Equipos para separar el aceite, agua y gas de la corriente de flujo.
Principales componentes del sistema de
produccin

Demostracin de la Ley de Darcy.

Dimensiones de la Permeabilidad
(
+ =
dz
dp kA
q
Experimento de la Ley de Darcy
Darcy, observ que el gasto que pasa a travs de
la arena es directamente proporcional al rea de
flujo y a la diferencia de alturas del fluido e
inversamente proporcional a la longitud del
empaque de arena.

(1)
Donde la constante de proporcionalidad est en
funcin de propiedades de la arena y del fluido
que pasa a travs de ella: Directamente
proporcional a la permeabilidad de la arena y a la
densidad del fluido e inversamente proporcional
a la viscosidad del fluido.

. .(2)
h
L
A
z
L
h
A q
A
A
~
L
h kA
q
A
A
=
g =
Sustituyendo (5) en (4):

..(6)
h
L
A
z
De la figura se observa que la altura h es:

.(3)
Obteniendo la derivada a partir de (3):

..(5)
z
g
p
h + =
Tomando el lmite cuando L0, en (2)

.(4)
dz
dh kA
q
=
1
1
+ =
dz
dp
dz
dh
+ =
dz
dp kA
q
Finalmente considerando que la diferencia de presiones en sentido contrario a la
direccin de flujo se considera el signo negativo:
(
+ =
dz
dp kA
q
Obteniendo la ecuacin de Darcy en su forma general.

Para el caso de flujo horizontal, el efecto de gravedad es cero:
dx
dp kA
q
=
Dimensiones de la Permeabilidad:

Anlisis Dimensional:
p A
L q
k
A
=

Gasto:
Viscosidad:
rea:
Longitud:
Presin:
(
T
L
q
3
(
TL
M
| | L L
| |
2
L A
(
A
L T
M
p
2
Dimensiones de la Permeabilidad:

| |
| |
2
2 3
2 5
2
2
3
L
T ML
T ML
L T
M
L
L
TL
M
T
L
k = =
(
=
p A
L q
k
A
=

Dimensiones de longitud al cuadrado
Ecuacin de Darcy en Coordenadas Radiales
Partiendo de la ecuacin obtenida para flujo horizontal
Geometra de Flujo:
dx
dp kA
q
=
Introduciendo el rea de flujo en la ecuacin de Darcy:
rea de Flujo
2r
r
h
h
dr
dp rh k
q
t 2
=
Ecuacin de Darcy en Coordenadas Radiales
Separando para Integrar y considerando las cadas de presin en direccin del flujo:
Despejando el gasto:
} }
=
pe
pwf
re
rw
dp dr
r h k
q 1
2t
( )
wf e
w
e
p p
r
r
h k
q
=
|
|
.
|
\
|
ln
2t
( )
wf e
w
e
p p
r
r
kh
q
|
|
.
|
\
|
=
ln
2
t
Ecuacin de Darcy en su forma ms simple (Estado Estacionario)
Periodos de Flujo en un Pozo
Existen 3 Periodos de flujo en un pozo:
1) Transitorio (Yacimiento Infinito)

2) Pseudoestacionario (Frontera Cerrada)

3) Estacionario (Frontera Presin Constante)
iable
dt
dp
var =
cte
dt
dp
=
0 =
dt
dp
Periodos de Flujo en un Pozo
4605
4610
4615
4620
4625
4630
4635
4640
4645
0 20 40 60 80 100
P
w
f
Tiempo
Periodos de Flujo
iable
dt
dp
var =
cte
dt
dp
=
0 =
dt
dp
II.-COMPORTAMIENTO DE AFLUENCIA
Solucin de la Ecuacin de Difusin
Coordenadas Cilndricas:
) ( t un en todo masa de cambio sale que masa entra que Masa A =
Balance de Masa en las tres direcciones
Ecuacin de Continuidad en Coordenadas Cilndricas:
( ) ( ) ( ) ( ) |
u
u
t
v
z
v
r
v r
r r
z r
c
c
=
c
c
+
c
c
+
c
c 1 1
Ecuacin de Estado y Momento necesarias
|
.
|
\
|
c
c
= =
r
p k
v
A
q
r
r
dp
d
c
1
=
T
p
c
|
|
.
|
\
|
c
c
=
u
u
1
Ecuacin de Difusin en Coordenadas Cilndricas:
Caractersticas:

- Ecuacin diferencial Parcial
- 2 Orden
- 1er Grado
- Lineal
- Homognea
- 3 Variables Independientes
t
p
k
c
z
p p
r r
p
r r
p
t
c
c
=
c
c
+
c
c
+
c
c
+
c
c |
u
2
2
2
2
2 2
2
1 1
Consideraciones para Deduccin:

- Medio Homogneo e Isotrpico
- Fluido ligeramente compresible
- Flujo laminar
- Gradientes de presin pequeos
- Efectos de gravedad despreciables

Ecuacin de Difusin considerando flujo radial nicamente:
t
p
k
c
r
p
r r
p
t
c
c
=
c
c
+
c
c | 1
2
2
t
p
k
c
r
p
r
r r
t
c
c
=
|
.
|
\
|
c
c
c
c | 1
Solucin en Estado Estacionario:
t
p
k
c
r
p
r
r r
t
c
c
=
|
.
|
\
|
c
c
c
c | 1
0 =
dt
dp
0
1
=
|
.
|
\
|
c
c
c
c
r
p
r
r r
0
1
=
|
.
|
\
|
dr
dp
r
dr
d
r
0 =
|
.
|
\
|
dr
dp
r
dr
d
1
C
dr
dp
r =
( )
2 1
ln C r C p + =
..(1)
( )
2 1
ln C r C p + =
Para evaluar las constantes se utilizan las condiciones de frontera
( )
2 1
ln C r C p
w wf
+ =
( )
2 1
ln C r C p
e e
+ =
Restando las Ecuaciones (3)-(2):
|
|
.
|
\
|
=
w
e
wf e
r
r
C p p ln
1
|
|
.
|
\
|
=
w
e
wf e
r
r
p p
C
ln
1
..(1)
..(2)
..(3)
Sustituyendo C1 en (2)
Sustituyendo C1 y C2 en (1):
( )
2
ln
ln
C r
r
r
p p
p
w
w
e
wf e
wf
+
(
(
(
(
(
|
|
.
|
\
|
=
( )
w
w
e
wf e
wf
r
r
r
p p
p C ln
ln
2
(
(
(
(
(
|
|
.
|
\
|
=
( ) ( )
w
w
e
wf e
wf
w
e
wf e
r
r
r
p p
p r
r
r
p p
p ln
ln
ln
ln
(
(
(
(
(
|
|
.
|
\
|
+
|
|
.
|
\
|
=
|
|
.
|
\
|
|
|
.
|
\
|
+ =
w
w
e
wf e
wf
r
r
r
r
p p
p p ln
ln
..(4)
De la ley de Darcy:
dr
dp khr
q
t 2
=
Obteniendo la derivada a partir de (4):
|
.
|
\
|
|
|
.
|
\
|
=
r
r
r
p p
dr
dp
w
e
wf e
1
ln
Substituyendo en la Ley de Darcy:
( )
|
|
.
|
\
|

=
w
e
wf e
r
r
p p kh
q
ln
2
t
Ecuacin de Darcy en Estado Estacionario
Solucin en Estado Pseudoestacionario:
t
p
k
c
r
p
r
r r
t
c
c
=
|
.
|
\
|
c
c
c
c | 1
cte
dt
dp
=
( )
(
|
|
.
|
\
|

=
4
3
ln
2
w
e
wf
r
r
B
p p kh
q
t
Solucin en Estado Pseudoestacionario unidades de Campo:
( )
(
|
|
.
|
\
|
=
4
3
ln 2 . 141
w
e
wf
r
r
B
p p kh
q
| |
| |
| |
| | psi p
ft h
mD k
bpd q
| |
| |
| |
| | ft r
ft r
cp
stbl bl B
w
e
/
Inclusin del Dao en Ecuacin de Darcy para Flujo Pseudoestacionario:
( )
(
+
|
|
.
|
\
|
=
s
r
r
B
p p kh
q
w
e
wf
4
3
ln 2 . 141
s
p
qB
kh
s A =
2 . 141
Van Everdingen and Hurst
Curva de Afluencia a partir de la Ley de Darcy
0
500
1000
1500
2000
2500
0 5000 10000 15000 20000 25000
P
w
f

(
p
s
i
)
Qo (bpd)
Curvas de Afluencia Darcy en Estado Pseudoestacionario
Dao
0
5
10
-1
Diferentes Mtodos de Curvas de Afluencia
Los diferentes mtodos de Afluencia o IPR son los siguientes:

I. Mtodo de Afluencia de Darcy*
II. Mtodo de Afluencia por Vogel*
III. Mtodo de Afluencia por IP (ndice de Productividad)*
IV. Mtodo de Afluencia de IPR Generalizada*
V. Mtodo de Afluencia Composite*
VI. Mtodo de Afluencia de Fetkovich
VII. Mtodo de Afluencia de Jones*
VIII. Mtodo de Afluencia de Standing
IX. Mtodo de Afluencia Transitorio
X. Mtodo de Afluencia Pozos Fracturados Hidrulicamente
XI. Mtodo de Afluencia Pozos Horizontales*
XII. Mtodos de Afluencia Pozos de Gas
Mtodo de Afluencia de
Vogel
Mtodo de Afluencia de Vogel
Vogel (1968) observo que el comportamiento de afluencia o curvas
IPR (Inflow Performance Relationship) en un pozo con entrada de
gas en solucin, tenan un comportamiento diferente a la curva de
afluencia generada por el modelo de Darcy. La cual es una recta
(producto de considerar una sola fase y fluido ligeramente
compresible).

Artculo SPE 001476:
Resultados por programa de clculo (Simulador) para diferentes
fluidos y propiedades de la formacin.
Variables Adimensionales
Ejemplo:
-Pws=2000 psi
-Pwf=1500 psi
-Qo=100 bpd
Pwf (psi) Pwf/Pws Qo/Qomax Qo (bpd)
1500 0.75 0.39 100
(Qomax=256.4)
1000 0.5 0.69 177
500 0.25 0.9 231
0 0 1.0 256.4
Ecuacin de Ajuste de la Curva Tipo de Vogel
2
max
8 . 0 2 . 0 1
|
|
.
|
\
|
|
|
.
|
\
|
=
ws
wf
ws
wf
o
o
p
p
p
p
q
q
Forma general de la ecuacin de Vogel
( )
2
max
1 1
|
|
.
|
\
|

|
|
.
|
\
|
=
ws
wf
ws
wf
o
o
p
p
c
p
p
c
q
q
Metodologa de Aplicacin de la ecuacin de Vogel
2
max
8 . 0 2 . 0 1
|
|
.
|
\
|
|
|
.
|
\
|
=
ws
wf
ws
wf
o
o
p
p
p
p
q
q
1) Datos de Prueba: pws, pwf, qo
2) Determinar qomax a partir de los datos de la prueba
3) Determinar qo para diferentes valores de pwf desde pws hasta
cero.
4) Graficar los resultados pwf vs qo
Ejercicio1: Aplicar la ecuacin de Vogel para generar la curva de
afluencia para la siguiente prueba:

Pws= 2000 psi
Pwf= 1500 psi
qo= 100 bpd

Para los valores de Pwf: 2000, 1500, 1000, 500 y 0 psi
2
max
8 . 0 2 . 0 1
|
|
.
|
\
|
|
|
.
|
\
|
=
ws
wf
ws
wf
o
o
p
p
p
p
q
q
Mtodo de Afluencia
por IP
Mtodo de Afluencia por IP
ndice de Productividad:
) (
wf ws
o
p p
q
IP
=
Medida de la facilidad de flujo de los fluidos a travs del medio poroso, se define
como el gasto de produccin del pozo que puede obtener ante una cada de
presin en el yacimiento:
Este mtodo de Afluencia es el ms simple y se utiliza cuando no se cuenta con
pruebas de produccin en el pozo y solo se tiene conocimiento o idea del Indice
de productividad del pozo.
Forma de Lnea Recta:
IP
q
p p
o
ws wf
=
Pwf
qo
IP
m
1
=
Como se Utiliza este mtodo cuando existe gas en solucin,
presin de fondo fluyendo menor a la presin de burbujeo
(Pwf< Pb)?.

Se considera que el ndice de productividad representa un valor
caracterstico del pozo para una condicin cercana a la presin del
yacimiento es decir: p es pequeo, de tal manera que el valor del IP
representa el inverso de la pendiente de la curva de afluencia en dicho
punto.
Como se Utiliza este mtodo cuando existe gas en solucin,
presin de fondo fluyendo menor a la presin de burbujeo
(Pwf< Pb)?.
Pwf
qo
IP
m
1
=
Para obtener la variacin del gasto los diferentes valores de
pwf , lo cual es igual al IP, se recurre a la ecuacin de Vogel
para obtener:
wf
o
dp
dq
IP =
Partiendo de la ecuacin de Vogel:
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
max
8 . 0 2 . 0 1
ws
wf
ws
wf
o o
p
p
p
p
q q
(
|
|
.
|
\
|
=
ws ws
wf
ws
o
wf
o
p p
p
p
q
dp
dq 1
2 8 . 0
2 . 0
max
Bajo la consideracin antes planteada de que la cada de presin es muy
pequea
Por lo cual:
(
=
ws ws
o
wf
o
p p
q
dp
dq 6 . 1 2 . 0
max
1 ~
ws
wf
p
p
(
=
ws
o
wf
o
p
q
dp
dq 8 . 1
max
Dado que la pendiente o derivada es negativa el signo negativo de la
ecuacin resultante se elimina:
ws
o
wf
o
p
q
dp
dq
max
8 . 1
=
Dado que:
Finalmente el gasto mximo ser:
ws
o
p
q
IP
max
8 . 1
=
wf
o
dp
dq
IP =
8 . 1
max
ws
o
p IP
q
=
La ecuacin obtenida es la relacin entre el modelo de afluencia de IP
con el modelo de Vogel que considera el gas en solucin.
Se debe observar que la derivacin anterior se generaliza como:
Donde el gasto mximo del modelo de afluencia utilizado puede ser:

Darcy
IP
Pozo Horizontal
Pozo Fracturado
Fetkovich
Etc.
( )
( )
8 . 1
mod
max
max
afluencia de elo
o
Vogel
o
q
q =
Procedimiento:
1) Datos: Definir el valor de IP a utilizar y la Pws del yacimiento.

2) Determinar el gasto mximo.

3) Aplicar el mtodo de Vogel determinando qo para diferentes valores
de pwf desde pws hasta cero.

4) Graficar los resultados pwf vs qo.
Ejercicio 4:
Se realiz una toma de informacin en el pozo C-2 y como resultado de
la interpretacin de la prueba se obtuvieron los siguientes resultados.

Se desea obtener el ndice de Productividad del pozo y su curva de
Afluencia IPR.

k (mD) 84
rw (ft) 0.35
re (ft) 2500
h (ft) 100
pws (psi) 2356
vis (cp) @ pws 1.30
Bo @ pws 1.05
Dao (s) 5
Pb (psi) 2500
Ejercicio 4:
(
=
(
+
|
|
.
|
\
|
=
A
=
psi
bpd
s
r
r
B
kh
p
q
IP
w
e
32 . 3
4
3
ln 2 . 141
bpd
p IP
q
wS
o
4347
8 . 1
max
=
=
0
500
1000
1500
2000
2500
0 1000 2000 3000 4000 5000
P
w
f

(
p
s
i
)
Qo (bpd)
Curva de Afluencia
Ejercicio 4:
0
500
1000
1500
2000
2500
0 2000 4000 6000 8000
P
w
f

(
p
s
i
)
Qo (bpd)
Curva de Afluencia
IPR Darcy-Vogel
IPR Darcy
Mtodo de IPR
Generalizada
Mtodo de IPR Generalizado
Este mtodo es aplicable tanto para la condicin del yacimiento bajosaturado
como saturado.

Se usa para generar la curva IPR de pozos donde la presin esttica del
yacimiento pws se encuentra por arriba de la presin de burbujeo.

Es evidente que con el tiempo de explotacin del yacimiento, la pws disminuye y
eventualmente en algn tiempo la presin ser menor a la presin de burbujeo.

Bajo lo anterior la Curva IPR que se construya tendr dos partes:

a) IPR Parte Bajosaturada (pwf > pb)

b) IPR Parte Saturada (pwf < pb)

IPR Parte Bajosaturada

Fase: Aceite IPR Lineal

Ecuacin que define el comportamiento: Indice de Productividad

Rango de Presin:

) (
wf ws
o
p p
q
IP
=
) (
wf ws o
p p IP q =
ws wf b
p p p < <
IPR Parte Bajosaturada

Forma Grfica:

Pws
qo
Pb
qob
) (
b ws ob
p p IP q =
IPR Parte Saturada

Fase: Aceite+Gas IPR Curva

Ecuacin que define el comportamiento: Vogel

Rango de Presin:

(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
max
8 . 0 2 . 0 1
b
wf
b
wf
o o
p
p
p
p
q q
b wf
p p < < 0
IPR Parte Saturada

Forma Grfica:

Pws
qo
Pb
qob qomax
8 . 1
max
b
ob o
p IP
q q

+ =
IPR Generalizada:

Suma de ambas Partes: Bajosaturada y Saturada

Pws
qo
Pb
qob qomax
( )
(
(
|
|
.
|
\
|
|
|
.
|
\
|
+ =
2
max
8 . 0 2 . 0 1 ) (
b
wf
b
wf
ob o b ws o
p
p
p
p
q q p p IP q
Su aplicacin depende de Pws y Pwf
Para la aplicacin de la ecuacin generalizada existen tres diferentes casos:
a) Pws>Pb y Pwf>Pb
b) Pws>Pb y Pwf<Pb
c) Pws<Pb

Pws
qo
Pb
qob qomax
Pwf
Pws

Pws>Pb y Pwf>Pb

Pws
qo
Pb
qob qomax
Pwf
Pws
) (
wf ws o
p p IP q =

Pws>Pb y Pwf<Pb

Pws
qo
Pb
qob qomax
Pwf
Pws
( )
(
(
|
|
.
|
\
|
|
|
.
|
\
|
+ =
2
max
8 . 0 2 . 0 1 ) (
b
wf
b
wf
ob o b ws o
p
p
p
p
q q p p IP q

Pws<Pb

Pws
qo
Pb
qob qomax
Pwf
Pws
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
max
8 . 0 2 . 0 1
ws
wf
ws
wf
o o
p
p
p
p
q q
Ejercicio 5:
Se realiz una toma de informacin con ULA al pozo C-3007, obteniendo la
siguiente informacin:

Si se sabe que la presin de burbujeo del yacimiento es Pb (kg/cm2)= 152.69 ,
obtener la curva IPR del pozo.
pws (Kg/cm2) 261.44
pwf (Kg/cm2) 212.50
Qo (bpd) 8526
Ejercicio 5:
Pwf (Kg/cm2) Qo (bpd)
0.00 33720
10.89 33449
21.79 33057
32.68 32546
43.57 31914
54.47 31161
65.36 30288
76.25 29295
87.15 28182
98.04 26948
108.93 25594
119.83 24119
130.72 22524
141.61 20809
152.51 18974
163.40 17076
174.29 15179
185.19 13282
196.08 11384
206.97 9487
217.87 7589
228.76 5692
239.65 3795
250.55 1897
261.44 0
0
50
100
150
200
250
300
0 5000 10000 15000 20000 25000 30000 35000 40000
P
w
f

(
k
g
/
c
m
2
)
Qo (bpd)
IPR (Generalizada)
Mtodo de IPR
Composite
Mtodo de IPR Composite
El modelo de Pozo Composite es aquel en el cual se considera la produccin de
aceite de un pozo con una fraccin de agua producida.

Se conceptualiza como si fuera un pozo con dos formaciones productoras en
donde una aporta aceite y la otra agua. El manejo de este modelo considera
aspectos que deben tenerse claros para comprender dicho comportamiento.

Partiendo del esquema grfico siguiente:

C.A.A
Aceite
Agua
Qo Qo
Qw Qw
Se debe tener en cuenta que este modelo considera el aporte de las fases dentro
del intervalo disparado, es decir, que el contacto agua-aceite ya est presente y
no est planteado para canalizaciones o conificaciones de agua.

El gasto de lquidos ser la suma de la produccin de aceite y agua:

Considerando la premisa de que la produccin de ambas fases es de la misma
formacin, el ndice de productividad es el mismo para calcular el flujo de aceite
y de agua. Adems considerando un yacimiento saturado.

El clculo de aceite lo obtenemos mediante la ecuacin de Vogel.

w o L
q q q + =
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
max
8 . 0 2 . 0 1
ws
wf
ws
wf
o o
p
p
p
p
q q
El gasto mximo se puede obtener de:

Sustituyendo en la ecuacin de Vogel

Para calcular la produccin de agua se utiliza la definicin de ndice de
productividad:

8 . 1
max
ws
o
p IP
q

=
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
8 . 0 2 . 0 1
8 . 1
ws
wf
ws
wf
ws
o
p
p
p
p
p IP
q
( )
wf ws w
P P IP q =
Dado que existe produccin de agua y aceite, los gastos a calcular debern ser
afectados por su correspondiente fraccin:

Agua:

Aceite:

El gasto total ser:

( )
wf ws w w
P P f IP q =
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
8 . 0 2 . 0 1
8 . 1
ws
wf
ws
wf
o
ws
o
p
p
p
p
f
p IP
q
( )
wf ws w
ws
wf
ws
wf
o
ws
L
P P f IP
p
p
p
p
f
p IP
q +
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
8 . 0 2 . 0 1
8 . 1
Despejando el ndice de productividad; que es el objetivo del modelo de pozo
composite:

De esta manera el IP puede ser calculado partiendo de una prueba de produccin
que contiene:

- Pws
- Pwf
- qL
- fw

( )
wf ws w
ws
wf
ws
wf
o
ws
L
P P f
P
P
P
P
f
P
q
IP
+
(
(
|
|
.
|
\
|
|
|
.
|
\
|
|
.
|
\
|
=
2
8 . 0 2 . 0 1
8 . 1
Se debe tener en cuenta que esta ecuacin aplica solo para el punto de la
prueba.

Para la generacin de la Curva IPR, se utiliza el mtodo planteado por
PETROBRAS, en donde se considera un comportamiento como el mostrado en la
siguiente grfica:

0.00
20.00
40.00
60.00
80.00
100.00
120.00
0 100000 200000 300000 400000 500000 600000
P
w
f
(
K
g
/
c
m
^
2
)
Ql (bpd)
IPR Composite
Qo (bpd) Qw (bpd) IPR TOTAL
El anlisis de la curva IPR se divide en tres secciones:

Parte bajosaturada

Parte donde

Parte donde

Parte bajosaturada

Para la parte bajosaturada, como no existe liberacin del gas en solucin el
comportamiento de la curva IPR es una lnea recta:

) (
wf ws l
p p IP q =
max o l b
q q q < <
max o l
q q >

Parte donde

En esta parte se debe tener en cuenta que la referencia de qomax es el potencial
mximo del pozo para su produccin de aceite sin considerar la fraccin de agua,
es decir para una fraccin de aceite igual a uno.

Se utiliza la ecuacin obtenida anteriormente para el ndice de productividad

Y la ecuacin para determinar la presin de fondo fluyendo para cada valor de
gasto de lquido
max o l b
q q q < <
( )
wf ws w
ws
wf
ws
wf
o
ws
L
P P f
P
P
P
P
f
P
q
IP
+
(
(
|
|
.
|
\
|
|
|
.
|
\
|
|
.
|
\
|
=
2
8 . 0 2 . 0 1
8 . 1

Cabe aclarar que esta ecuacin no es el valor de pwf despejado de la ecuacin que
define el ndice de productividad, sino que se obtiene al igualar la produccin de
aceite a la del lquido para determinar pwf al aceite y de igual manera se iguala la
produccin de agua a la de lquido para obtener pwf al agua.

|
.
|
\
|
+
|
|
.
|
\
|
|
|
.
|
\
|
=
IP
q
p f
q
q
f p p
l
ws w
o
l
o ws wf
1 80 81 125 . 0
max
IPR Composite
0.00
20.00
40.00
60.00
80.00
100.00
120.00
0 100000 200000 300000 400000 500000 600000
Ql (bpd)
P
w
f
(
K
g
/c
m
^
2
)
Qo (bpd) Qw(bpd) IPRTOTAL
De tal manera que la generacin de la IPR hasta esta etapa quedara de la
siguiente manera:

A partir de este punto, significa que ya no hay capacidad de aportacin de aceite,
por lo cual el resto de la aportacin es de agua con un comportamiento lineal,
pero con un nuevo valor de ndice de productividad (IP), el cual es calculado
como se muestra en la siguiente parte.
IPR Composite
0.00
20.00
40.00
60.00
80.00
100.00
120.00
0 100000 200000 300000 400000 500000 600000
Ql (bpd)
P
w
f

(
K
g
/
c
m
^
2
)
Qo (bpd) Qw(bpd) IPR TOTAL

Parte donde

Como ya se mencion anteriormente, cuando , ya no puede existir
aportacin de aceite, por lo que el resto de la aportacin es de agua, con un
comportamiento lineal, pero con un nuevo valor de ndice de productividad (IP),
dicho ndice es calculado con el valor de la pendiente de la curva en una seccin
final donde , tomando este punto y uno muy cercano que se establece
como

max o l
q q >
max o l
q q >
max o l
q q =
max
999 . 0
o l
q q =
IPR Composite
0.00
20.00
40.00
60.00
80.00
100.00
120.00
0 100000 200000 300000 400000 500000 600000
Ql (bpd)
P
w
f

(
K
g
/
c
m
^
2
)
qo max 0.999*qo
max
De la definicin de ndice de Productividad:

Una vez obtenida la expresin anterior, se puede definir que la presin de fondo
fluyendo , para los valores de gastos , se calculan con la ecuacin
resultante de evaluar para y restando el valor de correspondiente para
el valor de , utilizando el nuevo valor de .

p
q
IP
o
A
=
max
001 . 0
'
|
.
|
\
|
+
|
|
.
|
\
|
|
|
.
|
\
|

=
IP
q
f
q
q
f p
q
IP
o
w
o
o
o ws
o
max
max
max
max
001 . 0
1
999 . 0
80 81 125 . 0
001 . 0
'
max o l
q q >
wf
p
max o
q
max o l
q q
' IP
Si

restando el valor de correspondiente para el valor de , utilizando el
nuevo valor de .

max o l
q q =
|
.
|
\
|
+
|
|
.
|
\
|
|
|
.
|
\
|
=
IP
q
p f
q
q
f p p
o
ws w
o
o
o ws wf
max
max
max
1 80 81 125 . 0
|
.
|
\
|
=
IP
q
p f p
o
ws w wf
max
' IP
max o l
q q
'
max max
IP
q q
IP
q
p f p
o l o
ws w wf

|
.
|
\
|
=
Ahora la grfica de IPR completa queda de la siguiente forma

IPR Composite
0.00
20.00
40.00
60.00
80.00
100.00
120.00
0 100000 200000 300000 400000 500000 600000
Ql (bpd)
P
w
f

(
K
g
/
c
m
^
2
)

Para evaluar el gasto mximo de lquido del pozo se utiliza la ltima ecuacin
obtenida:

El gasto mximo de lquido qLmax se obtendr cuando pwf=0 ; por lo tanto:

'
max max
IP
q q
IP
q
p f p
o l o
ws w wf

|
.
|
\
|
=
'
0
max max max
IP
q q
IP
q
p f
o l o
ws w

|
.
|
\
|
=
|
.
|
\
|
+ =
IP
q
p f IP q q
o
ws w o l
max
max max
'
Procedimiento:
1) Datos: Pws, Pwf ,qL ,fw.

2) Determinar ndice de Productividad:

3) Determinar el gasto mximo de aceite qomax para 100% aceite.

4) Calcular ndice de Productividad prima

( )
wf ws w
ws
wf
ws
wf
o
ws
L
P P f
P
P
P
P
f
P
q
IP
+
(
(
|
|
.
|
\
|
|
|
.
|
\
|
|
.
|
\
|
=
2
8 . 0 2 . 0 1
8 . 1
|
.
|
\
|
+
|
|
.
|
\
|
|
|
.
|
\
|

=
IP
q
f
q
q
f p
q
IP
o
w
o
o
o ws
o
max
max
max
max
001 . 0
1
999 . 0
80 81 125 . 0
001 . 0
'
8 . 1
max
ws
o
p IP
q

=
Procedimiento:
5) Determinar el gasto mximo de lquido qLmax .

6) Establecer valores del gasto de lquido para el rango de valores de
cero hasta qLmax .

7) Determinar Pwf para cada valor de gasto de lquido establecido.

|
.
|
\
|
+ =
IP
q
p f IP q q
o
ws w o l
max
max max
'
|
.
|
\
|
+
|
|
.
|
\
|
|
|
.
|
\
|
=
IP
q
p f
q
q
f p p
l
ws w
o
l
o ws wf
1 80 81 125 . 0
max '
max max
IP
q q
IP
q
p f p
o l o
ws w wf

|
.
|
\
|
=
Ejercicio 6:
Se realiz una toma de informacin en un pozo con alta produccin de agua, los
resultados se muestran en la tabla.

Si se sabe que la presin de burbujeo del yacimiento es Pb (kg/cm2)= 152.69 ,
obtener la curva IPR del pozo.
pws (kg/cm2) 113.12
pwf (kg/cm2) 112.73
ql (bpd) 2965
fw (fraccin) 0.8
Ejercicio 6:
fo (fraccin) 0.2
IP (bpd/psi) 534.69
qomax (bpd) 477925
IP' (bpd/psi) 208.42
qlmax (bpd) 597156
0
20
40
60
80
100
120
0 100000 200000 300000 400000 500000 600000 700000
P
w
f

(
K
g
/
c
m
^
2
)
Ql (bpd)
IPR Composite
Ejercicio 6:
0
20
40
60
80
100
120
0 100000 200000 300000 400000 500000 600000 700000
P
w
f

(
K
g
/
c
m
^
2
)
Ql (bpd)
IPR Composite
80% Agua
0% Agua
Mtodo de IPR Jones
Mtodo de IPR de Jones
Este modelo de comportamiento de afluencia considera el parmetro o la
variable del intervalo disparado dentro del intervalo productor, por lo que el flujo
radial hacia el pozo puede llegar a convertirse en flujo esfrico dependiendo de la
relacin . Basado en lo anterior este modelo considera 2 trminos dentro
de la ecuacin, uno para el flujo laminar del fluido (flujo Darcyano) y otro para
considerar el flujo turbulento (flujo no-Darcyano).

Donde es el factor de tortuosidad;

) / ( h h
p
( ) | |
w p
o o o w e o o o
wf ws
r h
q B x
kh
s r r B q
P P
2
2 2
14
10 3 . 2 75 . 0 ln 2 . 141 |
+
+
=
|
201 . 1
10
10 33 . 2
k
x
= |
Si en la ecuacin anterior se considera la siguiente agrupacin:

Trmino de flujo laminar

Trmino de flujo turbulento

Entonces;

Ecuacin cuadrtica.
( ) | |
kh
s r r B
A
w e o o
+
=
75 . 0 ln 2 . 141
w p
o o
r h
B x
B
2
2
14
10 3 . 2 |
=
( ) 0
2
= +
wf ws o o
P P Aq Bq
B
p p B A A
qo
wf ws
2
) ( 4
2
+
=
Con la solucin de la ecuacin cuadrtica se pueden obtener los gastos a
diferentes presiones de fondo fluyendo y construir la curva IPR.

Sin embargo una de las mayores aplicaciones del mtodo de Jones es evaluar la
longitud ptima de disparos de un pozo, a fin de obtener la mayor produccin
con la menor longitud de disparos.

Evaluando lo anterior en trminos del qomax para pwf=0
B
Bp A A
q
ws
o
2
4
2
max
+ +
=
Para determinar la longitud ptima de disparos se genera una tabla como la
siguiente:

B
Bp A A
q
ws
o
2
4
2
max
+ +
=
hp (ft) B (psi/bpd^2) qomax (bpd)
1 B1 qomax1
3 B2 qomax2
5 B3 qomax3
10 B4 qomax4

h Bh qomaxh
w p
o o
r h
B x
B
2
2
14
10 3 . 2 |
=
( ) | |
kh
s r r B
A
w e o o
+
=
75 . 0 ln 2 . 141
Al graficar los resultados se obtiene:

hp
qomax
hp ptimo
Ejercicio 8:
Se encuentra interviniendo un pozo con una plataforma semisumergible, debido
a los altos costos de renta de la plataforma se desea optimizar el nmero de
corridas de disparos en el pozo sin afectar la produccin del mismo.

Se ha solicitado al grupo de Productividad de Pozos definir la longitud ptima de
disparos.
La informacin disponible de la formacin y fluidos es la siguiente:
Datos
k (mD) 422
rw (ft) 0.35
re (ft) 715
h (ft) 656
s 10
pws (psi) 1648
vis (cp) @ pws 2.77795
Bo @ pws 1.25811
o (lb/ft3) 49.58
Ejercicio 8:
hp (ft) B (psi/bpd^2) Qomax(bpd)
3
9.3E-06 11807
5
3.3E-06 18167
10
8.3E-07 29941
15
3.7E-07 37492
20
2.1E-07 42377
30
9.3E-08 47798
50
3.3E-08 51863
100
8.3E-09 54041
150
3.7E-09 54486
200
2.1E-09 54645
300
9.3E-10 54760
400
5.2E-10 54801
500
3.3E-10 54819
600
2.3E-10 54830
0
10000
20000
30000
40000
50000
60000
0 100 200 300 400 500 600 700
q
o
m
a
x
(
b
p
d
)
hp (ft)
Sensibilidad hp
hp ptimo = 100 ft
Mtodo de IPR de Pozos
Horizontales
Mtodo de IPR de Pozos Horizontales (Estacionario)
Una ecuacin de Flujo en estado estacionario para determinar el
comportamiento de pozos horizontales que fue determinada por Joshi (1991),
bajo la definicin de un rea de drene elptica del pozo la cual esta definida por la
longitud de la seccin horizontal y el radio de drene.
La ecuacin Obtenida por Joshi es la siguiente:

( )
(
(
(
(
(
|
|
.
|
\
|
|
.
|
\
|
+
(
(
(
(
(
|
.
|
\
|
|
.
|
\
|
+
=
w
o o
wf ws h
o
r
h
L
h
L
L
a a
B
p p h k
q
2
ln
2
2
ln 2 . 141
2
2
| |

Donde:

(Geometra del Pozo)

y

(Anisotropa)

Unidades

4
2
4
1
2
1
2
|
.
|
\
|
+ +
|
.
|
\
|
=
L
r L
a
e
v
h
k
k
= |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| | Drene de Radio ft r
pozo del Radio ft r
pozo del horizontal Longitud ft L
aceite del idad Vis cp
aceite del volumen de Factor stbl bl B
fluyendo fondo de esin psi p
esttica fondo de esin psi p
formacin la de Espesor ft h
Vertical dad Permeabili mD k
Horizontal dad Permeabili mD k
produccin de Gasto bpd Qo
e
w
o
o
wf
ws
v
h
;
;
;
cos ;
; /
Pr ;
Pr ;
;
;
;
;

Como puede observarse es de forma similar a la ecuacin de Darcy para pozos
verticales, la ecuacin anterior considera que no hay dao en la vecindad de la
formacin. La forma del rea de drene y la ubicacin del pozo en el mismo
tendr que corresponder al siguiente diagrama:

Area de Drene Elptic a - J os hi
-400
-300
-200
-100
0
100
200
300
400
-1500 -1000 -500 0 500 1000 1500
L
a=re
b
Al considerar el dao en la Ecuacin de Joshi

( )
(
(
(
(
(
|
|
|
|
|
.
|
\
|
|
|
.
|
\
|
+
(
(
(
(
(
|
|
.
|
\
|
|
.
|
\
|
+
(
(
(
(
(
|
.
|
\
|
|
.
|
\
|
+
=
w
e
w
o o
wf ws h
o
r
r
s
r
h
L
h
L
L
a a
B
p p h k
q
ln
1
2
ln
2
2
ln 2 . 141
2
2
| |

Evaluacin del ndice de Productividad del Pozo:

Evaluacin del gasto mximo del Pozo:

(
(
(
(
(
|
|
|
|
|
.
|
\
|
|
|
.
|
\
|
+
(
(
(
(
(
|
|
.
|
\
|
|
.
|
\
|
+
(
(
(
(
(
|
.
|
\
|
|
.
|
\
|
+
=
w
e
w
o o
h
r
r
s
r
h
L
h
L
L
a a
B
h k
IP
ln
1
2
ln
2
2
ln 2 . 141
2
2
| |
8 . 1
max
ws
p IP
q

=
Caractersticas del Mtodo:

- Aplicable para flujo Estacionario.

- Evaluacin del comportamiento del pozo en el centro del rea de flujo.

- No permite realizar evaluacin del efecto de la posicin del pozo dentro del
rea de drene.

- Se debe garantizar un rea de drene elptica (No considera efectos de
frontera)
Procedimiento:
1) Datos: L, re, kh , kv , h, rw ,vis, Bo, Pws , s, b

2) Determinar a y

3) Determinar el ndice de Productividad del pozo
4
2
4
1
2
1
2
|
.
|
\
|
+ +
|
.
|
\
|
=
L
r L
a
e
v
h
k
k
= |
(
(
(
(
(
|
|
|
|
|
.
|
\
|
|
|
.
|
\
|
+
(
(
(
(
(
|
|
.
|
\
|
|
.
|
\
|
+
(
(
(
(
(
|
.
|
\
|
|
.
|
\
|
+
=
w
e
w
o o
h
r
r
s
r
h
L
h
L
L
a a
B
h k
IP
ln
1
2
ln
2
2
ln 2 . 141
2
2
| |

Procedimiento:
4) Determinar qomax :

5) Establecer rango de pwf y determinar sus gastos correspondientes
con Vogel.

6) Construir curva IPR.

8 . 1
max
ws
p IP
q
=
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
max
8 . 0 2 . 0 1
ws
wf
ws
wf
o o
p
p
p
p
q q
Ejercicio 11:
Determinar la Curva IPR para el pozo Horizontal cuya informacin se muestra en
la tabla inferior. Se sabe que el yacimiento presenta una entrada de agua
importante, manteniendo la presin del yacimiento.
Datos
L (ft) 2625
re (ft) 984
kh (mD) 1100
kv/kh (adim) 0.1
h (ft) 394
rw (ft) 0.354
vis(cp) 2.7
Bo (bls/stbl) 1.16
pws (psi) 1636
S 5
Longitud Zona (ft) 8202
Ejercicio 11:
Clculos
kv (mD) 110
3.16
a 1469
IP (bpd/psi) 149

Qomax(bpd) 135499
0
200
400
600
800
1000
1200
1400
1600
1800
0 20000 40000 60000 80000 100000 120000 140000 160000
P
w
f

(
p
s
i
)
Qo (bpd)
IPR Pozo Horizontal Estado Estacionario - Joshi
Aspectos a Evaluar en
Pozos Horizontales
Mtodo de IPR de Pozos Horizontales (Pseudo-estacionario)
- Incremento del ndice de Productividad

Flujo Estacionario:

Flujo Pseudo-estacionario:

(
(
(
|
.
|
\
|
+ + +
|
|
|
.
|
\
|
=
s
L
b
s C
r
A
B
k k b
IP
R H
w
o o
z x
75 . 0 ln ln 2 . 141
2
1
(
(
(
(
(
|
|
|
|
|
.
|
\
|
|
|
.
|
\
|
+
(
(
(
(
(
|
|
.
|
\
|
|
.
|
\
|
+
(
(
(
(
(
|
.
|
\
|
|
.
|
\
|
+
=
w
e
w
o o
r
r
s
r
h
L
h
L
L
a a
B
kh
IP
ln
1
2
ln
2
2
ln 2 . 141
2
2
| |

- Evaluar prdidas de Presin por friccin en la seccin horizontal

Perfil de Aportacin Unitaria
0
100
200
300
400
500
600
700
0 100 200 300 400 500 600 700 800 900
Seccin Horizontal de Taln a Punta (Izquierda- Derecha) (mts)
Q
n

(
b
p
d
/
m
)
Dimetro 6.625
0
100
200
300
400
500
600
700
0 100 200 300 400 500 600 700 800 900
Q
n

(
b
p
d
/
m
)
Dimetro 6.184
0
100
200
300
400
500
600
700
0 100 200 300 400 500 600 700 800 900
Q
n

(
b
p
d
/
m
)
Dimetro 4.892
0
100
200
300
400
500
600
700
0 100 200 300 400 500 600 700 800 900
Q
n

(
b
p
d
/
m
)
Dimetro 3.958
- Determinar longitud horizontal ptima

1000
1200
1400
1600
1800
2000
2200
2400
2600
2800
3000
0 100 200 300 400 500 600 700 800 900 1000
Q
l

(
b
p
d
)
Longitud Horizontal (mts)
Comportamiento de Horizontalidad
- Regular y Homogeneizar el flujo a lo largo de la seccin horizontal (ICD)

200 md
T.R. 30
800 md T.R. 13 3/8
6 1/8 P.T @ 3482.83 md
2890 mV
Liner 7 3252.74 md, 2890 mv
Estado Mecnico
Propuesto
sssv
A
P
A
R
E
J
O

D
E

P
R
O
D
U
C
C
I
O
N

5

V.T. 5 1/2 150 M.

Liner 4 con sistema
Equalizer
Colgador 5 x 4 2900 md
F
r
a
c
t
u
r
a
m
i
e
n
t
o

3
3
1
4

3
3
8
7

m
d
E
m
p
a
c
a
d
o
r

3
3
1
0

m
d
E
m
p
a
c
a
d
o
r

3
3
9
5

m
d
E
m
p
a
c
a
d
o
r

3
4
8
2

m
d
IPR Futuro
IPR Futuro Mtodo de Eickmer
Este mtodo establece una ecuacin cbica que relaciona los gastos mximos
presentes y futuros en funcin del cambio de la presin del yacimiento.

Principal consideracin:

3
max max
|
|
.
|
\
|
=
actual ws
futura ws
actual futuro
p
p
q q
cte p = A
Aplicable para evaluar el Comportamiento IPR de un pozo en funcin de las
propiedades de un pozo correlacin.
Seleccin del Pozo Vecino cercano,
produciendo en la misma formacin y con
informacin reciente.
cte p = A
Procedimiento:
1) Datos originales del Pozo: qo, Pws , Pwf , Pws futura

2) Determinar qomax actual

3) Determinar qomax futuro
2 max
8 . 0 2 . 0 1
|
|
.
|
\
|
|
|
.
|
\
|
=
ws
wf
ws
wf
o
actual o
p
p
p
p
q
q
3
max max
|
|
.
|
\
|
=
actual ws
futura ws
actual o futuro o
p
p
q q
Procedimiento:
4) Determinar los gastos para diferentes Pwf , a partir de Vogel

5) Construir curva IPR
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
max
8 . 0 2 . 0 1
wsfutura
wf
wsfutura
wf
futuro o o
p
p
p
p
q q
Ejercicio 14:
Determinar el IP y la Curva IPR para un pozo en funcin de la informacin de un
pozo vecino cercano que servir para correlacionar la informacin de afluencia.
Si la presin estimada actual es de Pws = 113 kg/cm2.

Clculo de IPR Futuro
DATOS DEL POZO DE CORRELACION
Pozo Correlacin
C-XXXX
Pws 121.60 (kg/cm)
Pwf 117.02 (kg/cm)
Qo 7367 (BPD)
Ejercicio 14:
CONDICIONES FUTURAS
Pwsf 113.00 (kg/cm)
Pwf 108.42 (kg/cm)
Qmaxf
(BPD)
88685 (bpd)
IP 99.32 (BPD/psia)
Qof 6354 (BPD)
0
20
40
60
80
100
120
140
0 20000 40000 60000 80000 100000 120000
P
r
e
s
i
n

(
K
g
/
c
m
2
)
Qo (BPD)
Curva de IPR
IPR Original
IPR Futura
IPR Futuro Otros Mtodos
Fetkovich

Requiere prueba Isocronal

Standing

Requiere conocer el comportamiento de
Permeabilidades relativas con la presin

(
(
(
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
p
o o
ro
f
o o
ro
p f
B
k
B
k
IP IP
( )
n
wf ws
ws
ws
o o
P P
P
P
J q
2 2
1
1
2
1 2
'
|
|
.
|
\
|
=
(
(
|
|
.
|
\
|
|
|
.
|
\
|
=
2
8 . 0 2 . 0 1
8 . 1
ws
wf
ws
wf ws f
o
p
p
p
p p IP
q
III.- FACTOR DE DAO Y SU RELACIN CON
COMPORTAMIENTO DE AFLUENCIA
Definicin de Dao
Es una restriccin adicional al flujo que se manifiesta como una cada de presin
adicional en el trayecto del fluido del yacimiento al pozo.
Esta restriccin se da en la vecindad del pozo principalmente ya que es donde se
encuentran los fluidos de invasin.
Dao a la Formacin
Factores que provocan dao:
a) Invasin de Fluidos
b) Disparos
c) Penetracin Parcial
d) Desviacin
e) Pseudo-daos
I. Fracturamiento
II. Empacamiento de grava
III. Liner
IV. Cedazos
V. Etc.

Dao Van Everdingen y Hurst
Invasin de Fluidos:

Originado por el filtrado de los fluidos de perforacin y terminacin de los
pozos, generalmente este tipo de dao es el mas comn y el de mayor impacto
sobre la produccin. Dicho dao puede ser removido mediante tratamientos
como Estimulaciones y Limpiezas a la formacin.
Su radio de alcance es de algunos pies dependiendo de la permeabilidad de la
formacin o si existen fracturas.

Definicin por Van Everdingen y Hurst (1953). Pelcula Infinitesimal

s
p
qB
kh
s A =
2 . 141
Dao a la Formacin
Invasin de Fluidos:

s
p A
Dao a la Formacin Hawkins
Invasin de Fluidos:

Definicin de Dao por Hawkins

k k s
p p p
s
A A = A
rw
rs,,kS
re,k
0
0
0
< >
= =
> <
s k k
s k k
s k k
s
s
s
kh
r
r
qB
h k
r
r
qB
kh
s qB
w
s
s
w
s
|
|
.
|
\
|
|
|
.
|
\
|
=
ln 2 . 141 ln 2 . 141
2 . 141

|
|
.
|
\
|
=
k k h
r
r
qB
kh
s qB
s
w
s
1 1
ln 2 . 141
2 . 141

|
|
.
|
\
|
= 1 ln
s w
s
k
k
r
r
s
Dao por Disparos
Dao por Disparos:

Originado por mala eficiencia al disparar el Liner de Produccin para
comunicar el pozo con el yacimiento.

Principal objetivo de un Disparo?

Obtener una produccin del pozo igual o mayor a la produccin que tendra el
pozo si produjera en agujero descubierto. Lo anterior se logra al incrementar el
radio efectivo del pozo con el disparo, sin disparar todo el intervalo productor.

Factores que influyen:

a) Dimetro de Disparos
b) Penetracin de Disparos
c) Densidad de Disparos
d) Fase
e) Dimetro del pozo
f) Anisotropia (kz/kr)

w
r
p
l
p
h
u
p
r
Dao por Disparos

Formas de Evaluar el dao por disparos:

a) Nomograma; Hong (1975)

b) Analtica; Karakas y Tariq (1991)

Dao por Disparos
Nomograma; Hong (1975)

Dao por Disparos
Evaluar Dao por Disparos:

Densidad Disparos: 2 cargas/ft
Dimetro del Pozo: 6
Kz/Kr=1
Fase: 180
Penetracin: 4
Dimetro de Disparos:

Sp= 0.3
Dao por Disparos
Analtica; Karakas y Tariq (1991)

La determinacin del Dao por disparos fue desarrollado por Karakas y Tariq
en 1991, presentando una solucin semianaltica, dividiendo el efecto de dao
por disparos en 3 componentes:

Efecto de flujo en el plano horizontal ( )

Efecto de convergencia Vertical ( )

Efecto del agujero ( )

El dao total debido a los disparos ser:

H
s
v
s
wb
s
wb v H p
s s s s + + =
Dao por Disparos

Las variables requeridas para el clculo de los componentes de dao por
disparos son:

Radio del pozo ( )

Radio de los disparos ( )

Penetracin de los disparos ( )

Fase de las pistolas ( )

Densidad de cargas ( )

w
r
p
r
p
l
u
p
h / 1
Dao por Disparos
Clculo del Efecto de Flujo en el Plano Horizontal ( )

Lo cual significa que el efecto en el plano horizontal es una relacin del radio
de flujo que tendra el pozo en agujero descubierto, con respecto a un radio
equivalente generado en funcin de la fase de las pistolas. Esta definicin es
similar a la ecuacin de dao por invasin de fluidos a partir de radios
equivalentes.

H
s
( )
|
|
.
|
\
|
=
u
ln
w
w
H
r
r
s
Dao por Disparos
Para evaluar el radio equivalente en funcin de la fase de las pistolas se utilizan
las siguientes ecuaciones, determinadas por Karakas y Tariq:

El parmetro se determina a partir de la Tabla siguiente:

( )
( )
= +
=
=
0
0
4
u
u
u
u
para l r a
para
l
r
p w
p
w
u
a
fase
360 0.25 -2.091 0.0453 5.1313 1.8672 0.16 2.675
180 0.5 -2.025 0.0943 3.0373 1.8115 0.026 4.532
120 0.648 -2.018 0.0634 1.6136 1.777 0.066 5.32
90 0.726 -1.905 0.1038 1.5674 1.6935 0.0019 6.155
60 0.813 -1.898 0.1023 1.3654 1.649 0.0003 7.509
45 0.86 -1.788 0.2398 1.1915 1.6392 0.000046 8.791
u
a
1
a
2
a
1
b
2
b
1
c
2
c
Dao por Disparos
Al observar el comportamiento de la tabla, se determina que conforme
disminuye el ngulo de la fase, el factor aumenta, lo cual significa un radio
equivalente mayor y por consecuencia un dao por efecto de flujo horizontal
menor.

Clculo del Efecto de Convergencia Vertical ( )

Para el clculo de este efecto se determinan 2 variables adimensionales:

Donde son la permeabilidad horizontal y vertical respectivamente.

Dao por Disparos
v
s
v
h
p
p
D
k
k
l
h
h =
|
|
.
|
\
|
+ =
h
v
p
p
D
k
k
h
r
r 1
2
v h
k k ,
; es una relacin de la separacin entre las cargas, verticalmente, con
respecto a la penetracin de los disparos.

; es una relacin del radio de los disparos, con respecto a la separacin
entre las cargas.

El efecto de convergencia vertical se determina de la siguiente manera:

De donde a y b, se determinan con las siguientes ecuaciones:
Dao por Disparos
D
h
D
r
b
D
b
D
a
v
r h s
1
10
=
( )
2 1
log a r a a
D
+ =
2 1
b r b b
D
+ =
Donde son obtenidos de la tabla anterior.

El efecto de dao vertical, , generalmente es el de mayor contribucin,
principalmente para disparos con densidades de cargas pequeas.

Clculo del Efecto del Agujero ( )

Para determinar este efecto se requiere calcular una variable adimensional
, la cual es una relacin del radio del pozo con respecto a un radio equivalente
resultante de la penetracin de los disparos en la formacin, considerando una
densidad de cargas muy grande y un ngulo pequeo de fase. Siendo este el
efecto por el agujero.

Dao por Disparos
2 1 2 1
, , b y b a a
v
s
wb
s
wD
r
w p
w
wD
r l
r
r
+
=
El dao por el efecto del Agujero es:

Las constantes y son obtenidas de la tabla anterior.

Dao por Disparos
wD
r c
wb
e c s
2
1
=
1
c
2
c
Ejercicio 3.1

Determinar el Dao por Disparos a travs del Mtodo de Karakas y Tariq, para
un pozo con la siguiente informacin:

Dao por Disparos
lp (plg) Penetracin de los Disparos 8
dp (plg) Dimetro del Agujero de los disparos 0.25
rw (plg) Radio del Pozo 3.936
spf (cargas/pie) Cargas Explosivas por Pie 2
Fase ( ) Fase de las cargas 180
Kh (mD) Permeabilidad radial (a partir de la curva de
Incremento) 8000
Kz (mD) Permeabilidad vertical (Generalmente Kh/10) 800
Ejercicio 3.1

Resultados

Dao por Disparos
Dao por Disparos debido al efecto Radial (Sh) -0.416
Dao por Disparos debido al efecto Vertical (Sv) 5.426
Dao por Disparos debido al efecto del Agujero (Swb) 0.116

Suma del Dao por disparos 5.125
Dao por Penetracin Parcial y Desviacin
El Dr. Heber Cinco Ley (1975) determin una solucin semianaltica, en donde
present tablas de estos daos para diferentes combinaciones de penetracin
parcial y desviaciones, bajo la siguiente configuracin de Pozo-Yacimiento.

w
r
u
w
z
w
h
h

Para el clculo del dao por desviacin y Penetracin parcial, se requieren
evaluar los siguientes parmetros Adimensionales:

Con las variables anteriores, se entra a las siguientes tablas y se determinan los
valores de dao por desviacin y dao por penetracin parcial .

d c d c
s s s + =
+
w
D
r
h
h =
w
w
wD
r
z
z =
w
w
wD
r
h
h =
( )
d
s ( )
c
s
Existe una relacin Lineal entre Sc+d vs log(hD)

Ejercicio 3.2

Determinar el Dao por Penetracin Parcial y Desviacin para un pozo con la
siguiente configuracin, rw=3.936:

Desviacin ( ) Desviacin del Pozo 45
Cima (mV) Profundidad vertical de la
cima de la formacin 0
Base (mV) Profundidad vertical de la
base de la formacin 100
Cima disparos (mV) Profundidad vertical de la
cima de los disparos 0
Base disparos (mV) Profundidad vertical de la
base de los disparos 10
base disparos
cima disparos
cima
zw
h hw
base
Ejercicio 3.2

hD=1000
ZwD= 950
hwD= 101

S c+d = 41.521 -10.677 = 30.844

Ecuacin General de Dao
Ecuacin General de Dao:

ptotal = pradial +pc+d +pf + pp
rw
Flujo Radial
Penetracin Parcial
Invasin Disparos
Ecuacin General de Dao
ptotal = pradial +pc+d +pf + pp

El dao que puede ser removido es el de invasin de fluidos, mediante
Estimulacin o Limpieza de Intervalo productor:

p f d c total dao
p p p p A + A + A = A
+
p
w
f
w
d c t
s
kh
qB
s
kh
qB
s
kh
qB
s
kh
qB 2 . 141 2 . 141 2 . 141 2 . 141
+ + =
+
p
w
f
w
d c t
s
h
h
s
h
h
s s + + =
+
( )
p d c t
w
f
s s s
h
h
s =
+
Dao a la Formacin
Inclusin el Dao en la Ecuacin de Darcy:

( )
(
|
|
.
|
\
|
=
4
3
ln 2 . 141
w
e
wf
r
r
B
p p kh
q
( )
(
+
|
|
.
|
\
|
=
s
r
r
B
p p kh
q
w
e
wf
4
3
ln 2 . 141
Determinacin del Dao total
Este dao es calculado a partir del anlisis de una prueba de presin, el cual
generalmente se obtiene de la curva de incremento, por garantizar un gasto
constante en el pozo (qo=0).
Partiendo de la toma de informacin del pozo, se define la curva de
incremento que se va a analizar.

La tcnica de anlisis de curvas de incremento ms usual es la tcnica de
Horner, la cual consiste en graficar los datos de presin contra el logaritmo del
tiempo de Horner, el cual se obtiene a partir de la aplicacin de principio de
superposicin.

Es decir se grafica . Donde es el tiempo de Horner.

Del resultado de la grfica semilog, se obtiene una seccin de lnea recta, la
cual corresponde al comportamiento del yacimiento, como un yacimiento
infinito Homogneo.

|
|
.
|
\
|
A
A +
t
t t
vs p
p
ws
log
|
|
.
|
\
|
A
A +
t
t t
p
La ecuacin de define el comportamiento del yacimiento Homogneo infinito
para una curva de incremento es la siguiente:

De donde podemos observar que es la forma de una ecuacin de lnea recta, en
donde la pendiente es:

De donde puede ser despejada la permeabilidad de la formacin:

|
|
.
|
\
|
A
A +
+ =
t
t t
kh
qB
p p
p
wf ws
log
6 . 162
kh
qB
m
6 . 162
=
mh
qB
k
6 . 162
=
Partiendo de la aproximacin logartmica de la solucin lnea fuente:

El dao puede ser despejado:

( )
(
(
+
|
|
.
|
\
|
+ A + = s
cr
k
t
h k
qB
p p
w
i ws
8686 . 0 2274 . 3 log log
6 . 162
2
|
(
(
+
|
|
.
|
\
|
= A
= 2274 . 3 log
) 0 (
1513 . 1
2
1
w
wf hr
cr
k
m
t p p
s
|
Ejercicio 3.3
Determinar el Dao del Pozo a partir de la Siguiente Informacin
tiempo pws (psi) tiempo pws (psi) tiempo pws (psi) tiempo pws (psi)
25/07/2006 20:09:03 3850.459 25/07/2006 20:13:13 4212.998 25/07/2006 20:17:23 4484.915 25/07/2006 20:21:33 4597.726
25/07/2006 20:09:13 3865.7 25/07/2006 20:13:23 4226.335 25/07/2006 20:17:33 4492.597 25/07/2006 20:21:43 4599.221
25/07/2006 20:09:23 3881.548 25/07/2006 20:13:33 4239.234 25/07/2006 20:17:43 4500.231 25/07/2006 20:21:53 4600.575
25/07/2006 20:09:33 3896.582 25/07/2006 20:13:43 4252.052 25/07/2006 20:17:53 4507.478 25/07/2006 20:22:03 4601.79
25/07/2006 20:09:43 3912.219 25/07/2006 20:13:53 4264.577 25/07/2006 20:18:03 4514.416 25/07/2006 20:22:13 4602.875
25/07/2006 20:09:53 3927.272 25/07/2006 20:14:03 4277.151 25/07/2006 20:18:13 4521.189 25/07/2006 20:22:23 4603.863
25/07/2006 20:10:03 3942.614 25/07/2006 20:14:13 4289.742 25/07/2006 20:18:23 4527.446 25/07/2006 20:22:33 4604.716
25/07/2006 20:10:13 3957.654 25/07/2006 20:14:23 4301.678 25/07/2006 20:18:33 4533.646 25/07/2006 20:22:43 4605.499
25/07/2006 20:10:23 3972.631 25/07/2006 20:14:33 4313.874 25/07/2006 20:18:43 4539.36 25/07/2006 20:22:53 4606.168
25/07/2006 20:10:33 3987.591 25/07/2006 20:14:43 4325.431 25/07/2006 20:18:53 4544.905 25/07/2006 20:23:03 4606.782
25/07/2006 20:10:43 4002.293 25/07/2006 20:14:53 4337.264 25/07/2006 20:19:03 4550.072 25/07/2006 20:23:13 4607.312
25/07/2006 20:10:53 4017.248 25/07/2006 20:15:03 4348.533 25/07/2006 20:19:13 4554.979 25/07/2006 20:23:23 4607.789
25/07/2006 20:11:03 4031.688 25/07/2006 20:15:13 4359.915 25/07/2006 20:19:23 4559.678 25/07/2006 20:23:33 4608.199
25/07/2006 20:11:13 4046.587 25/07/2006 20:15:23 4370.749 25/07/2006 20:19:33 4564.005 25/07/2006 20:23:43 4608.577
25/07/2006 20:11:23 4060.706 25/07/2006 20:15:33 4381.636 25/07/2006 20:19:43 4568.203 25/07/2006 20:23:53 4608.909
25/07/2006 20:11:33 4075.282 25/07/2006 20:15:43 4392.16 25/07/2006 20:19:53 4571.961 25/07/2006 20:24:03 4609.199
25/07/2006 20:11:43 4089.394 25/07/2006 20:15:53 4402.595 25/07/2006 20:20:03 4575.596 25/07/2006 20:24:13 4609.459
25/07/2006 20:11:53 4103.648 25/07/2006 20:16:03 4412.645 25/07/2006 20:20:13 4578.898 25/07/2006 20:24:23 4609.688
25/07/2006 20:12:03 4117.885 25/07/2006 20:16:13 4422.544 25/07/2006 20:20:23 4582.031 25/07/2006 20:24:33 4609.897
25/07/2006 20:12:13 4131.671 25/07/2006 20:16:23 4432.277 25/07/2006 20:20:33 4584.878 25/07/2006 20:24:43 4610.081
25/07/2006 20:12:23 4145.839 25/07/2006 20:16:33 4441.629 25/07/2006 20:20:43 4587.498 25/07/2006 20:24:53 4610.254
25/07/2006 20:12:33 4159.307 25/07/2006 20:16:43 4450.924 25/07/2006 20:20:53 4589.981 25/07/2006 20:25:03 4610.411
25/07/2006 20:12:43 4173.188 25/07/2006 20:16:53 4459.64 25/07/2006 20:21:03 4592.164 25/07/2006 20:25:13 4610.555
25/07/2006 20:12:53 4186.422 25/07/2006 20:17:03 4468.485 25/07/2006 20:21:13 4594.225 25/07/2006 20:25:23 4610.681
25/07/2006 20:13:03 4199.973 25/07/2006 20:17:13 4476.719 25/07/2006 20:21:23 4596.043 25/07/2006 20:25:33 4610.806
Ejercicio 3.3
) ( 64 . 9
" 25 . 4
10 246 . 1
05 . 0
150
5014 . 1
2833 . 1
4539
1 5
cierre del antes produccin de tiempo hrs t
r
psi x c
mts h
cp
B
bpd q
p
w
t
o
o
=
=
=
=
=
=
=
=

|

Ejercicio 3.3
Grfica de Horner
3800
3900
4000
4100
4200
4300
4400
4500
4600
4700
1 10 100 1000 10000
(tp +dt)/dt
p
w
s

(
p
s
i
)
Grfica de Horner
4600
4605
4610
4615
4620
4625
4630
4635
4640
1 10 100
(tp +dt)/dt
p
w
s

(
p
s
i
)
Ejercicio 3.3

A partir de la pendiente obtenida de la lnea recta generada en la tendencia de
la curva, se puede calcular la permeabilidad de la formacin, mediante la
ecuacin antes mencionada:

1 4621.9
10 4616.3
|
|
.
|
\
|
A
A +
t
t t
p
) ( psi p
ws
( ) ciclo psi m / 6 . 5
1
6 . 5
) 10 log(
6 . 5
) 1 / 10 log(
9 . 4621 3 . 4616
=
=
mh
qB
k
6 . 162
=
( )( )( )
( )
mD k 516
492 6 . 5
5014 . 1 2833 . 1 4539 6 . 162
= =
Ejercicio 3.3
Grfica Semilog
dt vs pws
3800
3900
4000
4100
4200
4300
4400
4500
4600
4700
0.0010 0.0100 0.1000 1.0000 10.0000 dt (hrs)
p
w
s

(
p
s
i
)
1 4615.92
10 4620
t A
) ( psi p
ws
( ) hr psi m / 08 . 4
1
08 . 4
) 10 log(
08 . 4
) 1 / 10 log(
92 . 4615 4620
= = =
=
Ejercicio 3.3
La presin de fondo antes del cierre y la presin a 1 hora del cierre son las
siguientes:

Aplicando la ecuacin:

psi p
psi p
hr t ws
wf
92 . 4615
46 . 3850
1
=
=
= A
(
(
+
|
|
.
|
\
|
= A
= 2274 . 3 log
) 0 (
1513 . 1
2
1
w
wf hr
cr
k
m
t p p
s
|
( )( )( )
(
+
|
|
.
|
\
|
=

2274 . 3
12 / 25 . 4 10 246 . 1 5014 . 1 05 . 0
516
log
08 . 4
46 . 3850 92 . 4615
1513 . 1
2
5
x
s
6 . 208 = s
Ejercicio 3.3
Validacin del Resultado en Saphir:
Valor mnimo de Dao
Partiendo de la Ecuacin de Darcy en estado Pseudoestacionario con Dao,
obtenida anteriormente:

Para el caso de un pozo productor el gasto es positivo, lo cual implica que:

( )
(
+
|
|
.
|
\
|
=
s
r
r
B
p p kh
q
w
e
wf
4
3
ln 2 . 141
0
4
3
ln >=
(
+
|
|
.
|
\
|
s
r
r
w
e
Valor mnimo de Dao
El valor lmite ser cuando;

Despejando el Dao de esta expresin:

4
3
ln
min
+
|
|
.
|
\
|
=
w
e
r
r
s
0
4
3
ln =
(
+
|
|
.
|
\
|
s
r
r
w
e
Valor mnimo de Dao
Los radios de drene muy difcilmente pueden llegar a ser mayores a 1000 mts,
por lo cual se puede establecer que el dao es negativo no podr ser menor a
s=-8.5.

Lo anterior es una explicacin de la limitante de los valores de dao negativos
y el valor lmite depender de la relacin .
R e (mts )
-14
-12
-10
-8
-6
-4
-2
0
1 10 100 1,000 10,000 100,000
D
a
o

S

(
a
d
i
m
)
|
|
.
|
\
|
w
e
r
r
IV.- CURVAS DE DECLINACIN
Declinacin de la Produccin
La produccin de un pozo o un campo tiene la caracterstica de disminuir con el
tiempo, debido a la disminucin de la presin del yacimiento, dicho comportamiento
de declinacin de la produccin obedece a comportamientos matemticos, que
indican el tipo de declinacin:

Exponencial
Hiperblica
Armnica

Partiendo de la definicin del ritmo de declinacin: el cambio del gasto para un
periodo de tiempo determinado, respecto al gasto inicial.

t
q
q q
D
A
=
1
2 1
De forma general

Durante la vida productiva de un pozo a medida que el tiempo incrementa, disminuye
el gasto de produccin, por lo cual la ecuacin presenta valores negativos que
indican la declinacin. Por convencin a la ecuacin anterior se le introduce un signo
negativo para de esta manera indicar que valores positivos indican declinacin.

| | ao
q
t
q
D / 1
1
A
A
=
dt
dq
q
D
1
=
De la definicin de derivada de un logaritmo natural, tenemos que:
( ) | |
(
=
dy
dx
x dy
x d 1 ln
| |
dt
q d
D
) ln(
=
Arps (1954) defini que existen 3 diferentes formas de declinacin de la produccin,
antes mencionadas, Exponencial, Hiperblica y Armnica, las cuales estn definidas
por 3 parmetros:

- Gasto inicial,
- Factor o ritmo de declinacin,
- Grado de curvatura (variacin con el tiempo del inverso del factor de declinacin),
definido por la variable

El exponente que define el cambio respecto al tiempo del reciproco del ritmo de
declinacin, se establece de forma matemtica de la siguiente forma:

i
q
D
b
b
( )
dt
D d
b
/ 1
=
Al integrar la ecuacin anterior, desde un tiempo , hasta un tiempo , se tiene
lo siguiente:

0 = t
t
( )
} }
=
D
D
t
i
D d dt b / 1
0
D
D
t
i
D bt / 1
0
=
i
D D
bt
1 1
=
Despejando el valor del factor de declinacin:
1 +
=
i
i
btD
D
D
Sustituyendo en la ecuacin general del factor de declinacin

La cual es la ecuacin general de la declinacin de la produccin, donde el
exponente definir la forma de la curva de declinacin y por lo tanto el tipo
de declinacin:

; Declinacin Exponencial
; Declinacin Hiperblica
; Declinacin Armnica
| |
dt
q d
btD
D
i
i
) ln(
1
=
+
b
0 = b
1 0 < < b
1 = b
Declinacin Exponencial:

Esta declinacin se da cuando el exponente es igual a 0, por lo tanto la
ecuacin general de la declinacin ser la siguiente:

Establece la condicin de un Pozo produciendo a presin constante en un
yacimiento cerrado.

Al integrar de 0 a t;

| |
dt
q d
D
) ln(
=
b
| |
} }
=
q
q
t
q d Ddt
0
) ln(
0
| | ) ln( ) ln(
0
q q Dt =

t
q
q
t
q q
D
|
|
.
|
\
|
=
=
0
0
ln
) ln( ) ln(
t
q
q
D
|
|
.
|
\
|
=
0
ln
Dt
e q q
=
0
1,000
10,000
0 1 2 3 4 5 6 7 8 9 10
Q
o
(
b
p
d
)
tiempo (aos)

Para determinar el factor de Declinacin, D, se realiza una grfica de

t
q
q
D
|
|
.
|
\
|
=
0
ln
( ) t vs q log
1,000
2,000
3,000
4,000
5,000
6,000
7,000
8,000
9,000
10,000
11,000
0 1 2 3 4 5 6 7 8 9 10
Q
o
(
b
p
d
)
tiempo (aos)

Pronostico de Produccin:
Dt
e q q

=
0
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
1.6
1.8
11/2010 05/2011 11/2011 05/2012 11/2012 05/2013 11/2013
Q
o

(
M
B
P
D
)
Declinacin Exponencial
Historia Pronstico
t
q
q
D
|
|
.
|
\
|
=
0
ln
Dt
e q q

=
0
Resumen Declinaciones, Arps (1954)

Declinacin
Exponencial
Declinacin
Hiperblica
Declinacin
Armnica (b=1)
Gasto de Produccin,
q
Dt
e q q
=
0

( )
b
i
i
t bD
q
q
/ 1
1+
=
( ) t D
q
q
i
i
+
=
1

Factor de
Declinacin, D
t
q
q
D
|
|
.
|
\
|
=
0
ln

b
i
i
q
q
D D
|
|
.
|
\
|
=
|
|
.
|
\
|
=
i
i
q
q
D D

1,000
2,000
3,000
4,000
5,000
6,000
7,000
8,000
9,000
10,000
11,000
0 2 4 6 8 10 12
Q
o

(
b
p
d
)
tiempo (aos)
Curvas de Declinacin
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
b
1,000
10,000
0 2 4 6 8 10 12
Q
o

(
b
p
d
)
tiempo (aos)
Curvas de Declinacin
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
b
0.001
0.010
0.100
1.000
10.000
100.000
0.0001 0.001 0.01 0.1 1 10 100 1000
q
D
d
=

q
(
t
)
/
q
i
tDd=Di*t
Curva Tipo de Declinacin de Arps
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
Tiempos Cortos????
Declinacin Transitoria
b
Fetkovich:

Declinacin Transitoria:

Declinacin Transitoria
(Fetkovich)
Declinacin Pseudoestacionaria
(Arps)
Produccin Acumulada (Np):

periodos das n
q
MMBls N
n
i
i
p
= =
=1
6
10
) (
0.00
0.10
0.20
0.30
0.40
0.50
0.60
0.70
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
1.6
1.8
11/2010 05/2011 11/2011 05/2012 11/2012 05/2013 11/2013
N
p

(
M
M
b
l
s
)
Q
o

(
M
B
P
D
)
Relacin entre Declinacin Exponencial a Declinacin Lineal:

Igualando

Declinacin Lineal a Exponencial Declinacin Exponencial a Lineal

Dt
e q q
=
0
| | t D q q
% 0
1 =
| | t D q e q
Dt
% 0 0
1 =
| |
%
1 D e
D
=
| |
%
1 ln D D =
D
e D

=1
%
Relacin entre Declinacin Exponencial a Declinacin Lineal:

0
2,000
4,000
6,000
8,000
10,000
12,000
14,000
0.0 0.5 1.0 1.5 2.0 2.5
Q
o

(
b
p
d
)
Tiempo (aos)
Comparacin Declinacin Lineal vs Exponencial
Qo exp (bpd)
Qo lineal (bpd)
Ejercicio 1:

Se desea realizar un pronstico de produccin de un yacimiento a los 10 aos
de explotacin, en base al comportamiento de produccin de tres aos
indicado en la tabla inferior.

Calcular:

a) Verificar Tipo de Declinacin
b) Declinacin Exponencial
c) Declinacin Lineal (%)
d) Gasto (qo) del Yacimiento a los 10 aos de Explotacin
e) Produccin Acumulada (Np) durante los 10 aos
Tiempo
(aos) Qo (bpd)
0 2,164,012
0.25 2,052,041
0.5 2,005,899
0.75 1,944,727
1 1,984,086
1.25 1,890,846
1.5 1,767,937
1.75 1,705,712
2 1,606,503
2.25 1,588,357
2.5 1,554,197
2.75 1,410,148
3 1,219,762
Ejercicio 1:

0
500,000
1,000,000
1,500,000
2,000,000
2,500,000
0 0.5 1 1.5 2 2.5 3 3.5
q
o

(
b
p
d
)
t (aos)
Produccin Promedio Trimestral
100,000
1,000,000
10,000,000
0 0.5 1 1.5 2 2.5 3 3.5
q
o

(
b
p
d
)
t (aos)
Ejercicio 1:

t q
0 2164012
2.75 1410148
qo (bpd) 2164012
D (1/ao) 0.156
D % 14.42
0
500,000
1,000,000
1,500,000
2,000,000
2,500,000
0 2 4 6 8 10 12
q
o

(
b
p
d
)
t (aos)
Dt
e q q
=
0
t
q
q
D
|
|
.
|
\
|
=
0
ln
qo @ 10 aos =456,000 bpd
Ejercicio 1:

Np @ 10 aos = 3,872 MMbls
0
500
1000
1500
2000
2500
3000
3500
4000
0
500,000
1,000,000
1,500,000
2,000,000
2,500,000
0 1 2 3 4 5 6 7 8
q
o

(
b
p
d
)
t (aos)
V.- REGISTROS DE PRODUCCIN
Contenido

Registro Esttico de Presin por Estaciones

Registro Dinmico de Presin

Registro de Temperatura

Registro de aportacin de flujo (Molinete)

Interpretacin del Conjunto de registros

Definicin:

Medicin de la presin a diferentes profundidades del pozo mientras este permanece
cerrado, generando un perfil de presiones.

Aplicaciones:

- Determinar la presin de fondo esttica del pozo (Pws)
- Identificar los fluidos presentes en el pozo
- Evaluar el nivel de fluidos para el diseo de un sistema artificial de produccin
- Determinar presencia de agua en el intervalo disparado

Interpretacin de Registro Esttico:

GRADIENTE DE PRESION POR ESTACIONES CON POZO CERRADO
PROF. TVD PRESION
GRADIENTE
TEMPERATURA PRESION Y TEMP. EN CABEZA
mD mV psia Kg/cm
2
Kg/cm
2
/m F C psia Kg/cm
2
F
0 0 269.798 18.973 0.0000 77.87 25.48 267.534 18.814 74.76
500 500 279.798 19.676 0.0014 147.69 64.27 267.361 18.802 74.63
750 750 287.020 20.184 0.0020 155.94 68.85 267.304 18.798 74.67
1000 1000 293.630 20.649 0.0019 161.29 71.83 267.275 18.796 74.86
1253 1250 299.636 21.071 0.0017 168.76 75.98 266.405 18.735 75.01
1515 1500 363.278 25.547 0.0179 177.18 80.66 266.702 18.755 75.18
1780 1750 664.526 46.732 0.0847 185.76 85.42 266.593 18.748 75.43
2043 2000 961.523 67.618 0.0835 193.29 89.61 266.521 18.743 75.18
2297 2250 1253.327 88.138 0.0821 199.07 92.82 266.453 18.738 74.82
2547 2500 1541.840 108.428 0.0812 202.46 94.70 266.371 18.732 74.52
2695 2647 1711.491 120.358 0.0812 203.42 95.23 266.322 18.729 74.17
2730 2682 1759.567 123.739 0.0966 205.75 96.53 266.269 18.725 73.87

Nivel Medio del Intervalo Productor (NMIP):
Profundidad media del intervalo disparado
en el pozo.

NMIP = (2680+2710)/2 = 2695 md

ESTADO MECANICO
ACTUAL
150.00 md
T.R. 30
TR 11 7/8 y 9 5/8
1. Side track
Pez de 30 m. de long.
2112 2142 md.
Boca de pez a +/- 2786.54 md
(20 de agosto de 2005)
2807MD
161.00 md
158.00 md T.R. 16
T.R. 11 7/8
T.R. 13 3/8
sssv sssv
Vlvula de tormenta
de 7 a 146.06 md
Pws =120.36 kg/cm^2
GRADIENTE DE PRESION POR ESTACIONES CON
POZO CERRADO
PROF. TVD PRESION
GRADIENTE
mD mV psia Kg/cm
2
Kg/cm
2
/m
0 0 269.798 18.973
500 500 279.798 19.676
750 750 287.020 20.184
1000 1000 293.630 20.649
1253 1250 299.636 21.071
1515 1500 363.278 25.547
1780 1750 664.526 46.732
2043 2000 961.523 67.618
2297 2250 1253.327 88.138
2547 2500 1541.840 108.428
2695 2647 1711.491 120.358
2730 2682 1759.567 123.739

Gradiente: Cambio de la presin
medida respecto a la profundidad

1
1
= V
i i
i i
D D
p p
p
(
= V
m
cm kg
p
2
/
0014 . 0
0 500
973 . 18 676 . 19

GRADIENTE DE PRESION POR ESTACIONES CON
POZO CERRADO
PROF. TVD PRESION
GRADIENTE
mD mV psia Kg/cm
2
Kg/cm
2
/m
0 0 269.798 18.973 0.0000
500 500 279.798 19.676 0.0014
750 750 287.020 20.184 0.0020
1000 1000 293.630 20.649 0.0019
1253 1250 299.636 21.071 0.0017
1515 1500 363.278 25.547 0.0179
1780 1750 664.526 46.732 0.0847
2043 2000 961.523 67.618 0.0835
2297 2250 1253.327 88.138 0.0821
2547 2500 1541.840 108.428 0.0812
2695 2647 1711.491 120.358 0.0812
2730 2682 1759.567 123.739 0.0966
0
500
1000
1500
2000
2500
3000
0 20 40 60 80 100 120 140 160
0
500
1000
1500
2000
2500
3000
0 50 100 150 200 250 300 350
T
V
D

[
m
V
]
TEMPERATURA [C]
T
V
D

[
m
V
]
PRESION [Kg/cm
2
]
GRADIENTE ESTATICO POR ESTACIONES
PRESION
TEMPERATURA

0
500
1000
1500
2000
2500
3000
0 20 40 60 80 100 120 140 160
0
500
1000
1500
2000
2500
3000
0 50 100 150 200 250 300 350
T
V
D

[
m
V
]
TEMPERATURA [C]
T
V
D

[
m
V
]
PRESION [Kg/cm
2
]
PRESION
TEMPERATURA
Gas: Incremento pequeo de presin con la
profundidad.
Aceite: Incremento mayor de presin con la
profundidad.

0
500
1000
1500
2000
2500
3000
0 20 40 60 80 100 120 140 160
0
500
1000
1500
2000
2500
3000
0 50 100 150 200 250 300 350
T
V
D

[
m
V
]
TEMPERATURA [C]
T
V
D

[
m
V
]
PRESION [Kg/cm
2
]
PRESION
TEMPERATURA
Nivel de Fluidos dentro del Pozo:
Interseccin de las rectas.

0
500
1000
1500
2000
2500
3000
0 20 40 60 80 100 120 140 160
0
500
1000
1500
2000
2500
3000
0 50 100 150 200 250 300 350
T
V
D

[
m
V
]
TEMPERATURA [C]
T
V
D

[
m
V
]
PRESION [Kg/cm
2
]
PRESION
TEMPERATURA
Presencia de Agua en Intervalo
Productor.
Definicin:

Medicin de la presin a pozo fluyendo, bajando o subiendo la herramienta del
sensor generando un perfil de presiones.

Aplicaciones:

- Determinar las zonas de aportacin en la zona disparada.
- Identificar anomalas en la tubera.
- Evaluar el correcto funcionamiento de la vlvula de bombeo neumtico.
- Ajustar la correlacin de flujo multifsico del pozo para evaluar su comportamiento
de produccin a diferentes condiciones operativas.

E
S
T
A
D
O

M
E
C
A
N
I
C
O
A
C
T
U
A
L
1
5
0
.
0
0

m
d
T
.
R
.

3
0
T
R

1
1

7
/
8

y

9

5
/
8
.

S
i
d
e

t
r
a
c
k
P
e
z

d
e

3
0

m
.

d
e

l
o
n
g
.
2
1
1
2

2
1
4
2

m
d
.
B
o
c
a

d
e

p
e
z

a

+
/
-
2
7
8
6
.
5
4

m
d
(
2
0

d
e

a
g
o
s
t
o

d
e

2
0
0
5
)
2
8
0
7
M
D
1
6
1
.
0
0

m
d
1
5
8
.
0
0

m
d
T
.
R
.

1
6
T
.
R
.

1
1

7
/
8
T
.
R
.

1
3

3
/
8
s
s
s
v
s
s
s
v
V
l
v
u
l
a

d
e

t
o
r
m
e
n
t
a

d
e

7

a

1
4
6
.
0
6

m
d
Presin
P
r
o
f
u
n
d
i
d
a
d

Este registro da origen al registro gradio, amplificando el efecto de la entrada de fluidos.
Gradiente
Definicin:

Medicin de la temperatura a pozo fluyendo, bajando o subiendo la herramienta del
sensor generando un perfil de Temperatura.

Aplicaciones:

- Determinar las zonas de aportacin en la zona disparada.
- Identificar anomalas en la tubera.
- Determinar el tipo de fluido que entra al pozo (lquido o gas).

E
S
T
A
D
O

M
E
C
A
N
I
C
O
A
C
T
U
A
L
1
5
0
.
0
0

m
d
T
.
R
.

3
0
T
R

1
1

7
/
8

y

9

5
/
8
.

S
i
d
e

t
r
a
c
k
P
e
z

d
e

3
0

m
.

d
e

l
o
n
g
.
2
1
1
2

2
1
4
2

m
d
.
B
o
c
a

d
e

p
e
z

a

+
/
-
2
7
8
6
.
5
4

m
d
(
2
0

d
e

a
g
o
s
t
o

d
e

2
0
0
5
)
2
8
0
7
M
D
1
6
1
.
0
0

m
d
1
5
8
.
0
0

m
d
T
.
R
.

1
6
T
.
R
.

1
1

7
/
8
T
.
R
.

1
3

3
/
8
s
s
s
v
s
s
s
v
V
l
v
u
l
a

d
e

t
o
r
m
e
n
t
a

d
e

7

a

1
4
6
.
0
6

m
d
Temperatura
P
r
o
f
u
n
d
i
d
a
d

Aportacin de lquido a travs del intervalo productor.
E
S
T
A
D
O

M
E
C
A
N
I
C
O
A
C
T
U
A
L
1
5
0
.
0
0

m
d
T
.
R
.

3
0
T
R

1
1

7
/
8

y

9

5
/
8
.

S
i
d
e

t
r
a
c
k
P
e
z

d
e

3
0

m
.

d
e

l
o
n
g
.
2
1
1
2

2
1
4
2

m
d
.
B
o
c
a

d
e

p
e
z

a

+
/
-
2
7
8
6
.
5
4

m
d
(
2
0

d
e

a
g
o
s
t
o

d
e

2
0
0
5
)
2
8
0
7
M
D
1
6
1
.
0
0

m
d
1
5
8
.
0
0

m
d
T
.
R
.

1
6
T
.
R
.

1
1

7
/
8
T
.
R
.

1
3

3
/
8
s
s
s
v
s
s
s
v
V
l
v
u
l
a

d
e

t
o
r
m
e
n
t
a

d
e

7

a

1
4
6
.
0
6

m
d
Temperatura
P
r
o
f
u
n
d
i
d
a
d

Aportacin de Gas a travs del intervalo productor (Efecto Joule-Thomson).
Registro de Aportacin
Definicin:

Medicin de las revoluciones por segundo generadas en un molinete al pasar por
una zona donde existe flujo de fluidos. Su medicin continua respecto a la
profundidad genera un perfil de Aportacin.

Aplicaciones:

- Determinar las zonas de aportacin.
- Determinar la aportacin en gasto de las diferentes zonas identificadas.
- Identificar anomalas en la tubera del pozo.
- Apoyo para la realizacin del adecuado programa de estimulacin o tratamiento
en pozos

E
S
T
A
D
O

M
E
C
A
N
I
C
O
A
C
T
U
A
L
1
5
0
.
0
0

m
d
T
.
R
.

3
0
T
R

1
1

7
/
8

y

9

5
/
8
.

S
i
d
e

t
r
a
c
k
P
e
z

d
e

3
0

m
.

d
e

l
o
n
g
.
2
1
1
2

2
1
4
2

m
d
.
B
o
c
a

d
e

p
e
z

a

+
/
-
2
7
8
6
.
5
4

m
d
(
2
0

d
e

a
g
o
s
t
o

d
e

2
0
0
5
)
2
8
0
7
M
D
1
6
1
.
0
0

m
d
1
5
8
.
0
0

m
d
T
.
R
.

1
6
T
.
R
.

1
1

7
/
8
T
.
R
.

1
3

3
/
8
s
s
s
v
s
s
s
v
V
l
v
u
l
a

d
e

t
o
r
m
e
n
t
a

d
e

7

a

1
4
6
.
0
6

m
d
rps
0
P
r
o
f
u
n
d
i
d
a
d

Produccin
E
S
T
A
D
O

M
E
C
A
N
I
C
O
A
C
T
U
A
L
1
5
0
.
0
0

m
d
T
.
R
.

3
0
T
R

1
1

7
/
8

y

9

5
/
8
.

S
i
d
e

t
r
a
c
k
P
e
z

d
e

3
0

m
.

d
e

l
o
n
g
.
2
1
1
2

2
1
4
2

m
d
.
B
o
c
a

d
e

p
e
z

a

+
/
-
2
7
8
6
.
5
4

m
d
(
2
0

d
e

a
g
o
s
t
o

d
e

2
0
0
5
)
2
8
0
7
M
D
1
6
1
.
0
0

m
d
1
5
8
.
0
0

m
d
T
.
R
.

1
6
T
.
R
.

1
1

7
/
8
T
.
R
.

1
3

3
/
8
s
s
s
v
s
s
s
v
V
l
v
u
l
a

d
e

t
o
r
m
e
n
t
a

d
e

7

a

1
4
6
.
0
6

m
d
rps
0
P
r
o
f
u
n
d
i
d
a
d

Inyeccin
Ejemplo 1 Interpretacin de Registros de presin-produccin
E
S
T
A
D
O

M
E
C
A
N
I
C
O
A
C
T
U
A
L
1
5
0
.
0
0

m
d
T
.
R
.

3
0
T
R

1
1

7
/
8

y

9

5
/
8
.

S
i
d
e

t
r
a
c
k
P
e
z

d
e

3
0

m
.

d
e

l
o
n
g
.
2
1
1
2

2
1
4
2

m
d
.
B
o
c
a

d
e

p
e
z

a

+
/
-
2
7
8
6
.
5
4

m
d
(
2
0

d
e

a
g
o
s
t
o

d
e

2
0
0
5
)
2
8
0
7
M
D
1
6
1
.
0
0

m
d
1
5
8
.
0
0

m
d
T
.
R
.

1
6
T
.
R
.

1
1

7
/
8
T
.
R
.

1
3

3
/
8
s
s
s
v
s
s
s
v
V
l
v
u
l
a

d
e

t
o
r
m
e
n
t
a

d
e

7

a

1
4
6
.
0
6

m
d
rps
0
P
r
o
f
u
n
d
i
d
a
d

Presin
Grad Temp
Pozo productor de lquidos (Aceite-Agua) con todo el Intervalo aportando.
E
S
T
A
D
O

M
E
C
A
N
I
C
O
A
C
T
U
A
L
1
5
0
.
0
0

m
d
T
.
R
.

3
0
T
R

1
1

7
/
8

y

9

5
/
8
.

S
i
d
e

t
r
a
c
k
P
e
z

d
e

3
0

m
.

d
e

l
o
n
g
.
2
1
1
2

2
1
4
2

m
d
.
B
o
c
a

d
e

p
e
z

a

+
/
-
2
7
8
6
.
5
4

m
d
(
2
0

d
e

a
g
o
s
t
o

d
e

2
0
0
5
)
2
8
0
7
M
D
1
6
1
.
0
0

m
d
1
5
8
.
0
0

m
d
T
.
R
.

1
6
T
.
R
.

1
1

7
/
8
T
.
R
.

1
3

3
/
8
s
s
s
v
s
s
s
v
V
l
v
u
l
a

d
e

t
o
r
m
e
n
t
a

d
e

7

a

1
4
6
.
0
6

m
d
rps
0
P
r
o
f
u
n
d
i
d
a
d

Presin
Grad Temp
Pozo productor de Gas con la parte superior del Intervalo aportando.
Clculo de Produccin a partir de Molinete
Los molinetes consisten en una propela que gira debido al paso del fluido, dicho
movimiento se convierte en impulsos elctricos mediante una bobina. La corriente
elctrica generada se monitorea en superficie estableciendo una relacin directa
entre las revoluciones por segundo registradas y el gasto del pozo.

El calculo de la produccin de aceite en los registros PLT mediante el conteo de las
revoluciones de la propela del molinete est basado en el concepto bsico de la
definicin del gasto en funcin de la velocidad

Donde v es la velocidad con que gira el molinete y A el rea transversal de flujo del
interior de la tubera donde se encuentra el molinete.

Por definicin,

A v Q
o
=

4
2
d
A
t
=
La velocidad puede ser obtenida a partir de las revoluciones del molinete que estn
en funcin de la velocidad del flujo del aceite.

Para lo anterior se toman los datos de calibracin del molinete para establecer la
ecuacin que represente el comportamiento del molinete.

Para determinar el gasto de aceite tenemos que conocer una ecuacin que represente el
comportamiento del molinete. Por esta razn se realiza una calibracin previa a la
produccin para simular la produccin del pozo, paseando el molinete en una seccin de la
tubera tomando mediciones o diferentes velocidades del cable en forma descendente y
ascendente para simular una produccin o inyeccin, respectivamente.

Lo anterior esta basado en que el flujo de aceite es el que hace mover las aspas del molinete
cuando este se encuentra esttico, pero se tiene el mismo efecto si el fluido es el que se
encuentra esttico y el molinete es el que se mueve haciendo girar sus aspas dependiendo
de la velocidad del cable con la que se mueva el molinete.

Se puede concluir que la velocidad del cable refleja la velocidad de flujo del aceite en
condiciones de operacin del pozo.
Con los datos de la calibracin que son: , se realiza una
grafica que se ajustara a una lnea recta a la cual se le establece una ecuacin de la forma de
lnea recta

Esta ecuacin se aplicar para las condiciones de flujo del pozo. Cabe mencionar que la
ecuacin establecida para los datos de calibracin bajando (produccin) y subiendo
(inyeccin) son diferentes.

Una vez establecida la ecuacin:

La velocidad con la que se mueve el aceite puede ser despejada de la ecuacin anterior:

molinete cable
es revolucion vs v
b mv rps + =
b mv rps + =
m
b rps
v

=
Sustituyendo la velocidad y el rea en la ecuacin original:

Adicionalmente se maneja un factor de correccin, debido a que el molinete no cubre el
rea total de la tubera, lo anterior esta relacionado con el perfil de velocidades en un
tubera circular generado a partir de la ecuacin de Hagen-Poiseville.

Factor de Correccin
Comn= 0.86

|
|
.
|
\
|
|
.
|
\
|
=
4
2
d
m
b rps
Q
o
t

Adems como las revoluciones son medidas en el fondo, se debe afectar la ecuacin por el
factor de volumen del aceite, adems de utilizar una constante de conversin de unidades
para obtener el gasto de aceite en unidades de campo. La ecuacin anterior queda de la
siguiente forma:

o
c
o
B
F
d
m
b rps
Q
|
|
.
|
\
|
|
.
|
\
|
=
4
7809 . 1
2
t
| |
| |
| |
| |
| |
| |
unidades de conversin de factor
aceite del volumen de factor stbl bl B
correccin de factor a F
tubera la de erno dimetro p d
n calibraci la en ajuste de ecuacin la de pendiente
ft seg
m
n calibraci la en ajuste de ecuacin la de origen al ordenada seg b
molinete del uto por promedio es revolucion seg rps
aceite de gasto bpd Q
o
c
o
; 7809 . 1
; /
; dim
int ; lg
;
min
; / 1
min ; / 1
;
(

Obtencin del factor de conversin de unidades:

Obtencin de m y b:

| | | | 7809 . 1
1
min 1440
615 . 5
1
lg 12
1
lg
min
3 2 2
2
2
=
(
dia ft
bl
p
ft
p
ft
bpd Q
o
Subiendo (Inyeccin) Bajando (Produccin)
Velocidad (ft/min)
Revoluciones
(1/seg) Velocidad (ft/min)
Revoluciones
(1/seg)
-53.89 -3.577 52.54 2.939
-86.9 -6.285 82.57 5.71
-114.87 -8.607 116.85 8.501

-10
-8
-6
-4
-2
0
2
4
6
8
10
-150 -100 -50 0 50 100 150
R
e
v
o
l
u
c
i
o
n
e
s

d
e
l

M
o
l
i
n
e
t
e

(
r
p
s
)
Velocidad del Cable durante la calibracin (ft/seg)
Calibracin del Molinete
Pendiente (1/seg/ft/min) 0.082
Ordenada al Origen (1/seg) 0.868
Pendiente (1/seg/ft/min) 0.086
Ordenada al Origen (1/seg) -1.605
Ejercicio 1 Clculo de Produccin a partir de Molinete
Datos

Dimetro Interno de Tubera (plg) 2.992
Revoluciones promedio (1/seg) 8
Bo (bl/stb) 1.16
Factor de Correccin 0.86
Velocidad (ft/min)
Revoluciones
(1/seg) Velocidad (ft/min)
Revoluciones
(1/seg)
-53.89 -4 52.54 3
-86.9 -5 82.57 4
-114.87 -6 116.85 5
Datos de Calibracin del molinete
-8
-6
-4
-2
0
2
4
6
-150 -100 -50 0 50 100 150
R
e
v
o
l
u
c
i
o
n
e
s

d
e
l

M
o
l
i
n
e
t
e

(
r
p
s
)
Velocidad del Cable durante la calibracin (ft/seg)
Calibracin del Molinete
Pendiente (1/seg/ft/min) 0.033 Pendiente (1/seg/ft/min) 0.031
Ordenada al Origen (1/seg) -2.233 Ordenada al Origen (1/seg) 1.366
rea de la Tubera (plg^2) 7.031
Cte conversin de Unidades 1.781
m (1/seg/ft/min) 0.031
Velocidad del Aceite (ft/min) 213.32

Qo (bpd) 1980
o
c
o
o
B
F
d
m
b rps
Q
A v Q
|
|
.
|
\
|
=
=
4
7809 . 1
2
t
VI.- ANLISIS INTEGRAL DEL POZO
Contenido

Importancia de la Caracterizacin del Fluido y el
efecto de la Temperatura.

Flujo en el Yacimiento

Flujo en la tubera de produccin

Flujo en el Estrangulador

Metodologa del Anlisis Nodal

Importancia de la Caracterizacin del Fluido
Anlisis PVT:

El comportamiento de fases de una sustancia o mezcla queda definido por un
anlisis de Presin-Volumen-Temperatura (PVT).

Generalmente realizando pruebas en laboratorio del comportamiento del volumen
del fluido (Aceite y gas) a diferentes valores de Presin para una determinada
Temperatura (Yacimiento).

Propiedades Obtenidas:

- Factor de Volumen del Aceite (Bo)
- Factor de Volumen del Gas (Bg)
- Relacin de Solubilidad (Rs)
- Compresibilidad del Aceite (Co)
- Compresibilidad del Gas (Cg)
- Viscosidad del Aceite (
o
)
- Viscosidad del Gas (
g
)
- Factor de Desviacin o Compresibilidad del gas (z)
Anlisis PVT

El objetivo de establecer correlaciones PVT que representen el comportamiento del fluido es
poder generar clculos de las propiedades PVT a diferentes presiones y Temperaturas.
Presin (psi) Presin (kg/cm2) Rs (ft3/bl) Bo (bl/stb) Bg (ft3/scf)
3840 270.0 488 1.306
3500 246.1 488 1.306
3000 210.9 488 1.306
2500 175.8 488 1.306
2092 147.1 488 1.306
1800 126.6 440 1.285 0.00845
1600 112.5 401 1.27 0.00956
1400 98.4 365 1.255 0.01103
1200 84.4 328 1.24 0.01298
1000 70.3 290 1.225 0.01569
800 56.2 252 1.21 0.01985
600 42.2 212 1.194 0.02688
400 28.1 170 1.177 0.0405
200 14.1 118 1.154 0.08037
100 7.0 84 1.135 0.15312
0 0.0 0 1.053
Anlisis PVT (Relacin de Solubilidad, Rs)

o
o s
s
g
M
R
R
y

+
=
350
3 . 379
3 . 379 /
119118 . 0 36226 . 1 10526 . 3 043 . 5
2 3
+ + + =
g g g f
y y y p
0
100
200
300
400
500
600
0 1000 2000 3000 4000 5000
R
s

(
f
t
3
/
b
l
)
Presin (psi)
Ajuste Relacin de Solubilidad
Rs (ft3/bl)
0
100
200
300
400
500
600
0 1000 2000 3000 4000 5000
R
s

(
f
t
3
/
b
l
)
Presin (psi)
Rs (ft3/bl)
Rs cor(ft3/bl)
Parametros de Ajuste
Parametro1
1.05
Parametro2
75
2 1
P R P R
s corregida s
+ =
0
100
200
300
400
500
600
0 1000 2000 3000 4000 5000
R
s

(
f
t
3
/
b
l
)
Presin (psi)
Rs (ft3/bl)
Rs cor(ft3/bl)
( )
g
f
T p
p
460 +
=
Lasater
Anlisis PVT (Factor de Volumen del Aceite ,Bo)

Parametros de Ajuste
Parametro1
0.925
Parametro2
0.117551
1
1.05
1.1
1.15
1.2
1.25
1.3
1.35
0 1000 2000 3000 4000 5000
B
o

(
m
3
/
m
3
)
Presin (psi)
Ajuste Factor de Volumen del Aceite
Bo (bl/stb)
2 1
P B P B
o corregida o
+ =
| |
| |
| |
| |
| | bl bl B
F T
a
a
bl ft R
T R B
o
o
g
s
o
g
s o
/
dim
dim
/
25 . 1 00014 . 0 972 . 0
3
175 . 1
5 . 0
(
(
+
|
|
.
|
\
|
+ =
1
1.05
1.1
1.15
1.2
1.25
1.3
1.35
0 1000 2000 3000 4000 5000
B
o

(
m
3
/
m
3
)
Presin (psi)
Bo (bl/stb)
Bo (bl/stb)_Ajuste
1
1.05
1.1
1.15
1.2
1.25
1.3
1.35
0 1000 2000 3000 4000 5000
B
o

(
m
3
/
m
3
)
Presin (psi)
Bo (bl/stb)
Bo (bl/stb)_Ajuste
Aplicaciones del PVT en el Anlisis Nodal:

Comportamiento de Afluencia (Yacimiento)

( )
(
+
|
|
.
|
\
|
=
s
r
r
B
p p kh
q
w
e
wf
4
3
ln 2 . 141

Flujo Multifsico (Pozo-Superficie)

N ACELERACI ELEVACIN FRICCIN TOTAL
dL
dp
dL
dp
dL
dp
dL
dp
|
.
|
\
|
+
|
.
|
\
|
+
|
.
|
\
|
=
|
.
|
\
|
(
+ =
A
A
m
L
m
D x
M q f
h
p
5 11
2 2
10 9652 . 2 144
1
Hagedorn&Brown
Estado Mecnico
Propuesto
16
800md
20
200md
TR 11 7/8
2246md
LINER 9 5/8
2681md
BL 9 5/8 2176md
LINER 7
3188md
Agujero de 6 1/8
a 3588md
BL 7 2581md
Vlvulade Tormenta 5 @ 150 md
Estado Mecnico
Propuesto
16
800md
20
200md
TR 11 7/8
2246md
LINER 9 5/8
2681md
BL 9 5/8 2176md
LINER 7
3188md
Agujero de 6 1/8
a 3588md
BL 7 2581md
Comportamiento Bo
Estado Mecnico
Propuesto
16
800md
20
200md
TR 11 7/8
2246md
LINER 9 5/8
2681md
BL 9 5/8 2176md
LINER 7
3188md
Agujero de 6 1/8
a 3588md
BL 7 2581md
Comportamiento o
Flujo en el Yacimiento
El flujo del fluido en el Yacimiento queda definido por cualquier modelo de Afluencia;
Vogel, Darcy, Composite, Fetkovich, Fracturamiento Hidrulico, etc. Dependiendo de
las condiciones del yacimiento-pozo.

Modelo de Darcy:

(
+
|
|
.
|
\
|
=
s
r
r
B
kh
IP
w
e
4
3
ln 2 . 141
0
200
400
600
800
1000
1200
1400
1 10 100 1000
I
P

(
b
p
d
/
p
s
i
)
Viscosidad (cp)
Comportamiento de IP vs Viscosidad
Flujo en la Tubera de Produccin
Para determinar el flujo en la tubera, se hace uso de correlaciones de Flujo
Multifsico; entre las ms comunes se encuentran:

Duns and Ros

Hagedorn and Brown

Fancher Brown

Mukerjee Brill

Beggs and Brill

Orkszewski

Poettman and Carperter

Hagedorn and Brown:

Informacin bsica requerida (Pozo con Bombeo Neumtico):

Datos
API 21.4
g 0.922
w 1.026
Ql (bpd) 5000
WOR 0
Pws (Kg/cm2) 141.6
Tws (C) 100
Piny (mv) 2000
Qgi (mmpcd) 3
g lift 0.675
Pwh (Kg/cm2) 13
(
+ =
A
A
m
L
m
D x
M q f
h
p
5 11
2 2
10 9652 . 2 144
1
Hagedorn and Brown:
Perfil de Presin

Flujo a travs del Estrangulador
Existen diferentes correlaciones en la literatura para poder determinar el
comportamiento de flujo a travs de estranguladores las cuales se clasifican en dos
tipos:

a) Flujo subcrtico
b) Flujo crtico

El flujo crtico se da cuando existe una relacin de presiones entre la entrada y la
salida del estrangulador de aproximadamente 0.52, es decir y el
flujo subcrtico obedece a comportamientos donde dicha relacin tiene un valor
mayor a 0.52.

Cuando se habla flujo crtico, se indica que para el dimetro de estrangulador
definido se ha logrado un grado de saturacin tal que ante una variacin de presin
a la entrada no se tiene efecto alguno en la de salida, manteniendo la relacin de
0.52.

entrada salida
p p 52 . 0 =
Entre las correlaciones ms comunes existentes en la literatura se encuentran las de
Gilbert, Ros, Baxendell y Achong. Las cuales obedecen a la ecuacin:

( )
C
B
l
th
RGA q A
p
64
=
|
A B C
GILBERT 10 0.546 1.89
ROS 17.4 0.5 2
BAXENDELL 9.56 0.546 1.93
ACHONG 3.82 0.65 1.88
| |
| |
| |
| | lg
/
3
p
bl ft RGA
bpd q
psi p
l
th
|
Ejemplo:

Determinar la presin en cabeza del pozo requerida para un flujo de 5000 bpd de
aceite a travs de un estrangulador de 2, donde el RGA es de 350 ft^3/bl.

( )
C
B
l
th
RGA q A
p
64
=
|
p
th
(psi)
23
16.6
18.1
20.3
A B C
GILBERT 10 0.546 1.89
ROS 17.4 0.5 2
BAXENDELL 9.56 0.546 1.93
ACHONG 3.82 0.65 1.88
p
baj
(psi)
12.0
8.6
9.4
10.6
Metodologa de Anlisis Nodal
Anlisis Nodal: Punto de solucin entre la capacidad de Aportacin del Yacimiento y
la capacidad de Levantamiento del Pozo.

Metodologa General

a) Determinar la capacidad de Aportacin del Yacimiento (IPR)
b) Determinar la curva de levantamiento del Pozo (VLP)
c) Interseccin entre curvas ser la produccin del Pozo

Determinar la capacidad de Aportacin del Yacimiento (IPR)
0
20
40
60
80
100
120
140
160
0 5000 10000 15000 20000
P
r
e
s
i
n

(
k
g
/
c
m
2
)
Gasto (bpd)
Determinar la curva de levantamiento del Pozo (VLP)

Generar gradientes de presin para diferentes gastos manteniendo el resto de las
variables constantes
0
500
1000
1500
2000
2500
3000
3500
0 20 40 60 80 100 120 140 160
P
r
o
f
u
n
d
i
d
a
d

(
m
)
Presin (kg/cm2)
300
0
500
1000
1500
2000
2500
3000
3500
0 20 40 60 80 100 120 140 160
P
r
o
f
u
n
d
i
d
a
d

(
m
)
Presin (kg/cm2)
300
500
0
500
1000
1500
2000
2500
3000
3500
0 20 40 60 80 100 120 140 160
P
r
o
f
u
n
d
i
d
a
d

(
m
)
Presin (kg/cm2)
300
500
1000
0
500
1000
1500
2000
2500
3000
3500
0 20 40 60 80 100 120 140 160
P
r
o
f
u
n
d
i
d
a
d

(
m
)
Presin (kg/cm2)
300
500
1000
2000
0
500
1000
1500
2000
2500
3000
3500
0 20 40 60 80 100 120 140 160
P
r
o
f
u
n
d
i
d
a
d

(
m
)
Presin (kg/cm2)
300
500
1000
2000
5000
0
500
1000
1500
2000
2500
3000
3500
0 20 40 60 80 100 120 140 160
P
r
o
f
u
n
d
i
d
a
d

(
m
)
Presin (kg/cm2)
300
500
1000
2000
5000
7000
0
500
1000
1500
2000
2500
3000
3500
0 20 40 60 80 100 120 140 160
P
r
o
f
u
n
d
i
d
a
d

(
m
)
Presin (kg/cm2)
300
500
1000
2000
5000
7000
10000
Determinar la curva de levantamiento del Pozo (VLP)

Graficar los valores de Pwf requeridos para cada gasto.
qo (bpd) Pwf (kg/cm^2)
300 104.8
500 112.5
1000 113.9
2000 116.6
5000 124.9
7000 131.4
10000 140.1
0
20
40
60
80
100
120
140
160
0 2000 4000 6000 8000 10000 12000
P
r
e
s
i
n

(
k
g
/
c
m
2
)
Gasto (bpd)
Interseccin entre curvas IPR y VLP
Equilibrio entre lo que puede aportar el Yacimiento y lo que puede levantar el pozo.
El valor de la interseccin entre las curvas indica el gasto que se podr obtener del
pozo

Interseccin entre curvas IPR y VLP

Si se varia la condicin operativa (pwh)

Productividad de Pozos

Ing. Ricardo Posadas Mondragn

Taller de Productividad de Pozos Capítulo SPE

Cargado por

Copyright:

Formatos disponibles

Taller de Productividad de Pozos Capítulo SPE

Cargado por

Información del documento

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

Taller de Productividad de Pozos Capítulo SPE

Cargado por

Copyright:

Formatos disponibles

Ricardo Posadas Mondragn Productividad de Pozos

Captulo Estudiantil SPE, UNAM

Tomando el lmite cuando L0, en (2)

Ricardo Posadas Mondragn Productividad de Pozos

Mtodo de IPR de Pozos Horizontales (Estacionario)

Mtodo de IPR de Pozos Horizontales (Estacionario)

Mtodo de IPR de Pozos Horizontales (Estacionario)

Ricardo Posadas Mondragn Productividad de Pozos

V.T. 5 1/2 150 M.

Ricardo Posadas Mondragn Productividad de Pozos

Ricardo Posadas Mondragn Productividad de Pozos

También podría gustarte