CAPÍTULO 6

GENERALIZACION DE ANÁLISIS DE PRUEBAS DE POZOS

Las soluciones presentadas para yacimientos produciendo bajos diversas

condiciones de flujo son usualmente complejas y de difícil aplicación por las
dificultades numéricas que presentan. Una alternativa es la solución del problema
usando el principio de superposición. Con la aplicación del principio de
superposición, es posible producir tabulaciones de presiones adimensionales
determinadas en varias localizaciones, dentro del sistema de flujo rectangular
que tiene un pozo localizado en una determinada localización. El sistema
cuadrado prevé un bloque de construcción que puede ser usado para generar el
comportamiento de flujo, para cualquier forma rectangular, para una razón dada
en sus lados. Los yacimientos rectangulares pueden tener cualquier número de
pozos produciendo a cualquier tasa de flujo definida. Las condiciones de contorno
pueden ser de diversas naturalezas.

EL PROBLEMA

Considere un flujo de comprensibilidad constante, fluyendo en un yacimiento

homogéneo, de porosidad, permeabilidad y espesor constantes, produciendo a
una tasa constante. Bajo estas condiciones la ecuación de difusividad en términos
de presión describe el flujo a través del sistema.

1
En 1942 Van Everdingen y Hurst publicaron soluciones para el problema de
intrusión de agua dentro de un yacimiento cilíndrico.
Como el problema es idéntico matemáticamente al de un yacimiento produciendo
desde, un pozo situado en el centro de un yacimiento rectangular, la solución de
van Everdingen y Hurst1, puede ser aplicada al problema del pozo.

Si la ecuación de difusividad es:

6-1
∂ 2p 1 ∂p φµc ∂t
+ = (6-1)
∂r 2 r ∂r kt ∂r

Van Everdingen y Hurst1, demostraron que la solución al problema del pozo

produciendo a tasa de flujo constante (para yacimientos finitos e infinitos), podría
ser expresada en función de las variables adimensionales.

Así:

2πkh (p i − p)
pD (rD , t D ) = (6-2)
qµ

Donde:

rD = r / rw
(6-3)

kt
tD = (6-4)
φµcrD2

Si el radio del pozo se vuelve cero, rw → 0 , la solución de la línea fuente puede

ser usada en la ecuación (6-2) para sistemas infinitos.

1  φµcr 2  1  rD2 
pD (rD , tD ) = − Ei  −  = − Ei  −  (6-5)
2  4kt  2  4 tD 
donde -Ei (-x) es la integral exponencial. Si el valor del argumento es muy
pequeño ( x < 0.01) ,

1  φµcr 2  
pD (rD , tD ) = − ln  + 0.5772 (6-6)
2   4kt  

6-2
las ecuaciones (6-5) y (6-6) son aproximaciones a la ecuación (6-2) bajo ciertas
condiciones.

2
En 1954, Mathews, Brons y Hazebroek demostraron que las soluciones tales
como la ecuación (6-5) podían ser superpuestas para generar el comportamiento
de yacimientos cerrados para diversas formas geométricas, esto es, el
comportamiento de un pozo en un sistema cerrado en el contorno, podría ser
calculado añadiendo las presiones causadas por un arreglo infinito de pozos,
produciendo en un yacimiento infinito. Esta superposición puede ser representada
analíticamente como:
∞
pD ( xD , yD , tD ) = ∑ pD (aiD2 , tDA ) (6-7)
i =1

donde:
aiD = ai / A
ai = Distancia del pozo “i” al punto (xD , yD)
A = Área de drenaje de un yacimiento limitado (ó área por pozo)
y:

kt  rw2 
tDA = = tD   (6-8)
φµcA  A

NOTA: Las Figuras (6-1) y (6-2) son adaptadas de Earlougher y asociados 3. La

Figura 6.2 presenta el arreglo infinito de pozos usados para generar el
comportamiento de un pozo fluyendo a tasa de flujo constante en el centro de un
yacimiento del área cuadrada.
La Tabla 6.1 presenta los resultados en forma de presión adimensional como
función de tiempo adimensional para varios puntos del yacimiento cuadrado y
cerrado al flujo de su contorno exterior.
En el caso de las presiones en el pozo, los cálculos fueron efectuados para

A / rw = 2000 . Estos resultados pueden ser convertidos para otro valor de

6-3
 A / rw usando las ecuaciones (6-9) ó (6-10). Esto es, agregando y restando el

término logarítmico. ln ( rw2 / 4 AtDA ) , se obtiene :

   
1  1   1   1 
pD (rwD , tDA ) = − ln  − ln + +
( ) ( ) 2tDA   (
 ) DA 
2  2
4 A / rw2 tDA  4 A / rw2 4 2000 t
   

∞  2 
a
0.5772 + ∑ Ei  − 4 iDtDA 
i= 2

(6-9)

 
 1 1  1 
( 1
pD rwD , tDA , A / rw = 2x10 3 = − 
 1
)
 + − ln +

2

4 2x10 tDA 
3

2 2  4 At
(
DA 
)
∞   a 2  
0.5772 + ∑ Ei − iD  
  4 tDA  
(6-10)
i= 2

ó
 2
1  A / rw 
pD (rwD , tDA ) = − ln
( )
+ pD ( rwD , tDA A / rw = 2x10 3 )
2 
2 

2x10 3 
 ( )
(6-11)
(Ver Tabla 6.2)

Luego:
 A / r2 
(
pD (rwD , tDA ) = pD rwD , tDA , A / rw = 2x10 3 + ) 1
2
ln  w

 4 x10 6  (6-12)

(Ver Tabla 6.1)

6-4
El término ½ ln [(A / r ) / 4x10
2
w
6
] es un término de corrección.

Los resultados para puntos diferentes al pozo son independientes de (A/rw2) en

este tratamiento, excepto para tiempos muy pequeños, para las cuales la solución
de la línea fuente debe sustituirse por la solución exacta de radio del pozo finito.

Por conveniencia, las variables adimensionales de espacio (xD , yD) , usadas en la

Tabla 6.1, son definidas como fracción de la mitad del lado del cuadrado . El
origen se toma en el pozo, (ver Fig. 6.1). Esto es:

x 2x
xD = = (6-13)
L A

y 2y
yD = = (6-14)
L A

donde:

L=½ A.

Luego:

la razón

rw / L =0.001. (6-15)

La Figura 6.3 es un gráfico de presión adimensional para algunos de los puntos

considerados en la Figura 6.1. La Figura 6.3 muestra, que a tiempos
suficientemente grandes el sistema exhibe flujo pseudo-continuo, esto es, la
presión adimensional se convierte en una función lineal de tiempo adimensional.

6-5
Obsérvese que no todos los puntos del yacimiento alcanzan el Flujo semi -
continuo al mismo tiempo. Las presiones, a tiempos más grandes que las
mostradas en la Tabla 6.1, pueden ser obtenidas por extrapolación lineal.

Supóngase que se requiere el comportamiento de presión para un pozo en el

centro de un rectángulo con lados en la relación de 2:1. El comportamiento de este
yacimiento podría ser determinado usando la ecuación (6-9) para un arreglo
infinito de pozos (Ver Figura 6.2) pero con las distancia entre los pozos en la
dirección x. Este arreglo particular resulta de un segundo arreglo de pozos como
el mostrado en la Figura 6.2, de tal forma que, el pozo origen del segundo arreglo
este localizado a xD = 0, yD = 1.0. Esto es equivalente a sobreponer un sistema
cuadrado, con un pozo en el centro, con otro sistema cuadrado pero, con el
segundo pozo en el punto (0,1.0). El área del rectángulo 2:1 formado, A′ , es la
mitad del área, A, del cuadrado inicial. La caída de presión en cualquier punto
dentro del rectángulo 2:1, es simplemente la suma de las caídas de presión
causadas por cada uno de los sistemas cuadrados en el punto común.

La Figura 6.4 ilustra la superposición de los dos arreglos de cuadrados descritos

previamente. Los círculos abiertos indican las localizaciones de los pozos
pertenecientes al arreglo inicial, mientras que los círculos llenos indican las
localizaciones de los pozos del arreglo superpuesto. Es claro que la nueva unidad
de simetría es un rectángulo de 2:1 con un pozo en el centro. (Ver Anexo 6).

CÁLCULO DE PD (XD,YD,TDA ) PARA UN RECTÁNGULO 2:1 ( TABLA 6.2)

Una línea de transición bastante amplia sería detectada si la posición del pozo se
apartara del centro del área de drenaje.
Otros efectos pueden enmascarar la línea recta semilog inicial, tal como el efecto
de llene. Sin embargo, esto no será considerado en el análisis que a continuación
se tratará.

6-6
 1.0
Punto de localizaciones de x
presiones notadas con X’s
0.8

x x
0.6

x x x
0.4

x x x
x 0.2

pozo x
0.0
x x x x x
0.0 0.2 0.4 0.6 0.8 1.0

Figura 6.1. Octante de un cuadrado mostrando el pozo y las localizaciones de los

puntos de presión.
6 5 4 5 6

5 3 2 3
2A 5

4 2 1
2A 2 4

A A

5A
5 3 2 3 5

8A
6 5 4 5 6

FIGURA 6.2. Superposición de un arreglo infinito de pozos para generar el

comportamiento de presión de un yacimiento cuadrado con un pozo en el centro.

6-7
Tabla 6.1. Resultados en forma de presión adimensional como función de tiempo
adimensional para varios puntos del yacimiento cuadrado y cerrado al flujo de su
contorno exterior

PRESIÓN ADIMENSIONAL
Tiempo
adimensional
XD=0,000 XD=0,250 XD=0,250 XD=0,500 XD=0,500 XD=0,500 XD=0,750 XD=0,750
tDA YD=0,000 YD=0,000 YD=0,250 YD=0.000 YD=0,250 YD=0,500 YD=0,000 YD=0,250
0,0010 4,5516 0,0021 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0015 4,7543 0,0109 0,0004 0,0000 0,0000 0,0000 0,0000 0,0000
0,0020 4,8910 0,0261 0,0021 0,0000 0,0000 0,0000 0,0000 0,0000
0,0025 5,0097 0,0456 0,0056 0,0001 0,0000 0,0000 0,0000 0,0000
0,0030 5,1009 0,0675 0,0109 0,0004 0,0001 0,0000 0,0000 0,0000
0,0040 5,2447 0,1141 0,0261 0,0021 0,0007 0,0000 0,0000 0,0000
0,0050 5,3563 0,1607 0,0456 0,0056 0,0021 0,0001 0,0001 0,0000
0,0060 5,4474 0,2053 0,0675 0,0109 0,0047 0,0004 0,0002 0,0001
0,0070 5,5245 0,2475 0,0906 0,0178 0,0085 0,0011 0,0006 0,0003
0,0080 5,5913 0,2871 0,1141 0,0261 0,0135 0,0021 0,0012 0,0007
0,0090 5,6502 0,3243 0,1376 0,0354 0,0194 0,0036 0,0021 0,0013
0,0100 5,7029 0,3592 0,1607 0,0456 0,0261 0,0056 0,0034 0,0021
0,0150 5,9056 0,5063 0,2676 0,1023 0,0675 0,0218 0,0154 0,0109
0,0200 6,0494 0,6211 0,3592 0,1607 0,1142 0,0456 0,0350 0,0266
0,0250 6,1610 0,7147 0,4379 0,2164 0,1609 0,0735 0,0597 0,0474
0,0300 6,2522 0,7939 0,5065 0,2685 0,2061 0,1032 0,0876 0,0716
0,0400 6,3965 0,9232 0,6224 0,3628 0,2906 0,1650 0,1485 0,1263
0,0500 6,5099 1,0279 0,7192 0,4470 0,3685 0,2276 0,2125 0,1854
0,0600 6,6050 1,1178 0,8041 0,5242 0,4415 0,2904 0,2772 0,2466
0,0700 6,6888 1,1983 0,8815 0,5968 0,5112 0,3532 0,3418 0,3086
0,0800 6,7654 1,2728 0,9539 0,9539 0,6661 0,5786 0,4160 0,4061
0,0900 6,8374 1,3434 1,0231 0,7334 0,6446 0,4788 0,4700 0,4338
0,1000 6,9063 1,4114 1,0902 0,7992 0,7095 0,5417 0,5363 0,4965
0,1500 7,2311 1,7347 1,4119 1,1186 1,0274 0,8558 0,8482 0,8106
0,2000 7,5468 2,0501 1,7271 1,4335 1,3421 1,1700 1,1636 1,1248
0,2500 7,8611 2,3644 2,0414 1,7478 1,6563 1,4841 1,4778 1,4390
0,3000 8,1753 2,6786 2,3556 2,0620 1,9705 1,7983 1,7919 1,7531
0,4000 8,8036 3,3069 2,9839 2,6903 2,5988 2,4266 2,4202 2,3814
0,5000 9,4320 3,9352 3,6122 3,3186 3,2271 3,0549 3,0486 3,0098
0,6000 10,0603 4,5636 4,2406 3,9469 3,8555 3,6833 3,6769 3,6381
0,7000 10,6886 5,1919 4,8689 4,5752 4,4838 4,3116 4,3052 4,2664
0,8000 11,3169 5,8202 5,4972 5,2036 5,1121 4,9399 4,9335 4,8947
0,9000 11,9452 6,4485 6,1255 5,8319 5,7404 5,5682 5,5618 5,5230
1,0000 12,5735 7,0269 6,7538 6,4602 6,3687 6,1965 6,1902 6,1513
2,0000 18,8567 13,3600 13,0370 12,7433 12,6519 12,4797 12,4733 12,4345
4,0000 31,4230 25,9263 25,6033 25,3097 25,2182 25,0460 25,0397 25,0009
8,0000 46,5557 51,0590 50,7360 50,4423 50,3509 50,1787 50,1723 50,1335
10,0000 69,1220 63,6253 63,3023 63,0087 62,9172 62,7450 62,7386 62,6999

6-8
Tabla 6.1. Resultados en forma de presión adimensional como función de tiempo
adimensional para varios puntos del yacimiento cuadrado y cerrado al flujo de su
contorno exterior. (Continuación).

PRESIÓN ADIMENSIONAL
Tiempo
adimensional
XD=0,750 XD=0,750 XD=0,875 XD=0,875 XD=1,000 XD=1,000 XD=1,000 XD=1,000 XD=1,000
tDA YD=0,500 YD=0,750 YD=0,000 YD=0,500 YD=0,000 YD=0,250 YD=0,500 YD=0,750 YD=1,000
0,0010 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0015 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0020 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0025 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0030 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0040 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0050 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0060 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0070 0,0000 0,0000 0,0001 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0080 0,0001 0,0000 0,0002 0,0000 0,0000 0,0000 0,0000 0,0000 0,0000
0,0090 0,0003 0,0000 0,0004 0,0001 0,0001 0,0001 0,0000 0,0000 0,0000
0,0100 0,0005 0,0001 0,0008 0,0001 0,0003 0,0002 0,0000 0,0000 0,0000
0,0150 0,0040 0,0008 0,0055 0,0016 0,0031 0,0023 0,0009 0,0002 0,0001
0,0200 0,0121 0,0036 0,0164 0,0060 0,0111 0,0087 0,0042 0,0014 0,0005
0,0250 0,0245 0,0091 0,0329 0,0143 0,0249 0,0203 0,0112 0,0045 0,0023
0,0300 0,0404 0,0177 0,0539 0,0264 0,0436 0,0365 0,0219 0,0102 0,0062
0,0400 0,0805 0,0437 0,1050 0,0600 0,0913 0,0793 0,0532 0,0307 0,0223
0,0500 0,1281 0,0800 0,1628 0,1030 0,1469 0,1308 0,0947 0,0623 0,4998
0,0600 0,1807 0,1241 0,2237 0,1525 0,2065 0,1871 0,1431 0,1029 0,0872
0,0700 0,2366 0,1740 0,2859 0,2064 0,2678 0,2460 0,1962 0,1502 0,1321
0,0800 0,2948 0,2279 0,3486 0,2632 0,3299 0,3064 0,2525 0,2023 0,1826
0,0900 0,3146 0,2846 0,4114 0,3219 0,3925 0,3677 0,3109 0,2579 0,2369
0,1000 0,4153 0,3334 0,4744 0,3820 0,4551 0,4296 0,3708 0,3157 0,2939
0,1500 0,7257 0,6500 0,7888 0,6913 0,7692 0,7421 0,6797 0,6209 0,5976
0,2000 1,0393 0,9632 1,1030 1,0047 1,0834 1,0561 0,9931 0,9338 0,9103
0,2500 1,3534 1,2772 1,4172 1,3188 1,3975 1,3702 1,3071 1,2478 1,2243
0,3000 1,6676 1,5913 1,7313 1,6330 1,7117 1,6843 1,6213 1,5620 1,5384
0,4000 2,2959 2,2196 2,3597 2,2613 2,3400 2,3127 2,2496 2,1903 2,1667
0,5000 2,9242 2,8479 2,9880 2,8896 2,9683 2,9410 2,8779 2,8186 2,7950
0,6000 3,5525 3,4763 3,6163 3,5179 3,5966 3,5693 3,5062 3,4469 3,4233
0,7000 4,1808 4,1046 4,2446 4,1462 4,2249 4,1976 4,1346 3,0756 4,0517
0,8000 4,8092 4,7329 4,8729 4,7745 4,8533 4,8269 4,7629 4,7036 4,6800
0,9000 5,4375 5,3612 5,5012 5,4029 5,4816 5,4542 5,3912 5,3319 5,3083
1,0000 6,0658 5,9895 6,1296 6,0312 6,1099 6,0826 6,0195 5,9602 5,9366
2,0000 12,3490 12,2727 12,4127 12,3144 12,3930 12,3657 12,3027 12,2434 12,2198
4,0000 24,9153 24,8321 24,9791 24,8807 24,9594 24,9321 24,8690 24,8097 24,7861
8,0000 50,0480 49,9717 50,1117 50,0134 50,0921 50,0647 50,0017 49,9424 49,9188
10,0000 62,6143 62,5381 62,6781 62,5797 62,6584 62,6311 62,5680 62,5087 62,4851

6-9
 TABLA 6.2. CÁLCULO DE pD (xD, yD, tDA) para un rectángulo 2:1

PD (XD. YD. TDA) para cuadrados en las localizaciones pD (xD. ND. IDA) en rectángulo 2:1

(1) (2) (3) (4) (5) (6) (7) (8)

rw (0.5.0) tDA
tDA Iw (0.5.0) (1.0) (1.0.5) (2) + (3) + 2 (1)
(4) (5)
* ** ***
.0010 4.5516 .0000 .0000 .0000 4.5516 .0000 .0020
.0015 4.7543 .0000 .0000 .0000 4.7543 .0000 .0030
.002 4.8981 .0000 .0000 .0000 4.8981 .0000 .0040
.0025 5.0097 .0001 .0000 .0000 5.0097 .0001 .0050
.0030 5.1009 .0004 .0000 .0000 5.1009 .0004 .0060
.0040 5.2447 .0021 .0000 .0000 5.2447 .0021 .0080
.0050 5.3563 .0056 .0000 .0000 5.3563 .0056 .0100
.0060 5.4474 .0109 .0000 .0000 5.4474 .0109 .0120
.0070 5.5245 .0178 .0000 .0000 5.5245 .0178 .0140
.0080 5.5913 .0261 .0000 .0000 5.5913 .0261 .0160
.0090 5.6502 .0354 .0001 .0000 5.6503 .0354 .0180
.0100 5.7029 .0456 .0003 .0000 5.7032 .0456 .0200
.0150 5.9056 .1023 .0031 .0009 5.9087 .1032 .0300
.0200 6.0494 .1607 .0111 .0042 6.0605 .1649 .0400
.0250 6.1610 .2164 .0249 .0112 6.1859 .2276 .0500
.0300 6.2522 .2685 .0436 .0219 6.2958 .2904 .0600
.0400 6.3965 .3628 .0913 .0532 6.4878 .4160 .0800
.0500 6.5099 .4470 .1469 .0547 6.6568 .5417 .1000
.0600 6.6050 .5247 .2065 .1431 6.8115 .6673 .1200
.0700 6.6888 .5968 .2678 .1962 6.9566 .7930 .1400
.0800 6.7654 .6661 .3299 .2525 7.0953 .9186 .1600
.0900 6.8374 .7334 .3925 .3109 7.2299 1.0443 .1800
.1000 6.9063 .7992 .4551 .3708 7.3614 1.1700 .2000

6-10
 14
LIMITE EXTERIOR CERRADO

12

10
PD

LIMITE EXTERIOR A PRESIÓN CONSTANTE

2

0
102 103 104 105 106 107
tD

Figura 6.3. pD vs. tD para un pozo en el centro del cuadrado.

6-11
 Área original de cada cuadrado

(x,y)
Area = 1
A

Figuras 6.4. Superposición de dos arreglos cuadrados para formar un rectangular

2:1.

6-12
 10

PD
8

6 Comienzo aproximado del flujo seudo-continuo.

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
tDA

Figura 6.5. PD vs. tD para un pozo en el centro de un cuadrado finito produciendo a

tasa de flujo constante.

ECUACIÓN BÁSICA DE PRUEBAS DE RESTAURACIÓN DE PRESIÓN.

Si un pozo se cierra por un tiempo ∆t , después de haber estado produciendo por

un tiempo (t) la caída de presión en el pozo (pws), vendrá dada por la solución de
la ecuación (6-16):

p ws = pi −
. qµB
1412
kh [
pD ( t + ∆t)D − pD ( ∆tD ) ] (6-16)

El efecto de daño no aparece en la ecuación (6-14) ni en la ecuación (6-16). Es

necesario emplear la siguiente ecuación, que define la presión de flujo en el pozo,
antes de cerrarse, para determinar el efecto de daño (S).

6-13
kh (pi − p wf )
= pD ( tD ) + S (6-17)
. qµB
1412

Restando la ecuación (6-16) de la ecuación (6-17), se obtiene:

kh
141.2qµB
( [
pi − p wf ) − ( pi − p ws ) =] kh
( p − p wf ) =
. qµB ws
1412

p D ( t D ) + S − p D ( t + ∆t ) D + p D ( ∆t D ) (6-18)

la cual puede ser resuelta para S:

S=
kh
. qµB
1412 [
(pws − pwf ) + pD ( t + ∆t)D − pD (∆t)D − pD ( tD ) ] (6-19)

Si ∆t << t (esto no necesariamente supone flujo semi-continuo)

pD ( t + ∆t)D ≈ p ( tD ) (6-20)
D

y la ecuación (6-19) puede escribirse :

kh
S=
1412
. qµB
(pws − pwf ) − pD ( ∆tD ) (6-21)

Un concepto muy útil es el de presión promedia volumétrica, p . Un balance

volumétrico para un yacimiento cerrado de área A, es:

Volumen producido del yacimiento = Expansión de los fluidos inicialmente en el

yacimiento.

Entonces:

5.615qBt = Ahφc( pi − p ) (6-22)

6-14
que puede ser escrita :

kh 2π ( 0.00634) kt
. qµB
1412
( pi − p ) =
φµcA
= 2π tDA (6-23)

La presión promedia volumétrica, también puede definirse como el límite de la

presión de cierre pws, cuando el tiempo de cierre tiende a infinito:

p = lim (pws )
(6-24)
∆t→∞

También p es la presión volumétrica promedio a cualquier tiempo; si:

vt A
1 1
p=
vt ∫ p dv =
A ∫ p dA (6-25)
0 0

El método de Muskat7 y el método de Miller-Dyes-Hutchinson8, consideran la

diferencia de presión ( p -pws), que puede ser originada de las ecuaciones (6-15) y
(6-24).

kh kh
141.2qµB
( pi − p ws ) − ( pi − p ) =
. qµB
1412
( p − p ws ) = (6-26)

p D ( t + ∆t ) D − p D ( ∆t D ) − 2π t DA

Horner9 sugirió que un gráfico de pws contra el [ ]

ln ( t + ∆t) / ∆t debería ser una

recta y podría ser extrapolada a una falsa presión, p*. Añadiendo y restando ½ ln
( t + ∆t)D al miembro derecho de la ecuación (6-15). (Figura 6.9)

kh
( p − p ws ) = p D ( t + ∆t )D − p D (∆t D ) +
141.2qµB i
(6-27)

6-15
 1
1 / 2 ln ( t + ∆t ) D − ln ( t + ∆t ) D
2

La ecuación (6-27) puede ser usada para expresar pD ( ∆tD ) para tiempos de

cierres cortos ∆t , y para muchas áreas de drenaje:

kh
[
( p − p ws ) = p D ( t + ∆t )D − 1 / 2 ln (∆t D ) + 0.80907 +
141.2qµB i ]
1 / 2 ln ( t + ∆t ) D − 1 / 2 ln ( t + ∆t ) D (6-28)

La ecuación (6-28) puede ser arreglada en la forma siguiente:

kh  t + ∆t 
141.2qµB
( pi − p ws ) = 1 / 2 ln 
 ∆t 
 + p D ( t + ∆t ) D

[
−1 / 2 ln ( t + ∆t ) D + 0.80907 ] (6-29)

Note que el suscrito D puede ser eliminado de la razón de tiempos.

Los últimos dos tiempos de la ecuación (6-28) se eliminarían si pD ( t + ∆t)D

pudiese ser representada por la aproximación logarítmica de la ecuación (6-18).

Bajo estas condiciones la ecuación (6-28) indica que pws sería una función lineal
de ln ( t + ∆t) / ∆t , para todo tiempo de cierre. Esto implicaría que la presión

falsa de Horner9, (p*), sería igual a la presión inicial (pi), cuando la producción
acumulada es más pequeña y no suficiente para reducir la presión media p
significativamente.

Supongamos ahora que ∆t es muy pequeño comparado con t.

6-16
Esto es ∆t << t . Observe que esta condición es casi siempre posible para
cualquier tiempo de producción. No se ha supuesto que el tiempo de producción
es muy grande, tal que el estado de flujo semi-continuo haya sido alcanzado.
Bajo esta condición:

t + ∆t ≈ t

y la ecuación (6-29) puede escribirse :

kh  t + ∆t 
141.2qµB
( pi − p ws ) = 1 / 2 ln 
 ∆t 
+

(
p D ( t D ) − 1 / 2 ln t D + 0.80907 ) (6-30)

La ecuación (6-30) es una línea recta en papel semi-log y dicha ecuación puede
ser evaluada a ∆t = ∞ .

Con:

t + ∆t
=1 y p ws ≈ p *
∆t

Entonces:

kh
. qµB
1412
(pi − p *) = pD ( tD ) −
1
2 [( ln (t
D ) + 0.80907 )] (6-31)

p* podría ser igual a pi, sólo cuando el tiempo de producción sea bastante
pequeño, de tal forma que pD (tD) pueda ser representado por la aproximación
semilog.

Matthews - Brons - Hazebroek 2 definieron una diferencia de presión en términos

de (p* - p ). Esta función se usa para obtener p , a partir de la falsa presión p*, y

6-17
puede ser obtenida a partir de las ecuaciones (6-23) y (6-30). Restando (6-30) de
(6-24)

Así:

kh kh
. qµB
1412
( pi − p ) − ( pi − p *) ( p − p) = 2π t DA −
. qµB i
1412

[
p D ( t D ) + 1 / 2 ln t D + 0.80907 ] (6-32)

La ecuación (6-31) puede ser usada para hallar la presión adimensional Matthews
- Brons - Hazebrock2 que fue definida como dos veces la ecuación (6-31).

p DMBH =
kh
70.65qµB
(
( p * − p) = 4 π t DA − 2 p D ( t D ) + ln t D + 0.80907 ) (6-33)

Si ( tD ) es pequeño, tal que, la aproximación semi-log de la ecuación (6-18) se

aplica a los dos últimos términos de la ecuación (6-32), estos se cancelan y la
función pDMBH sería igual a 4 π (tDA). Si (tD) es suficientemente grande, que pudiese
ser reemplazado por la ecuación (6-20), la ecuación (6-31) se convertirá en:

kh 4A
( p * − p) = 2πt DA − 1 / 2 ln 2
− 2 πt DA +
141.2qµB γCA rw

4At DA
1 / 2 ln = 1 / 2 (ln CA t DA ) (6-34)
2
γrw

donde γ es igual a la exponencial de la constante de Euler:

γ = exp (0.5772...) = 1.781.....

Para tDA ≥ comienzo del período de flujo semi -continuo.

6-18
Luego:

kh
70.65qµB
(p * − p ) = ln (C A tDA ) (6-35)

para la misma condición de tDA que la ecuación (6-34).

Brons y Miller5, presentaron ecuaciones para estados semi - continuos, para áreas
rectangulares en términos de la diferencia de presión ( p -pwf). Usando las
ecuaciones (6-17) y (6-23):

kh 1  4A 
( p − pwf ) = ln   +S
 (6-36)
. qµB
1412 2  γC Arw2 

Ahora utilizaremos las ecuaciones desarrolladas para analizar los métodos

convencionales de Análisis de Pruebas de Pozos.

MÉTODO DE MUSKAT7 EXTENDIDO.

El método extendido de Muskat7 emplea un gráfico semi-log de ( p − p ws ) vs ∆t .

La base del método es la ecuación (6-26).

kh
. qµB
1413
( p − pws ) = pD ( t + ∆ t)D − pD ( ∆ t)D − 2πtDA (6-26)

El análisis previo del método de Muskat7 había supuesto que el pozo produjo
hasta la condición de flujo semi - continuo.

6-19
Esta suposición puede removerse y usar la ecuación (6-26) para cualquier tiempo
de producción como puede verse en la Figura 6.6 donde todas la líneas se
convierten en línea rectas.

La Figura 6.6 presenta el gráfico extendido de Muskat7 para diversos tiempos de

producción, tDA como parámetro. Se requiere obtener una línea recta cuando el
valor correcto de p ha sido seleccionado. Las pendientes de las líneas rectas son
iguales y están relacionadas a la difusividad, mientras que el intercepto está
relacionado al producto (kh), en una función de tiempo de producción, excepto
para tDA>0.1.

La Figura 6.7 presenta un gráfico del intercepto como función de tiempo de

producción. Es claro que con el gráfico de Muskat7, el producto (kh) debe
determinarse para el tiempo de producción mayor a aquel al cual se alcanza flujo
semi-continuo o hacer las correcciones necesarias mediante la Figura 6.6.

La línea de trazos en la Figura 6.6 indica un límite inferior o comienzo de la línea

recta. El tiempo parece aumentar al disminuir el tiempo de producción y está
comprendido entre ∆tDA de 0.05 y 0.06. La línea de trazos de la derecha
representa un límite superior para la línea recta, el valor de ∆tDA está
comprendido entre 0.02 y 0.1 y parece decrecer al disminuir el tiempo de
producción.

7
El método de Muskat contiene también la posibilidad teórica de estimar ( φcA )
de la pendiente de la línea recta. (La pendiente de todas las líneas rectas es
0.056-1 log 10 ~.)

La ecuación de la pendiente es:

−0.113 k
φcA = (Pendiente, log 10 ~ días 1) (6-37)
µ

6-20
El efecto de daño puede ser determinado de la ecuación (6-16) después de
sustituir la expresión apropiada para pD (tD).

10

10-1

10-2

10-3

10-4
0 0.4 0.8 1.2
∆tDA

Figura 6.6. Gráfico generalizado de Muskat para un pozo localizado en el centro

de un cuadrado.

0.5

0.4

0.3

0.2

0.1

10-4 10-3 10-2 tDA 10-1 1 10

Figura 6.7. Intercepto del gráfico de Muskat, como función adimensional de

producción, para un pozo localizado en el centro de un cuadrado.

6-21
 0

4
10-5 10-4 10-3 10-2 10-1 1
∆tDA

Figura 6.8. Gráfico generalizado de Millar – Dyes –Hutchinson (MDH), para un

pozo localizado en el centro de un pozo cuadrado.

VARIAS CONCLUSIONES PUEDEN SUMARIZARSE CON RESPECTO AL

MÉTODO DE MUSKAT.

1. Es importante hacer chequeo del tiempo adimensional para determinar si el

tiempo de cierre fue lo suficientemente grande. A bajos tiempos de cierre, las
líneas son cóncavas hacía arriba y con la aplicación del método a estos datos,
daría valores bajo de p .

2. El tiempo de cierre para determinar algunas variables tal como (kh), efecto de
daño y presión estática, es bastante alto en comparación con el tiempo requerido
para determinar las mismas propiedades con otros métodos.

3. Solo una porción muy corta de la curva se vuelve línea recta (un tercio de log)
y para tiempo de producción pequeño, pws es muy aproximada a p , lo cual
podría dificultar la aplicación del método.

6-22
MÉTODO DE MILLER - DYES - HUTCHINSON8.

La ecuación (6-26) es la base para el análisis del método:

kh ( p − p ws )
= pD ( t + ∆t)D − pD ( ∆tD ) − 2πtDA (6-26)
. qµB
1412

El método establece que un gráfico de la presión de cierre, contra el logaritmo del

tiempo de cierre debe producir una línea recta cuya pendiente es proporcional al
producto (kh).
La presentación original de este método supone que el pozo ha producido hasta el
período semi - continuo. Sin embargo la ecuación (6-26) puede usarse para
generalizar el método.

La Figura 6.8 presenta el gráfico de restauración de presión para diversos tiempo

de producción. Todas las curvas mostraron la misma pendiente igual a, 1.151, a
tiempos cortos y de cierre.
La pendiente de 1.151 es la razón 2.032/2. Así que todas las curvas deberían dar
el valor apropiado de (kh) si se selecciona la línea recta correctamente. Sin
embargo, note la línea horizontal en la Fig. 6.8: esta línea indica el tiempo máximo
de cierre en el tope de la línea recta. La duración de la línea recta, depende del
tiempo de cierre si se incrementa al aumentar el tiempo de producción. La
transición entre la línea recta apropiada y las condiciones de presión estática
parecen aumentar al disminuir el tiempo de producción. Note que a pequeños
tiempos de producción la duración de la línea recta apropiada es muy corta y
cualquier fenómeno que obscurezca la línea recta (efecto de llene y/o daño),
podría resultar en un bajo valor de la pendiente y en consecuencia, en un alto
valor para (kh).

El valor de la permeabilidad se determina de la fórmula:

6-23
 qµB
kh = 162.6 (6-38)
mh

Luego, es posible determinar tDA y ∆tDA para determinar con la Figura 6.8, si los
datos de restauración de presión caen debajo de la línea horizontal.

Si los datos de restauración de presión caen debajo de la línea recta horizontal , la

Figura 6.9 puede usarse para determinar la presión estática, p . El procedimiento
consiste en calcular tDA para un tiempo de restauración sobre la línea recta y
después leer la presión adimensional correspondiente de la Figura 6.9. Luego, p ,
puede determinarse de la presión de cierre, pws , que corresponde a tDA .
Para tiempos cortos de producción, este procedimientos podría dar valores altos
de p .
El Skin puede ser determinado con las ecuaciones (6-16) y (6-17).

Donde:
p − p wf k 
.  1hora
S = 1151 − log + 3.23  (6-39)
 m φµc2w 

MÉTODO DE HORNER9.

La ecuación (6-29) es la ecuación básica de interpretación:

kh 1  t + ∆t 
. qµB
1412
( pi − p ws ) = ln
2  ∆t 
 + p D ( t + ∆t ) D

−
1
2 [
ln ( t + ∆t ) D + 0.80907 ]
(6-29)

6-24
La ecuación (6-29) puede usarse para preparar la Figura 6.8. Como puede
observarse, en dicha figura todos los tiempos de producción tDA < 0.06, originan
la misma línea recta y tienen la pendiente de 1.151 y la extrapolación es pi (a una
razón de tiempos iguales a uno). p Será diferente para todos los diferentes
tiempos de producción; esto puede ser determinado de las líneas horizontales a
razón de tiempos cercanos a la unidad. El valor de p viene dado por el balance
volumétrico de la ecuación (6-23).

kh
. qµB
1412
(pi − p ) = 2πtDA (6-23)

Es interesante notar que p se alcanza para todos los casos cuando ∆tDA = 0.1
para un pozo en el centro del cuadrado cerrado. Esto puede ser visto por la línea
punteada de la Figura 6.9.

4
Ramey y Cobb demostraron que el gráfico de Horner9 producía una línea recta
de mayor longitud que en el gráfico de MDH para diversos tiempos de producción
y concluyeron que para el caso del yacimiento cuadrado con un pozo en el centro,
el método de Horner14 es superior al método MDH.

Es posible demostrar que:

 t + ∆t 
ln
 ∆t  p

ws = p
(
= 4 πt DA − 2 p D ( t D ) + ln t D + 0.80907 ) (6-40)

y si el tiempo de producción es muy corto , tal que, la aproximación logarítmica se

cumpla para PD,

 t + ∆t 
  = exp ( 4πtDA ) (6-41)
 ∆t  p
ws = p

6-25
y si tD es muy grande tal que la aproximación de flujo semi - continuo se cumpla
para PD.

 t + ∆t 
  = CA tDA (6-2)
 ∆t  p
ws = p

6-26
 LISTA DE FIGURAS
CAPÍTULO 6

Figura 6.1.- Octante de un cuadrado mostrando el pozo y las localizaciones de

los puntos de presión.

Figura 6.2.- Superposición de un arreglo infinito de pozos para generar el

comportamiento de presión de un yacimiento cuadrado con un pozo en el centro.

Figura 6.3.- PD vs tDA para un pozo en el centro de un yacimiento cuadrado

finito produciendo a tasa de flujo constante.

Figura 6.4.- pD vs tD para un pozo en el centro de cuadrado.

Figura 6.5.- Superposición de dos arreglos cuadrados para formar un rectangular

2:1.

Figura 6.6.- Gráfico generalizado de Muskat para un pozo localizado en el centro

de un cuadrado.

Figura 6.7.- Intercepto del gráfico de Muskat, como función del tiempo
adimensional de producción, para un pozo localizado en el centro de un cuadrado.

Figura 6.8.- Gráfico generalizado de Miller-Dyes-Hutchinson (MDH) para un pozo

localizado en el centro de un cuadrado.

Figura. 6.9.- Gráfico de Horner generalizado, para un pozo en el centro de un

cuadrado.

6-27
 REFERENCIAS CAPITULO 6

1. Van Everdingen, A.F. and Hurst, W.: “The Application of the Laplace
Transformation to Flow Problems in Reservoir, “Trans. AIME (1949) 186, 305 -
324.

2. Matthews, C.S., Brons, F., and Hazebroek , P.: “A Method for Determination of
Average Pressure in a Bounded Reservoir“ Trans., AIME (1954) 201, 182-191.

3. Earlougher , R.C. , Jr. , Ramey , H.J. , Jr. , Miller , F. G. , and Mueller , T.D. :
“Pressure Distributions in Rectangular Reservoirs , JPT. (Fab. 1968) 199 - 208.

4. Ramey , H.J. , Jr. and Cobb , W. M . : “A General Pressure Buildup Theory for
a well in a closed Drainage Area” JPT (Dec. 1971) 1493, 1505; Trans., AIME, 251.

5. Brons, F. and Miller, W.C.: “A Simple Method for Correcting Spot Pressure
Readings, “ JPT. (Aug. -1961) 803 - 805 .

6. Matthews, C.S.: “Analysis Pressure Build-up and Flow Test Data”, JPT, 1961.
13(9), 862-870.

7. Muskat, M.: “Use of Data on the Build-up of Bottom Hole Pressure”, Trans,
AIME (1937) 123, 44-48.

8. Miller, C.C. -Dyes A-Hutchinson, C.A., Jr.: “The Estimation of Permeability and
Reservoirs Pressure from Bottom-Hole Pressure Build-up Characteristics”, Trans.,
AIME (1950) 189, 91-104.

9. Horner, D.R.: “Pressure Build-up in Well”, Proc., third World Pet. Cong., E.J.
Bill, Leiden (1951) II, 503.

6-28