Introduccion Al Analisis de Pruebas de Presion PDF
Introduccion Al Analisis de Pruebas de Presion PDF
Introduccion Al Analisis de Pruebas de Presion PDF
DISEÑO DE PRUEBAS DE
PRESIÓN
tD
www.daprat.com
1
GIOVANNI DA PRAT
INTRODUCCIÓN AL ANÁLISIS
Y DISEÑO DE PRUEBAS
DE PRESIÓN
3
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 5
PRÓLOGO
El objetivo principal de este libro o manual es el de capacitar al profesional sin
experiencia en el tema de pruebas de presión, sobre todo en los conceptos
básicos que sustentan el análisis de los datos de presión, que se obtienen de
las pruebas que se realizan en los pozos exploratorios o de desarrollo.
La mayoría de los conceptos teóricos que se presentan en este libro, han sido
ya publicados en numerosos artículos; los más relevantes para cada capítulo
se citan en las referencias que se presentan en la Bibliografía y se recomienda
su lectura en caso de querer profundizar en un tema particular. Con el fin de
complementar los conceptos teóricos, se presentan las ecuaciones que definen
las soluciones básicas para cada método, así como ejercicios, en donde se
enfatiza el cálculo manual para arribar a la solución y se presenta su
comparación con los resultados obtenidos vía software comercial (Saphir).
ÍNDICE
9
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 11
CAPÍTULO I
INTRODUCCIÓN A LAS PRUEBAS DE PRESIÓN
1.1. Definición
Desde el punto de vista operacional y de campo, una prueba de presión
consiste en la adquisición de datos de presión, producción y muestra del fluido
a condiciones de pozo (aperturas y cierres) controladas. La completación del
pozo objeto de prueba puede ser temporal o permanente. A continuación se
presenta el ejemplo de una prueba tipo multi-tasa a fin de ilustrar el concepto
señalado.
Período de flujo a
diferentes
reductores
1/4
Período de
3/8 restauración
de presión
Presión de fondo
(registrada por sensor
de presión) 1/2
Tasas de q4 = 0
producción q3
(pozo cerrado a
q2
producción)
q1
Presión
DST (con taladro)
Muestra de fluido para análisis PVT
Muestreo
Exploratorio
Permeabilidad y daño
Prueba sin taladro
Potencial del pozo e índice de
Probador de Formación
productividad
Permeabilidad y daño
Restauración, multitasa
Presión actual y promedio
Interferencia
Productor Tipo de límites asociados con el área
Sensores de presión
de drenaje
permanentes
Monitoreo continuo de presión de
Gradientes de presión
fondo
Pozo Productor: obtener el valor de la presión actual del área de drenaje así
como evaluar la eficiencia de flujo del pozo. Determinar los límites del área de
drenaje. Probar comunicación hidráulica (interferencia) entre pozos.
El departamento de ingeniería de yacimientos y producción define la frecuencia
con que deben de hacerse las pruebas en estos pozos (ver Capítulo III para el
análisis de los datos).
Las pruebas de presión se aplican también en pozos que forman parte del
proceso de recuperación secundaria (inyectores de agua, por ejemplo). Como
se detalla en el Capítulo VII, el análisis de los datos de presión registrados en
el pozo inyector, durante el período de cierre de la inyección (período
denominado fall off), permite determinar el valor de la presión media del patrón
de inyección al momento de realizar la prueba; dicho valor resulta muy
importante y está asociado con la eficiencia del barrido de petróleo por agua.
Pc
qs
4-1/2” (tubería)
perfil de presiones
p = p(r,tp)
9-5/8” sensor de
presión
7”
P=Pi P = Pi (t = 0)
pf
p
AD
límite
qf h
rd
re
rw r
Fig. 1.2. Perfil de presión en el yacimiento. Pozo en producción a una tasa de qs. Pi es la
presión inicial (t=0) del yacimiento antes de la apertura del pozo a producción
Esto se deduce con base al perfil de presiones, ya que a partir del valor del
radio de drenaje rd la presión no muestra variación radial apreciable, es decir,
es prácticamente constante, por lo tanto no hay gradiente de presión; lo que se
traduce en que la tasa es prácticamente igual a cero a partir del radio de
drenaje rd.
La forma más simple de visualizar una prueba de presión y digamos una cuyos
objetivos de evaluación sean los de determinar la permeabilidad efectiva, daño
y Pactual o la presión del yacimiento luego de un acumulado de producción al
tiempo tp, sería la de registrar las presiones en el fondo del pozo, pf, así como
medir la tasa de producción durante el tiempo de producción y cierre del pozo
(durante el cierre la tasa es cero).
4850
pf = Pi
pf(t)
4650
4450
∆p (psi)
100
10
Presión derivativa
(∆t.∆p’)
1
1E-4 1E-3 0.01 0.1 1 10
Datos conocidos: petróleo como fluido productor, radio del pozo, rw = 0.3 ft,
viscosidad=0.4cp, espesor de la capa (h)=30ft, Factor volumétrico,
(Bo)=1.31B/STB, compresibilidad total, ct=12.510-6psi-1, Temperatura del
yacimiento = 212 °F, Pi = 5000 psia, Área de drenaje total = 650 Acres (re =
3000 ft).
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 19
21
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 23
CAPÍTULO II
TEORÍA DE FLUJO DE FLUIDOS EN EL YACIMIENTO
Vr
h
r r re
Fig. 2.1. Balance de masa (Flujo radial)
(r r )hVrt (1)
24 // GIOVANNI DA PRAT
1
rVr
r r t (4)
kr p
Vr
r (5)
1 rkr p
r r r t t t (6)
La densidad del fluido está relacionada con la compresibilidad del mismo. Dos
casos son de interés en ingeniería de yacimientos: liquido (petróleo) y gas.
c p po
o e
(7)
2 p 1 p c p
r r r k t
2
(8)
M p
(9)
RT z
26 // GIOVANNI DA PRAT
1 p p p
r
r r ( p) z ( p) r k t z ( p) (10)
Presión adimensional:
PD
kh
Pi P(r , t )
141.2qB (11)
0.000264kt
tD
cr 2 w (12)
r 2w
tDA tD (13)
A
Radio adimensional:
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 27
r
rD (14)
rw
Las unidades para los parámetros involucrados en estas definiciones son las
de campo:
2 pD 1 pD pD
rD rD rD tD
2 (15)
r P( r , t )
h re
k
pozo
1 rD 2
PD rD, tD i (16)
2 4tD
eu
i x du (17)
x u
1 tD
PDrD, tD
2 rD 2
ln 0.80907 (19)
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 29
tD
Válida para 100. Sin embargo, la diferencia es de sólo un 2% cuando
rD 2
tD
es mayor que 5. Esta solución define el régimen transiente de flujo
rD 2
denominado tipo infinito (infinite acting), comúnmente también denominado
“flujo radial”.
Desde el punto de vista práctico, nos interesa la solución de la presión a nivel
de pozo, que es donde comúnmente se mide la misma, es decir en r rw (
rD 1 ). En este caso la EC. 19 se puede expresar como:
1
PD ln tD 0.80907
2
(20)
Ejemplo1
Cálculo del valor de la presión en el pozo al cabo de 500 horas de producción.
Solución:
Substituyendo las definiciones de tiempo adimensional, presión adimensional y
radio adimensional en la EC. 20, se obtiene que la presión en el pozo viene
dada por:
162.6qB k
Pwf Pi log t log
crw2 3.2275 (21)
kh
2850
Solución para la presión en
función del tiempo obtenida
mediante software comercial
2750
[psia]
2550
20
[m3/D]
10
500 Hr.
0
0 100 200 300 400 500 600 70
Ejemplo 2
Cálculo del valor de la presión a 300 metros (984 ft) del pozo.
Se desea saber cuál es la presión a una distancia radial de 300 metros del
pozo, al cabo de 500 horas de producción:
Solución:
Substituyendo las definiciones de tiempo adimensional, presión adimensional y
radio adimensional en la EC. 19, se obtiene que:
162.6qB t k
Pr , t Pi log log
crw2 3.2275 (22)
r
2
kh
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 31
tD
El valor de 2
es igual a 11.6. Por lo tanto podemos usar la aproximación
rD
logarítmica con un 2% de error. Substituyendo valores, se obtiene que la
presión a 300 metros del pozo es de 2945 psia, al cabo de 500 horas de
producción.
3000
2980
[psia]
2960
2945,6
i
2940
20
[m3/D]
10
500 hr.
0
0 100 200 300 400 500 600 70
Área Circular
Pozo
re
Área Rectangular
L3
L4
L2
pozo
L1
1 A 2.2458
PD 2tDA ln 2 ln (23)
2 rw CA
cA
tpss tDA pss (24)
0.0002637 k
Definiendo m* como:
0.23395qB
m*
chA
70.60qB A 2.2458
Pwf Pi m*t ln 2 ln (26)
kh rw CA
Ejemplo 3
Cálculo del valor de la presión en el pozo a 240 horas producción.
34 // GIOVANNI DA PRAT
Asumamos los mismos valores para los parámetros de yacimiento y pozo que
los dados en el ejercicio 1:
Solución:
Primero debemos asegurarnos que el régimen de flujo es semi-estacionario: el
valor de tDA es de 1.68, por lo tanto mayor que 0.1. De aquí que se cumpla la
condición requerida en cuanto al tiempo de prueba.
La Fig. 2.6. muestra la solución para la presión en función del tiempo generada
mediante software asumiendo los mismos parámetros de diseño en cuanto a
yacimiento y pozo:
2900
Comienza conducta de
decaimiento lineal de la
presión (14.6 horas)
[psia]
2500
2387.3 psia
2100
20
[m3/D]
10
240 hrs.
0
0 100 200 300 400
550
Solución sistema con límite cerrado
450
350
P
(3000-pwf) Pozo
250
Re = 300 m
log t
0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6
0.23395qB
m*
chA
0.23395qB
hA
ct m*
1000
Solución presión
límite de no-flujo
log P
Solución presión
sistema infinito
log t
-2 -1 0 1 2
Los valores del cambio de presión usados para elaborar la Fig. 2.8. son los
mismos que los usados en la Fig. 2.7., sólo se cambiaron las escalas de las
gráficas. Esto es con la finalidad de visualizar la respuesta en esta escala.
31.62 0.1
30.8828 0.1
27.6 0.2
31.6 0.1
12.9851 0.7
Cuadrada
21.8369 0.3
4.5141 1.5
Rectangular 1:2
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 39
5.3790 0.8
0.2318 4
Rectangular: 1:4
Asumamos los mismos valores para los parámetros de yacimiento y pozo que
los del ejemplo anterior. En este caso el pozo se encuentra ubicado en un área
de drenaje mayor, tal como se presenta a continuación:
40 // GIOVANNI DA PRAT
pozo
701 m
1041 m
2800
2400
2200
re
2
re
PD ln (27)
rw
khPi Pw
q
re (28)
141.2 B ln
rw
Esta es la Ley de Darcy que aplica para el caso de flujo tipo estacionario.
Ejemplo 4
Cálculo de la presión en el pozo al cabo de 240 horas de producción.
42 // GIOVANNI DA PRAT
Usando los mismos datos del ejemplo 3, tenemos que tDA = 1.68 (mayor que
0.40). De aquí que se cumpla que el estado de flujo es estacionario.
141.2qB re
Pw Pi ln (29)
kh rw
2900
2800
[psia]
20
[m3/D]
10
240 hr.
0
0 100 200 300 400 500
Fig. 2.11. Solución gráfica generada usando software (Límite de Presión constante)
330
290
270
P= (Pi-Pwf)
Zona de transición:
Flujo radial a Estacionario
250
PD
1
ln tD 0.80907 (30)
2
Tenemos que:
dPD dPD
tD (31)
d ln tD dtD
dPD 1
tD , para todo tD (32)
dtD 2
1 A 2.2458
PD 2tDA ln 2 ln (33)
2 rw CA
dPD
tDA 2tDA (34)
dtDA
Otro ejemplo importante son los casos donde tenemos límites tipo
mantenimiento de presión. Es decir, de tener régimen de flujo en estado
estacionario durante una prueba, la presión adimensional viene dada por
re
PD ln y, por lo tanto, la presión derivativa es igual a cero.
rw
Ejemplo 5
Calcular el valor de la presión derivativa en estado de sistema infinito.
dPD 1
Tenemos la relación: tD
dtD 2
d P 1 141.2qB
t (35)
dt 2 kh
d P
Substituyendo valores tenemos que t 16.92 psia
dt
Limite de Sello
Limite de Mantenimiento de Presión
1000 Sistema Infinito
Límite Sellante
Cambio de Presión
100
16.89 psi
Presión derivativa
Presión constante
10
kt
rinv 0.029
c (36)
Ejemplo 7
Cálculo del radio de investigación al cabo de 100 horas de flujo.
Con el fin de enfatizar el hecho de que la EC 36 usada para el cálculo del radio
de investigación es válida sólo durante el régimen de flujo radial y que no
existen efectos de límites, presentamos a continuación un ejemplo de
ilustración.
Esto es debido a que al cabo de 100 horas, el estado de flujo es del tipo semi-
estacionario y el concepto de radio de investigación ya no aplica. La conducta
de presiones y su derivativa en función del tiempo de fluencias generadas
mediante software se muestran a continuación:
48 // GIOVANNI DA PRAT
100
10 Re = 200
Claro que la caída de presión medida en el pozo es diferente según sea la tasa
de producción, como se aprecia en la siguiente gráfica:
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 49
2800
q = 10 m3/d q = 30 m3/d
(62.8 BPD) q = 1 m3/d (188.7 BPD)
[psia]
(6.28 BPD)
2400
q = 100 m3/d
(628 BPD)
2000
20
[m3/D]
10
0
0 10 20 30 40 50 60 70 80 90 100
Fig. 2.15. Solución para la presión de pozo considerando varias tasas de producción
Las pruebas de presión, por lo general se efectúan en este medio, de aquí que
debamos generar soluciones o modelos que consideren estos factores, que
son comunes en la mayoría de los yacimientos.
50 // GIOVANNI DA PRAT
Pozo
objeto
de prueba
Infinito
Almacenamiento de pozo y daño Homogéneo
Una falla
Horizontal Doble porosidad
Circulo
Vertical Doble permeabilidad
Fallas paralelas
Hidráulicamente fracturado Radial compuesto
Fallas intersectantes
Penetración Parcial Linear compuesto
Rectángulo
re
kf
(fractura)
pozo k
kf w xf
FCD (yacimiento)
kx f
1000
Flujo lineal
(fractura)
100
Flujo radial
(formación
10
Matriz
fractura
1000
flujo sólo
en Período de transición
fracturas flujo de matriz a
100
fracturas
efecto de
límite
10
1E-3 0.01 0.1 1 10 100
zona de transición
Zona 1 Zona 2
h
(k, , c, )1 (k, , c, )2
Ri
Pozo
Fig. 2.21. Geometría sistema radial compuesto. Ejemplo: zona 1 gas, zona 2 agua.
54 // GIOVANNI DA PRAT
1000
(k/µ) zona 1
100
(k/µ) zona 2
10
100
10
1
1E-3 0.01 0.1 1 10 100
100
10
1
1E-3 0.01 0.1 1 10 100
10
0.1
Fig. 2.25. Interferencia entre Pozos. Respuesta del cambio de la presión en pozo observador.
200
100
Fallas
-100
Pozo Inyector
Pozo Productor
-200
(objeto de prueba)
10
0.1
1E-3 0.01 0.1 1 10 100
Bibliografía
59
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 61
CAPÍTULO III
ANÁLISIS DE DATOS DE PRESIÓN
3.1. Introducción
El objetivo de este capítulo es presentar los métodos básicos que permiten el
análisis de los datos de presión adquiridos durante las pruebas de pozo. Dado
que las presiones se adquieren en el pozo y que durante el proceso de
perforación y completación del mismo pueden alterarse las propiedades del
yacimiento en la zona cercana al pozo (por ejemplo, la permeabilidad efectiva),
es necesario incorporar las condiciones del pozo en las soluciones obtenidas a
la ecuación de difusividad, presentadas en el capítulo anterior.
recorrido desde el área de drenaje del yacimiento, hasta que finalmente llega al
separador de prueba.
Pc P c = Pc-Ps
Caída de presión en
reductor
Entorno completación del
pozo: perforaciones,
cemento, daño ó Presión en
estimulación (permeabilidad choke Separador
de esta zona puede ser (Ps)
diferente a la del yacimiento).
Pp
Caída de presión
en la tubería
re
Pc Pf kd Pr
Caída de presión a través Caída de presión en el
de la completación rw yacimiento
rd
3.2. Daño
En esta sección la intención no es cubrir los aspectos relacionados con el
origen del daño, ni su tratamiento en el sentido de reducirlo, sólo
mencionaremos que, durante la perforación, es posible que se alteren las
propiedades de la formación, sobre todo la permeabilidad, en la zona cercana
al pozo, lo que puede repercutir en la obtención de una tasa de producción
menor a la esperada, según el diseño o las expectativas.
Pf
Yacimiento
P (daño -)
Pf (sin daño)
Kr
P (daño +) rd (permeabilidad
del yacimiento
Pf (con daño) kd
Zona de Daño
Fig. 3.2. Idealización del concepto de daño
pD s
kh
pi pwf (1)
141.2qB
141 . 1 qB
ps s (2)
kh
k r
s 1 ln d (3)
kd rw
rwa rw e s (4)
khPi Pw
q
r
141.2B ln e S
rw
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 65
-3
100
Solución sin daño
-5
Valor de la presión derivativa
durante flujo tipo radial infinito es igual
para todos lo casos
10
Como veremos en este capítulo, el valor del factor de daño (S) se obtiene del
análisis de los datos de presión. Cabe mencionar que hay varios factores que
inducen al daño, no sólo la alteración de la permeabilidad en la zona cercana al
pozo. Por ejemplo, el hecho que el pozo atraviese la formación de forma
desviada (slant well) induce un daño negativo. Sin embargo, el valor del daño
que se obtiene del análisis de los datos de presión es el daño total, resultado
de la suma de varios daños parciales.
N
s si
i 1
El daño es uno de los parámetros que se obtiene del análisis de los datos de la
prueba y cuyo valor permite tomar acciones que pueden conducir al
mejoramiento de la productividad del pozo.
p p wf p s
EF (5)
p p wf
pD
1
ln t D 0 . 80907 (7)
2
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 67
pD S
1
ln t D 0 . 80907 2 S
2
Sin embargo esta solución considera que el flujo del yacimiento es hasta la
cara de la arena, es decir, no incluye el hecho práctico de que existe un pozo,
el cual almacena un volumen del fluido y que sirve para la conducción del
mismo hasta la superficie. Igualmente se asume que la tasa es constante a
partir del inicio y durante todo el período de fluencia.
qs
qf/qs efecto de
almacenamiento
de pozo mínimo
Zona de Efecto de o despreciable
almacenamiento (tasa de fondo
de pozo constante)
(tasa de fondo
variable)
0 Tiempo
diTiempo,
i hrs.
l
qf yacimiento
V
C (8)
P
En donde:
C = Vw c (9)
Vu
C (10)
g
144 g c
5.6146 C
CD (11)
2chrw
2
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 69
dp tD
D
tD CD (13)
d (ln )
CD
t D 60 3 . 5 s C D (14)
t
200000 12000 S C
kh (15)
C=0
C=0.01
C=0.1 Cambio de Presión
C=0.001
0
0.0001
100
pD
1
ln t D / C D 0.80907
2
Presión Derivativa
10 0.01
45 °
Tendencias
pD t D (expresión
C = 0.1 Bbl/psi CD adimensional)
1
0.01 0.1 1 10 100 1000 10000
Fig. 3.5. Efecto de almacenamiento de pozo. No hay factor de daño (S) presente
C= 0
Ci/Cf = 10
Ci/Cf = 0.1
Cf = 10 Ci
10
C=0
Cf = 0.1 Ci
1
1E-3 0.01 0.1 1 10 100
Fig. 3.6. Efecto de almacenamiento de pozo variable durante la prueba. No hay efecto de daño
(f: valor final, i: valor inicial)
Ejemplo 1
Duración del período de almacenamiento de pozo.
Solución:
Substituyendo valores en la EC. 15, tenemos que el tiempo de fluencia a partir
del cual el efecto es despreciable es de 2.54 hrs.
100 Solución
C = 0 y S =5
10
45 °
2.54
1
1E-3 0.01 0.1 1 10 10
5000
Pi = 5000 psia
(valor de la presión del yacimiento justo antes
de iniciar el período de fluencia)
4900
[psia]
4800
100
0
0 10 20 30 40 50
pD
1
ln t D 0 . 80907
2
Tal como se vio en el capítulo anterior la presión derivativa vine dada por:
dPD dPD
tD
d ln tD dtD
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 75
pD s
kh
p i p wf
141 . 2 qB
dPD 1
tD
dtD 2
0.000264kt r
tD (16) rD (17)
crw 2 rw
162.6qB k
p wf pi log t log 3.2275 0.86859 s
(18)
crw
2
kh
162.6qB k
p logt log
2
3.2275 0.86859s (19)
kh crw
162.6qB
m (20)
kh
76 // GIOVANNI DA PRAT
k
p1hr pi m log 3.2275 0.86859s
(21)
crw
2
tD
p D
C D
d P
t 0.434m
dt
Con base a esta expresión se infiere que la presión derivativa es una constante
y no depende del factor de daño durante el período de flujo radial.
La EC. 22, sugiere que si se grafica la presión fluyente en función del logaritmo
del tiempo de fluencia (gráfica comúnmente denominada tipo Semi-log),
debemos observar una tendencia tipo línea recta para la respuesta de los datos
de la presión. De aquí que, en la práctica, se puede elaborar una gráfica con
los datos de presión registrados durante la prueba de fluencia, y se pueden
determinar los parámetros de yacimiento así como el valor del daño del pozo.
Usando la EC. 20 y el valor obtenido para la pendiente m de los datos de
campo se puede calcular el valor del producto kh:
162.6qB
kh (23)
m
p pi k
s 1.1513 1hr log 3.2275
2 (24)
m crw
Zona de
almacenamiento
100
de pozo puro
Zona de
10
transición a
flujo radial
1
1E-4 1E-3 0.01 0.1 1 10
Fig. 3.9. Gráfica tipo Log-Log. Presión y presión derivativa. Ejemplo simulado prueba de
fluencia.
4950
162.6qB
pendiente ( m)
kh
4850
4750
p wf m log t p1hr
4650
4550
p 1 hr
4450
-3 -2 -1 0 1
Ejemplo 2
Los valores para los parámetros de yacimiento y pozo que se usaron a fin de
generar la solución para la presión de fluencia mediante software, y cuyos
resultados en forma gráfica se presentan en las Figs. 3.9. y 3.10. son los
siguientes:
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 79
P = 5000 psi, k = 33.3 mD, h = 30 ft, B = 1.31 BY/BN, = 10%, rw = 0.3 ft, =
0.48 cp, Factor de daño, S = 5. C = 0.01 bbl/psi, constante durante toda la
prueba. El período de fluencia tiene una duración de 24 horas y la tasa de
producción es de 500 BPD.
1 qB
Sabemos que p t , en donde el valor de p correspondiente
24 C
a t se lee sobre la línea recta que define el efecto de almacenamiento de pozo
(45°) en escala Log-Log.
De aquí que:
qB t
C
24 p
= (500)(1.31)(0.003841)
(24)(10.3066)
= 0.010 Bbl/psi
162.6qB 162.6(500)(1.31)(.48)
kh 959 mD. ft
m 53.29
kh 959
k 32mD
h 30
4489.78 5000 32
s 1.1513 log 5
3.2275
2
53.29 (.10)(.48)(1.2510 (.3)
Los valores obtenidos, tanto para la permeabilidad como para el daño, no son
exactamente iguales a los usados en el diseño de la prueba (k = 33.3 mD) y S
= 5. Sin embargo están dentro de un error menor del 10%, lo cual es aceptable
desde el punto de vista práctico, como veremos en la próxima sección. En la
obtención de los valores de estos parámetros juega, sin duda, un papel crítico
la elección de la correcta línea recta, puesto que esta define el valor de la
pendiente m y el valor de la intersección P1hr. Lo ideal es elegir la línea recta,
considerando para su trazado los puntos que estén menos afectados por efecto
de almacenamiento de pozo, en nuestro caso, a partir de las 10 horas. Para
concluir el análisis, y si se quiere, como un proceso de control de calidad, es
necesario simular toda la historia de presiones de la prueba usando los valores
obtenidos de los parámetros de yacimiento y pozo.
Pi = 5000 psia
4850
4650
4450
200
24 horas a una tasa constante e igual 500 BPD
0
0 4 8 12 16 20 24
Ejemplo 3
Definición de la correcta línea recta. Tiempo insuficiente de prueba.
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 81
En las siguientes figuras se comparan las dos soluciones (con y sin efecto de
almacenamiento de pozo y daño), a fin de ilustrar estos efectos.
C = 0, S = 0
C = 0.1 y S = 20 Solución no afectada por
almacenamiento de pozo y daño
4800
Pwf
4300
Solución considerando efectos de
almacenamiento de pozo
cambiante y daño de 20
3800
Posible línea
recta?
(Ver Fig. 3.12)
-3 -2 -1 0 1
Log t
Cabe mencionar que aunque en este caso se trata de una prueba de fluencia,
la importancia de la identificación de la correcta línea recta se presenta en las
pruebas de restauración de presión, en donde un objetivo importante es la
determinación de la presión media del área de drenaje (que se obtiene del valor
de la presión al extrapolar la línea recta hasta tiempo infinito) como veremos
más adelante.
C = 0, S = 0
C = 0.1 y S = 20
Solución
C = 0.1, Ci/Cf =0.01
1000 S = 20
Solución
C = 0, S = 0
Ver estos
100 puntos
en la Fig.
3.12)
45°
primera
10
estabilización
Valor de la derivativa
(flujo radial)
Este ejemplo sirve también para ilustrar que de ser una prueba monitoreada en
tiempo real, que significa que la presión se registra y analiza a medida que se
desarrolla la prueba (ver Capítulo VI) se extendería su duración hasta 100
horas de fluencia con base a la tendencia a esperar a partir de unas 25 horas
de fluencia para la presión en función del tiempo según inspección visual de la
Fig. 3.13., a fin de lograr la definición de la estabilización en la presión
derivativa y, por consiguiente, una mejor definición en cuanto al cálculo de la
pendiente de la línea recta.
Para finalizar esta sección vale mencionar que las técnicas de análisis de datos
de presión basadas en el uso de curvas tipo como veremos en la próxima
sección, permiten la obtención de los parámetros de yacimiento aún
considerando efectos de almacenamiento de pozo y daño que podrían ser una
limitante en la aplicación del método tipo semi-logarítmico presentado en esta
sección.
PD AP tD B t (25)
Modelo teórico
PD vs. tD
CDe2S
(curva tipo)
1000
100
PD
P 10
100 (psi)
Datos reales
1
10
1E-3 0.01 0.1 1 10 100
tD 0.1 1 10 100
t (hr)
Los valores obtenidos, y que definen el punto de ajuste o “match”, se usan para
el cálculo de los parámetros de yacimiento y pozo, así como parámetros
asociados con la geología del área de drenaje afectada al pozo durante la
prueba.
kh
A
141.2qB
0.000264k
B
crw 2
P 0.000264k t
kh 141.2qB D c
p M rw 2 tD M
kh 0.000264k
PD p tD t
141.2qB crw 2
tD kh 0.8936C 0.8936 2 S
0.000295 t CD CDe 2 S Ce
CD C chrw 2 chrw 2
86 // GIOVANNI DA PRAT
Ejercicio 4
A continuación, con el fin de ilustrar el proceso de análisis de datos mediante
ajuste por curvas tipo, se presenta un ejemplo simulado, pero que representa
situaciones de campo.
100
Régimen de flujo
radial
10
1
1E-4 1E-3 0.01 0.1 1 10
Fig. 3.15. Ejemplo de aplicación de análisis vía curva tipo. Gráfica de diagnóstico Log-Log.
El objetivo del análisis es obtener los parámetros básicos del yacimiento, como
identificar que tipo de modelo o curva tipo puede aplicarse para el ajuste de los
datos.
Una vez completada esta tarea, que puede involucrar el probar con diferentes
modelos existentes y que presenten soluciones teóricas similares en cuanto al
comportamiento de la presión y su derivativa, se procede con el ajuste. Es
decir: obtener los valores para el punto de ajuste, y con esta información
obtener los parámetros de yacimiento como el producto permeabilidad espesor,
el factor de daño y la descripción geológica del área de drenaje afectada
durante la prueba.
Análisis y diagnóstico
4980 tipo Semi-log
4960
4920
-3 -2 -1 0 1
En la Fig. 3.16. se puede apreciar la existencia de dos líneas rectas, como era
de esperarse según el diagnóstico tipo Log-log. El valor de la pendiente m1 es
cercano al doble del valor de la pendiente m2. Basados en la ecuación que
relaciona el valor de la pendiente con el producto permeabilidad espesor, (EC.
23) podemos concluir que hay una discontinuidad lateral en la permeabilidad.
Dado que en este caso m1 es mayor que m2, se deduce que la permeabilidad
lateral incrementa a más del doble a partir de cierta distancia radial (a
determinar del análisis) del pozo.
Vale mencionar en este punto que la elección del posible modelo de análisis se
basa sólo en la inspección visual de las gráficas de diagnóstico, y que, si bien
el diagnóstico ayuda en la elección del modelo, pueden igualmente existir otras
posibilidades o modelos de ajuste. En la práctica debemos de integrar
información tanto la geología sísmica como la completación del pozo a fin de
realizar interpretación de calidad y consistente con la realidad.
Estos modelos han sido publicados en la literatura y las Figs. 2.17. a 2.27. son
ejemplos de algunos de los más comunes. Los softwares comerciales disponen
de la mayoría de los modelos publicados a la fecha. Dado que nuestro ejemplo
es simulado a fines de ilustrar el proceso de ajuste, sabemos que el modelo
tipo radial compuesto es el que se adecua a la respuesta de los datos de
presión.
100
10
1
1E-4 1E-3 0.01 0.1 1 10
Fig. 3.17. Ajuste de los datos mediante modelo tipo radial compuesto
Resultados:
El software efectúa en forma automática mediante regresión el ajuste de las
curvas real y teórica y provee el punto de ajuste: para t = 1 hora, tD = 1210 y
para ∆P= 1 psia, PD = 0.118.
Con los valores obtenidos para el punto de ajuste se procede con el cálculo de
los parámetros y los resultados se muestran a continuación:
kh tD kh
PD p 0.000295 t
141.2qB CD C
p
kh 141.2qB D 141.2(188.8)(1.30689)(0.487563)(.118) 2004mD ft
p M
0.000295kh t 0.000295(2004) Bbl
C (0.0008264) 0.001
tD 0.487563 psi
CD M
Del análisis tipo semi-log podemos efectuar el cálculo del factor de daño y la
movibilidad M, puesto que se definieron con resolución las líneas rectas:
5000
4800
Pq4 (t) =?
[psia]
4600
Presión de
fluencia (psi)
4400
q4 = 900BPD
[S TB /D]
Tasa
q1 q1
q2 q3
q2
Tasa
P(t) =?
0 t1 t
Tiempo, t
(1) solución para el caso en que el pozo que produce a un tasa q1 desde
tiempo 0 hasta el tiempo t, y (2) solución que corresponde a un segundo pozo,
que produce a un tasa de (q2-q1) a partir del tiempo t1 y hasta el tiempo t, de
esta manera la tasa neta es de q1+ (q2-q1)=q2.
p i p wf
qj q log t t b '
m
' j 1
j 1 (26)
qN qN
k
b ' m ' log 3 . 2275 0 . 86859 S
cr 2 (27)
w
m’ se define como:
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 93
162.6 B
m' (28)
kh
p i p wf
qj q log t t
j 1
versus
j 1 (29)
qN qN
Los datos de presión deben de exhibir una tendencia lineal. Se asume que el
efecto de almacenamiento de pozo no causa distorsión. El producto kh y el
factor de daño se pueden determinar usando la expresión que involucra la
pendiente m’ así como la intersección de la línea recta definida por b’.
Ejemplo 5
A continuación se presenta ejemplo simulado de una prueba multitasa a fin de
ilustrar el método de cálculo usado en el análisis de los datos.
pi pwf (t ) 3 q q
m'
j j 1
logt t j 1 b'
qN j 1 qN
q
q j 1 log t t j 1 1 . 158
3
1
j
943 . 5 j 1
k
b' m ' log 3.2275 0 0.576
2
crw
162.6B
m' 0.1037
kh
4900
4700
4500
4300
N=3
500
N=2
N=1
0
0 5 10 15 20 25 30 35 40
100
10
1E-3 0.01 0.1 1 10
96 // GIOVANNI DA PRAT
Fig. 3.21. Gráfica tipo Log-Log: presión y derivativa. 3 er período de flujo (N=3)
4530 Pwf
4480
4430
4380
4330
1 3
943.5 j 1
qj qj1logt t j1
-0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.2
Fig. 3.22. Gráfica tipo semi-log. La sumatoria incluye todos los períodos de flujo.
Cabe aclarar, sin embargo, que son muchas las situaciones de campo en
donde no se realizan pruebas de restauración, debido, en parte, a que obliga al
cierre del pozo cuya producción puede estar comprometida (diferir la
producción).
5000
Período de producción
pws (t) = ?
previo al cierre
4900
[ps ia]
250 t
[S TB /D]
125
tp
0
0 10 20 30 40 50 60 70
Tasa = q
q
Período de cierre
t
0
Tasa = 0
-q
Tasa = - q
tp
Fig. 3.24. Ilustración principio de superposición
98 // GIOVANNI DA PRAT
141 . 2 qB
p i p ws ( t ) p D tp t D p D t D (30)
kh
141 . 2 qB
pi p ws t 0
k
p t D p D (31)
pD
1
ln t D 0.80907 (32)
2
162 . 6 qB
m (34)
kh
4920
162 . 6 qB
m
kh
efecto de
4820
almacenamiento
de pozo
1 2 3 4
t D 50 C D e 0 . 14 S (35)
En términos dimensionales:
0 . 14 S
170000 Ce
t (36)
kh
100 // GIOVANNI DA PRAT
Se asume que el régimen de flujo es del tipo radial infinito. Esto también quiere
decir que el yacimiento es de extensión infinita lo cual no se da en la práctica y
de aquí que se prefiere usar el termino P* en lugar de Pi como veremos más
adelante.
5000
Período de producción
(tasa variable)
4800
Período
de cierre
[ps ia]
4600
4400
[S TB /D]
500
q1 q3 q4 = 0
q2
0
0 10 20 30 40 50
N qj t t N 1 t
p ws ( t ) p i m log N (39)
qN t t t
j 1 N j
Ejemplo 6
Cálculo del valor de la presión al cabo de 30 horas de cierre.
102 // GIOVANNI DA PRAT
12 314.5 IPR
12 628.9 IPR
12 943.5 IPR
36 0 cierre
162.6(943.5)(1.3)(0.487) 3 q j 36 t 2 30
pws ( 30 ) 5000 log 4990 psia
(33.3)(30) j 1 943.5 36 t j 30
4985 psia
4900
[psia]
4700
4500
[STB/D]
500
66 hrs.
0
0 10 20 30 40 50 60 70
Como era de esperar, las soluciones son iguales desde el punto de vista
práctico.
Ejemplo 7
En este caso se trata de un pozo en donde el tipo de prueba que se realizó fue
multitasa con período de cierre final.
Período de
producción Período de cierre final
(análisis de los datos)
6200
Pi = 7000 psia
5200
3000
2000
1000
0
30/03/2005 31/03/2005 01/04/2005 02/04/2005 03/04/2005 04/04/2005 05/04/2005
1000
100
10
1E-3 0.01 0.1 1 10
6800
6400
6000
5600
5200
4800
7000
Línea recta
Pendiente (m)
6600
6200
162 . 6 qB
m
kh
k = kh/h = 9557/150 = 64 mD
6692.13 4851.5 64
s 1.1513 log 6
3.2275 4.42
2
210.9 (.16)(3.159)8.10 (0.35)
kt
rinv 0.029
c
Substituyendo valores:
(64)(72)
rinv 0.029 978 ft
(.16)(3.159)(8.10 6 )
p p wf p s
EF
p p wf
Substituyendo valores:
Una vez obtenidos los valores para los parámetros de yacimiento y pozo, que
en este caso usamos técnicas tipo semi-log, debemos simular la historia
completa de la prueba a fin de reproducir las presiones. A continuación se
presentan los resultados:
6200
5200
3000
2000
1000
0
30/03/2005 31/03/2005 01/04/2005 02/04/2005 03/04/2005 04/04/2005 05/04/2005
P* = 2933 psia
Presión actual
Régimen de flujo radial de la capa
2650 (sistema infinito) =2731 psia
Período de Estado
transición estacionario
Área de drenaje
2550 limitada o
mantenimiento
de presión
-5 -4 -3 -2 -1
Presión
Inicial
4200
[psia]
Primer cierre
Cierre Final
3200
Período de Producción
[STB/D]
500
0
0 40 80 120
Los dos cierres tienen la misma duración (24 horas) y tal que se logra
condiciones de estabilización de la presión en cada cierre. Se asume que el
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 111
4991
P actual = 4998.9psia
4986
4981
-2.5 -2 -1.5 -1 -0.5
3400
P* = 3613psia
3200
-5 -4 -3 -2 -1
tp t 0 .0002637 ktpC A
C A t pDA (42)
t p cA
Aplicando esta fórmula a nuestro ejercicio, dado que el régimen de flujo semi
estacionario se alcanzó durante el período de producción previo al cierre final,
como se muestra mediante la gráfica de diagnóstico tipo log-log aplicado al
período de fluencia previo al último cierre:
Flujo Radial
1000
100
Flujo Semi
estacionario
10
1E-3 0.01 0.1 1 10 100
Para el caso de un área tipo circular, tenemos que CA = 31.62; de aquí que
aplicando la ecuación anterior:
tp t 0 . 0002637 ( 33 . 3 )( 100 )( 31 . 62 )
C A t pDA 93 . 76
t p (. 1)( 0 . 487 )( 0 . 0000125 ) ( 393 . 6 ) 2
P = 3422 psia
3400
3350
P = 3424 psia
(Vía software)
3300
3250
3200
3150
0.5 1 1.5 2 2.5 3 3.5 4 4.5
Debimos mencionar que la duración del período de cierre podía haber sido de
tan solo 1 hora de duración y no de 24 a efectos del cálculo de la presión
media, por lo cual, y dado de que no se consideran efectos de almacenamiento
de pozo, la identificación o tendencia de línea recta está bien definida durante
la primera hora del cierre. La idea de este ejemplo fue estar prevenidos en
cuanto al significado de la presión extrapolada en una gráfica tipo semi-log, así
como presentar método simple para el cálculo de la presión actual de la capa.
Pozo
h k, Falla
(sellante)
k, h
P Y P’
X
-L L
p
Esta consideración permite que se cumpla la condición 0 , para
x x 0
todo valor de Y, en X = 0, lo cual es una condición de no-flujo.
141 .2 qB 2L
p p D (t D ,1) S p D (t D , ) (44)
kh rw
p D (t D ,1)
1
ln t D 0 .80907 (45)
2
2L 1 tD
(46)
p D (tD , ) ln 2
0 . 80907
rw 2 2L
rW
141 . 2 qB
2
1 2L
p 2 p D S ln (47)
kh 2 rw
2L
2
p wf 2 ( m log t p 1 hr ) pi m 0 . 86859 S log (48)
rw
162.6qB
en donde m es la pendiente de la línea recta en gráfica tipo
kh
semi-log y S es el factor de daño (no se consideran efectos de almacenamiento
de pozo).
De aquí que una gráfica tipo semi-log de la presión de fluencia, en función del
tiempo y de existir una falla tipo sellante que se pueda definir durante el
período de prueba (es decir, el área de investigación de la prueba es mayor
que la distancia a la falla) y en el caso de un yacimiento tipo infinito, la
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 117
Ejemplo 8
Gráficos de diagnóstico de falla sellante y determinación de la distancia a la
falla “L”.
48 816
144 0
production #1
build-up #1 (ref)
100
10
production #1
140 build-up #1 (ref)
120
100
80
60
40
-4 -3 -2 -1 0 1
4940
m = -14.7 psi
4920
4900
m = -29.3 psi
4880
t = tx = 0.34 hrs
4860
4990
k = 692 mD
S = 4.42
4970
m = -29.4 psi
4950
m =- 14.7 psi
4930
4910
(tp+t)/t = 135
2
1 2L
pD S 2 pD S ln
2 rw
Al sustituir pD por
1
ln t x D 0 .80907
2
Se obtiene que el valor de la distancia L viene dado por:
ktx
L 0 . 01217 (49)
c
(213)(0.34)
L 0.01217 63.8ft (19.2m)
(0.15)(0.6)(0.00003)
tD 1 tp t
pD 2 ln t
(50)
x
2
2L
rw
0 . 0002637 ktp
L
4 c arg p D (51)
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 121
Período
5450
multitasa
Cierre Final
Restauración de
presión
5250 (150 horas)
Cierre
inicial
5050
3500
1500
1000
cierre inicial
100
(12 horas)
10
cierre final
(150 horas)
1
1E-4 1E-3 0.01 0.1 1 10 100
Vale mencionar que la prueba se realizó con taladro por lo que los costos son
considerables; sin embargo, el valor de la información en este caso justificó la
duración de la prueba.
300
200
Primer cierre
Resultados análisis semi-log;
A este efecto y luego de revisar las respuestas que se obtienen según las
soluciones dadas por las curvas tipo (ver Capítulo II y sección anterior),
podríamos pensar en yacimiento con falla cercana y límite tipo sello o de
mantenimiento de presión. Igualmente, podríamos pensar en un yacimiento tipo
radial compuesto y límite de mantenimiento de presión o cerrado. Con esto se
quiere decir que, independiente del tipo de yacimiento, el efecto de límite es el
mismo y hay dos posibilidades: tipo volumétrico o mantenimiento de presión
(acuífero activo). A continuación se presenta el ajuste de los datos obtenidos
mediante estos dos modelos:
1000
100
10
1
1E-4 1E-3 0.01 0.1 1 10 100
Fig. 3.48. Ajuste de los datos del último cierre mediante dos modelos posibles
P i = 5500 psia
k.h = 57600 md.ft
k = 720md
Skin =4.83
R i = 115 ft
M =2.22
D = 4.82
Re =1540 ft
zona de transición
Zona 1 Zona 2
h
k, , c, k, , c,
Ri
Pozo
Fig. 3. 49. Modelo de yacimiento aplicado en el ajuste de los datos
Los resultados obtenidos nos permiten inferir que a una distancia Ri del pozo
igual 115 ft hay zona de transición que se caracteriza por un cambio en la
movilidad “M” definida como (k/) zona 1/(k/) zona 2. El valor obtenido para M
de 2.22 indica que (k/) es mayor en la zona 1 que en la zona 2 por más del
doble. Así mismo, la constante de difusividad “D” definida como (k/c) zona
1/(k/c) zona 2, cuyo valor obtenido del análisis es de 4.82, es indicativo de
que el valor de (k/c) en la zona 1 es mayor que la zona 2 (por más de 4
veces).
126 // GIOVANNI DA PRAT
El área de drenaje está confinada a una geometría tipo círculo cuyo radio es e
1540 ft.
P* = 5497 psia
5440
5400 efecto de
limite
5360
pendiente “m”
= 36 psi/ciclo
5320
5450
5350
Simulación de toda la
5250
prueba usando Modelo
de interpretación
5150
3500
1500
Como se aprecia en la Fig. 3.51. la solución obtenida del modelo usado para el
ajuste de los datos del último cierre permite reproducir la historia de presiones
durante toda la prueba.
5000
Presión de fondo Efecto del cambio de
(pozo observador) tasas en el pozo activo
4998
4996
p (máx.)
4994 = 5.34 psi
4992
0
0 40 80 120 160 200 240
1 rD 2
PD rD, tD i
2 4tD
1 tD
PD rD, tD ln 2 0.80907
2 rD
tD tD
Válida para 2
100. Sin embargo la diferencia es de sólo un 2% cuando 2
rD rD
es mayor que 5.
tD
PD rD, tD vs. 2
rD
kh p * p (r , t ) tD 0.0002637kt
pD (52) (53)
141.2qB rD
2
ct r 2
141.2qB p D 0.0002637k
k (54) c t (55)
h p M t /r 2
r 2 D D
t M
Ejemplo
Consideremos dos pozos a fin de ilustrar el efecto de interferencia (asumimos
que existe tal efecto) de manera simple. El pozo observador en este caso es un
productor, que fluye a una tasa constante de 300 BPD. El pozo activo, el cual
estaba cerrado durante 90 horas, se abre a producción por un período de 48
horas, a fin de crear el efecto de interferencia y seguidamente se cierra.
5000
Presión registrada en
el fondo del pozo
r = 1500 ft
observador (caudal es
de 300 BPD) Pozo Pozo
activo observador
4980
∆t
(∆P, ∆t) ∆P
4960 Datos usados en
el análisis
Abierto
Historia de tasas a producción
1000 (pozo activo) (48 horas)
cerrado 2000 BPD cerrado
0
0 40 80 120 160 200 240
Aplicamos análisis tipo ajuste por curva tipo durante el período de cierre. A
continuación se muestra los resultados del ajuste:
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 131
Curva tipo
tD
PD rD, tD vs 2
rD
10
Datos de presión
registrados en el
pozo observador
∆P Derivativa (curva tipo)
Punto de ajuste:
T M = 0.0639 hrs-1
PM= 0.0556 psi-1
0.1
0.1 1
∆t 10 100 100
0.0002637k 0.0002637(331)
c t 1.24 10 6 ( psi) 1
t /r
2 2
(1500) (0.4875)(0.0639)
r 2 D D
t M
Faltó mencionar que como premisa debemos de tener completados los dos
pozos en el mismo yacimiento. Cuando se opta por la prueba de interferencia,
es porque existen incertidumbres en cuanto a la conectividad hidráulica de los
pozos, aún estando completados en la misma formación.
Método Semi-log
Al igual que para el análisis de pruebas explicado en secciones anteriores, en
el caso de interferencia se pueden aplicar métodos de diagnóstico tipo semi-
log, siempre y cuando se logre el régimen de flujo radial infinito, lo que depende
del tiempo de duración de la prueba.
tD
La EC. 56 es válida siempre y cuando 2
100
rD
162.6qB k
m (57) p1hr pi m log 3.2275
(58)
ct r
2
kh
k p p1hr
ct anti log i 3.2275 (59)
r
2
m
En la EC. 59, el factor de daño no aparece porque se asume que para efectos
de este análisis el pozo observador no produce. Sin embargo y relacionado con
el efecto de almacenamiento de pozo y daño en el pozo activo, se recomienda
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 133
(y esto para efectos de diseño, ver Capítulo VI) que el tiempo de la prueba, o
expresado en forma adimensional, el factor tD/rD**2, cumpla con la relación:
0.86
tD C
230 15s D2
(60)
rD rD
Presión
Pozo activo
Pulsos
pozo
activo
1, 2, 3, número de pulsos.
tp: tiempo de duración del período de cierre del pulso.
tL1: duración entre el final del pulso y el pico de presión causado por el mismo.
134 // GIOVANNI DA PRAT
p1: magnitud del cambio de presión, causado por el pulso: distancia vertical
entre las dos líneas mostradas en la figura (la línea superior es paralela a la
inferior).
2950
Presión de fondo
2850
medida en el pozo
observador
Pozo
2750
en Pozo observador cerrado
producción (Restauración de presión)
2650
Pulsos
50
pozo 1 2 3 4 5 6
activo
0
95 105 115 125 135 145
2850
2750
Presión pozo
observador
2650 Inicio secuencia de
pulsos en pozo
activo
2550
1000 Pulsos
Pozo activo 1 2 3 4 5 6
0
70 80 90 100 110 120 130 140
Almacenamiento de
pozo Flujo zona
y Radial de
zona de transición a transición
flujo radial Flujo
(originada radial
100 por la falla) (falla)
10
1E-4 1E-3 0.01 0.1 1 10
Fig. 3.58. Gráfica tipo Log-Log. Presión y Presión derivativa vs. Log del tiempo de prueba
(35 hrs. de duración). Falla localizada a 100 ft del pozo.
100
10
1E-4 1E-3 0.01 0.1 1 10
Fig. 3.59. Simulación usando cuatro modelos diferentes. Todos los modelos convergen a la
misma solución
Similar situación tenemos con los otros modelos: la solución que considera el
modelo de doble porosidad asume que el yacimiento es naturalmente
fracturado, es decir, contentor de fracturas naturales, que no es el caso. El
modelo de doble permeabilidad asume que existe un gran contraste de
permeabilidades en el yacimiento: zonas de alta y zonas de baja permeabilidad
(tampoco es nuestro caso). Finalmente, el modelo de pozo ubicado en bloque
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 139
puede que sea compatible con la solución, puesto que uno de los lados está a
100 ft del pozo, que simularía ser la falla. Sin embargo, sabemos que la
solución no considera límites (sistema infinito).
De aquí que en nuestro caso se logró el proceso de ajuste aplicando los cuatro
modelos y que cada uno de ellos sea una posible solución desde el punto de
vista matemático y no necesariamente físico.
En donde:
140 // GIOVANNI DA PRAT
q sf (t )
q D (t D ) Siendo q ref Bq (62)
q ref
2kh
PwD (t D )
P Pwf (t )
q ref i
(63)
dp D
P' D (t D ) (64)
dt D
Pwf (t n )
m
in
q D (t i 1) q D (t i )log(t n t i ) b (65)
q D (t n ) i 1 q D (t n )
tD
1
PsD (t D ) ( ) PwD (t D )d (67) en donde (t D ) L _
1
s q (s)
0
D
4800
4400
4000
3600
1000
500
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65
horas STB/D
5 0
5 500
6 1000
5 400
5 1000
4 300
5 700
20 1200
5 500
2 (cierre) 0
100
10
0.01 0.1 1 10
Sin embargo el análisis del periodo de cierre es posible, asumiendo claro esta
de que de la tendencia que se observa para la caída de presión, es la que se
espera continúe hasta alcanzar flujo radial. Igualmente debemos asumir el tipo
de yacimiento, por ejemplo homogéneo, así como de que no existen
discontinuidades cercanas al pozo. En otras palabras podemos derivar del
análisis del periodo de cierre, de tan solo 2 horas de duración, los parámetros
de yacimiento básicos como lo son permeabilidad y daño, en el entendido de
que los resultados deben de considerase cualitativos y objeto de discusión.
100
10
0.1 1 10 100
la de alcanzar flujo radial, como primer régimen de flujo, luego de los efectos de
almacenamiento de pozo.
Esto implica de que si bien podemos realizar análisis, aunque limitado, usando
los datos de presión registrados durante el periodo corto de cierre, la respuesta
obtenida de la presión vía deconvolucion, si permite obtener los parámetros de
yacimiento y pozo, como lo son la permeabilidad y daño, con un alto grado de
certidumbre, por cuanto se tiene una buena definición de flujo radial.
A continuación se presenta a simulación de los datos de presión obtenidos vía
deconvolucion:
100
10
0.01 0.1 1 10
4800
4000
3600
500
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65
Bibliografía
4. “A new set of type curves simplifies well test analysis”. Bourdet et al.
World Oil publication (1983).
5. Well Test Analysis with Changing Wellbore Storage. P.S. Hegeman, D.L.
Hallford, and J.A. Joseph, SPEFE (September 1993) 201; Trans. AIME,
295.
10. Burns William A. “New single well test for determining vertical
permeability” J. Pet. Tech. (Junio 1969).
11. Prats Michael: “a method for determining the net vertical permeability
near a well from in situ measurements” J. Pet. Tech. (Mayo 1970).
14. SPE 71574: "Deconvolution of Well Test Data as a Nonlinear Total Least
Squares Problem", Thomas von Schroeter, Florian Hollaender, Alain C.
Gringarten, Imperial College (2001).
15. SPE 77688: "Analysis of Well Test Data From Permanent Downhole
Gauges by Deconvolution", Thomas von Schroeter, Florian Hollaender,
Alain C. Gringarten, Imperial College (2002).
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 147
18. Da Prat, G: “Well test analysis for fractured reservoir evaluation. Ch.1
Elsevier Publishing. Developments in Petroleum science, 27 (1990)
19. Van Everdingen, A.F. and Hurst W. The application of the Laplace
transformation to Flow problems in Reservoirs. Transactions AIME,
186:305-324. SPE 949305-G
Capítulo IV
149
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 151
CAPÍTULO IV
PRUEBAS DE PRESIÓN EN YACIMIENTOS DE GAS
4.1. Introducción
Aunque, como veremos en este capítulo, la metodología de análisis de pruebas
de presión transiente en pozos de gas sigue los mismos lineamientos que se
presentan en el Capítulo III, es necesario dedicar un capítulo aparte para el
caso de gas, a fin de enfatizar el proceso de linearización de la ecuación de
difusividad, ya que a diferencia del caso de petróleo no podemos considerar al
gas como fluido poco compresible sino dependiente de la presión, al igual que
la viscosidad. Otra diferencia importante en el caso de pruebas en pozos de
gas es el factor de daño, como veremos en este capítulo, ya que el daño total
puede contener una componente que es dependiente de la tasa, efecto
denominado “turbulencia” o “flujo tipo no-darcy”.
S’ = 37
3000
S’ = 23
S’ = 30
[ps ia]
2000
20
23
[M M s c f/D ]
10 q = 15 MMscf/D
8
0
0 10 20 30 40 50 60 70 80 90
La Fig. 4.2. muestra la dependencia del daño (se incrementa su valor) con el
valor de la tasa para la prueba que se muestra en la Fig. 4.1.
S’ (daño efectivo)
30
S’ = 30
26
Tasa (MMpcd)
S’ = 23
22
8 10 12 14 16 18 20 22
1 rkr p
(1)
r r r t t t
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 153
µ es independiente de la presión
ct es pequeña y constante.
M p
(2)
RT z
1 p p p
r
r r ( p ) z ( p ) r k t z ( p )
(3)
2 p 2 1 p 2 ci p 2
(4)
r 2 r r k t
154 // GIOVANNI DA PRAT
1.4 (p)
0.03 (p)
(cp
0.02
.06
cg
(1/psi)
.04
cg
.02
p
p
m( p ) 2 dp
pr
( p ) ( p ) (5)
siendo pr una presión arbitraria, aunque en general se toma el menor valor del
intervalo de variación de presiones durante la prueba de pozo, ya que solo
interesa la diferencia de seudo-presiones. Las unidades de seudo presión son
psi2/cp.
2 p 1 p c p
(6)
2
r r r k t
Podemos decir que son similares; sin embargo, la EC. 5 aún es no lineal, por lo
cual tanto la viscosidad como la compresibilidad dependen de la presión.
t
1
Ip= dt (7)
0
ct
mD pD
kh m( p i ) m( p wf ) t D (t D ) pss
válido para (8)
1424 qT
0.000264kt
tD (9)
( ct )i rw 2
pi p pi p p wf p
m( p) m( pi ) m( pwf ) dp dp dp (10)
pw f 0 0
158 // GIOVANNI DA PRAT
m( pwf ) m( p*)
1637qT
log tD 0.3513 0.87 s 0.87 Dq (11)
kh
s1 = s + D q1
s 2 = s + D q2
0.1
610 5 k s h
D (13)
rw h 2 perf
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 159
Ejemplo de cálculo
Supongamos que h = 30 ft, rw= 0.3 ft, = 0.025 cp, = 0.65, k = 33.3 md, el
valor de D usando la EC. 13, es igual a:
Ejercicio1
Compruebe con los datos de este ejemplo, que si el espesor cañoneado es
igual al 40% del valor del espesor de la formación, y que si la permeabilidad en
la zona vecina al pozo es de 5 md, debido a daño, el factor o constante de
turbulencia D es igual a 0.00092 (MSCFD)-1.
2
4.4. Rango de validez de las aproximaciones (p y p )
Si graficamos la seudo-presión en función de la presión para un rango de
presiones desde 0 hasta 6000 psi, tal como se muestra en la siguiente figura
(Fig. 4.6.), se aprecia que dependiendo del rango de presión se pueden
obtener expresiones simples para el valor de la seudo presión.
Sin embargo, en la zona de transición (rango de presión entre 1400 psia a 5000
psia) es necesario evaluar la integral de forma numérica, es decir, no hay
expresión analítica. Para presiones menores que 1400 psia, aproximadamente,
la seudo presión puede aproximarse con el cuadrado de la presión como
método de solución, tal como en el caso de un gas ideal.
160 // GIOVANNI DA PRAT
4.E+05
4.E+05
3.E+05
p
m( p) 2 ( P *i Pwf )
3.E+05
2(p/m uZ)
p
p
2.E+05 m( p) 2 dp
pr
( p) ( p)
2.E+05
P *i Pwf
2 2
1.E+05
m( p)
5.E+04
0.E+00
0 1000 2000 3000 4000 5000 6000
Presión (psia)
p 2 p 02
p
2
m( p )
Z p pdp i Z i (14)
0
p
pi
2p
Z p0
m( p ) dp p p (15)
i Zi
o
es decir, a altas presiones los pozos de gas se comportan como un fluido poco
compresible y los datos de presión pueden usarse directamente en el análisis.
kh( pr pwf )
2 2
q( MSCFD) (16)
0.472re
1424T ln( s'
rw
2kh p( p r pwf )
q ( MSCFD) (17)
0.472re
1424T ln( s '
rw
En las ECs. 16 y 17, y p son los valores promedio del producto viscosidad-
factor de desviación del gas y la presión media entre p r , que es la presión
actual del yacimiento y la pwf , respectivamente. A fines prácticos, es el
valor de calculado a la presión media p . La presión media se puede
calcular mediante la expresión:
p r p 2 wf p pwf
2
q ( MSCFD )
kh m( pr ) m( pwf )
0.472re (18)
1424T ln s '
rw
Ejemplo de cálculo:
Con la finalidad de visualizar la utilidad de la función m(p), como herramienta
de cálculo, deduciremos la tasa de producción esperada usando las tres
ecuaciones (ECs. 16, 17 y 18).
Mediante uso de software se realizó simulación (ver Capítulo VI) de una prueba
de presión típica para pozos de gas y de una duración tal que que se logren las
condiciones de régimen de flujo semi-estacionario durante el período de
fluencia previo al cierre final que es un requisita para la aplicación de las
ecuaciones.
DURACIÓN TASA
(hr) (Mscf/D)
5 0
8 5000
8 0
8 10000
8 0
8 20000
8 0
16 10000
32 0
5000
p r 4946 psia
pi 5000 psia
4800
p wf 4736 psia
4600
10000
0
0 10 20 30 40 50 60 70 80 90 100
Compruebe que:
16 9988
17 9991
18 10000
Los resultados obtenidos, sobre todo con las ECs. 16 y 17, son aproximados;
es decir, ya que debemos usar valores promedios para la viscosidad, la presión
y el factor de desviación, es posible que, dependiendo del valor obtenido,
tengamos ligeras diferencias.
dp q
2
q
(19)
dr k A A
kh( p p 2 wf )
2
q ( MSCFD )
r (20)
1424 T ln( d ) s Dq
rw
re rd
rw
rd r 0.000264 kt
ln pD (t D ) 2t D ( w ) (22) tD (23)
rw re ( ct )i rw 2
166 // GIOVANNI DA PRAT
2
1 r
pD (t D ) ln t D 0.80907 válida cuando 100 t D e
1
(24)
2 4 rw
2 2
r 3 r 1 re
pD (t D ) ln e 2t D w que es válida cuando t D (25)
rw 4 re 4 rw
Como se explicó en el Capítulo II, la solución dada por la EC. 24, representa
régimen de flujo radial en yacimiento infinito, y la solución dada por la EC. 25
representa régimen de flujo tipo semi-estacionario, que es válido para
yacimientos volumétricos o de límites sellantes.
ln t D 0.80907
rd 1
ln (26)
rw 2
rd 1
log log t D 0.1757 (27)
rw 2
La EC. 27 indica que, durante el régimen de flujo tipo radial infinito, el radio de
drenaje aumenta con el tiempo y su valor es igualmente dependiente del radio
del pozo. Por ejemplo, el valor del radio de drenaje para un tiempo tD = 100 es
de 14.98 rw, para tD = 1000, rd = 47.39 rw. Una vez alcanzado el régimen de
flujo tipo semi estacionario, el radio de drenaje viene dado por:
rd 0.472re
ln ln (28)
rw rw
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 167
P (media) = 2300 psia, rw = 0.5 ft, q = 500 MSCFD, h = 10 ft, = 0.1, = 0.7,
Tsc = 60°F, Psc = 14.67 psia, T (yacimiento) = 130°F. Re = 300 ft (sistema
cerrado). La prueba consiste de un período de fluencia de 24 horas de
duración. Asuma que no hay daño, ni efecto de turbulencia.
1E+7
45 °
1E+6
1 2 3 4 5 6
1E+5
1E-4 1E-3 0.01 0.1 1 10
a) Demuestre que el radio de drenaje (rd) durante los tiempos elegidos (puntos)
viene dado por los siguientes valores:
Ayuda…
La relación tasa de gas, y cambio de presión, viene dada por la EC. 21.
Substituyendo los valores para las condiciones Standard y asumiendo s = 0 y D
= 0 (esto sólo para simplificar el ejercicio, por cuanto sabemos que el factor de
turbulencia es necesario de considerar tratándose de gas) tenemos la siguiente
relación:
rd
2.041 10 7 m( p ) ln
rw
rd
1.5 t D
rw
1.02 10 4 m( p)
q
rd
ln
rw
q 1.8 10 5 m( p )
d) Efectúe gráfica tipo Log-Log del radio de drenaje, en función del tiempo. Se
deben obtener los siguientes resultados:
100
rd (ft)
10
1
0.001 0.01 0.1 1 10 100
tiempo (hrs)
Fig. 4.10. Gráfica tipo Log-Log. Radio de drenaje vs. tiempo
e) Repita el mismo ejercicio usando una tasa de 750 MSCFD. Compruebe que
se obtiene el mismo resultado en cuanto a los valores obtenidos para el radio
de drenaje.
m( pwf ) m( p*)
1637 qT
log t D 0.3513 0.87 s 0.87 Dq
kh
m( pi ) m( p1hr ) k
s´ s Dq 1.151 log 3.23 (29)
b ( ct )i r w
2
170 // GIOVANNI DA PRAT
qT
b 1637 (30)
kh
2000
Datos de yacimiento y
fluido:
rw = 0.5 ft
Presión (psia) h = 10 ft
=10%
1000
P = 2300 psia
= 0.7
T = 130°F+ 460 =
590°R
3000
Tiempo (hrs.)
production #1
production #2 (ref) Segunda fluencia
Primera fluencia
1E+8
1E+7
Flujo Radial
Permite el cálculo de el daño efectivo s’ así
como el producto kh, durante cada período de
flujo, mediante ajuste por curva tipo o análisis
tipo semi-log (pendiente claramente definida)
1E+6
1E-4 1E-3 0.01 0.1 1
4E+8
production #1
production #2 (ref)
Segunda fluencia
3E+8
2E+8
Primera fluencia
1E+8
s’ = 1.37
s' s Dq
1.2
s’ (para q = 0) = s = 0.03
1
0.8 ds'
D 0.000429(MSCFD)1
dq
s’= 0.685
1600 1800 2000 2200 2400 2600 2800 3000 3200
Skin vs Rate
Fig. 4.14. Daño efectivo, s’, en función de tasa
q( MSCFD)
kh m( pr ) m( pwf )
0.472re
1424T ln s'
rw
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 173
q( MSCFD)
48 m( pr ) m( pwf )
0.472re
1424(590) ln 0.03 (0.000429)q
.5
Si asumimos un radio externo re = 2000 ft, el tiempo mínimo a partir del cual se
obtienen condiciones de flujo estabilizado se puede calcular mediante las
relaciones:
0.000264kt
0.472re 1.5r w t D tD
ctr 2 w
Substituyendo valores:
2
0.472( 2000) (.10)(0.0183)(.000439)(0.5) 2
t 249.5 ( hrs ) 10.4 ( dias )
(1.5)(.5) (0.000264)(4.83)
q( MSCFD)
48 m( pr ) m( pwf )
0.472(2000)
1424(590) ln 0.03 (0.000429)q
.5
q ( MSCFD )
5.71 10 5 m( p ) m( pwf )
7.57 4.3 10 q
4
1E+8
1000
C and N - m(p) Vertical Well IPR (Bottom Hole Pressures - Flow after flow) Plot - Test Design 6: m(Pavg)-m(Pf) [psi2/cp] vs Q [Mscf/D]
q 1.32 104 4.15 108 m( pwf ) 0.8586
q 0.001746 (2300) 2 p 2 wf 0.9339
2000
1500
1000
500
4.5.2.1. P* y P media
A fin de introducir el concepto de P* y presión media para gas, consideraremos
el caso en que la duración de la prueba permite identificar el límite del
yacimiento y que sea del tipo sellante. En la siguiente gráfica se presenta la
historia de la prueba, la cual consiste en dos períodos de fluencia y un cierre
final.
2300
2200
período
2100
de
2000
cierre
1900
500
0
0 10 20 30 40 50 60 70 80 90
1E+8
45 °
1E+7
1E+8
1E+7
1E+6
1E-4 1E-3 0.01 0.1 1 10
i ct i A
t pss (t DA ) pss
0.000264k
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 177
(0.1)(0.0183)(0.000439)( )(300) 2
tpss (0.1) 17.8 horas
0.000264(4.83)
3.9E+8
3.7E+8
3.5E+8
3.3E+8
P* = 2269.86 psia
3.1E+8
m(p*)= 4.06 10 (exp 8) psi**2/cp
-5 -4 -3 -2 -1
qT ( ct ) i re 2
m( p ws ) m( p ) 1637 log (31)
kh 0.00266kt
178 // GIOVANNI DA PRAT
( ct )i re 2
t (32) se tiene que m( pws ) m( p )
0.00266 k
(0.1)(0.0183)(0.000439)(300) 2
t 5.66h 5.66
(0.00266)(4.8) hrs.
-4 -3 -2 -1 0 1
2400
1400
Presión de fondo (psia)
2000
Tasa (Mpcnd)
Rate [Mscf/D]
1000
0 10 20 30 40 50 60 70
Pi = 2500 psia, rw = 0.3 ft, porosidad = 15%, T res = 135°F, = 0.67, So (daño
tipo mecánico = 0), D = 0.001 (MSCFD)-1.
production #3
production #4
production #5
production #6 (ref)
1E+8
1E+7
1E+6
1E-4 1E-3 0.01 0.1 1 1
Podemos efectuar una gráfica del daño efectivo obtenido del análisis semi-log
de cada período en función de las tasas registradas durante la prueba.
production #3
production #4
production #5
3.5E+8 production #6 (ref)
2.5E+8
1.5E+8
5E+7
Los resultados (relevantes) obtenidos del análisis para cada período de flujo
son:
Skin=2.45 Skin=1.46
Skin=0.958 Skin=1.95
1.6
Resultados:
1.2
D = dS/dQ=9.9168E-4[Mscf/D]**-1
Skin vs Rate
Fig. 4.25. Variación del daño con la tasa
1E+8
1E+7
1E+6
1E-4 1E-3 0.01 0.1 1 10
1E+8
1E+7
1E+6
1E-4 1E-3 0.01 0.1 1 10
2400
Solución
daño variable
1900
1400
Solución
Daño variable (debido a efecto de turbulencia)
daño constante
Daño constante (no se considera turbulencia
2000
1000
0
0 10 20 30 40 50 60 70
2.4
1.6
Solución correcta
Daño efectivo: variable
s’ ~ 0 + (0.001) q
production #6
1.2 production #3
production #4
production #5
build-up #2 (ref)
Skin vs Rate
q C( p p 2 f )n
2
(32)
1E+8
n=1
1E+7
1E+6
AOFP
(pf = 14.7 psi)
1E+5
C and N Vertical Well IPR (Bottom Hole Pressures - Modified isochronal) Plot - Test Design: Ps²-Pf² [[psi]**2] vs Q [Mscf/D]
Por esta razón, si se usa este tipo de curva, se deben efectuar pruebas
periódicamente en el pozo, para reajustar el valor de C y n. De aquí que, como
veremos en la próxima sección, se prefieran las pruebas tipo isocronales.
4.6.1. Pruebas tipo flujo tras flujo (Flow After flow test)
Esta prueba es muy común (tanto en gas como en petróleo: prueba multi-tasa)
y se aplica con el fin de obtener un valor para el índice de productividad del
pozo (IPR), el cual puede ser de carácter transiente o estabilizado.
186 // GIOVANNI DA PRAT
2300
Pozo cerrado
Pf1
Presión = ACTUAL
2100 Pf 2
1
1900
2 Pf3
3
1700
Pf4
q3 4
q2 q4
1000
q1
0
12 16 20 24 28 32 36 40 44 48
Se registra, mediante sensor de fondo, las presiones de fluencia así como las
tasas correspondientes, mediante separador en superficie. La tasa debe ser
constante durante todo el período de flujo particular.
1
2
2000
Pf 4
1500
(psia)
1000
500
q (Mscf/D)
Fig. 4.32. Gráfica tipo IPR. Prueba tipo flujo tras flujo
C = 7.92743E-4[Mscf/D]/[psi]**2N
N = 0.99
q 7.92743 10 4 ( 2300) 2 pwf
2 0.99
production #3
production #4
production #5
production #2 (ref)
1E+8
1E+7
1E+6
Fig. 4.33. Gráfica de diagnóstico tipo Log-Log (todos los períodos de flujo)
2300
Ps Ps Ps
2200
Pf1
2100 Pf2
2000
Pf3
Prueba isocronal
1000
0
0 100 200 300 400 500
asumiendo la misma duración para los cierres que en el caso anterior, se logra
alcanzar la estabilización de la presión.
Nuevamente vale decir que en la práctica puede que no se logre este objetivo
durante el período de flujo extendido, sobre todo si la medición de presiones no
se realiza en la modalidad tiempo real, que permite verificar el régimen de flujo
existente mediante el análisis de los datos sobre la marcha de la prueba.
Radio de drenaje
efectivo
Radio de drenaje
aparente
Este desarrollo permitió ofrecer una alternativa a las curvas tipo contra-presión,
debido a que las mismas pueden presentar cambios en las constantes C y n. A
continuación se presentan las gráficas de comparación de la prueba isocronal y
tipo contra-presión efectuada en un mismo pozo:
Fig. 4.36. Pruebas tipo isocronal. Datos de campo (n = 0.867 para todas las curvas).
Fig. 4.37. Prueba tipo contra-presión efectuado en el mismo pozo (Fig. 4.36.)
Radio de drenaje
aparente (Ra)
Radio de drenaje
aparente
(valor igual al primer flujo)
Radio de drenaje
Aparente
(valor igual al primero y
segundo flujo)
Presión de fluencia al
finalizar el tercer
período de flujo
Ps1 Ps2
Ps3
2200
2000
Pf1
1800
Pf2
Presión de fondo
(psia) Pf3
1600
Pext
Tasas
(Mscf/D) 3 Flujo extendido
1000
2
1
0
0 10 20 30 40 50 60 70 80 90 100
1E+9
AOFP
1E+7
Log q (Mscf/D)
P(media)=2232.8 psia
AOFP=3732.4Mscf/D
C (trans.)=7.96635E-6[Mscf/D]/[psi2/cp]**n
n =1.0
isocronas =3
C (ext.) = 6.9492E-6[Mscf/D]/[psi2/cp]**n
P ext.=1467.5 psia
q ext=2000 Mscf/D
1E+8
45°
1E+7
1E+6
1E-4 1E-3 0.01 0.1 1 10
q 6.9492 10 6 3.947 108 m( pwf ) 1
Contra-presión
Isocronales e isocronales Contra-presión e
Pruebas válidas Flujo tras flujo e
modificadas isocronales
isocronales.
Coeficiente C:
prueba contra- Disminuye con la duración
Constante Constante
presión ó del período de flujo
isocronal
3000
2900
2800
10
0
0 40 80 120 160 200 240
1E+5
10000
1000
1E-3 0.01 0.1 1 10 100
1E+5
10000
1000
1E-3 0.01 0.1 1 10
De aquí que el análisis de los datos, en este caso, esté orientado a definir los
parámetros básicos de yacimiento y límite del área de drenaje así como la
obtención de la relación de disponibilidad del gas.
1E+5
10000
1000
1E+5
10000
1000
1E-3 0.01 0.1 1 10 100
Fig. 4.48. Ajuste de datos para el período de cierre final, usando modelo de yacimiento
homogéneo continuo, y de área de drenaje tipo barra
3000
2900
2800
10
0
0 40 80 120 160 200 240
800
N
400
0 O pozo E
-400
S
-800
1E+7
recta
(flujo extendido)
1E+6
(Isócronas de 8
horas)
1E+5
1 10 100
C and N Vertical Well IPR (Bottom Hole Pressures - Modified isochronal) Plot - Test Design 3: Ps²-Pf² [[psi]**2] vs Q [MMscf/D]
Resultados:
AOFP = 56.9617MMscf/D
C (trans.) =1.19838E-4[MMscf/D]/[psi]**2N
n=0.846
C (ext.) =7.41908E-5[MMscf/D]/[psi]**2N
P ext=2821.11 psia
q ext.=9.18181MMscf/D
En forma de ecuación:
q 1.2 10 4 pavg 2 p wf
2 0.84
El valor obtenido para n (0.846) refleja de que hay efecto de turbulencia, tal
como se podía anticipar, del análisis de los períodos de fluencia y cierre, en el
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 205
cual se obtuvo un daño por efecto de turbulencia con un valor para D de 0.283
(MMSCFD)-1.
Bibliografía
9. Railroad Commission of Texas. Back pressure tests for natural gas wells
(Abril, 1972).
PRUEBAS EXTENDIDAS
207
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 209
CAPÍTULO V
PRUEBAS EXTENDIDAS
5.1. Introducción
Las pruebas extendidas en su concepto más básico consisten en una serie de
aperturas y cierres controlados del pozo a realizar durante un período de
tiempo que depende del alcance de los objetivos de evaluación, aunque
basados en la experiencia, el período de tiempo total puede variar entre días y
meses. Durante la ejecución de la prueba se registra en forma continua tanto la
presión como otros parámetros de producción, por ejemplo, la tasa de
producción, tipo de fluidos que se producen, corte de agua, relación gas-
petróleo, entre otros.
El análisis de los datos de presión, como se detalló en los Capítulos III y IV,
permite determinar el área de drenaje dinámica y confirmar límites, acuíferos
activos y sellos a una escala que puede ser estructural. Igual de importante es
la determinación del daño y el seguimiento de su valor durante el periodo de
prueba.
cierre
Período de fluencia
Inicio del
efecto de límite
45°
Flujo radial
Fig. 5.2. Diagnóstico del cambio de la presión y su derivativa durante el período de fluencia
Hay casos en los que durante una prueba de restauración de presión, el tiempo
necesario para lograr condiciones de estabilización puede ser extremadamente
grande, dado que la presión muestra un incremento que, aunque pequeño o
212 // GIOVANNI DA PRAT
La inspección visual de las Figs. 5.3. y 5.4. indica que la presión de fondo no
presenta declinación con respecto a su valor inicial durante el período de cierre
final, como consecuencia del acumulado de producción. Igualmente, es
importante aclarar que de ser el sistema tipo volumétrico se logra alcanzar una
estabilización de la presión durante el periodo de cierre, pero la diferencia está
en que el valor de la presión estabilizada es menor que la presión inicial, como
se explicó en el caso A.
En este caso, la prueba extendida sirve para corroborar los resultados de las
pruebas iniciales, pero abarcando un radio de investigación mayor (escala del
área de drenaje).
5000 Pi P1 P2
P3 P5
P4
4000
3000
puntos
de control
0
0 40 80 120 160 200 240 280
Ejemplo 1
En la Fig. 5.6. se presentan los resultados en forma gráfica de una prueba
extendida, realizada en un pozo ubicado en un yacimiento o área de drenaje
contentora de gas-seco. La motivación a realizar la prueba extendida se debe a
que, durante las pruebas iniciales, se presentó tendencia de declinación de la
presión, lo que podría indicar que se trata de un área de drenaje limitada, por lo
que se sugirió realizar una prueba extendida a fin de probar la comercialidad
del pozo o el área de drenaje asociada al mismo.
P = 7296.7psia
8500 Gas Acumulado
1566 MMpcd
Z = 1.18974
7500 P(inicial)
9000 psia
Gas Acumulado P = 8653.6 psia
0 MMpcd Gas Acumulado
Z = 1.32165 290 MMpcd
6500
Z = 1.29434
40
20
0
06/10/2004 26/10/2004 15/11/2004 05/12/2004 25/12/2004
build-up #3
production #3
production #4 (ref) Períodos
build-up #4
de fluencia
1E+9
1E+8
Períodos
de cierre
1E+7
1E-4 1E-3 0.01 0.1 1 10
Período de transición a
régimen de flujo semi-
estacionario
1E+9
45°
1E+8 Flujo
Semi
estacionario
Flujo radial
6.6 días
1E+7
1E-4 1E-3 0.01 0.1 1 10
Fig. 5.8. Gráfica tipo Log-Log Período de flujo previo al cierre final
Como se aprecia en la Fig. 5.8., luego del período de régimen de flujo radial, la
presión derivativa muestra tendencia al aumento y a partir de 6.6 días
(inspección visual) los datos de presión presentan tendencia de caída tipo lineal
con el tiempo de prueba, es decir, dp/dt es constante, lo cual indica que se
alcanzó el régimen de flujo semi-estacionario, condición que se logra sólo de
ser el área de drenaje tipo volumétrica.
El cálculo de las reservas se puede obtener, ya sea, por análisis tipo cartesiano
aplicado al período de flujo semi-estacionario, por cuanto se define línea recta
cuya pendiente está relacionada con el volumen de drenaje, o por análisis
mediante ajuste por curva tipo, como se explicó en capítulos anteriores.
1E+9
Modelo de Ajuste
1E+8
1E+7
1E-4 1E-3 0.01 0.1 1 10
800
Volumen de reserva
400 15054 MMpcn 264m
0 588 m 62m
202m
-400
-800
-2400 -2000 -1600 -1200 -800 -400 0 400
7000
6000
5000
P/Z
4000
3000
2000
1000
0
0 5000 10000 15000 20000
223
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 225
CAPÍTULO VI
DISEÑO E IMPLEMENTACION DE LAS PRUEBAS
6.1. Introducción
El diseño de las pruebas de presión es una parte importante del proceso de
evaluación del pozo y se puede decir que está muy vinculada con el éxito de la
prueba, en cuanto al logro de los objetivos de evaluación propuestos.
Por ejemplo, la prueba tipo DST (Drill Steam Test) es muy común en pozos
exploratorios, y se toma ventaja de que el taladro está en sitio, lo cual provee
una logística de seguridad óptima. Sin embargo la prueba puede ser de corta
duración (lo que se traduce a que puede que no se investigue toda el área de
drenaje), por lo cual, y aunque el costo de la prueba per se no sea sino el 30%
del costo total, este último puede ser considerable, como es el caso de los
pozos costa fuera (promedio 250,000US$/día). De aquí que los resultados en
cuanto al potencial del pozo se tornen importantes y se sugiera prueba de
mayor duración una vez completado el pozo.
Presión de fluencia
4330
no consistente con
tasa de producción
4310
Respuesta de la presión
(tendencia) no es consistente
4290 durante un período de cierre
(fugas, comunicación o mal
funcionamiento del sensor de
4270 presión?)
40
20
0
168 172 176 180 184 188 192 196 200
1000
100
10
0.01 0.1 1 10 100
Fig. 6.2. Diagnóstico tipo Log-Log. Período de cierre. Sensor afectado por ruido.
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 227
3000
Pressure [psia]
2900
Rate [Mscf/D]
10000
0 10 20 30 40 50
Fig. 6.3. Prueba en pozo de gas. Efecto de comunicación de 2 intervalos vía mala
cementación.
Diseño de la Prueba
Duración prueba y procedimiento
de toma de muestra
Puede que el pozo atraviese varios intervalos, como por ejemplo en la cuenca
del golfo San Jorge, Argentina, donde los intervalos pueden ser contentores de
fluido diferentes (agua, petróleo o gas). Es decir que el diseño debe hacerse
considerando estos aspectos. El resultado final, de la evaluación del intervalo o
capa mediante la prueba debe ser conclusivo; es decir, mediante la prueba y
sus resultados debemos concluir si el intervalo o la zona de interés es
comercial o no, y, de presentarse dudas, sugerir prueba extendida (ver Capítulo
V) antes de comprometer la producción del mismo.
Intervalo
a
probar
pozo
0 1
Km.
Presión: 4200 psia Espesor bruto: = 74 ft, Espesor neto (h) = 47 ft, Porosidad =
18%, Sw = 35%, Sg = 65%, So = 0%, T = 196°F, = 0.57, rw = 0.34 ft
Profundidad = 8800 ft.
4”
8 ½”
En cuanto a los beneficios de probar con taladro en sitio, como en este caso, la
sarta tipo DST tiene una ventaje inherente: el cierre es de fondo. De aquí que la
duración de los períodos de cierre y los efectos de segregación de fases se
minimizan. Como veremos en el caso de un pozo de desarrollo, y de no
disponer de cierre de fondo, debemos estimar el período de duración del efecto
de almacenamiento de pozo, de manera tal que permita el diagnóstico del
régimen de flujo radial.
q(MSCFD)
kh m( pi) ) m( pw f ) k
1
Dado que tenemos un estimado para la presión inicial, así como la gravedad
especifica del gas, se puede construir la función m(p), mediante el método visto
en el Capítulo IV.
q ( MSCFD )
(50)(47) 1.024 10 9 m( p w f ) 50
1
log t log 2
3.23
1638(656) (0.18)(0.02451)(0.000124 )(0.34)
q( MSCFD)
0.00219 1.024 109 m( pwf )
log t 5.67
q ( MSCFD ) 3.10 10 4 1.024 10 9 m ( pwf )
A continuación se presentan resultados asumiendo varios valores para la
presión de fondo fluyente.
234 // GIOVANNI DA PRAT
Tabla 6.1.
4000
Presión, psia
3000
2000
AOF
1000
0
0 50 100 150 200 250 300 350
Tasa (MMscfD)
Fig. 6.8. IPR (transiente)
f f T q
2
0.0375 g sen L
s
T (3)
gq
N Re 20.09
D (4)
En nuestro caso podemos asumir varios valores para la presión en cabezal, por
ejemplo: 2000, 2500, 3000, y 3500 psia. Par cada uno de estos valores
calculamos la presión de fondo (en este caso a 8800 ft (ver Fig. 6.7.). A
continuación se presentan los resultados obtenidos aplicando la EC. 2:
236 // GIOVANNI DA PRAT
4500
4000
3000
3500
2500
3000
Presión, psia
2000
2500
IPR (Yacimiento)
2000 1500
1500 1000
1000
P(cabezal)
500
31 MMscfd 45 MMscfd
0
0 20 40 60 80 100 120 140
Tasas, (MMscfD)
( 200000 12000 s )C
Período de Fluencia t kh (5)
170000 C e0.14 s
Período de Cierre: t (6)
kh
(200000)(0.024)
t 0.05 hrs
2350
.02451
Aunque del análisis petrofísico, así como del probador de formación, se obtiene
información valiosa para el diseño, la mayor incertidumbre en el caso de pozos
exploratorios está en que no conocemos a priori o antes de realizar la prueba la
descripción del área dinámica de drenaje, es decir, si se trata de un yacimiento
homogéneo con discontinuidades o naturalmente fracturado, o capa de
permeabilidad variable, por mencionar algunos ejemplos. De no disponer de
ninguna información de geología, lo aconsejable es efectuar diseño
considerando un yacimiento homogéneo.
4150
4050
3950
20000
0
0 10 20 30 40 50
production #4
build-up #4 (ref)
1E+8
1E+7
1E+6
1E-4 1E-3 0.01 0.1 1 10
production #4
build-up #4 (ref)
1.1E+8
9E+7
Período de restauración
7E+7
5E+7
Período de fluencia
3E+7
1E+7
-4 -3 -2 -1 0
rd 1
log log t D 0.1757 (7)
rw 2
0.000264kt
en donde tD (8)
( ct )i rw 2
rd 1 0.000264(50)(6)
log log 0.1757
rw 2 (0.18)(0.0245)(0.00012397)(0.34) 2
de aquí que:
rd 570 ft
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 241
Cabe mencionar que, tratándose de gas, por lo general son dos los objetivos
que deben complementarse: obtener la IPR del pozo mediante la prueba de
producción (o comúnmente denominada contra-presión) y extraer información
de yacimiento, en cuanto a caracterizar el área de drenaje del pozo que se
logra o deriva del análisis de los datos de presión.
En este caso, por ejemplo, según el modelo de sísmica hay una falla cercana al
pozo, que dista cerca de los 200 m (656 ft) del mismo. Tal como está propuesto
el diseño, al finalizar el período de cierre, se investigó un área tipo circular de
348 m (1141 ft) de radio, con lo cual, de existir falla localizada en este radio, se
verificaría.
1E+8
2
hrs
sin falla
1E+7
Efecto de falla
1E+6
1E-4 1E-3 0.01 0.1 1 10
1E+8
Efecto falla
sellante (puro)
1E+7
Zona de transición
1E+6
1E-4 1E-3 0.01 0.1 1 10 100
Ejercicio 1
Elabore diseño de prueba usando el mismo procedimiento, pero considerando
que la permeabilidad del yacimiento es de 375 md, y que el pozo presenta un
daño total (no hay efecto de turbulencia) igual a 21. Asuma una tubería de
producción nominal de 3 ½ pulgadas. Las propiedades del fluido así como la
presión son las mismas. Estime el tiempo de prueba necesario o mínimo que
garantice la verificación de la existencia de falla sellante ubicada a 200 m (656
ft) radialmente al pozo.
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 243
Las pruebas pueden ser rutinarias, es decir, formar parte de una campaña de
pruebas (de frecuencia anual, por ejemplo), a fin de dar seguimiento a un
proceso de recuperación secundaria (ensayos tipo fall-off). O a ser realizadas
en el caso de presentarse un comportamiento anómalo en la conducta de
producción o presión del pozo: declinación inesperada, por ejemplo. En otros
casos se quiere corroborar la continuidad lateral entre dos o más pozos, por lo
que proponen pruebas tipo interferencia.
En la mayoría de los casos, tanto los parámetros básicos del yacimiento y pozo
(presión, permeabilidad efectiva, daño del pozo, por ejemplo) como la
composición del fluido son conocidos, y el pozo puede que ya tenga un historial
de pruebas. De aquí que los objetivos de evaluación sean concretos y las
alternativas en cuanto al tipo de pruebas y el tiempo de duración de las mismas
depende del objetivo de evaluación. El procedimiento empleado en el diseño de
las pruebas es similar al caso del pozo exploratorio, aunque se puede definir
con mayor resolución el tiempo de prueba por cuanto los parámetros de
yacimiento son conocidos.
A fin de ilustrar este punto, aplicamos el concepto usando datos que son
representativos de casos de campo.
C Vw c f
244 // GIOVANNI DA PRAT
C=0.001,S = 0
C=0.01,S=0
C=0.01, S=10
C=1.92E-4, S = 0.5 (Ensayo DST)
C=0.0011, S= 5 (Cierre en Superficie)
C=0 S=0
1000
C=0.01 y S= 10
Efecto de
falla?
100
Cierre en
Cierre de fondo Superficie
C (variable) S=5
Flujo radial C=0 y S=0 S = 0.5
(inicial)
10
1E-3 0.01 0.1 1 10
Por ejemplo, en el caso de que el valor de C sea igual a 0.01 bbl/psi y el factor
de daño igual a 10, aunque se tenga una línea recta semi-logarítmica, a partir
de 1 hora del ensayo la misma no corresponde al régimen de flujo radial infinito
inicial y de aquí que el cálculo de la permeabilidad resulte incorrecto (valor
menor). Esto afecta igualmente el diseño de prueba tipo larga duración, puesto
que se usarían los valores de los parámetros de yacimiento obtenidos de la
prueba para el diseño, y por ello el tiempo de prueba determinado según
diseño estaría en error.
DISEÑO
Sísmica Diseño
Y Geología
EJECUCIÓN
Completación
y Petrofísica, Núcleos
perforación registros
Pc
SISTEMA DE
CIRCULACIÓN y CONTROL
SISTEMA DE
CIERRE EN FONDO
Ph SISTEMA DE ADQUISICION
(DATOS DE PRESION)
EMPACADURA
yacimiento
Pf
Pc < P f < Ph
Inicio de la prueba
Una vez posicionada la sarta de prueba en profundidad y asentada la
empacadura, se procede a abrir la válvula de circulación para desplazar el
248 // GIOVANNI DA PRAT
Válvula de
circulación
Válvula
de prueba
Sensores
de presión
YACIMIENTO
Final de la prueba
Durante los períodos de flujo, mediante los equipos de superficie (separadores,
tanques de calibración, toma muestras), es posible medir las tasas de
producción, presión y temperatura en cabezal, así como la toma de muestras
de fluido.
valvula de
circulación
(Abierta)
valvula de
prueba
(cerrada)
Sensores
Ph de
presión
empacadura
Pf
Presión (anular)
Pressure [p sia]
8000
4000
0
Presión en superficie (14.7 psia)
0 50 100 150 200 250 300 350 400
A B C D
10
Sin válvula de
cierre de fondo
1
1E-4 1E-3 0.01 0.1 1 10
Cable eléctrico
Tubería de
producción
El sensor de presión transmite
información vía cable eléctrico
anexado a la tubería.
Sensor de presión
capa A
capa B
8 1/2” 7” @ 17690’
Separador trifásico
Tanque de
(Ps)
Calibración
Choke Manifold
Reductores
petróleo para prueba
P=Pc- Ps
Presión en
agua Cabezal (Pc)
gas
Quemador
Datos de Presión
El fluido proveniente del pozo pasa primero a través del choke manifold, el cual
permite el uso de varios reductores (choke) según sea la prueba. A la salida del
manifold el fluido pasa al separador, en donde se efectúan las mediciones de
tasas (gas, petróleo y agua, así como la toma de la muestra).
254 // GIOVANNI DA PRAT
1000
100
10
0.1
1E-5 1E-4 1E-3 0.01 0.1 1 10 100 100
Fig. 6.25. Sensor tipo cristal de cuarzo. Resolución 0.02 psia, drift = 0, ruido = 0
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 255
100
10
0.1
1E-4 1E-3 0.01 0.1 1 10 100
Fig. 6.26. Sensor: tipo strain, resolución = 0.5 psi, drift = 1.E-7/hr, ruido = 0
1000
100
10
0.1
1E-5 1E-4 1E-3 0.01 0.1 1 10 100 100
Fig. 6.27. Sensor tipo strain: resolución = 0.5 psi, drift = 1.E-7 /hr, ruido = 0.5 psia
256 // GIOVANNI DA PRAT
1000
100
10
0.1
1E-5 1E-4 1E-3 0.01 0.1 1 10 100 1000
Fig. 6.28. Sensor tipo cuarzo. Resolución = 0.02, Drift = 1.E-5 hr, ruido = 0.
1000
100
10
0.1
1E-5 1E-4 1E-3 0.01 0.1 1 10 100 1
Fig. 6.29. Sensor tipo mecánico. Resolución = 5 psi, drift = 1.E-7/hr, ruido = 0.
presión, según las especificaciones del sensor. Dado que son ejemplos
simulados, conocemos cuál es la respuesta según el parámetro o
especificación del sensor. Sin embargo, en la práctica puede que no sea tan
simple de inferir y, por lo general, se tiene la tendencia (no sugerida) de
“suavizar” (como se verá más adelante) los puntos de presión, con la finalidad
de minimizar el efecto de baja resolución del sensor o ruido.
P2
t2
P1
t1
X-L X X+L
258 // GIOVANNI DA PRAT
p1 p
t2 2 t1
t t2
p ' 1
t1 t2
10
0.1
1E-4 1E-3 0.01 0.1 1 10
Fig. 6.31. Gráfica de diagnóstico. Sensor de resolución = 0.3 psia, tipo strain, ruido de 0.05
psia y 1.E-7 hr-1 de drift. Factor de smoothing (L) = 0
baja resolución del sensor, ruido y drift. En estos casos es aconsejable usar
otras gráficas de diagnóstico, como la semi-log y cartesiana, a fin de
diagnosticar el régimen de flujo. Es decir, la derivativa es muy sensible y, de no
estar suavizada, presenta la tendencia que se muestra.
5000
4980
4960
4940
4920
-3 -2 -1 0 1
5000
4980
4960
4940
4920
0 5 10 15 20 25 30 35 40 45
Fig. 6.33. Gráfico tipo cartesiano: tendencia tipo línea recta de los datos de presión
260 // GIOVANNI DA PRAT
Como se aprecia de las Figs. 6.32. y 6.33., el régimen de flujo tipo semi-
estacionario es más resolutivo de identificar. De hecho, el cálculo del volumen
de drenaje se obtiene de la pendiente de la Fig. 6.33., y se corresponde con el
valor usado en la simulación. En conclusión, debemos tener presente que de
no tener un algoritmo de diferenciación para el cálculo de la derivativa que
considere smoothing, tenemos que recurrir a otras gráficas especializadas a fin
de identificar el régimen de flujo.
10
0.1
1E-4 1E-3 0.01 0.1 1 10
Fig. 6.34. Presión derivativa usando factor de smoothing L de 0.5 y Modelo usado en la
generación del ejemplo
Características A B
Diámetro 1,5” 3/4”
Longitud 48” 16”
Batería Litio 300 ºF Litio 300 ºF
Frecuencia de muestreo 1sec mínimo 3 sec mínimo
Presión (máx.) 10.000 psi 15.000 psi
Temperatura (máx.) 302°F 302 ºF
Exactitud 0,02% Escala Total 0,022% Escala Total
Resolución 0,01 psi 0,06 psi
Capacidad memoria 102.400 set de datos 260.000 set de datos
5500
Pressure [psia]
5300
5100
Efecto de valor incorrecto de la
temperatura de yacimiento usado
en el cálculo de la presión.
4900
Podría interpretarse como fugas
de presión durante un cierre
44 46 48 50 52 54 56 58 60 62 64 66
276
Anomalía de temperatura
274 asociada con el sensor y no
con la formación
272
270
268
266
262
45 46 47 48 49 50 51 52 53 54 55 56
sino porque también la muestra del fluido puede estar contaminada, con lodo
presente en el anular.
Para finalizar, cabe mencionar que tanto el tipo de yacimiento como el fluido
productor influyen en el éxito de la prueba.
6.7. Resumen
En este capítulo se presentaron los aspectos conceptuales involucrados con el
diseño de las pruebas y sobre todo el procedimiento de cálculo tanto de las
tasas de producción como la presión de fondo a esperar durante una prueba. El
éxito del diseño se basa en poder contar con información de geología,
petrofísica y sísmica en el caso de pozos exploratorios, ya que definen el tipo
de prueba y el fluido a esperar. Se complementa la sección de diseño con una
breve introducción a las consideraciones y a las contingencias a esperar
durante la implementación de las pruebas. Además, se presenta la
configuración de prueba o sarta y el procedimiento operacional involucrado
durante la ejecución de la prueba.
Bibliografía
267
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 269
CAPÍTULO VII
PRUEBAS EN POZOS INYECTORES DE AGUA
7.1. Introducción
El valor agregado de la información que se obtiene de las pruebas, que se
realizan en los pozos inyectores y productores en campos que están en la
etapa de recuperación secundaria, es muy importante como veremos en este
capítulo. En muchos campos se llevan a cabo campañas anuales de pruebas
(tipo fall off) con la finalidad de monitorear, entre otros parámetros, la presión
media del patrón de inyección y verificar la eficiencia del proceso de inyección
con el tiempo.
Dirección de flujo
Petróleo (Sro)
petróleo
Agua
Agua
Frente de
agua
(transición)
kw
w Sor
M (1)
ko
o Siw
Productor Productor
Petróleo
(banco)
Agua
(banco)
Pozo Inyector
Zona de transición
Productor Productor
Fig. 7.2. Patrón tipo “five spot”. Idealización de los bancos antes del llenado
600
200
Pozo Inyector
0
P = 2868 psia 2072 psia
200
2151 psia
400
-600 -400 -200 0 200 400 600 800 1000
(Sor)
Pozo Productor
Saturación
Petróleo antes
Petróleo antes del barrido del barrido
(So)
Agua inyectada
Agua connata
rw
Fig. 7.4. Perfil de saturaciones
k, , , ct , h
(petróleo)
k, μ, cta, , h
Pe
R
agua
re
6400
Presión de fondo
(inyección)
Fall off
6200
(24 horas)
Período de inyección
(240 horas)
6000
Tasa de
-5000
inyección
(-13560 BPD)
-10000
2 p pD 1 p pD a p pD
Válida en RD rD (2)
r 2 D rD rD p t D
paD p pD en rD RD (3)
paD p pD
M en rD RD (4)
rD rD
k k
a (5) p (6)
ct banco de abgua ct banco de apetroleo
Pe (constante) = 5820 psia, qiny = -13560 BPD, durante 250 horas, h = 168 ft,
rw= 0.3 ft, re = 5000 ft, kh = 13100 md.ft, R = 250ft, M = 2, y D = 2. D es igual a
la razón de difusividades definida como:
k
a ct agua
D (7)
p k
ct petróleo
M=2
D=2
1000
100
banco de Zona de banco de petróleo
agua transición
Inicio efecto de
10
mantenimiento
de presión
6300
Presión de fondo
(pozo inyector)
6100
Período de inyección Fall off
(48 horas)
5900
100
10
1
1E-3 0.01 0.1 1 10
6070
kh = 13100 md.ft
6020
5970
kh = 16500 md.ft
P* = 5820 psia
5920
5870
A fin de ilustrar esto último, se presenta una gráfica tipo semi-log usando los
datos de este mismo ejemplo, pero asumiendo un límite de mantenimiento de
presión.
5960
5920
5880
5800
2100
1700
1300
-200
-400
Fig. 7.14. Datos de presión obtenidos durante el período de cierre (fall off)
1000
100
10
1E-3 0.01 0.1 1 10
2100
1900
1700
1500
1300
-5 -4 -3 -2 -1
1000
100
10
1E-3 0.01 0.1 1 10
1000
100
10
1E-3 0.01 0.1 1 10
mpDMBH (t pDA )
p p* (9)
2.3025
Productor Productor
Isobárico: valor
igual al 50% de
caída de presión
Inyector
Productor Productor
1
p pwf ps pozo inyector
1 n
pwf ps productori (10)
2 n i 1
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 287
Ejemplo
Asumamos un patrón de inyección tal como el que se muestra en la Fig. 7.19.
La distancia entre los pozos productores es de 300 metros para este ejemplo.
Asumimos las siguientes condiciones de caudales de inyección y producción:
rw=0.35ft, h = 28 ft, porosidad = 0.17, B = 1 BR/BN, = 0.698 cp, ct= 6.62 10-6
psi-1. El período de inyección previo al cierre tuvo una duración de 300 horas.
El período de fall off fue de 100 horas de duración (lo cual no es común para
estas pruebas y se presenta sólo a efectos de ilustrar el cálculo de la presión
media). Se asumió un kh = 500 md.ft y una presión de reservorio inicial de
1000 psia. Esta es la presión de yacimiento antes de iniciar el proceso de
inyección y producción.
1500
1100
-250
-500
-750
1000
100
10
1E-3 0.01 0.1 1 10 100
1900
1700
m = -182 psi/ciclo
P* = 833.8 psia
1500
1300
1100
-5 -4 -3 -2 -1
800
1000
1000
400
1000
1100
1150 1050
1000
-400
1000
1000
1000
1000
1000
1000
-800
-1200 -800 -400 0 400 800 1200
Cabe recordar que este análisis es válido sólo para el caso de valores de la
razón de movibilidad, M, y razón de difusividades, D unitarias.
Presión de fondo
Registro continuo de datos de
4800
presión durante cada período de
inyección
4400
t: duración período de inyección
(constante para cada tasa)
4000
t
-500 tasas de Inyección
-1000
1 2 3 4 5 6 7 8 9
La duración del período de inyección, para cada caudal debería ser constante,
1 ó ½ hora de duración es bastante común. La idea es incrementar la presión
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 291
B re
pw pe q 141.2 ln s (11)
kh rw
El estado estacionario implica que tanto re como pe, en la EC. 11, permanecen
constantes, y de aquí que la presión de inyección sea proporcional al caudal y
la tendencia tipo línea recta (se asume régimen de flujo radial para todos los
períodos).
292 // GIOVANNI DA PRAT
162.6 B k
pw pi q log t 3.23 0.87 s (12)
kh ct rw 2
1
4.064qB 2
pw pi t (13)
x f h ct k
Basados en la EC. 13, y asumiendo tanto la longitud mitad de fractura (xf) como
el mismo intervalo de tiempo t, la presión de inyección es proporcional a la
tasa, y de aquí la tendencia de los puntos a definir una línea recta en gráfica,
como la que se muestra en la Fig. 7.25.
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 293
Bibliografía
6. Use of pressure fall off tests to locate the burning front in the Miga field-
Eastern Venezuela. Da Prat et al, SPE 13667 (1984).
295
CAPÍTULO VIII
PRUEBAS EN POZOS HIDRAULICAMENTE FRACTURADOS
8.1. Introducción
En este capítulo se presentan la metodología de análisis de datos presión,
adquiridos en pruebas efectuadas en pozos, en los cuales se realizo trabajo de
estimulación, del tipo fractura hidráulica. Vale mencionar que estas pruebas no
son muy comunes, por cuanto no todos los pozos se fracturan y aunque se
fracturen, el ensayo tipo “post-frac” en muchos casos no se realiza, por cuanto
se asume que de lograrse el incremento esperado en la producción del pozo
según el diseño de la fractura, garantiza el éxito de la misma por una parte, y
en otros casos, el tiempo de prueba puede ser considerable a fin de lograr
alcanzar régimen de flujo radial.
Por ejemplo en yacimiento tipo tight gas, debido a que la permeabilidad efectiva
puede ser tan pequeña como 0.01 md, el tiempo de cierre puede tener una
duración de hasta meses si queremos alcanzar el régimen de flujo radial, es
decir, y como veremos en este capítulo, el régimen de flujo puede ser el
asociado con el plano de la fractura (tipo lineal o bilineal) durante casi todo el
periodo de cierre o fluencia.
kf w
FCD (1)
kx f
297
298 // GIOVANNI DA PRAT
re
kf xf
rw
rwD (2)
xf
En donde, r’w, viene dado como
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 299
s f
r 'w rwe (3)
El radio de pozo efectivo se puede considerar como el radio del pozo, cuya
área se modificó, tal que es igual al área creada por la fractura de longitud xf.
Se puede decir que, el lograr obtener el mayor aumento del radio efectivo, es
objetivo del diseño de una fractura eficiente.
Se puede demostrar que cuando FCD es menor que 0.5 el radio efectivo r 'w es
igual a:
0.28k f w
r 'w (4)
k
En el caso de que FCD sea mayor que 20, el radio efectivo es igual a:
xf
r 'w (5)
2
Por otro lado, si asumimos una longitud de fractura digamos, xf = 100 ft, y una
permeabilidad de 10 mD, y conductividad de fractura, Fc = 2000md.ft (FCD=2)
Implica que de mantener la longitud de fractura constante, así como le valor de
FCD, en el caso de que la permeabilidad del reservorio sea de 50 md, es
300 // GIOVANNI DA PRAT
1000
100
10
Podemos preguntarnos a partir de qué tiempo se logra flujo radial tipo infinito,
de existir el mismo: con base a la regla del ciclo y medio deberíamos esperar
flujo radial a partir de 1 hora, sin embargo en ningún momento, la derivativa es
constante a partir de 1 hora.
a. Conductividad infinita
302 // GIOVANNI DA PRAT
b. Flujo uniforme
k f =
xf
pD t Df (6)
dpD 1
pD (7)
d ln(t Df ) 2
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 303
kh 0.000264k
pD p (8) t Df t (9)
141.2qB ct x 2 f
4.064qB
pw f Pi t (10)
hx f ct k
4.064qB
mFL (11)
hx f ct k
2
4.064qB
kx f
2
(12)
ct mFL h
Luego de este régimen de flujo, se sucede una transición hacia flujo radial
infinito (denominado también seudo-radial), una vez alcanzado el mismo, de
existir limites y los mismos se pueden detectar durante el tiempo de ensayo, se
presentan los regímenes de flujo asociados con los limites, por ejemplo estado
semi-estacionario de ser el área de drenaje volumétrica.
Durante el periodo de flujo semi-radial infinito, la solución a la presión en forma
adimensional viene dada por (con un error menor que el 1 %):
304 // GIOVANNI DA PRAT
pD
1
ln t Dxf 2.2 (13) valido para t D xf 2 (14)
2
1E+9
1E+8
1E+7
1E+6
1E+5
1E-5 1E-4 1E-3 0.01 0.1 1 10 100 1000 10000
La EC 14, permite el cálculo del tiempo mínimo, a fin de alcanzar flujo radial
infinito. Asumiendo tDxf = 2, tenemos que usando la EC 9:
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 305
ct 2
t (hrs ) 2 x f (15)
0 . 000264 k
En el caso del ejemplo simulado, tenemos que el tiempo mínimo es de 673 días
(error menor que el 1%).
Como se aprecia por la EC .15, el tiempo es directamente proporcional a la
distancia de la fractura al cuadrado. A continuación se muestra la solución al
mismo caso, pero una longitud de factura igual a 100 ft:
xf=100
1E+8
xf = 100 ft
1E+7
1E+6
xf = 1000 ft
1E+5
1E-4 1E-3 0.01 0.1 1 10 100 1000
Fig.8.5. Efecto de longitud de fractura, en la duración de los regímenes de flujo lineal y semi-
radial infinito
En este ejemplo se conocen todos los parámetros de reservorio y fluido. Un
parámetro importante de diseño, es el de la permeabilidad del reservorio. En la
práctica debemos de conocer este valor a fin de efectuar el diseño óptimo de
fractura. Basados en la EC.12, y de tener la prueba una duración que tan solo
2
se identifique flujo lineal, solo se puede estimar el producto kx f , es decir, si
queremos conocer la longitud de la fractura es necesario conocer la
permeabilidad del reservorio.
1E+9
1E+8
1E+7
1E+6
1E+5
1E-5 1E-4 1E-3 0.01 0.1 1 10 100 1000 10000
Fig.8.6. Ajuste de los datos de presión usando solución que considera solo daño (rango
negativo) y se asume el mismo valor para kh en un reservorio de infinita extensión
xf
r 'w rwe s (16), de aquí que x f 2rwe s (17)
2
Substituyendo el valor de s = -7.39, en la EC. 17, tenemos que Xf = 1062 ft, el
cual no difiere mucho del valor usado para la simulación del ejemplo (xf = 1000
ft). Este ejemplo fue solo para efectos de ilustración del concepto histórico en
cuanto al radio de pozo equivalente.
5000
4500
[p s ia ]
4000
[S T B /D ]
100
0
0 10 20 30 40 50 60 70
Fig.8.7. Historia de la prueba (se asume a efectos de este ejercicio de que el caudal es
constante durante la fluencia, lo cual sabemos que no es cierto ya que declina acorde con
régimen de fractura, y en la práctica seria la historia de caudales a usar para el análisis de los
datos del periodo de cierre)
Pi = 5000 psi, rw = 0.3 ft , h = 100 ft, = 10%, B = 1.30 RB/STB, = 0.48 cp, ct
= 1.25(10-5)psi-1. El fluido es petróleo. Durante el periodo de producción, el
caudal fue de 200 BPD.
10000
1000
100
10
1
1E-5 1E-4 1E-3 0.01 0.1 1 10 100
Análisis de datos:
En este caso, mediante análisis de grafica tipo cartesiana, se puede obtener el
2
producto kx f .
En la siguiente figura se presenta la gráfica tipo cartesiana, de la presión en
función de la raíz cuadrada del tiempo de fluencia.
5000
4200
3800
0 1 2 3 4 5 6 7
2
El valor del producto kx f ., obtenido mediante el programa se basa en el uso
de la EC 12. Podemos corroborar este valor, substituyendo valores en la
EC.12.
kx f 1056mdft 2
2
0.100.0000125 201.5100
p ws p * 4.064
qB
hx f
ct k
t p t t (18)
100
10
1
1E-4 1E-3 0.01 0.1 1 10
4300
4100
3900
Resultados Saphir
mLF= -202 psi/(hrs)**0.5
2
3700 kx f = 1070 mDft2
Usando la información, así como datos del ejercicio anterior, compruebe de que
la duración del periodo de fluencia, tal que se logre alcanzar por lo menos el
tiempo de inicio del flujo radial tipo reservorio infinito, es de 280 días.
circulo sellado
constant p
Área: circulo
(volumterico)
1000
Área: Circulo
(Presión
100 constante)
10
0.01 0.1 1 10 100 1000 10000 1E+5
La fractura tipo flujo uniforme, asume que el caudal que entra al plano de la
fractura por unidad de área, y por unidad de tiempo es constante, a diferencia
de la condición que se asume en la fractura de tipo conductividad infinita, en
donde la presión es constante en la cara de la fractura.
Se puede demostrar que la solución a la fractura de conductividad infinita, es la
misma solución que se obtiene para el caso de fractura de caudal uniforme, si
se calcula a una distancia medida desde el origen, igual a x = 0.732xf.
1E+9
Modelo
Caudal en fractura
Constante
1E+8
1E+7
Modelo
Conductividad
1E+6
Infinita
1E+5
1E-3 0.01 0.1 1 10 100 1000 10000 1E+5
Como se aprecia en la Fig. 8.13, la diferencia entre una solución y la otra, está
principalmente en la duración del periodo de transición entre flujo lineal y radial,
que para el caso de la fractura tipo conductividad infinita es mayor.
De aquí que, y a efectos de entendimiento de este modelo, consideramos que
lo explicado en la sección anterior (fractura de conductividad infinita) se aplica
igualmente al caso de fractura vertical tipo caudal de fractura constante (flujo
uniforme).
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 313
re
w
kf
(fractura)
Xf k
kf w (reservorio)
FCD
kx f
qB
p 44.1 4
t
h k f wf ct4
100
10
kf w
FCD k f Dwf D 2
kx f
1
0.01 0.1 1 10 100 1000
2.45 dpD 1
pD (19) pD (20)
k f w f
4 t Dxf
d ln(t Dxf ) 4
D
qB
p 44.1 4
t (21)
h k f w f ct 4
qB
mBF 44.1 (22)
h k f w f 4 ct
qB 1
k f w f 1944.8 2
(23)
BL
hm ct k
Valor de la pendiente,
mBF= 69 psi/(hr)**1/4
250
150
50
1 2 3 4 5
4
Fig.8.16. Gráfico cartesiano del cambio de presión en función de t
316 // GIOVANNI DA PRAT
(100)(1.65)(0.85) 1
k f w f 1944.8 2651md . ft
2
FCD= 240
Xf =500ft
10000 kh=5 md.ft
Régimen de
flujo lineal
1000
Régimen de
flujo Bilineal
100
10
1
1E-3 0.01 0.1 1 10 100 1000 100
Conductividad Infinita
Caudal Uniforme
Conductividad
Infinita
Flujo Uniforme:
Pendiente = 1/2
1000
Conductividad Finita
Pendiente = 1/4
100
10
0.01 0.1 1 10 100 1000 1000
En cuanto al tiempo de inicio del régimen de flujo radial tipo infinito, podemos
usar la expresión dada por la EC 15, válida para cualquier de los modelos
usados en el ajuste. Lo mismo aplica en cuanto al estudio de límites, es decir
podemos usar la solución dada por la EC. 19.
pws p * 44.1
qB
t p t
0.25
t
0.25
(24)
h k f w f 4 ct k
2950
2850
[psia]
2750
[STB/D]
50
0
0 200 400 600 800 1000 1200 1400
250
200
150
Pendiente (mBF):
100
69.6 psi/(hrs)**25
50
2940
2890
Resultados
kh = 205 mD.ft
k= 6.84 mD
2840 S = -5.34
2790
2740
2690
-5 -4 -3 -2 -1
Fig.8.21. Gráfica tipo Horner aplicada a periodo final de flujo radial tipo infinito
Los valores obtenidos del análisis son prácticamente iguales a los usados en la
simulación. Usando el valor de la permeabilidad se puede estimar la longitud de
la fractura.
Sin embargo, la determinación de la longitud de fractura xf, mediante concepto
de radio de pozo equivalente dado por la EC 5, en este caso no aplica, por
cuanto la conductividad adimensional de la fractura no es grande (mayor que
20 por ejemplo). Haciendo uso de correlaciones (Ver Ref.1), se tiene que para
una conductividad adimensional (FCD) de 2, como es nuestro caso, el radio de
pozo efectivo es igual a 0.3, de aquí que:
320 // GIOVANNI DA PRAT
El daño puede deberse a varios factores, al igual que en el caso de pozos sin
fractura, como se explicó en él capitulo tres. Los factores más importantes que
pueden ocasionar daño en fracturas son
1. Que la permeabilidad de la fractura en la zona cercana al pozo, kfd, no
sea igual a la del resto de la fractura, kf. Este efecto se conoce como
fractura “estrangulada” (choked fracture).
Zona de daño
(alrededor del
pozo)
xd
kfd kf
Reservorio de permeabilidad k xf
El efecto de daño, y dependiendo del tipo, puede causar una caída de presión
adicional a la solución de la presión esperada sin daño, tal como se explicó en
él capitulo 3, ó puede reducir la productividad esperada del pozo, con base al
diseño de la fractura.
A continuación se presenta la solución a la presión considerando varios valores
de daño y su efecto:
1E+8
Daño de fractura, sf = 0.5
Daño de fractura, sf = 0
1E+7
1E+6
1E-3 0.01 0.1 1 10 100
1/2
1/4
100
10
1
0.01 0.1 1 10 100
1E+9
1E+8
1E+7
8.4.1. Caso 1
Los datos de presión son simulados, pero se basan en parámetros de
reservorio y caudales obtenidos del análisis de un ensayo de restauración de
presión, realizado en un pozo que fue hidráulicamente fracturado, es decir,
usando los valores de los parámetros de reservorio, fluido, conductividad y
longitud de fractura, que se presentan en el ejemplo de campo C (Ver. Ref.6),
se construyó la historia de la prueba, en cuanto a la secuencia de eventos: flujo
y cierre y las presiones.
del pozo durante 300 horas, a fin de evaluar la condición de productividad del
pozo, principalmente reflejada en el incremento en la producción, y verificar si
los valores obtenidos para la longitud y conductividad de fractura son los
esperados según diseño.
2800
2400
2000
200
cierre
Periodo de producción
[S T B /D ]
100 post-fractura
0
0 200 400 600 800 1000 1200 1400 1600 1800 2000 2200
Método de análisis de datos (se asume que los datos de presión tienen la
calidad de medición, según las especificaciones del sensor de presión usado):
Pendiente = ½
1000
Pendiente = ¼
100
Pendiente = 0
10
En este caso tenemos información del yacimiento, que se puede obtener del
periodo de flujo radial. Podemos calcular el daño de pozo, y su valor está
relacionado con la longitud de la fractura mediante el concepto de radio de
326 // GIOVANNI DA PRAT
2600
S = -4.69
kh=103 md.ft
2400 k=2.11 mD
2200
2000
-5 -4 -3 -2 -1
1000
Conductividad Finita
100
100
1
0.01 0.1 1 10 100
radioequivalente
conductividad finita
2700
2300
1900
200
[STB/D]
100
Modelo Conductividad finita
0
0 200 400 600 800 1000 1200 1400 1600 1800 2000 2200
8.4.2. Caso 2
Al igual que el caso anterior, en este caso se presenta un ejemplo de prueba
inspirado en datos de campo. Se trata de un pozo terminado en yacimiento de
gas de muy baja permeabilidad. A continuación se presenta la historia del
ensayo:
3800
3300
[psia]
2800
2300
[Mscf/D]
1000
0
0 50 100 150 200 250 300 350 400
Tendencia de
Pendiente de 1/2
1E+8
1E+7
1E+6
1E-3 0.01 0.1 1 10 100
3400
3200
3000
kx f 1510mdft 2
2
2800
2600
6 7 8 9 10 11 12 13 14
xf
Ya que se trata de una fractura de conductividad infinita, se tiene que r 'w ,
2
en donde
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 331
Flujo Uniforme
conductividad infinita
Conductividad Finita
1E+8
1E+7
1E+6
1E-3 0.01 0.1 1 10 100
Los valores obtenidos para los parámetros de reservorio y fractura, usando los
tres modelos, se presentan a continuación:
Modelo k(mD) Xf (ft) Fc (mD.ft)
Flujo Uniforme 0.0496 222 infinita
Conductividad Infinita 0.0159 374 infinita
Conductividad Finita 0.019 325 3390 (infinita)
Basados en inspección visual de la figura 8.35., pareciera que el modelo que
definitivamente logra un mejor ajuste de todos los datos, es el de fractura tipo
flujo uniforme.
3500
[psia]
3000
2500
[Mscf/D]
1000
0
0 50 100 150 200 250 300 350 400
3470
3450
Conductividad Infinita
3430 Conductividad Finita
c
Fig.8.37. Simulación usando los tres modelos (ultimas 28 horas del cierre)
Determine mediante la fórmula pertinente el tiempo mínimo necesario de
prueba, tal que se muestre el inicio del régimen de flujo radial tipo infinito.
Sugiere una prueba de esta duración?
Bibliografía
Apéndice A
APÉNDICE A
MODELOS DE YACIMIENTO Y POZO: EJEMPLOS (CURVAS
TIPO)
En esta sección se presentan las respuestas esperadas para el cambio de
presión y su derivativa, de varios modelos de yacimiento y pozo y las cuáles se
aplican en el proceso de ajuste por curvas tipo. No es la intención presentar
todas las soluciones de los modelos existentes, sino las que consideramos que
son de uso común en la interpretación de datos de presión. Sin embargo, se
sugiere al lector que consulte las referencias para información adicional y otras
soluciones. La familiarización visual de estas soluciones es de gran ayuda en la
identificación inicial del modelo de yacimiento aplicable para el análisis de los
datos.
Cada figura muestra la respuesta del cambio de presión y la derivativa
correspondiente en función de los parámetros de yacimiento y pozo que
caracterizan el modelo. Se asume que la respuesta corresponde a la de un
pozo productor de petróleo o gas fluyendo a una tasa constante. La idea es la
de presentar sólo la respuesta esperada en cuanto a la forma o tendencias de
las curvas, según la variación de los parámetros de yacimiento y pozo que
caracterizan cada modelo. La generación de la respuesta o solución se obtuvo
mediante la opción de diseño del software (en este caso Saphir). Se asumieron
parámetros de yacimiento así como de pozo, a fin de generar las soluciones
que aquí se presentan.
200
10000 100
50
40
30
20
10
5
1000
S=0
100
10
=1 10-6
1000
100
=0.9
=0.1
10
=0.05
=0.01
1
0.01 0.1 1 10 100
1000
100
Flujo transiente
10
Flujo semi-estacionario
1
1E-3 0.01 0.1 1 10 100
= .01
1000
100
-4
= 1.10
-6
= 1.10 = 1.10
-8
10
1
1E-3 0.01 0.1 1 10 100
1000 -6
C = 10 Bbl/psi
C=0.001Bbl/psi
100
10
C = 0.01 Bbl/psi
1000
100
10
1
1E-3 0.01 0.1 1 10 100
10000
1000
100
= 0.1
= 0.5
10
= 0.99
1000
100
10
Datos de presión registrados en
un yacimiento de doble porosidad
-6
( = 0.01 y = 10 )
1
1E-3 0.01 0.1 1 10 100 1000
1000
100
10
Ajuste de los datos de presión
registrados en un yacimiento de doble
porosidad aplicando modelo de doble
permeabilidad:
= 0.01 y = 10-6 y = 1
1
1E-3 0.01 0.1 1 10 100 1000
Fig. 5. La figura muestra que el modelo de doble permeabilidad (ver Fig. 4.)
puede ajustar los datos de presión registrados en un yacimiento de doble
porosidad (Fig.2). Este se debe a que los dos modelos son idénticos desde el
punto de vista matemático si el parámetro que caracteriza al modelo de doble
permeabilidad es igual a 1. Sistema infinito.
342 // GIOVANNI DA PRAT
Xf = 500 ft
1000
Xf = 2000 ft
100
10
1
0.01 0.1 1 10 100 1000
FC = 500 md.ft
100
10000 md.ft
10
0.01 0.1 1 10 100 1000
1
1000
.5
.1
S=0
100
10
0.01 0.1 1 10 100
Efecto de daño (S). Fractura de conductividad Finita.
C = 0.01 Bbl/psi
Xf = 300 ft
1000 ft
1000
2000 ft
100
10
Fig.7. Pozo con fractura hidráulica tipo conductividad finita: Efecto de daño y
Modelo de fractura hidráulica tipo FLUJO UNIFORME (sensibilidad a la
longitud de la fractura). Se asume almacenamiento de pozo, C = 0.01 Bbl/psi.
Sistema infinito.
344 // GIOVANNI DA PRAT
1000
L = 500 ft
100
L = 50 ft
L = 5000 ft
10
L = 1000 ft
1
0.01 0.1 1 10 100 1000
1000
FALLA sellante
100
10
FALLA no conductiva
1000
100
10 Falla no conductiva
L1
L1 = L2 = 500 ft
L2
1
Falla no conductiva
10000
1000
L1 = 50 ft
500 ft
1000 ft
100 5000 ft
L2
yacimiento
L1
10000
ANGULO = 30
ANGULO = 60
ANGULO = 90
1000
= 30ª
= 60ª
= 90ª
100
10
SELLO SELLO
kz
h
kr
ACUIFERO ACTIVO
Acuífero Activo
10000
1000
kz/kr = 0.0005
kzkr = 0.005
100
10
0.01 0.1 1 10 100 1000
sello
sello
1E+5
10000
L = 1500 ft
L = 800 ft
1000
L = 200 ft
100
sello
Kz
h
h/2 L
Kr
sello
1000
L = 300 ft
1000 ft
3000 ft
100
10
sellante
Capa 1
h1 = 20 ft Pi1 = 5000 psia, k = 5mD, S = 0, = 10 %
sellante
h2 = 5 ft Capa 2 Pi2 = 4500 psia, k = 20 mD, S = 0, = 15%
sellante
10000
1000
100
Fig.14. Sistema de dos capas (no hay flujo cruzado a nivel de formación solo
en el pozo. Las capas no tienen el mismo valor para la presión inicial.
Sensibilidad a la presión inicial de la capa inferior. No hay discontinuidades, ni
limites presentes es decir un sistema lateralmente infinito. Note que la solución
es diferente a la del modelo de doble permeabilidad (two layers) la cual asume
flujo cruzado a nivel de formación (Ver Fig.4)
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 351
k h
2000 ft
Kh = 50 md.ft
Kh = 500 md.ft
1E+7
Kh = 5000 md.ft
1E+6
1E+5
pozo
k1, 1 k1 , 1 k2 , 2
L
1E+9
1E+8
M = 50
1E+7
M= 5
M = 0.5
1E+6
Zona cañoneada
hw
h
Kz
Kr
h = 30
h=5
h = 10
h = 20
1000 h=3
3 ft
5 ft
10 ft Influencia de
flujo tipo
100
20 ft esférico
30 ft
Espesor cañoneado
10
(hw)
Bibliografía
357
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 359
APÉNDICE B
SOLUCIÓN DE LA PRESIÓN TRANSIENTE A ESPERAR
DURANTE PERÍODO DE CIERRE (EJEMPLOS)
Pozo horizontal
INTRODUCCIÓN AL ANÁLISIS Y DISEÑO DE PRUEBAS DE PRESIÓN // 365