Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

ANALIPPRES1

Descargar como ppt, pdf o txt
Descargar como ppt, pdf o txt
Está en la página 1de 603

ANALISIS DE PRUEBAS DE

PRESION.

DR. HEBER CINCO LEY.

ANALISIS MODERNO DE PRUEBAS DE


PRESION Y DATOS DE PRODUCCION
Dr. Hber Cinco Ley
Contenido
1. INTRODUCCION
- Antecedentes
- Pruebas de presin y caracterizacin de yacimientos
2. FLUJO DE FLUIDOS EN YACIMIENTOS
- Naturaleza del flujo en yacimientos
- Modelos bsicos de flujo y ecuaciones
. Ecuacin de difusin y condiciones de frontera
. Geometras de flujo en yacimientos
- Ecuaciones de flujo y grficas
. Flujo lineal, radial, esfrico y bilineal
. Flujo estacionario
. Flujo pseudo-estacionario
. Declinacin exponencial de flujo
APP001A

- Variables adimensionales
. Definicin y caractersticas
. Presin, tiempo, distancia y flujo
- Principio de superposicin
. Superposicin en espacio
. Superposicin en tiempo
3. EFECTO DEL POZO Y DE SUS VECINDADES
- Factor de dao y de pseudodao
. Invasin de fluidos
. Penetracin parcial, disparos y desviacin
. Flujo no-darciano
. fracturas hidrulicas
- Efectos de almacenamiento
. Expansin de fluidos
. Movimiento de nivel de lquido
. Comportamiento de presin en el pozo
- Efectos de inercia y de sgregacin de fluidos

APP001B

4. DIAGNOSTICO DE REGIMEN DE FLUJO


- Funciones de presin y de derivada
. Cambio de presin
. Funciones de primera y segunda derivada
- Grfica doble logartmica
5. ANALISIS DE AJUSTE DE CURVA TIPO
- Curvas tipo para un modelo de flujo
. Definicin de una curva tipo
. Curva tipo doble logartmica
- Ajuste de curva tipo
. Seleccin de curvas
. Estimacin de parmetros
6. METODOLOGIA GENERAL PARA ANALISIS DE PRUEBAS
DE PRUEBAS DE PRESION
- Interpretacin general de pruebas
. Suficiencia y consistencia de datos

APP001C

. Definicin del tipo de prueba


. Estrategia de anlisis
. Normalizacin de datos
. Diagnstico de flujo
. Seleccin del modelo de flujo
. Anlisis de curva tipo
. Grficas especializadas
. Estimacin de parmetros
. Validacin del modelo de flujo
. Informe
- Combinacin de informacin
7. PRUEBAS DE DECREMENTO DE PRESION
- Introduccin
- Interpretacin
. Normalizacin de datos
+ Normalizacin
+ Convolucin

APP001D

+ Deconvolucin
. Grficas de anlisis
- Pruebas multiflujo
. Propsito
. Interpretacin
- Pruebas de lmite de yacimiento
. Propsito
. Interpretacin
. Limitaciones
- Diseo y conduccin de una prueba
. Objetivo, duracin y condiciones
. Aspectos prcticos
8. PRUEBA DE INCREMENTO DE PRESION
- Introduccin
- Conceptos y ecuaciones
. Ecuaciones para la presin de cierre
. Radio de investigacin y rea de drene

APP001E

- Interpretacin
. Normalizacin de datos
+ Normalizacin del cambio de presin
+ Convolucin
+ Deconvolucin
+ Funcin de impulso
. Grficas de anlisis
+ Diagnstico de flujo
+ Estimacin de parmetros
. Presin inicial, promedio y dinmica
- Diseo y conduccin de una prueba
9. PRUEBAS DE INTERFERENCIA
- Introduccin
- Interpretacin
. Modelos de flujo
. Ajuste de curva tipo

APP001F

. Grficas especializadas
. Estimacin de parmetros
- Diseo y conduccin de una prueba
10.PRUEBAS ESPECIALES
- Introduccin
- Pruebas de formacin (DST)
. Tipos
. Interpretacin
. Diseo y conduccin
. Aspectos prcticos
- Multiprueba de formacin
. Conduccin
. Interpretacin
. Aspectos prcticos

APP001G

11.APLICACIONES DE PRUEBAS DE PRESION


- Introduccin
- Caracterizacin de yacimientos
- Evolucin de la presin en un yacimiento
- Evolucin de un proyecto de inyeccin
- Evaluacin de la estimulacin de un pozo
12.APLICACION DE LA COMPUTADORA EN
PRUEBAS DE PRESION
- Introduccin
- Adquisicin de datos
- Diseo e interpretacin de pruebas

APP001H

Optimizacin de Explotacin de
un Campo
Caracterizacin
Simulacin
Esquema Optimo

Caracterizacin de un
Yacimiento

Definicin:
Detectar y evaluar los elementos que
constituyen y afectan el comportamien
comportamie
de un yacimiento.
Tipos:
. Esttica
. Dinmica

Caracterizacin Esttica

Definicin:
Deteccin y evaluacin de los elemen
que constituyen un yacimiento.
Herramientas:
Datos Geofsicos
Datos Geolgicos
Registros de Pozos
Datos de laboratorio

Caracterizacin Dinmica
Definicin:
Deteccin y evaluacin de los
elementos
que afectan el comportamiento de
un
yacimiento.
Herramientas:
. Pruebas de presin
. Datos de produccin
. Registro de flujo y
temperatura
. Pruebas de trazadores

Caracterizacin Dinmica
Caracterizacin Esttica
. Pruebas de presin
. Datos de produccin
. Registro de flujo
. Pruebas de trazadores
. Registros de
temperatura
Modelo Dinmico del
Yacimiento

Elementos que Afectan el


Comportamiento de un
Yacimiento
*
*
*
*
*
*
*
*

Permeabilidad, Porosidad y Anisotropa


Fuerzas Capilares y Mojabilidad
Estratificacin
Fallas Geolgicas
Discordancias
Acuamientos
Fracturamiento
Compartamentalizacin

CASO 1 EVALUACION DE LA CAPACIDAD DE FLUJO

FIGURA 1

CASO 1 EVALUACION DE UN FRACTURAMIENTO HIDRAULICO

FIGURA 2

CASO 1
RESULTADOS
Prefrac

Posfrac

K = 0.115 md

K = 0.14 md

S = 1.8

xf = 664 pies
FCD = 22
kfbf = 2045 mdpie

CASO 2

DETECCION DE UNA FALLA CONDUCTIVA

-1
1

1/4
1

FIGURA 3

CASO 2 FALLA CONDUCTIVA


RESULTADOS

df

FCD

CASO 3 DETECCION DE CASQUETE DE GAS

FIGURA 9

CASO 3 DETECCION DE CASQUETE DE GAS

Gas

Petrleo

Frontera a presin
constante

CASO 4 DETECCION DE CONTACTO AGUA - PETROLEO

1/2
1
FLUJO LINEAL

-1/2
1
FRONTERA A
PRESION
CONSTANTE

FIGURA 10

CASO 4 DETECCION DE CONTACTO AGUA - PETROLEO

MODELO CONCEPTUAL

C A/P

Mecayucan - Modelo Dinmico Conceptual

Casquete

Acufero

MECAYUCAN 51
10

PERMEABILIDAD VERTICAL (MD)

3PUNTOS
0.1

5PUNTOS
Khor

0.01

0.001
1

13

17

21

25

29

33

37

41

45

49

53

57

61

MUESTRAS

65

69

73

77

81

85

89

93

97 101 105 109

MODELOS DE FLUJO PARA YACIMIENTOS


NATURALMENTE FRACTURADOS
Homogneo
Anisotropa
Zonas Mltiples
Canal Dominante (Fracturas, Fallas y
Cavernas)
Doble Permeabilidad
Doble (Mltiple) Porosidad

MODELOS DE FLUJO PARA YNFS (PRUEBAS DE PRESION)

SAL.
111

114

115
105

103

6000

408

429
125 5500

121
6400

109

107

101 B

120

127
129

428
448

123
145

5500

301 A

119

117

149

5200

426
422

446

147

447

468

6400

169

62

6000

488
6400

58

56

14 D

NOMENCLATURA.

5500

8
26

5200

43

44

67
65

63
62

6400

83

49
69

6000

45

22 A

24

6000

6400

PENETRACION PARCIAL.

47

42

RADIAL HOMOGENEO.

FLUJO RADIAL COMPUESTO.

27

23 A

DOBLE POROSIDAD.

FLUJO LINEAL O BILINEAL .

29

3
2A

18

438

12

25

6000
6400

7
5

14

16

15

5800

13 A

38
439

32

34

36 A

6000

33

54

189

459

53

466

5500

167

444

5500

89

348

AGUA FRIA 847


MODELOS DE FLUJO

Zona invadida por agua


de inyeccin
Fractura (Porcin
abierta)
Fractura (Porcin
cerrada)

AGUA FRIA 847


RESULTADOS DEL
ANALISIS
PERIODO DE CIERRE

MODELO
K (md)

PARAMETROS
S

(Xf)

Lrad (Pies)

w
1
1.9

1.8
3

2
1.94
3
2.4
4
2.8
5
1.8
6
1.8

FCD
7

-1.8

(3.96)

2.5

RADIAL
4

COMPUESTO

-1.45 (2.78)

2.8

2.2

-1.7

(3.58)

15

1.6

-3.8

(29.32)

85

2.15

-3.45 (20.64)

160

2.15

-3.83 (30.20)

310

3
1.5
1.3
FRACTURA VERTICAL DE
CONDUCTIVIDAD FINITA
1.3
CON ZONA DE
PERMEABILIDAD
Kda (md) REDUCIDA
bd (pies)

K (md)
2.3

Xf (Pies)

140

Sf

AGUA FRIA 847


RESULTADOS DEL
ANALISIS
ZONA DE DAO

Xf = 350
pies

ZONA DE
DAO

bd =18.46
pies

METODO BASADO EN
DEFINICION DE Sf

METODO BASADO EN FINAL


DE FLUJO LINEAL

Sf = 0.48

Telf= 0.35 hrs


6
psi-1

ct = 6x10-

Kd = 0.25 md

= 0.12

K = 1.7 md

Kd = 0.25 md
350 pies

Xf =

= 0.375 cp

bd = 18.44 pies

bd = 18.49 pies

EXPLORACION Y PRODUCCION
SUBDIRECCION REGION NORTE
DISTRITO POZA RICA

POZO COYOTES No. 184


HISTORIA DE PRODUCC1N

250

ACEITE; RGA; AGUA

200

01 ENERO 1998
NP= 80 176 BLS
GP= 119.55 MMPC
WP= 1214 BLS

150

100

50

AOS; MESES

ACEITE (BPD)
RGA (M3/M3)/10
AGUA (BPD)

COYOTES 184
1.0000

1/QO (1/BPD)

0.1000

Serie1

1/2
0.0100

FLUJO LINEAL

0.0010
1.0

10.0
TIEMPO (MESES)

100.0

MODELO DE FLUJO LINEAL

EXPLORACION Y PRODUCCION
SUBDIRECCION REGION NORTE
DISTRITO POZA RICA

POZO COYOTES No. 427


HISTORIA DE PRODUCC1N

300

250

01 ENERO 1998
NP= 346 750 BLS
GP= 409.54 MMPC
WP= 0 BLS

ACEITE; RGA; AGUA

200

150

100

50

AOS; MESES

ACEITE (BPD)
RGA (M3/M3)/10
AGUA (BPD)

COYOTES 427

1/QO (1/BPD)

0.1000

0.0100

Serie1

1/4

FLUJO BILINEAL

0.0010
0

10
TIEMPO (MESES)

100

1,000

MODELO DE FLUJO BILINEAL EN UN CANAL

Caracterizacin Dinmica de
Yacimientos
Metodologa:
. Control de Calidad de la Informacin
. Sincronizacin de Datos de Presin y
Produccin
. Correccin de Datos de Presin y
Produccin
. Diagnstico de Geometras de Flujo
. Estimacin de Parmetros del
Yacimiento
. Clculo de Volumen de Drene
. Deteccin de Interferencia entre Pozos
. Integracin del Modelo de Flujo

CASO 6 HISTORIA DE PRODUCCION Y PRESIONES MEDIDAS

FIGURA 5

CASO 6 SIMULACION DE PRUEBAS

FIGURA 7

ASO 4 PRUEBAS DE INCREMENTO


CASO 4 GRAFICA SEMILOGARITMICA
DE PRUEBAS DE INCREMENTO

FIGURA 6

35

1800

30

1600

25

1400

20

1200

15

1000

10

800

600

400

-5

200

-10
0

50000

100000

150000
TIEMPO (HRS)

200000

250000

0
300000

QL (BPD)

FACTOR DE DANO

ARENQUE 13A

S
QL

CASO 4 RESULTADOS

MODELO DE FLUJO :
RADIAL HOMOGENEO
AREA DE DRENE RECTANGULAR
( EMPUJE HIDRAULICO )
PERMEABILIDAD

K = 7.3 MD

DAO DEL POZO

S = -3.5 (VARIABLE)

PRESION INICIAL Pi = 8338 LB/PLG2

11
,

00
0

PI
ES

CACTUS 1
MODELO DE FLUJO

DOBLE PERMEABILIDAD

k1 = 9 md
S = -4.3
=ct h 1 / ct h t = 0.3
= (k h)1 / (k h)t = 0.28
= 3.7x10-7
pi = 6426 psi
El pozo siente los efectos de
interferencia de los pozos
vecinos.

YACIMIENTO
COMPARTAMENTALIZADO
BELLOTA 94

FALLA
CONDUCTIVA
SEGUNDO
COMPARTIMENTO

NOVILLERO 14

DELTA P / Q (PSI/MMPCD)

1000

EFECTOS DE
FRONTERA
DELP/Q

100

DELPC/Q

1/2
1

10
100

FLUJO LINEAL

1000

10000
TIEMPO (HORAS)

100000

NOVILLERO 14
GRAFICA DE FLUJO LINEAL
450

400

350

(Pi - Pwf)/Q (PSI/HR1/2)

300

250
DELP/Q
DELPC/Q

EFECTOS DE
FRONTERA

200

150

100

50

0
0

20

40

60
RAIZ (T) (HORAS1/2)

80

100

120

116
PIES

4500
PIES
252
PIES

121
PIES

L = 4621
pies

POZO FRACTURADO EN UN YACIMIENTO CON ARENAS


MULTIPLES CASO ARCOS 10

L-18

L-20-21-22

L-24
L-25

L-26

SIMULACION DEL COMPORTAMIENTO DEL POZO ARCOS 10


DATOS UTILIZADOS EN EL AJUSTE
ARENA
L-18
L-20-22
L-24
L-25
L-26

K (MD)
0.58
0.98
0.28
0.27
0.08

H (PIES) POROSIDAD
32.8
0.18
29.52
0.17
32.8
0.2
39.37
0.19
75.46
0.19

Sw
0.25
0.3
0.2
0.32
0.21

Xf (PIES)
600
600
650
620
550

FCD
20
16
30
33
52

Sf
0
0
0
0.007
0

D (1/MPCD)
1.30E-05
1.30E-05
2.00E-06
1.40E-06
1.40E-06

CASO 18
ARENA
L-18
L20-22
L-24
L-25
L-26

Pi (PSI)
7144.5
7106.73
9174.36
8962.07
8508.23

L1 (PIES)
300
300
300
300
300

AREA DE DRENE
L2 (PIES) L3 (PIES) L4 (PIES)
4500
600
4500
4500
600
4500
4500
2300
3000
4500
2300
3000
4500
2300
3000

K = 0.135 md
FCD = 60

Xf =450 pies

3100 pies

K = 0.135 md
FCD = 60

Xf =450 pies

3100 pies

K = 0.135 md
FCD = 60

Xf =450 pies
2000 pies

K = 0.135 md
FCD = 60

Xf =450 pies

3100 pies

K = 0.135 md
FCD = 60

Xf =450 pies
2000 pies

CAMPO ARCOS
1000.00

A-11

800.00

A-42

A-6

A-51

600.00

A-20
400.00

A-36
200.00

A-33D

Y (M)

A-10
0.00

A-10

A-13
A-52D

-200.00

A-55
-400.00

A-34

A-25

-600.00

A-81

A-75D
-800.00

-1000.00
-1000.00

-800.00

-600.00

-400.00

-200.00

0.00

200.00

400.00

600.00

800.00

1000.00

1200.00

K = 0.135 md
FCD = 60
Xf =450 pies

A-20

A-10
A-13

4000 pies

QG (MMPCD)

CULEBRA 600
18

9000

16

8000

14

7000

12

6000

10

5000

QGREAL
QGCE400
PWFREAL

4000

3000

2000

1000

0
0

200

400

600

800

1000
TIEMPO (DIAS)

1200

1400

1600

1800

0
2000

PWFCALCE400

CULEBRA 600
18

16

14

12

QG (MMPCD)

QGREAL
10

QGCE400
QGE600
QGE800
QGE1000

QGE1200

0
0

200

400

600

800

1000
TIEMPO (DIAS)

1200

1400

1600

1800

2000

POZO
FRACTURADO
K = 1 MD
Xf = 2000 PIES

JUJO 523
300

250

RS RGA (M3/M3)

200

RS

150

RGA

100

50

0
0

50

100

150

200

250
PRESION (KG/CM2)

300

350

400

450

500

2500
0
2700
0
2900
0

SAL

132
114

3100
0

201-A

111
101-B
121

301-A

143-A

115
105

103
125

123

147

145

109

107

127

129

169

408

429 428

53

466

468

52
54

488

189

426

446 422-A 444


447

448

149

167

3300
0

117

56
459

58

34
12

36

3500
0

38
439

16

3700
0

13-A
5

22-A

26

29
27

23-A

24

49
47
67
89

65

42
63

POZOS
ANALIZADOS

62

3900
0
4100
66000.
0
00

69

45
43

44

348

25

2-A
4

15

14-D

18
438

33
32

83

J.S. KIMMERIDGIANO

68000.
00

70000.
00

72000.
00

74000.
00

76000.
00

78000.
00

80000.
00

CONCLUSIONES
* La caracterizacin dinmica detecta los elementos y
evala los
parmetros que afectan el comportamiento de un
yacimiento y
determina cmo los fluidos se mueven bajo condiciones
de
explotacin.
* Este proceso se realiza analizando informacin tomada
bajo
condiciones
de flujo (dinmicas) en el medio tal como
datos
de
produccin, presin, trazadores, temperatura, flujo, etc.
* El problema de unicidad se resuelve combinando
informacin de
varias fuentes.

I. INTRODUCCION
Objetivo
Sealar:
* Importancia de las pruebas de
de presin en la caracterizacin
de yacimientos.
* Describir los tipos de pruebas de
presin, sus ventajas y desventajas.
* Analizar el desarrollo histrico
de las pruebas de presin.
APP005

OPTIMIZACION DE LA EXPLOTACION
DE UN YACIMIENTO

Caracterizacin

Simulacin de Comportamiento

Esquema Optimo de Explotacin

APP006

GEOFISICA

GEOLOGIA

CARACTERIZACION

DATOS DE POZO

LABORATORIO

APP007

GEOLOGOS
Ambiente de
depositacin
Petrografa
Paleontologa

VOLUMEN
ORIGINAL
RESERVAS
PRODUCCION
ACUIFERO

INGENIEROS
PETROLEROS
Anlisis de:
Registros
Muestras
Pruebas

GEOFISICOS
Interpretacin de
datos ssmicos

Fase de evaluacin
APP008

GEOLOGOS
Correlacin
Continuidad
Mapas
Secciones

Plan de explotacin
Localizacin de
Pozos y plataformas

INGENIEROS
PETROLEROS
Caracterizacin
Simulacin
Aspecto econmico

GEOFISICOS
Interpretacin
Continuidad
Fallas
Acufero

Fase de planeacin
APP009

P(t)

PRUEBA
DE
PRESION

ELEMENTO DE
PRESION

APP010

P VS t
q VS t
X

INFORMACION
ADICIONAL

MODELOS DE
INTERPRETACION
ESTADO DEL POZO

APP011

DATOS PARA ANALISIS DE


UNA PRUEBA DE PRESION
- PRESION VS TIEMPO
- PRODUCCION VS TIEMPO
- GOR, WOR
- TEMPERATURA VS TIEMPO
- CONDICIONES MECANICAS DEL POZO
- ANALISIS PVT DE LOS FLUIDOS
- REGISTRO DE FLUJO
- MUESTRAS DE ROCA
- DATOS GEOLOGICOS
- DATOS GEOFISICOS
- INFORMACION DE OTROS POZOS
APP012

PRUEBA DE PRESION
Medicin continua de la presin
de fondo y del caudal (gasto)
en un pozo

Estmulo

Yacimiento

Respuesta

APP013

Pruebas de Presin
Un Pozo
Estmulo
Respuesta

Yacimiento
?

Dos Pozos
Estmulo
Respuesta

Yacimiento
?

Varios Pozos
Estmulo

Yacimiento
?

Respuesta

Respuesta
Respuesta
Respuesta
Respuesta

1
2
3
4 n
APP014

PRUEBAS DE PRESION
TIPO

GASTO
q

PRESION
P
wf

1.- DECREMENTO
to
q

Pw

2.- INCREMENTO

tP
q

t0 t1 t2
q
0

4.- INYECCION

t
t

tp

Pw

tp
-q

t0 t1 t 2

t
Pw

-q

5.- ABATIMIENTO

P
wf

tP

tP

3.- Q VARIABLE

to

t
APP015

PRUEBAS DE PRESION
(CONT.)

TIPO

GASTO

PRESION
P
wf

P
wf

6.- PRESION CONSTANTE


0

P
wf

7.- POTENCIAL
t1 t

t1 t 2 t3 t 4

-q

t
t
P
w

t
q1

q2

t1

t
P
w

10.- MULTIPRUEBA
DE FORMACION

P
wf

8.- INTERFERENCIA
VERTICAL

9.- PRUEBA DE
FORMACION

t1

t2

t
APP016

PRUEBAS DE PRESION
(CONT.)
TIPO

GASTO

PRESION
P
wf

11.- PRUEBAS DE
ESCALERA
0

t 1 t2 t
3 t 4

Active Well

t 1-t 2 -t 3-t

Observation Wells

12.- INTERFERENCIA
t0

t
Active Well

13.- PULSOS

to

P
wf

t
Observation Wells

APP017

Pruebas de Decremento de Presin


(Abatimiento)
Drawdown Test

p wf

t0
p wf

p
i

t
APP018

Pruebas de Caudal Mltiple


pw f
q

vs t

q
q

t1

t1

pwf

t1

Dos gastos

t2 t

wf

t1

t2

Gastos Mltiples

APP019

Pruebas de Incremento de Presin


(Build up Test)

Pws vs t

pws

-p
pwf

pws
t

tp

APP020

Pruebas de Inyeccin
piny vs t

iny

q iny
t0
0

-q

piny

t0

APP021

Pruebas de Cierre ( Pozo Inyeccin )


(FALL-OFF TEST)

ws

q=0

q
0

t iny
-q
tiny

ws

t
APP022

Pruebas de Interferencia (Horizontal)


(INTERFERENCE TEST)

vs t

0
Zona de Estudio

p
p = f(t)

Activo

Observacin

t
APP023

Pruebas de Interferencia Vertical


Un Pozo

Seccin Activa

t0

=
=
=
=

=
=

X X

=
=

Registrador de
Presin

Seccin de Observacin

t
APP024

Pruebas de Interferencia Vertical


Dos Pozos
Pvs t

Pozo Activo

X X

t0
X X

Pozo de Observacin

t
APP025

Pruebas de Pulsos de Presin


Horizontal
p

q vs t

vs t

Pozo Activo

t
p
XX

XX

Pozo de
Observacin
APP026

Pruebas de Pulsos de Presin


Vertical

Pozo Activo

t
=
=

=
=

= X X=
=
=
Registro de Presin

Pozo de
Observacin
APP027

Pruebas de Formacin
(Drillstem Test)
Presin Atmosfrica
o del Colchn de Fluidos

Vlvula

X X

Elemento de Presin

t
APP028

Pruebas a Presin Constante

q vs t
p
wf

Pwf= cte.

t0

t0

X X

APP029

INTERFERENCIA VERTICAL (PRATS)


pw

APP030

Pruebas RFT- MDT

APP028

ANALISIS DE PRUEBAS
DE PRESION
DESARROLLOS
Perodo
1950-70

1970-80

Mtodo
Lnea recta
(Horner)
(MDH)

Caractersticas
Yacimiento
homogneo

Curva tipo
(Ramey)

Efecto de Pozo
y sus vecindades

1980-85

Curva tipo
con Parmetros

Pozo fracturado
Doble Porosidad

1984-90

Derivada
(Bourdet)

Yacimiento
heterogneo

1990-

Anlisis con
Computadora

Integracin
de Informacin
APP031

Produccin de un Pozo

q- o =

ct

koh
Bo

- Factor de
Resistencia
de Drene

req =

q- o =

wf

rw

- wf )

ln (

req
rw

- Factor de Dao
APP032

PRODUCCION DE UN POZO DE ACEITE

qo

k h ( p pwf )
= ---------------------------141.2 Bo o ( ln ( re/rw) + s )

Vlido para comportamiento de


mediano y largo plazo

CAUSAS DE BAJA PRODUCTIVIDAD

Baja capacidad de flujo de la formacin

qo

k h

( p pwf )
= ---------------------------141.2 Bo o ( ln ( re/rw) + s )

Remedio:

Fracturamiento hidrulico

CAUSAS DE BAJA PRODUCTIVIDAD

Baja energa disponible

qo

k h ( p pwf )
= ---------------------------141.2 Bo o ( ln ( re/rw) + s )

Remedios:

Sistemas artificiales
Mantenimiento de presin

CAUSAS DE BAJA PRODUCTIVIDAD

Alta viscosidad del fluido

qo

k h ( p pwf )
= ---------------------------141.2 Bo o ( ln ( re/rw) + s )

Remedio:

Recuperacin trmica
Reductores de viscosidad

CAUSAS DE BAJA PRODUCTIVIDAD

Dao a la formacin

qo

k h ( p pwf )
= ---------------------------141.2 Bo o ( ln ( re/rw) + s )

Remedio:

Estimulacin del pozo

CAUSAS DE BAJA PRODUCTIVIDAD

Area de drene irregular

qo

k h ( p pwf )
= ---------------------------141.2 Bo o ( ln ( re/rw) + s )
= 0.472

Remedio:

= 16.79

1
10

CAUSAS DE BAJA PRODUCTIVIDAD


DE UN POZO

* Baja capacidad de flujo kh


* Baja presin media del yacimiento
* Alto factor de dao
* Alta viscosidad del fluido
* Baja eficiencia de drene

APP033

PERFILES DE PRESION

Pw

t=0
t1 t
2
t3

rw

r
APP035

PRUEBAS DE PRESION
OBJETIVOS
* ESTIMAR LOS PARAMETROS DEL YACIMIENTO
* CALCULAR LA PRESION PROMEDIO DEL AREA DE DRENE
* DETECTAR LAS HETEROGENEDADES DEL YACIMIENTO
* HALLAR EL GRADO DE COMUNICACIN ENTRE ZONAS
DEL YACIMIENTO
* DETERMINAR EL ESTADO DE UN POZO (DAADO)
* ESTIMAR EL VOLUMEN POROSO DEL YACIMIENTO

APP036

PRUEBAS DE PRESION
OBJETIVOS
* ESTIMAR LAS CARACTERISTICAS DE UNA FRACTURA
QUE INTERSECTA AL POZO
* ESTIMAR LOS PARAMETROS DE DOBLE POROSIDAD DE
UNA FORMACION
* DETERMINAR LAS CONDICIONES DE ENTRADA DE AGUA
* CONFIRMAR LA PRESENCIA DE UN CASQUETE DE GAS
* ESTABLECER EL GRADO DE COMUNICACION DE VARIOS
YACIMIENTOS A TRAVES DE UN ACUIFERO COMUN
* ESTIMAR EL COEFICIENTE DE ALTA VELOCIDAD EN
POZOS DE GAS
APP037

PRUEBAS DE PRESION
OBJETIVOS
* ESTIMAR LOS FACTORES DE PSEUDO DAO (PENETRACION
PARCIAL, PERFORACION S, DESVIACION, FRACTURA, ETC.)
* ESTIMAR EL AVANCE DEL FRENTE DE DESPLAZAMIENTO EN
PROCESOS DE INYECCION.

APP038

COMENTARIOS
* Las pruebas de presin constituyen una
herramienta poderosa para la caracterizacin
de yacimientos.
* Existen diversos tipos de pruebas con
objetivos diferentes.
* La interpretacin confiable de una prueba se
logra mediante la combinacin de informacin
de diversas fuentes.
APP038A

ELEMENTOS QUE CONTROLAN


EL FLUJO DE FLUIDOS EN UN
YACIMIENTO

NIVEL
Microscpico

Macroscpico

Megascpico

- Distribucin de
Tamao de Poro

- Estratificacin

- Geometra del
Yacimiento

- Geometra de
Poro

- Variacin de la
Permeabilidad

- Espacio poroso
sin salida

- Distribucin de
Fracturas

- Sistemas de
Fracturas y
Fallas

- Microfracturas

APP039

GEOMETRIAS DE FLUJO

LINEAL

RADIAL

ESFERICO

APP040

GEOMETRIAS DE FLUJO

Flujo hacia un pozo


totalmente penetrante.

Flujo hacia un pozo


parcialmente penetrante

APP041

ECUACION DE DIFUSION
2

P=

Ct
k

P
t

ECUACION DIFERENCIAL EN
DERIVADAS PARCIALES LINEAL

CONDICIONES INICIALES Y DE FRONTERA

SOLUCION

p = f ( r, rw , q, k, , , c t , h, p i , t )
APP044

CONDICIONES INICIALES
p ( x, y, z,..., t=0 )= p i

CONDICIONES DE FRONTERA

Especificar:
- Produccin ( Caudal )
o
- Presin
APP045

GASTO CONSTANTE

k
T

q= Constante

Area A
Frontera

q = - k (A p)
n Frontera

Ecuacin de Darcy

p =
n

q
kA

Frontera

Frontera
APP046

GASTO CONSTANTE
k
q

p
Frontera

p i, t = 0
t1

t2

t3

p = cte
n
s

APP047

FRONTERA A PRESION
CONSTANTE
q(t)

Frontera
p = cte

p i, t = 0
t1

Frontera
p0

t2

t3

q = f (t)

APP048

FRONTERAS IMPERMEABLES
( Gasto constante = 0 )

NO FLUJO

p
CONDICION DE FRONTERA

=0
Frontera
APP049

YACIMIENTO INFINITO

Lim

p ( s, t ) = p

APP050

Compresibilidad
C= - 1
v

v
p

Compresibilidad
de la formacin

Cf = 1

Compresibilidad
Total

c t = c f+ so co + s g c g+ swc w

Definicin

Roca
Agua
Gas

-1

( psi )

-1
2

kg cm

Aceite
APP051

PARAMETROS DEL YACIMIENTO


DIFUSIVIDAD HIDRAULICA

K
Ct

T= K h

TRANSMISIBILIDAD

S=

CAPACIDAD DE ALMACENAMIENTO

Ct h

T
S
APP052

PERMEABILIDAD
6

K ( md ) 10 10 10 10

Acuferos

Bueno

Suelos

Rocas

Arena
Limpia

PERMEABLE

Grava
Limpia

PERMEABILIDAD

-2

-4

1 10 10

SEMIPERMEABLE

IMPERMEABLE

Pobre

No
Existe

Arena
muy
fina y
Arcillo

Limo

Yacimientos
Petroleros

Granito

Yacimientos
de Baja Permeabilidad
APP053

PROMEDIOS DE PERMEABILIDAD

k1
k1

k2

k1 > k 2

k2

k1
k2

k1

k2

k1

k2

Aritmtico

k1+ k 2

-k = -k =
A

Armnico

-kH=

2
1+ 1
k1 k2

Geomtrico

k1

k1

-k = k-G =
-kG =

k1 k 2

k k2... k
1

APP054

Ejemplo
k = 100 md
1

k =

k1 + k2
2

kH =

2
1
k1

1
k2

k = 5md.
2

= 100 + 5
2

1+ 1

100

= 52.5 md.

= .21 = 9.54 md.

k G = k 1 x k2 = 100 x 5 = 22.4 md.

APP055

FLUJO LINEAL
b

k, , , c

x=L

x=0

Frontera
externa

Frontera
interna
(Pozo)

p = ( x, t, p , k, , , ct, L, h, b, q )
i

APP057

DISTRIBUCION DE PRESION

p ( x, y, z, t ) = f ( x, y, z, t, k, , , c t, p ,... )
i

Es imposible presentar las soluciones


en forma grfica para diversos valores
de las variables independientes.

10 Parmetros
10 Valores para
cada parmetro

10

10

Casos

APP058

FLUJO RADIAL

re

k ct

rw

p = f ( p , k, c ,
t

, h, re, rw, q, t, r )
APP059

FLUJO ESFERICO

re

k ct

rw

p = f ( p , k, c ,
t

, re , rw, q, t, r )
APP060

VARIABLES ADIMENSIONALES

Objetivo: Eliminar la presencia


de variables del yacimiento en
la solucin
Definicin: Combinacin de
variables para formar grupos sin
dimensiones
Caractersticas:
Las variables adimensionales son
directamente proporcionales a las
variables reales.
APP061

VARIABLES ADIMENSIONALES

TIPOS
. Presin ( Cambio )
. Tiempo
. Distancia
. Gasto

APP062

VARIABLES ADIMENSIONALES

Caida de Presin ( Cambio )

Lineal:

pDL =

Radial:

pD =

Esfrico: pDsph

kbh p
qB L
L
kh p
qB
k rw p
=
sph q B
APP063

VARIABLES ADIMENSIONALES

Tiempo
Lineal:

tDL =

Radial:

tD =

Esfrico:

k t
2
ct L
k t
c t rw2

Espacio
Lineal:
Radial:
Esfrico:

xD= x

rD = r
rw
APP064

VARIABLES ADIMENSIONALES
p wf = cte.

q = f (t)

GASTO

Lineal:

q DL =

Radial:

qD =

Esfrico: qDsph=

L q (t)

kbh

q (t)

kh
sph

p0
p0

q (t)

k r w p0

p0 = p i - pwf = cte.
APP065

Sistema de Unidades
Variable

p
k
L,b,h,r

presin
permeabilidad

viscosidad

gasto
porosidad

c t compresibilidad
t

tiempo

Ingls
psi
md

kg/cm2
md

pie

cp

cp

STB/D
fraccin
psi

sph

-1

m3/D
fraccin
( kg/cm )-1
horas

horas
2.637X10

Mtrico

-4

3.489X10

141.2

19.03

887.2

119.58

70.6

9.52

-4

APP066

Ejemplo
k = 55 md

h = 95 pies

q = 600STB/D

B = 1.2
= 0.8cp.

rw= 0.25pie

-6

c t = 12X10 psi

= 0.11

pD = ?
tD = ?
Solucin

Flujo Radial
-1

p = 1 psi
t = 1 horas

Sistema Ingls
= 141.2

= 2.637X10

-4

APP067

pD = kh

p
qB

55 X 95 X 1
141.2X600X1.2X0.8

pD = 0.064
tD =

-4

kt

ct

r2

t D = 219,750 =

2.637x10 x55x1
-6

0.11x0.8x12x10 x(.25) 2

2.2x10

APP068

Ecuaciones de Flujo y Grficas

Geometra de Flujo
. Lineal
. Radial
. Esfrico
. Bilineal

Rgimen

. Transitorio
. Pseudoestacionario
. Estacionario
APP069

Condiciones de Frontera

. Pozo
+ Gasto constante
+ Presin constante
. Fronteras
+ Infinitas
+ Finitas cerradas
+ Presin constante
Condiciones Iniciales
. Presin inicial uniforme
p(s,t=0) = p

i
APP070

Flujo Lineal en un Yacimiento Infinito Hacia


un Pozo que Produce a Gasto Constante

b
h

x=0 (Pozo)

k, ct
x

oo

p = f ( x, t ) = ?
APP071

PRESION EN CUALQUIER PUNTO


(FLUJO LINEAL)

[2

kt
ct

qB

kbh

2
x
c
t
-(
)
e 4 k t - x erfc

1/2

c t x2
4 kt

p (x, t) =

erfc(x) - Funcin Error


Complementaria
APP072

PRESION EN EL POZO (X=0)

Lq B
bh

t
k ct

1/2

p (t) =

Area de flujo A

2.518 q B
(
p (t) =
W
k
A

ct

1/2

ct

Sistema Mtrico

8.128 q B
(
P (t) =
W
k
A

1/2

Sistema Ingls

1/2

1/2

APP073

FLUJO LINEAL (POZO)

p = m lf t

1/2

mlf

1
0

t 1/2

APP074

ESTIMACION DEL AREA DE FLUJO

Sistema Ingls

ct

Sistema Mtrico

2.518 q B
A=
mlf

(k

1/2

ct

(k

1/2

8.128 q B
A=
mlf

APP075

PRESION EN PUNTOS DE
OBSERVACION (X>0)
=2

t DL/ xD

4 t DL / x

2
D

- erfc (

1
2 tDL / x2

DL

xD

DL

Log p /x

-(

DL

p (x , t )

Log tDL / x D

APP076

APLICACION DE LA GRAFICA

p=?

x, t

2.

t=?

x,

DL

Log p /x

1.

1
Log TDL / xD2
APP077

APLICACIONES DEL MODELO DE FLUJO LINEAL

. Pozo Fracturado
. Arenas Lenticulares
. Pozos Horizontales
. Canales
. Yacimientos Fracturados

POZO FRACTURADO

A = 4 xfh
h

Pozo

Area de Fractura
Af = xf h f =

xf
Fractura

16.25 q B
4 (k

ct ) m lf

ARENAS LENTICULARES

Flujo Lineal

k2
k1
k1 >>>

k2
Area de Flujo
A=

16.25 q B
(k

1/2

c t )2 m lf

CANALES

b
h
Area de Flujo
A=bh=

8.12 q B
(kct )1/2 mlf

POZOS HORIZONTALES

h
Flujo Lineal

AREA DE FLUJO
A = 2 h Lw
L w- Longitud del intervalo abierto

YACIMIENTOS FRACTURADOS

Flujo Lineal
Fracturas

Flujo en un Yacimiento Lineal Finito

q = cte

k, c

h
x=0 (Pozo)

Frontera
Impermeable

x
x=L

p = f ( x, t ) = ?
APP084

Flujo Lineal en un Yacimiento Finito Cerrado


Comportamiento de Presin

. Tiempos Pequeos (tDL 0.25)

( p(x,t) )finito

t teia

= ( p(x,t) )infinito

0.25 ct L2
k

APP085

Flujo Lineal en un Yacimiento Finito Cerrado


Comportamiento de Presin

. Tiempos Largos (tDL 2.5)

p = mpss t + b*
m pss =

L q B
b h L ct
t t pss =

2.5 ct L2
k

tpss - Comienzo del flujo pseudoestacionario


APP086

Flujo Lineal en un Yacimiento Finito Cerrado


Comportamiento de Presin

p i , t=0

t1
t2
t3

t4

t5

p
= cte
t

Flujo pseudo-estacionario
x=0

x=L
APP087

GRAFICA DE FLUJO PSEUDO-ESTACIONARIO

p
w

m pss

t pss

0
V p= b h L

qB

c t m pss

APP088

Flujo en un Yacimiento Lineal Finito


con presin constante en la frontera

b
q = cte
x=0 (Pozo)

k, c

Presin
Constante

x
x=L

p = f ( x, t ) = ?
APP089

Flujo Lineal en un Yacimiento Finito


con Presin Constante en la Frontera
Comportamiento de Presin

. Tiempos Pequeos (t
(

p(x,t))Finito = (

0.25

t t eia =

DL

0.25)

p(x,t)) Infinito
2

ct L
k

t eia - Final de comportamiento


de Yacimiento Infinito
APP090*

Flujo Lineal en un Yacimiento Finito Cerrado


Comportamiento de Presin

. Tiempos Largos (t DL 2.5)

p(x,t)) Finito = Constante


Flujo Estacionario

APP091*

Flujo Lineal en un Yacimiento Finito


con Frontera a Presin Constante
Comportamiento de Presin

pi , t=0
t1
t2
t t 3
Flujo estacionario

x=0

x=L
APP092*

APLICACIONES

GAS

FALLA
PERMEABLE

APP093*

Flujo Lineal en un Yacimiento Infinito Hacia


un Pozo que Produce a Presin Constante

b
h

q(t)
Presin
Constante

k, ct

x=0 (Pozo)

oo

q=f(t)=?
APP094

Flujo Lineal en un Yacimiento Infinito Hacia


un Pozo que Produce a Presin Constante

q(t) =

bh
L

1/q(t) =

pw

ct k

1
t 1/2

bh

1/2

ct k

APP096

GRAFICA DE DECLINACION DEL GASTO


FLUJO LINEAL

q
m qlf
1

1/ t

APP097

POZO PRODUCIENDO CON PRESION


DE FONDO FLUYENDO CONSTANTE

q
pi
pwf

q(t)
p wf

t
APP098

FLUJO LINEAL HACIA UN POZO CON


PRESION DE FONDO CONTANTE EN
UN YACIMIENTO CERRADO
p =cte
wf

k c

q(t)
h
x=0

x=L

APP099

COMPORTAMIENTO DEL GASTO


Tiempos Pequeos t DL 0.25

q(t)|

finito

= q(t)|

infinito

Tiempos Grandes t DL 2.5


2

q DL = 2 e

t DL

APP100

DECLINACION EXPONENCIAL
2 k b h pw
q(t) =
BL

Log q(t) = Log

2kbh

kt
2

ctL
p

BL

k t
2.303 x 4

c t L2
APP101

GRAFICA DE DECLINACION EXPONENCIAL

Log q
b*qlf
1

m*qlf
tDL = 2.5

t
APP102

ESTIMACION DE PARAMETROS

Volumen Poroso Drenado

Vp =

L B b*qlf
18.424 ct pw m*qlf

APP103

FLUJO RADIAL

k c

rw
h

APP104

FLUJO RADIAL
Fuente Lineal

EXACTO

APROXIMADO
APP105

SOLUCION DE LINEA FUENTE


pD (r D , t D) =

1
2

E1 (

1
4t / r
D

2
D

E 1 (x) - Integral Exponencial


E1 (x) =
x

e- u du
u

APP106

SOLUCION DE LINEA FUENTE

Log p

Log t D /r D2
APP107

SOLUCION DE LINEA FUENTE

Log pD

rD = 1
2

20

Solucin
de Lnea
Fuente
25

Log t D /r D2
APP108

VALIDEZ DE LA SOLUCION
DE LA LINEA FUENTE
r 20 cualquier tD
D

rD = 1

para

t D 25

(pozo)

APP109

APROXIMACION LOGARITMICA

Para t D /r D2 5
p(rD , t D )
D

1
Ln(t D /rD 2 ) + 0.80907
2

APP110

APROXIMACION SEMILOGARITMICA

Pozo
1.151 q B
pw =

kh
Log t + Log ( k/

c t rw2 ) + 0.351

APP111

GRAFICA SEMILOGARITMICA

POZO
p

m=

1.151 q B
kh

1
Log t
APP112

GRAFICA SEMILOGARITMICA

Pozo de Observacin
m

p
1
m=

t D /r D2 = 5

1.151 q B
kh

Log t
APP113

YACIMIENTO FINITO CERRADO


q = constante
q

Frontera Impermeable
APP114

YACIMIENTO FINITO CERRADO


q = constante

Tiempos Cortos t t eia

p) Finito = (

p)

Infinito

eia - End of infinite acting


(Final del comportamiento
de yacimiento infinito)

APP115

YACIMIENTO FINITO CERRADO


q = constante
Tiempos Largos t t pss

pw =

qB
2kh

qB t
ct h A

Ln( A2 ) + Ln( 2.2458 ) + 2 s


CA
rw

pss - Beginning of pseudo-steady state


(Comienzo del flujo pseudoestacionario)
CA - Factor de eficiencia de drene
APP116

FLUJO PSEUDO-ESTACIONARIO
Efectos de fronteras estn presentes
en el comportamiento del pozo

t t pss
APP117

FACTORES DE FORMA
CA

Forma

60 0

1/3

t DAeia

t DApss

31.62

0.10

0.1

31.6

0.10

0.1

27.6

0.09

0.2

27.1

0.09

0,2

21.9

0.08

0.4

0.098

0.015

0.9

30.8828

0.09

0.1

APP116A

FACTORES DE FORMA
Forma

CA

t DAeia

t DApss

12.9851

0.03

0.7

4.5132

0.025

0.6

3.3351

0.01

0.7

21.8369

0.025

0.3

10.8374

0.025

0.4

4.5141

0.06

1.5

APP116B

FACTORES DE FORMA
CA

t DAeia

t DApss

2.0769

0.02

1.7

3.1573

0.005

0.4

0.5813

0.02

2.0

0.1109

0.005

3.0

Forma
1
2

APP116C

FACTORES DE FORMA
Forma

CA

t DAeia

t DApss

5.3790

0.01

0.8

2.6896

0.01

0.8

0.2318

0.03

4.0

0.1155

0.01

1.0

2.3606

0.025

1.0

5
APP116D

FLUJO PSEUDO-ESTACIONARIO

t=0

p
Flujo Pseudoestacionario

p
= constante
t
r

Frontera
cerrada
APP118

FLUJO PSEUDO-ESTACIONARIO
La presin declina de manera
uniforme en el yacimiento
CA , t eia y t pss dependen de la forma
y tamao del rea de drene y de la
posicin del pozo
La presin vara linealmente con
el tiempo

APP119

GRAFICA DE FLUJO PSEUDO-ESTACIONARIO

m*
1

b*

t pss

t
APP120

ESTIMACION DE PARAMETROS

Volumen poroso de drene

Vp =

hA =

qB
c t m*

Factor de eficiencia de drene

C A = f ( b*, m, s )
APP121

ESTIMACION DE PARAMETROS
m* b*
CA

eia

t eiaDA

t pss
t pssDA

V
Tabla
Forma de Area de Drene
Posicin del Pozo
APP122

PRODUCCION DE UN POZO A PRESION CONSTANTE


EN UN YACIMIENTO INFINITO (FLUJO RADIAL)

Para tD 8 x 10 4
1/q =

1.151
kh

Log t + Log ( k/

B
po

ct rw2 ) + 0.351

APP123

ANALISIS DE DATOS DE PRODUCCION

1/q

m
1
kh =

1.151
m

B
po

log t
APP124

PRODUCCION DE UN POZO A PRESION CONSTANTE


EN UN YACIMIENTO CERRADO (2-D)

q(t)

Tiempos pequeos t t eia


(q) Finito (q) Infinito
APP125

PRODUCCION DE UN POZO A PRESION CONSTANTE


EN UN YACIMIENTO CERRADO (2-D)

Tiempos Largos t t pss


q D=

2
A)
Ln( 2.2458
r w2 CA

Log q(t) = Log

e-

t DA

Ln( 2.2458 A )
r w2 CA

2kh

pw

BLn(2.2458 A )
r w2 CA
4

2.303 A c t

kt

Ln( 2.2458 A )
r w2 CA
APP126

GRAFICA DE DECLINACION EXPONENCIAL

bq
1

Log q

mq
t pss
t
APP127

ESTIMACION DE PARAMETROS

2
k
h
p
2.2458 A e
CA =
bq B
2
rw
A=

B
ct h

bq

pw m q

APP128

ESTIMACION DE PARAMETROS
Sistema Ingls

CA =

2.2458 A
r
A=

2
w

kh

pw

70.6 bq B

bq
2.438 B
ct h p w m q

APP129

FLUJO ESFERICO

Pozo

k ct

APP131

APROXIMACION PARA FLUJO ESFERICO

rw

Pozo de radio finito

Punto fuente
APP132

SOLUCION DE PUNTO FUENTE

PDsph =

p(r,t) =

1
rD

qB
kr

sph

rD
Erfc (
)
1/2
2 t D

Erfc

r
2

c t 1/2
)
kt

APP133

Log p Dsph rD

SOLUCION DE PUNTO FUENTE

Log t D / rD2
APP134

FLUJO ESFERICO
q

APLICACIONES

APP135

FLUJO ESFERICO
Presin en el pozo

wDsph

pw =

= 1-

qB
k rw

sph

1
(tD )1/2
-

sph q B

)1/2

( c t )1/2
k3/2 t1/2

3/2

rw - Radio de la esfera que representa


al pozo
APP136

GRAFICA DE FLUJO ESFERICO

b sph
1

msph

1/t

1/2
APP137

ESTIMACION DE PARAMETROS

k = (-

sph

qB

rw =

3/2

( ct )1/2 ) 2/3

)1/2 m sph

qB
k b sph

sph

APP138

FLUJO ESFERICO
Comentarios
El flujo esfrico en un yacimiento
infinito tiende hacia flujo estacionario
a tiempos grandes.
El pozo acta como una esfera

APP139

FLUJO BILINEAL
Flujo compresible lineal

k ct

kkff bf

Flujo incompresible lineal

kf >>> k
APP140

FLUJO BILINEAL
El flujo bilineal existe cuando :
Dos flujos lineales se superponen
El flujo en el medio de alta permea
bilidad es incompresible.
El flujo en el medio de baja permea
bilidad es compresible.
No existen efectos de frontera.

APP141

ECUACIONES DE FLUJO BILINEAL


pwD =

FCD =

kf b f
k xf

2.45
1/2
FCD

Conductividad adimensional
de la fractura
1/4

pw =

1/4

t Dxf

2.45

qB

1/2

h (k f bf ) (

t1/4
1/4

c t k)

APP143

GRAFICA DE FLUJO BILINEAL

m bf

p
1

t 1/4
APP144

APLICACIONES DEL MODELO


DE FLUJO BILINEAL
. Pozo Fracturado
. Pozos Horizontales
. Canales

APP142

POZO HIDRAULICAMENTE FRACTURADO

k c

k f bf

( k f b f) =

qB

h m bf (

1/4

ct k) 2

APP145

POZOS HORIZONTALES

k 1 h1
(

( k1 h1 ) = (

qB
L w m bf (

1/4

c t k)2

c t k) 2

APP146

POZO EN UN CANAL

(k 1 h 1 ) = (

qB
b mbf (

ct k)1/4
2

APP147

PRINCIPIO DE SUPERPOSICION
Ecuacin de Difusin
2

p =

ct

p
t

p = f(x, y, z, ..., t)
Ecuacin Diferencial en Derivadas Parciales
Ecuacin de Segundo Orden
Ecuacin Lineal
APP148

PRINCIPIO DE SUPERPOSICION
Si una EDDP tiene n soluciones independientes una
combinacin lineal de ellas es tambin una
solucin.
Si

p1 =

F1 ( x, ..., t)

Solucin 2

p2 =

F2 ( x, ..., t)

pn =

Fn ( x, ..., t)

Solucin 1

Solucin n

p =

n
i=1

Ci

F i ( x, ..., t )

Tambin es una Solucin.


APP149

PRINCIPIO DE SUPERPOSICION

Modos :

Espacio
Tiempo

Metodologa :
Superponer cadas de
de presin causadas
por distintos pozos

APP150

PRINCIPIO DE SUPERPOSICION
Funcin Influencia

p1 (t)

Cambio de presin en el yacimiento (pozo)


causada por produccin a gasto unitario.
La respuesta de presin correspondiente
a un pozo que produce a gasto constante
est dada por :
p (t) = q

p1 (t)
APP151

PRINCIPIO DE SUPERPOSICION
Superposicin en espacio
Consideremos n pozos
produciendo en un
yacimiento

La cada de presin en el
pozo j est dada por :

1 q1
2 q2
j

3 q3
4 q4

pj =
n qn

n
i=1

qi

1 i,j

pj = ?

APP152

PRINCIPIO DE SUPERPOSICION
Superposicin en tiempo
Consideremos un pozo
produciendo a gasto
variable

El gasto se puede aproximar


por escalones que representan el inicio de produccin
de pozos ficticios con gasto
q i - q i-1 en el tiempo t i

q2
q

p(t) = ?

qn

t- t n

q1

t- t 3
t- t 2
t- t 1

t
t

t1

t2

t3

tn t t
APP153

PRINCIPIO DE SUPERPOSICION
Superposicin en tiempo
La respuesta de presin a un tiempo t es la
suma de los efectos correspondiente a cada
pozo ficticio

p (t) =

n
i=1

( qi -qi-1 )

p (t-t )
1

APP154

PRINCIPIO DE SUPERPOSICION
Superposicin en tiempo
Si se considera una variacin contnua
del gasto se tiene :
p (t) =

t
q( )
0

p1 (t-

)d

* Integral de Duhamel
* Integral de Convolucin
* Integral de Superposicin
* Integral de Faltung
APP155

III. EFECTOS DEL POZO Y DE SUS


VECINDADES
Objetivo
Analizar
* Los efectos del pozo en el comportamiento
de presin
* Los efectos de las vecindades de un pozo.

APP156A

EFECTOS DEL POZO Y DE SUS VECINDADES


Los datos de presin medidos en un
pozo pueden estar afectados por :
Efectos de dao
- Invasin de fludos
- Penetracin parcial
- Desviacin del pozo
- Disparos (perforaciones)
- Alta velocidad de fludos
Efectos del pozo
- Almacenamiento
- Inercia
- Segregacin de fludos
APP156

EFECTOS DEL POZO Y DE SUS VECINDADES

Segregacin

Inercia
Almacenamiento

Disparos

Penetracin parcial

Dao
APP157

DAO POR INVASION


k

ks
rw
ks

rs
k

(p )
wf

id

(pwf )

p)

dao

real

rw

rs

Sin dao
Con dao
APP158

FACTOR DE DAO
Suposicin
(

FLUJO RADIAL

p) dao = Cada extra de presin

PD =

kh

qB

k h ( p) Dao
=s
qB
( p) Dao =

Factor de Dao

qB s
kh
APP159

FACTOR DE DAO
S =

p) dao
k h(
141.2 q B

El factor de dao representa la cada


extra de presin expresada en forma
adimensional.
S=0
S>0
S<0

no hay dao
Si hay dao
Estimulacin

Vlido para flujo radial


APP160

FACTOR DE DAO
Relacin de Hawkins
S = ln ( r s ) ( k

rw

ks

1)

rw k s k
rs

APP161

DAO EN EL POZO
RADIO EFECTIVO DEL POZO

rw' = rw e-S
EFICIENCIA DE FLUJO

qideal

ln
=

ln

req
rw

req
rw + S
APP162

Estimulacin
k
ks

(p wf ) real
(

(p wf ) id

rw

p )dao < 0

rs
ks > k
APP163

FLUJO RADIAL

k,
, c t

Zona Daada
k
-3.2275
ct r w2
(

[ Log t + Log(

pw = 162.6 q B
kh
+ ( ) dao

APP164

FLUJO RADIAL
pw
1
(

p )

w 1 hr

t = 1 hora

log t

p w ) 1hr
m

S = 1.151 [ (

162.6 q B
m
- Log

k
c t rw2

+ 3.2275

kh =

APP165

FLUJO RADIAL
p

Pozo Daado

p)dao

Pozo sin dao

Log t
kh =

162.6 q B
m
APP166

Ejemplo
Datos

Diagnstico

Flujo Radial

Grfica Semilog

m, ( pw ) 1 hr

75 psi
m=
Ciclo

( pw ) 1 hr = 120 psi

q = 1200 STB / D

B = 1.25

= 1.1 cp

h =190 pies

= 0.13

C t = 15 x 10

-6

psi

-1

r w = 0.29 pies
APP167

CALCULOS
kh =

162.6 q B

162.6 x 1200 x 1.25 x 1.1

75

kh = 3577 md pie
k o = 3577= 3577 = 18 md
190
h

18
120
S = 1.151 [ 75 - Log [ 0.13x1.1x15x10 x(0.29) 2 +3.2275

S = - 3.67
r'w = rw e-S = 0.29 e - (-3.67) = 11.38 pies

( p)dao =
=

141.2 q B

kh

S = m / 1.151 S

141.2 x 1200 x 1.25 x 1.1


(-3.67)= -239 psi
3577

APP168

CALCULOS
q
q ideal

Ln (
=
Ln (

req / r w )
req / r w ) + s

req = 1500 pies

r w = 0.29 pies

S = - 3.67
Tabla CA = 4.5132
=

7.055 / C A =

7.055 / 4.51

= 1.25
Ln (

req / rw ) = Ln (1.25 x 1500 / 0.29)

q
8.77
= 1.71
=
q ideal 8.77 - 3.67
APP169

DISTRIBUCION DE PRESION
ALREDEDOR DE UN POZO DAADO
q

t=0

rw

r
APP170

CAIDAS EXTRAS DE PRESION


1 - ZONA DAADA
2 - PENETRACION PARCIAL
3 - DESVIACION
4 - DISPAROS ( Perforaciones)
5 - FLUJO DE ALTA VEL.
Cada cada extra de presin se puede expresar
por un factor de dao
S - Factor de dao ( invasin )
S p , S , S disp. - Factores de pseudo dao

APP171

FACTOR DE PSEUDODAO DE
PENETRACION PARCIAL
Causa: Convergencia de lneas de flujo hacia zona
disparada

Z1
hw

h
kv
kh

ESTIMACION
- Mtodo de Papatzacos
APP172

Papatzacos
Sp = ( h - h w ) Ln ( h
2 rw
hw

k r / kz )

+ h / h w Ln { hw / h A - 1 }
2 + hw / h B - 1
A=

4h
4z 1 + h w

4h
B = 4z + 3h
w
1

kr = k h

kz = k v

APP173

Ejemplo
k v = kh

h = 600 pies

r w = 0.3 pies

hw=

90 pies

Terminacin en la parte superior de la formacin.


Z1 = 0
Solucin

Sp = ?

A=

4 x 600
4h
= 26.6
4Z 1 + h w = 4(0) +90

B=

4h
= 4 x 600 = 8.8
4Z 1 + 3h 4(0) + 3x 90

600 - 90 ) Ln {
x 600
Sp = (
2 x .3
90

kh / k v

600 ) Ln { 90 / 600 [ 26.6 - 1 ] 1/2 }


+ (
90
2 + 90 / 600 8.8 - 1
APP174

SOLUCION
S p = + 31.4
q
Ln ( req / rw )
=
qid
Ln ( req / rw ) + s
q
=
q id

7
7 + 31.4

0.18

APP175

Ejemplo 2
Datos

h = 350 pies

h w = 30 pies

r w = 0.29 pies Z1 = 25 pies


kr / k z = ? = 1, 2, 5, 10
Sp = ?
Solucin

h - h w = 350 - 30 = 10.66
hw
30
x 350 = 1895.8
h
=
2 x 0.29
2 rw
h / h w = 350/30 = 11.66
APP176

SOLUCION
h w /h
2 + h w /h
A=
B=

0.0857

4h
4z 1 + h

2 + 0.0857
=

= 0.041

4 x 350
4 x 25 + 30

= 10.76

4 x 350
4h
=
4 z 1 + 3h
4 x 25 + 3 x 30

S p = 10.66 Ln

( 1895.8

kr
kz

= 7.368

+ 11.66 Ln ( 0.041

10.76-1
7.368-1

S p = 10.66 Ln ( 1895.8 )

kr
kz

- 34.75
APP177

SOLUCION
kr
kz
1

+ 45.7

+ 49.4

+ 54.3

10

Sp

+ 58

APP178

Pozo Desviado

rw

Pozo totalmente penetrante


S = -(

41

2.06

- (

56

1.865

Log (

h )
100 rw

APP179

Ejemplo
Datos

= 24

rw = 0.29 pies
Solucin

h = 150 pies
S =?

150
24 )2.06 ( 24 )1.865
(
S = log
100 x 0.29
41
56
S = - 0.33 - 0.2059 x 0.71 = - 0.476

= 75

S = - ( 75 )2.06 - ( 75 )1.865 x 0.71


56
41
= - 3.47 - 1.72 x 0.71 = - 4.69
APP180

POZO DESVIADO PARCIALMENTE


PENETRANTE

hw

h
rw

+p

= ?

zw

METODO DE PAPATZACOS
z 1 = h - hw Cos(

) / 2 - zw
APP181A

FACTOR DE PSEUDODAO POR DISPAROS

ks

rs

ap

dp

rw

APP186

FACTOR DE PSEUDODAO DE DISPAROS


K z /k r

0.25
0.5
0.75
h r 1.
15 1.5
12 2.

8
6
5
4
3

90
120
6"
d = 12"

15
10
8
6
4
3
2 Sdisp
1
.5
0

180

ap Pulgs.
1
2
4

8 12

16

20

-.5

APP187

FACTOR DE DAO DE DISPAROS


E INVASION

ks

S d+disp
APP190

FACTORES DE PSEUDO DAO

S d+disp

Invasin + disparos

Desviacin + penetracin parcial

+p

APP191

DAO TOTAL
Vlido para el perodo
de flujo pseudoradial

Stotal = S

+p

+ (h/h w ) Sd+disp

Sd+disp = (hw /h) (S total - S


Prueba

+p

Correlaciones
APP192

DAO POR ALTA VELOCIDAD


Pozos de gas
S Total = S + D q
STotal
S

q
APP193

POZO HIDRAULICAMENTE FRACTURADO

K f bf

Xf

APP194

POZO HIDRAULICAMENTE FRACTURADO

K f b f - Conductividad de la Fractura
F CD , (K f b f) D ,... CONDUCTIVIDAD ADIMENSIONAL
FCD =
FCD 300

K fb f
KXf

S f = f (FCD , X f / rw )

CAIDA DE PRESION DENTRO


DE LA FRACTURA ES
DESPRECIABLE

FRACTURA DE CONDUCTIVIDAD INFINITA

APP194A

FRACTURA DE CONDUCTIVIDAD INFINITA.

r w' =

Xf
2

rw e - S f

S f = Ln ( 2 rw / x f )
Vlido para flujo pseudoradial

APP194B

Ejemplo
X f = 60 Pies
Sf=?

r w = 0.25 Pies
r w' = ?

= Ln 2 x 0.25 = - 4.78
60

r 'w = x f / 2 = 60 / 2 = 30 pies

APP194C

Almacenamiento

q
q
q sf
0

tp

APP195

ALMACENAMIENTO
q = qw + q sf

qsf

III

II

q
qw
q sf

APP196

COMPORTAMIENTO DE PRESION
p

nt
e
i
am

en
c
a
lm
to
A
n
n
e
i
i
S
am
n
e
ac
m
l
A
n
o
C
III
II

t
I

Periodo totalmente dominado por el


almacenamiento

II

Periodo de Transicin

III

Periodo libre de almacenamiento


APP197

Coeficiente de Almacenamiento C
C = Vw c
Volumen de fluido que hay que aadir o remover del pozo
para modificar la presin de fondo en una unidad.

L3
F/L2 ]

C [ Bbl/psi ]
C [ m3 / kg / cm2 ]
APP198

Periodo Dominado por Almacenamiento


pw

qBt
24 C

C = Coeficiente de Almacenamiento
pw

mws
1

ws = wellbore storage
0

C=

qB
24 m ws
APP199

De los datos del pozo


C = c Vw
c - Compresibilidad promedio de fluido
dentro del pozo
Vw - Volumen del pozo

VW
x

APP200

ALMACENAMIENTO CAUSADO POR


MOVIMIENTO DE NIVEL DE LIQUIDO
q

APP201

COEFICIENTE DE ALMACENAMIENTO CAUSADO


POR MOVIMIENTO DE NIVEL DE LIQUIDO
C= (

Vu =

Vu

g )
144 gc

Vol. de espacio anular por unidad


de longitud bbl / pie

= densidad lb / pie 3

g = aceleracin de la gravedad pie / seg

g = constante de conversin de unidades (32.17)


APP202

Flujo Radial
C - Coeficiente de Almacenamiento
CD - Coeficiente de Almacenamiento Adimensional
CD = 2
CD =

Periodo I

5.6146 C

c t h r w2
C

PD =

c t h r w2

(Ingls)

(Mtrico)

tD
CD
APP203

EFECTO DE ALMACENAMIENTO Y DAO


pw

( p)dao

C1

S=0

S>0

C2

C 2> C

Final del almacenamiento

t
APP204

Flujo Radial
Sin Almacenamiento
Sin dao

pwD = f ( tD )
Pozo con almacenamiento y dao

pwD = f ( tD , cD , s )

APP205

Final del Almacenamiento


Flujo Radial
Ramey

t ewsD = (60 + 3.5 S) C D


t ews =

Chen & Brigham

(200,000 + 12,000 S) C
(kh / )

(Ingls)

ewsD

t ews

= 50 C

e 0.14 S

0.14 S
170,000
C
e
=
(kh / )

(Ingls)

APP206

Flujo Radial
p

t ews

log t
ews - end of wellbore storage
APP207

Ejemplo
S = 10
k = 20 md
h = 150 pies
Flujo radial
t ews =

-2

C = 10 bbl/psi
= 1 cp.

(200,000+12000 S) C
( kh / )

-2
(200,000+12000x10
)
x10
t ews =
( 20 x 150 / 1 )

1.06 horas
APP208

t ews =

170,000 C e 0.14 S
(kh / )
-2

0.14x10

tews = 170,000x10 e
(20x150/1)
=

2.29 horas

Ramey

Chen-Brigham

s = 10

1.06

2.29

s = 20

1.46

9.32

APP209

Ejemplo
S = 10
t ews = ?

-2

C = 10

bbl / psi

2
3
4
k = .1, 1, 10, 10, 10 md

= 1 cp.
h = 150 pies
-2 0.14x10
170,000
x
10
e
tews =
k x 150 / 1
= 11.33 x 4.055 / k = 45.94
k
k (md)
.1
1
10 2
3
10
10 4

t ews (hrs)
459.4
45.9
0.45
0.045
0.0045
APP210

Efectos de inercia

APP211

Efectos de segregacin de gas en el pozo

Humping

t
APP212

GRAFICA ESPECIALIZADA

p
m
1

f(t)
APP213

GRAFICAS ESPECIALIZADAS
1. Lineal

p vs t

2. Radial

p vs Log t

3. Esfrico

p vs 1/ t

4. Pseudoestacionario

p vs t

5. Estacionario

p = cte

6. Almacenamiento

p vs t

7. Bilineal

p vs t 1/4
APP214

FLUJO LINEAL
pw =
A

16.25 q B
t1/2 pdao
b h (ct k)1/2

pw

m lf
1

m lf
1

pdao
0

t 1/2
bh=

qB
16.25
(ct k)1/2 mlf

APP215

FLUJO RADIAL
162.6 q B
kh

p =
w

Log t + Log

k
c r w2
t

-3.2275 + 0 .87 S
p

o
p

p) 1

t=1
1

( p)dao

Log t
APP216

ESTIMACION DE PARAMETROS
q
kh = 162.6 B
m

(pw)1hr

S = 1.151 (

m
(

- Log

m S
1.151

p)dano =
r w' = r W e

q
q

ideal

= ln (

r eq

rw

k
c r

+ 3.2275 )

t w

) / { ln (

r eq

rw

) +S }

APP217

FLUJO ESFERICO
sph

P =

qB

krwsph

sph

qB

1/ 2

3/ 2

1/ 2

( ct )

k 3/ 2 ( t ) 1/ 2

+prest

b sph

pw

1
m sph

p1

p
prest.

1/ t
APP218

ESTIMACION DE PARAMETROS
k=

(-

r wsph =

sph

sph

3/ 2

c t ) 1/ 2

qB
( m sph

2/ 3

qB

k b sph

Radio de esfera ideal


-1

{{

( rwsph )ideal

2
0.25 + ( r w ) k z
hw k r
= h w { ln{
- 0.5+ 0.25 + ( r w )2 k z
hw k r

0.5+

Rodrguez Nieto-Carter
APP219*

FLUJO BILINEAL
p

44.1 q B
1/4 + p
dano
=
t
1/ 2
1/
4
h f ( k f b f) (
ct k )
p
W

0
p

p)

m bf

m bf

dao

t 1/ 4
APP220

ESTIMACION DE PARAMETROS

h f ( k fb f )

(K
kf bf)

1/ 2

1/ 2

44.1 q B
mbf (

c f k )1/ 4

44.1q B
=
1/4
h m (
)
c
k
f
f
bf

APP221

FLUJO PSEUDOESTACIONARIO

0.23395 q B

pw =

m
2.303

hA

ln (

ct

A
2.2458
)
(
) + 2S
+
ln
2
CA
rw

APP222

GRAFICA DE FLUJO
PSEUDOESTACIONARIO

pw

m*

tpss

b*
0

t
APP223

ESTIMACION DE PARAMETROS
Vp

0.23395 q B
=
ct
m*

C A = 5.456

m
m*

t eia

2.303 (

p1hr - b* )

m
t pss

Forma del rea de drene


Posicin del pozo
APP224

ALMACENAMIENTO
pw =

qB t
24 C

pw
m ws

p
1

pc

t
0

t
APP225*

ESTIMACION DE PARAMETROS
qB
24 m

C =

ws

t prueba = t +
(

p )
w

p w

t
p

APP226

IV. DIAGNOSTICO DE REGIMEN DE FLUJO

Objetivo:
Analizar los mtodos para detectar los
diversos regmenes de flujo presentes
en una prueba de presin

APP227

DIAGNOSTICO DE FLUJO
Comportamiento de presin
. La geometra y el rgimen de flujo
definen la funcin del tiempo que
controla el cambio de presin
. Dados los datos de presin se requiere
hallar la geometra y el rgimen de flujo
que dominan la prueba.

APP228

DIAGNOSTICO DE FLUJO
Datos de presin
Diagnstico de flujo
Aplicacin de grficas especializadas
Conformacin del modelo de flujo

APP229

DIAGNOSTICO DE FLUJO
Herramienta: Funcin de derivada
Bourdet
t

p'

tiempo transcurrido durante


la prueba

p' derivada de cambio de presin


durante la prueba

APP230

DIAGNOSTICO DE FLUJO
Bourdet
Funcin de derivada
* La funcin de derivada
es la derivada con respecto al logaritmo natural del tiempo.
* La funcin de derivada
es proporcional a la
pendiente semilogartmica.

p' = d

p / d Ln t

Ln t

APP230A

Forma General
t

p = c t

Tipo de flujo
Almacenamiento
Pseudoestacionario
Lineal
Bilineal
Radial
Esfrico

n
1
1

0
-

APP231

DIAGNOSTICO DE FLUJO
t
Log t

p = c t

p' = Log c + n Log t

Log
n

p'

Log t
APP232

DIAGNOSTICO DE FLUJO
Almacenamiento
Pseudoestacionario
1

Lineal

Log t

Bilineal

Radial

Esfrico

Log t
APP233

TRAZO DE PENDIENTES
102
1

-
10
1

p
(psi)
1

10-1
10-2

10 -1

10

102

t (hrs)
APP234

Almacenamiento

Log

p
1

Log t

p'
1

Log t
APP234A

Flujo Lineal

Log

p
'

Log t d p'
dt

p ) dao

Log 2

Log t
APP235

Flujo Radial

Log

p
'

p
Log t d '
dt

f(s)

kh

Log t
APP236

Flujo Esfrico

log

log t

l
p

dt

1
-

Log t
APP237

Flujo Bilineal

dao (+)

Log

Log t d
dt

'
p'

Log 4

1
(-)

Log t
APP238

Flujo Pseudo-estacionario

Log
Log t

p'
1

Log t
APP239

Pozo Parcialmente Penetrante

Radial
Esfrico
Pseudo
Radial

APP240

Pozo parcialmente penetrante

Log
Radial

p
t

p'
1

Almacenamiento

t ews

ter

t bsph

Esfrico

-
1

t esph

t bpr
Pseudo
radial

log t
APP241

DIAGNOSTICO DE FLUJO
Comentarios
* El comportamiento de presin de un
pozo puede exhibir varios tipos de flujo.
* Existen perodos de transicin entre
perodos que pueden ser expresados
por funciones simples de tiempo
(radial, lineal, bilineal, esfrico, pseudo
estacionario, etc.)

APP241A

Funcin de 2a. Derivada


t

p" | = c t

0
-

Tipo de flujo
Almac. y Pseudoest.
Lineal
Bilineal
Radial
Esfrico

APP242

Grfica de diagnstico de flujo


con la segunda derivada
Lineal

t|

p" |
l

Bilineal

Radial
1

-
Esfrico
Log t
APP243

COMENTARIOS
* El diagnstico de flujo es una etapa indispensable
para lograr un anlisis confiable de una prueba de
presin.
* Las funciones de primera y de segunda derivada
constituyen herramientas confiables de diagnstico.
* Los tipos de flujo que afectan a una prueba ocurren
a diversos tiempos, por consiguiente no se
translapan.
APP244

V. AJUSTE DE CURVA TIPO

Objetivos:
* Presentar el mtodo de ajuste de curva
tipo, sus aplicaciones y limitaciones.
* Examinar las diversas curvas tipo
disponibles.

APPV01

Ajuste de Curva Tipo


Comentarios
* Las grficas especializadas nicamente se
aplican a la porcin de los datos de una prueba
que pueden se representados por una funcin
simple de tiempo.
* El comportamiento de algunos modelos de flujo
est dado por funciones complejas del tiempo.
* Es necesario un mtodo para analizar la totalidad
de los datos de una prueba simultneamente
incluyendo los perodos de transicin.
APPV02

Curva Tipo
Definicin:
Grfica que representa el comportamiento de presin
en un pozo o en un punto de observacin expresada
en trminos de
variables adimensionales;
generalmente se usan escalas logartmicas.

Log F1 (p D )

Log F2 (t D)
APPV03

Ajuste de Curva Tipo

Log

Log

kh
qB

Log pD
ct r 2
Log

ct r2

2
t
/r
Log D D

Log t

Las curvas tienen la misma forma


APPV05

Ajuste de Curva Tipo


Procedimiento
1. Seleccionar la curva tipo
2. Graficar datos de la prueba en un papel
semitransparente usando la escala de la
curva tipo.
3. Ajustar datos a la curva tipo deslizando
la hoja con datos sobre la curva tipo .
4. Seleccionar un punto de ajuste.
5. Estimar los parmetros usando el punto de
ajuste y las definiciones de las variables
adimensionales que representan los ejes de
la curva tipo.
APPV06

Ajuste de Curva Tipo


Prueba de interferencia
Paso 1
Seleccionar la curva tipo
- Flujo radial en medio homogneo
- Flujo lineal en medio homogneo
- Flujo esfrico en medio homogneo
- Flujo radial afectado por falla
- Flujo radial en medio de doble porosidad
APPV07

Ajuste de Curva Tipo


Prueba de interferencia
Paso 2

Graficar datos en papel


semitransparente

10

10
p (psi)

SOLUCION DE
LINEA FUENTE

.1

t (hrs)

100

100 t /r
D

2
D

1000

APPV08

Ajuste de Curva Tipo


Prueba de interferencia
Paso 3

Ajustar datos con la curva tipo

10
p

10
p (psi)

1
SOLUCION DE
LINEA FUENTE

.1
.1

100

t (hrs)
.1

10

100 tD/r D2 1000


APPV09

Ajuste de Curva Tipo


Prueba de interferencia
Paso 4

Seleccionar el punto de ajuste

10
p

10
p (psi)

1
M

SOLUCION DE
LINEA FUENTE

.1
.1

Punto de ajuste

100

t (hrs)
.1

10

100 t /r 2 1000
D D
APPV10

Ajuste de Curva Tipo


Prueba de interferencia
Paso 5

Estimar parmetros

Datos del punto de ajuste:

p) M (t ) M ( p D ) M ( t D / rD2 ) M

Definicin de variables adimensionales:

kh p
p =
qB
D

t D / r D=

kt

ct r 2

Estimacin de parmetros :
kh=

qB
(

(p ) M
D

p )M

ct =

k (t) M
r 2 (tD / r D2 )M
APPV11

Curvas Tipo
CARACTERISTICAS DE UNA BUENA CURVA TIPO

* La curva debe poseer una forma con


curvatura caracterstica.
* En caso de una familia de curvas,
stas deber emerger de o converger a
una curva comn.

APPV11A

Curvas Tipo

Log p

Pozo con Almacenamiento y Dao


(Flujo Radial)

2s

Log t D p'

CD e

Log t D/C

APPV12

Curvas Tipo
1. Flujo radial con almacenamiento y dao
2. Prueba de interferencia
(Flujo lineal, radial y esfrico)
3. Prueba de 1 pulso (Flujo Radial)
4. Prueba de 1 pulso (Flujo Lineal)
5. Prueba de 1 pulso (Flujo Esfrico)
6. Yacimiento de doble porosidad
(Pozo, Modelo de Flujo Transitorio)
7. Yacimiento de doble porosidad
(Pozo, Modelo de Flujo Pseudoestacionario)
8. Yacimiento de doble porosidad
(Interferencia, Modelo de Flujo Transitorio)
9. Yacimiento de doble porosidad
(Interferencia, Modelo de Flujo Pseudoestacionario)
10. Pozo Hidrulicamente fracturado
11. Pozo cercano a una falla
APPV12A

Curvas Tipo

2s

p'

CDe

Log t Dp'

Log p

Pozo con Almacenamiento y Dao


(Flujo Radial)

Log t D/C

t (hrs)
APPV14

Curvas Tipo
Flujo Lineal, Radial y Esfrico
Punto de Observacin
Lineal

Esfrico

Log F (p )

Radial

Log F 2 (tD )

APPV16

Curvas Tipo
Prueba de un solo Pulso
Flujo Lineal
t pD /xD2

Log F (p )

Lineal

Log F2(t D )
APPV17

Curvas Tipo
Prueba de un solo pulso
(Flujo Radial)

Log F1 (p )

tpD /r 2D

Log F2(t ) D
APPV18

Curvas Tipo
Prueba de un solo pulso
(Flujo Esfrico)

Log F 1 (p D)

t pD /r D2

Log F2 (t D )

APPV19

Curvas Tipo

Log p

Pozo en Yacimiento Naturalmente


Fracturado (Flujo Transitorio)

Log t p'

CD e2S

CD /(1-

'

)2

Log t D / CD
APPV20

Curvas Tipo

Log p

Pozo en Yacimiento Naturalmente


Fracturado (Flujo Pseudoestacionario)

Log tD p'D

CD e2S

e-2S

CD /(1 -

CD /(1-

Log t D /CD
APPV21

Curvas Tipo
Interferencia en Yacimiento Naturalmente
Fracturado (Flujo Transitorio)

Log p

rD

Log t D / rD 2
APPV22

Curvas Tipo

Log p

Interferencia en Yacimiento Naturalmente


Fracturado (Flujo Pseudoestacionario)

r2
D

Log tD / rD2
APPV23

Curvas Tipo
Pozo Hidrulicamente Fracturado
(Fracturas Largas)
Log t D p'D F
Log p F
D CD
CD

FCD

FCD
Log t

Dxf CD
APPV24

Curvas Tipo
Pozo Hidrulicamente Fracturado
(Fracturas Cortas)

Log p

FCD

Log tD p'

FCD

Log t

Drw'
APPV25

Curvas Tipo
Flujo Bilineal con Almacenamiento
y Dao

Log F (p ')

Log F (p )

F4 (Sf )

F4 (Sf )
Log F2 ( tDxf )
APPV26

Curvas Tipo
Pozo cercano a una Falla Conductiva
Falla impermeable

Log t p'

Sf

Falla a
Presin
Constante

-1
1
Log t

FCD

1/4

Ddf

APPV27

VI. METODOLOGIA GENERAL


DE INTERPRETACION
Objetivo:
* Presentar y discutir una metodologa
general para analizar pruebas de
presin, de tal forma que se produzcan resultados confiables.

APPVI01

HERRAMIENTAS DISPONIBLES PARA


LA INTERPRETACION
* Grficas especializadas
* Curvas tipo
* Grficas de diagnstico de flujo

APPVI02

METODOLOGIA GENERAL
DE INTERPRETACION
DIAGNOSTICO DE PRUEBA
FILTRADO DE DATOS
NORMALIZACION
DIAGNOSTICO DE FLUJOS
CONFORMACION DEL MODELO
AJUSTE DE CURVA TIPO
GRAFICOS ESPECIALIZADOS
ESTIMACION DE PARAMETROS
VALIDACION DE MODELO
INFORME

APPVI03

DIAGNOSTICO DE PRUEBA

* Consistencia de datos
* Tipo de prueba
* Condiciones de prueba
* Estrategia de interpretacin

APPVI04

CONSISTENCIA DE DATOS

q
p

Tiempo
APPVI05

CONDICIONES DE PRUEBA
Mediciones

q3

q4
q1

q2
t1

t2

t3

t4

t5

Tiempo
APPVI06

TIPO DE PRUEBA Y CONDICIONES

* Prueba de incremento
* Flujo variable antes de la prueba
* Los datos de presin no solo dependen del gasto estabilizado mantenido antes de la prueba.

APPVI07

ESTRATEGIA DE INTERPRETACION

* Determinar qu perodos de flujo


afectan drsticamente a la prueba
* Diagnosticar las posibilidades de
interpretacin considerando perodos
de flujo y tiempos de duracin
* En caso de conocer el modelo de flujo
aplicable utilizar la historia de flujos
completa.
APPVI08

EVALUACION DEL EFECTO DE LA HISTORIA


DE FLUJO EN UNA PRUEBA

Datos de Flujo
y Presin

Modelo(s)
de Flujo

Superposicin en tiempo

Efecto Relativo de
cada Perodo de Flujo
APPVI09

EVALUACION DEL EFECTO DE LA HISTORIA


DE FLUJO EN UNA PRUEBA

t1 q 1
t2 q 2
t3 q 3

p (t) =
ws

(qi-qi-1)

p1 (t-ti)

i=1

tn q n
Efecto sobre la variacin
del Cambio de Presin:

p'(t) =
ws

n
i=1

(qi-qi-1)

p1'(t-ti)
APPVI10

EVALUACION DEL EFECTO DE LA HISTORIA


DE FLUJO EN UNA PRUEBA
Mediciones

q3

Efecto de historia de
produccin en prueba
de incremento

q4
q1

t1

p'(t)
=
ws

q2
t2

t3

(qi-qi-1)

t4

t5

Tiempo

p1'(t-ti)

i=1
APPVI11

EVALUACION DEL EFECTO DE LA HISTORIA


DE FLUJO EN UNA PRUEBA

p'(t) =
ws

(qi-qi-1)

p1'(t-ti)

i=1

Flujo radial

p1'(t) =
p'(t) =
ws

m1

1
t

2.303
m1

2.303

i=1

1
(q i-qi-1)
(t-ti)
APPVI12

EVALUACION DEL EFECTO DE LA HISTORIA


DE FLUJO EN UNA PRUEBA

Ejemplo

Estimar el efecto de cada perodo


de flujo en una prueba de incremento

t1 = 0 hrs
t2 = 24 hrs
t3 = 96 hrs
t4 = 192 hrs
t5 = 216 hrs

para

q1 = 650 STB/D
q2 = 0
q3 = 1350 STB/D
q4 = 1200 STB/D
q5 = 0

Flujo radial

t = 1, 2, 12, 24 hrs.

APPVI13

EVALUACION DEL EFECTO DE LA HISTORIA


DE FLUJO EN UNA PRUEBA

Solucin

p'(t) =
ws

m1

2.303

i=1

qi-qi-1

1
(qi-qi-1)
(t-ti)

t - ti

217

2.99 -3.37 11.15 -6.00 -1200.

218

2.98

-3.35 11.06 -5.77

-600.

12

228

2.85 -3.18 10.22 -4.17

-100.

24

240

2.71 -3.01

9.37 -3.12

-50.
APPVI14

FILTRADO DE DATOS

Los datos de presin medidos en un


pozo estn afectados por:
* Ruido generado en el yacimiento
* Ruido causado por la herramienta
* Efectos de tendencia de presin
* Efectos de marea.

APPVI15

FILTRADO DE DATOS
Ruido en mediciones
Datos originales

Suavizacin

Anlisis de datos
APPVI16

FILTRADO DE DATOS

p vs t
Suavizar datos

Calcular derivada

Calcular derivada
(Suavizada)

p' vs t
APPVI17

FILTRADO DE DATOS
Suavizacin de datos

t
APPVI18

FILTRADO DE DATOS
Suavizacin de datos

Ventana de
Suavizacin

t
psuav = ?

t
APPVI19

FILTRADO DE DATOS
Suavizacin de datos

Promedio Mvil

p (t) =
suav

t+

t /2

p( ) d
t-

t /2

APPVI20

FILTRADO DE DATOS
Suavizacin de datos
Caso Discreto

Ventana de
Suavizacin

psuav

i-2 i-1

i+1 i+2

Ventana de suavizacin N puntos

t
APPVI21

FILTRADO DE DATOS
Suavizacin de datos
Caso Discreto
i+(N-1)/2

psuav =

1
N

pj
j= i+(N-1)/2

N Ventana de suavizacin
(Impar)
APPVI22

FILTRADO DE DATOS
Suavizacin de datos
* Prueba de un solo pozo
Usar escala logartmica de tiempo
* Prueba multipozos
Usar escala normal del tiempo

APPVI23

FILTRADO DE DATOS
Estimacin de la derivada suavizada
Ventana de
Diferenciacin

p'suav = ?

p'j+1/2 =

pj+1

j+1/2

j+1

- pj
t

APPVI24

FILTRADO DE DATOS
Efecto de la Tendencia de Presin

1
pmed

mtend
pi
pprueba

(p)prueba

to

t
APPVI25

FILTRADO DE DATOS
Determinacin de la Tendencia de Presin
* Medicin
* Anlisis de datos
Correccin de la Respuesta de Presin

pprueba = pi - pmed (t) - mtend (t-to )


APPVI26

FILTRADO DE DATOS
Efecto de Mareas

pmed

to

Hora y Fecha

Inicio de prueba

t
APPVI27

FILTRADO DE DATOS
Efecto de Mareas

pmed
pcorr
hmarea
Hora y Fecha

pcorr = pmed - M hmarea


APPVI28

FILTRADO DE DATOS

COMENTARIOS
* El filtrado de datos se requiere en datos
de pruebas en yacimientos de alta permeabilidad y en pruebas de interferencia.
* El anlisis de datos no filtrados puede
producir resultados errneos.

APPVI29

NORMALIZACION

COMENTARIOS

* Las tcnicas de anlisis y de diagnstico


son aplicables a pruebas realizadas con
un cambio de gasto (caudal) constante.
* En la prctica las pruebas se realizan
bajo condiciones de gasto variable.

APPVI30

NORMALIZACION

NORMALIZACION
Estimacin de la respuesta de presin
correspondiente a un gasto constante
(unitario).

TECNICAS
* Deconvolucin
* Convolucin.
APPVI31

NORMALIZACION

DECONVOLUCION
* Calcula la respuesta de presin para
un gasto (caudal) base.
* No supone modelo de flujo.
CONVOLUCION
* Supone un modelo de flujo.
* Superpone los efectos de cambios
del gasto (caudal).

APPVI32

DIAGNOSTICO DE FLUJO

HERRAMIENTAS

Pruebas de un solo pozo


Funcin de primera derivada.
* Funcin de segunda
derivada.
Derivadas especficas

Pruebas multipozos
* Ajuste de curva tipo.
APPVI33

CONFORMACION DEL MODELO

Tipos de Flujo

Patrones de Comportamiento

Modelo Integral de Flujo


APPVI34

GRAFICAS ESPECIALIZADAS

p = bflujo + mflujo f(t)


Flujo

p
mflujo
1
bflujo

f(t)

1/2
t
Lineal
t1/4
Bilineal
Log t
Radial
Esfrico
t-1/2
Almacen.
t
Pseudoest. t
Pres. Cte.
t-1

f(t)
APPVI35

AJUSTE DE CURVA TIPO

2s

p'

CD e
t

Log tD p'D

Log p

Pozo con Almacenamiento y Dao


(Flujo Radial)

Log tD /C D

t (hrs)
APPVI36

ESTIMACION DE PARAMETROS

Grficas
Especializadas

Ajuste de
Curva Tipo

Valor de Parmetros

APPVI37

VALIDACION DEL MODELO

Principio
Comparar respuesta de presin calculada
con la respuesta de presin medida.
Parmetros
Estimados

Condiciones de
la prueba.

Modelo seleccionado

Respuesta de presin calculada


APPVI38

VALIDACION DEL MODELO

Log
pm
pc

pm
pc

t pm'
t pc'
t

Log t
APPVI39

INFORME DE RESULTADOS

Contenido:
* Resumen de prueba y resultados
* Bitcora de la prueba
* Diagrama del pozo con herramienta
* Datos medidos
* Secuencia de anlisis
* Grficas de interpretacin
* Validacin de modelo(s)
* Conclusiones y recomendaciones.
APPVI40

VII. PRUEBAS DE DECREMENTO


(Drawdown Test)
Objetivo:
Presentar, discutir y analizar los
mtodos de interpretacin de datos
de los diversos tipos de pruebas de
decremento de presin.

APPVII00

PRUEBA DE DECREMENTO
(Drawdown Test)
Medicin contnua de la presin de fondo
de un pozo durante un perodo de flujo

pwf

APPVII01

PRUEBA DE DECREMENTO
Ventajas:
* Estimacin de la capacidad de
flujo del pozo.
* Anlisis simple de datos.
Desventajas
* Variacin del gasto (caudal)
durante la prueba.
Tipos:
* Decremento sencillo
* Prueba multiflujo
* Lmite de yacimiento
APPVII02

PRUEBA DE DECREMENTO

Anlisis
* Diagnstico de prueba
* Suavizacin de datos
* Normalizacin
* Diagnstico de flujo
* Conformacin del modelo
* Grficas especializadas
* Ajuste de curva tipo
* Estimacin de parmetros
* Validacin de modelo(s)
* Informe de resultados.
APPVII03

PRUEBA DE DECREMENTO

Comentarios
* Las tcnicas de anlisis de datos
consideran un gasto (caudal) constante
* Es necesario transformar la informacin
tomada durante la prueba para estimar
la respuesta de presin correspondiente
a un gasto base constante.

APPVII04

PRUEBA DE DECREMENTO
( Drawdown Test )
pi

pwf
q(t)

pwf ( t )

t
pwf vs t

Normalizacin

p1 vs t

q vs t
APPVII05

PRUEBA DE DECREMENTO
Respuesta de Presin
pwf ( t ) =

t
q'( )

p1( t- ) d

0
pwf ( t ) =

t
q( )

p1'( t- ) d

0
medido

medido

?
APPVII06

NORMALIZACION

Estimacin de la respuesta de presin


correspondiente a un gasto (caudal)
constante.

Mtodos:
- "Normalizacin"
- Convolucin
- Deconvolucin
APPVII07

"Normalizacin"
Datos:

pwf ( t )

p1 ( t )

q(t)

pwf ( t )
q(t)

( Gladfelter )

Mtodo aproximado
Produce resultados aceptables
para variaciones suaves de q.
pqbase ( t )

pwf ( t )
q(t)

x qbase
APPVII08

NORMALIZACION
Convolucin
Supone la forma de la funcin
influencia . ( seleccionar modelo )
t
q'( )

pwf ( t ) =

p1( t-

)d

0
Para N Periodos de flujo
n

pwf ( t ) =

( qi . qi.1 ) x

p1 ( t - t i )

i=1
APPVII09

CONVOLUCION
Flujo Lineal
1 /2
16.25 B
t
A(k
c t ) 1 /2

p1 ( t ) =

pwf ( t ) =

pwt ( t )
qn

16.25 B
A(k
c t ) 1 /2
= 16.25 B
Ak
c t ) 1 /2
+

p1,dao

( q i - qi -1) ( t - ti )

/2

i=1
n

i=1

/2

(qi - qi.1) ( t - ti )
qn

APPVII10

Grfica de Convolucin
Flujo Lineal
Pwf

qn
ml f 1 =

n
i=1

( qi - qi-1 )
qn

16.25 B
A(k

ct )1/2

( t - ti )1/2
APPVII11

CONVOLUCION
Flujo Radial
p1 (t) = m1 ( log t + log
n

k
ctrw

k
- 3.2275 + 0.87 S
2
ctrw

= m1

( qi-qi-1 )

i=1

+ m1 { log

qn

log ( t-ti )

k
- 3.2275 + 0.87 S
2
ctrw

+ qn m1 { log

APPVII12

i=1

qn

-3.2275 + 0.87 S)

( qi-qi-1 ) log ( t-ti ) x m1

pwf (t) =

pwf (t)

Grfica de Convolucin
Flujo Radial
Pwf

162.6 B
m1 =
kh

qn

1
b
0

n
i=1

qi - qi-1
qn

Log ( t - ti )
APPVII13

ESTIMACION DE PARAMETROS
162.6 B
kh =
ml

k
b
3.2275
+
S = 1.151 { m - Log
ct r w2
l
Limitaciones
. Vlido para el modelo de
flujo radial.
APPVII13A

CONVOLUCION
Comentarios:
* Aplicable a datos dominados por el
flujo seleccionado para generar la
grfica de convolucin.
* El nmero de trminos, n, en la
serie de convolucin es tal que
tn+1 t > tn.
* Los parmetros del modelo se estiman de la pendiente de la lnea
recta y de la ordenada al origen.
APPVII14

DECONVOLUCION
Definicin:
Estimacin de la respuesta de presin
correspondiente a gasto constante
(unitario o base).
t

Respuesta
de presin
medida

pwf vs t
q vs t

pwf (t) =

q( )

0
Deconvolucin

No supone modelo de flujo.

p1 ( t- ) d

p1 vs t

APPVII15

DECONVOLUCION
p1 ( t )

Estimacin de

(Kuchuc & Ayestarn)


ti

pwf ( tn ) =

q' ( tn- )
i=1

ti-1

p1( ti-1/2 )

=
i=1

tn-ti-1

ti
q' ( tn- ) d
ti-1

p1 ( ) d
ti
q' ( tn- ) d
ti-1

q' ( ') d '


tn-ti

APPVII16

DECONVOLUCION
n

p1 ( ti-1/2 ) x [ q ( tn-ti-1 ) -q (tn - ti )]

pwf ( tn ) =
i=1

n=1

pwf ( t1 ) =
p1 ( t1/2 ) =

pwf ( t1 )
q ( t1 )

n>1

n-1

pwf ( tn ) p1 ( tn-1/2 ) =

p1 ( t1/2 ) q ( t1 )

p1 ( ti-1/2 ) [ q ( tn-ti-1 ) -q (tn - ti )]


i=1

q ( t1 )
APPVII17

FUNCION INFLUENCIA

Caractersticas:
p1

p1 0

p1'

p'1 > 0
p"1 < 0

t
APPVII18

DECONVOLUCION

Respuesta de presin para qbase


pqbase(t)

= qbase

p1(t)

Aplicar metodologa de anlisis a


p1(t)

pqbase(t)

APPVII19

Diagnstico de Flujo
pqbase vs t
Diagnstico
Tipo(s) de Flujo
Conformar el modelo
de flujo
Grficas especficas

APPVII20

DISEO
+ OBJETIVO
+ DURACION
. t ews
. t bh
}
. t eh

tiempo de efectos de
heterogeneidades

+ CONDICIONES
. Caudal (gasto)
. Herramienta
- Resolucin
- Limitaciones

APVII21A

PRUEBA MULTIFLUJO
Objetivos
+ Anlisis Nodal
+ Efectos de Alta velocidad
( turbulencia )
q4
p
i
q3
pwf
q2
q

q1

t1

t2

t3

t4

t5

t
APPVII22

PRUEBA DE LIMITE DE YACIMIENTO

Objetivos:
+ Volumen Poroso de drene
+ Forma del rea de drene
+ Posicin del pozo dentro de rea de
drene.
Principio:
Alcanzar durante la prueba los efectos
de las fronteras del rea de drene.
(Flujo Pseudo-estacionario).
APPVII23

PRUEBA DE LIMITE DE YACIMIENTO


Datos: pwf vs t

Diagnstico de Flujo
Log
pwf
t

Almacenamiento

tpss

p'wf

1
1

tews
Radial

1
teia

PseudoEstacionario

Log t
APPVII24

PRUEBA DE LIMITE DE YACIMIENTO

Grfica de Flujo Radial


Almacenamiento

tews

( pwf)1hr

teia

m
1
Radial
t=1hr
kh = 162.6 q B

tpss

/m

S = 1.151 [ (pi-pwf1hr)/m - log(kh/

PseudoEstacionario

Log t
ctrw2) + 3.2275

pwf

APPVII25

PRUEBA DE LIMITE DE YACIMIENTO


Grfica de Flujo Pseudoestacionario

pi
pwf
b

Almacenamiento

tews
teia

tpss

Radial

PseudoEstacionario

m*
1

t
Vp = 0.23395 q B / ct m*

CA = 5.456 (m/m*) e

- 2.303 (b*-

pw1hr)/m
APPVII26

PRUEBA DE LIMITE DE YACIMIENTO

Vp
teia

teiaDA

tpss

tpssDA

CA

- Forma del rea de drene


- Posicin del pozo dentro
del rea de drene.

APPVII27

PRUEBA DE LIMITE DE YACIMIENTO

Diseo
-Condiciones de prueba
*Usar mximo gasto permisible
*Mantener gasto constante
*Medir gasto contnuamente de
preferencia en el fondo del pozo
*Medir presin en la cabeza del pozo
*Usar herramienta de alta resolucin
-Duracin
*Duracin mnima de 2 tpss
APPVII28

VIII PRUEBAS DE INCREMENTO

Objetivo:
Presentar, discutir y analizar los
diversos mtodos de interpretacin
Para pruebas de incremento de Presin

APVIII01

PRUEBA DE INCREMENTO
PRESSURE BUILDUP TEST
Definicin:
Medicin contnua de la presin de cierre
de un pozo despus de un perodo de Flujo
pi
pw

pwf

pws (

t)

tp

t (tiempo de cierre)
t

APVIII02

PRUEBA DE INCREMENTO
Objetivo:
Estimar parmetros del yacimiento
Estimar el factor de dao del pozo
Determinar la presin media del rea de
drene
Ventajas:
Mediciones suaves de presin
Gasto (caudal) constante (q=o)
Desventajas
Produccin diferida de hidrocarburos
Anlisis de datos complejo.
APVIII03

PRUEBA DE INCREMENTO
Ecuaciones de comportamiento
pws = q

p1 (tp +

t)-q p1 (

t)

pi

pws

pw
Por cada medicin de
Presin existen dos
incgnitas.
Medicin
pws ( t)

Incgnitas
p1 (tp +
t)
p1 ( t)

pwf

pws (

t)

q
p
t
tp

t
APVIII04

PRUEBA DE INCREMENTO
Respuesta de Presin
pws ( t)
pw

pws
q(t)

t
tp

pws (

t)=

tp
q, (
0

t
)

p 1 ( tp +

t-

)d
APVIII06

PRUEBA DE INCREMENTO
Interpretacin

pws vs t
q vs t

Normalizacn

p1 vs t

Mtodos de pruebas
de Decremento
APVIII07

PRUEBA DE INCREMENTO
pi
pw
pws
q(t)
q

Q( tp )
tp

tp
2 tp

2 tp

3 tp

pws depende del ltimo gasto antes del cierre


pws depende de q ( t ) y de tp
APVIII08
pws depende de Q ( tp )

PRUEBAS DE INCREMENTO
Normalizacin
* Convolucin
- Modelo preseleccionado
- Historia de flujo
* Deconvolucin
- No requiere un modelo preseleccionado
- Historia de flujo
* Impulso
- No requiere un modelo preseleccionado
- Historia de flujo
APVIII09

PRUEBAS DE INCREMENTO

Convolucin:
n

pws ( t) =

(qi - qi-1) p1 (tp + t - ti)


i=1

Si se considera un gasto base qbase ( qltimo)


n

pws ( t) =
i=1

(qi - qi-1)
qbase

pqbase (t + t - t )
p
i
APVIII10

PRUEBAS DE INCREMENTO
Convolucin
Flujo Radial

pws = pi n+1

i=1

1.151

qbase

kh

qi - qi-1
qbase

pws ( t)

pi

pws ( t)

pw
x
q3
q2

Log (tp +

t - ti)
q

qn

q1
t
t1 t2 t3

tn tp
tn+1

t
APVIII11

PRUEBAS DE INCREMENTO
Grfica de Convolucin (Flujo Radial)
pws

pi
1
mqbase=

n+1

i=1

qi - qi-1
qbase

1.151

Log (tp +

qbase

kh

t - ti)

APVIII12

PRUEBAS DE INCREMENTO
Convolucin para Flujo Radial
Suposicin
pws ( t) = pi -

pi

q = constante
1.151

qB

Log

kh

tp +

t
t

Grfica de Horner

pws

m=

Log

1.151

qB
kh

tp +

t
t

APVIII13

PRUEBAS DE INCREMENTO

Estimacin del Factor de Dao


S = 1.151

pws (

t = 1 ) - pwf
m
- Log

+ Log
k
ctrw2

tp + 1
tp
- 0.35137

Si tp >> 1
S = 1.151

pws (

t = 1 ) - pwf
m

- Log

k
ctrw2

- 0.35137
APVIII14

PRUEBA DE INCREMENTO
Convolucin

Modelo de Flujo j
p = A + m f (t)
1
j j

p ( t) = p - m
ws
j
i

n+1
i=1

( q i - q i -1 ) f j ( t p +

t - ti )

pi
pws

n+1
i=1

mj

( q i - q i -1 ) f j ( t p +

t - ti )

APVIII15

Prueba de Incremento
Tiempo de Superposicin
. Flujo Radial
tsup =

n+ l
i= l

q - qi i
1 ) Log (t +
( q
p
base

t - ti )

. Flujo j (Lineal, bilineal, esfrico, etc.)


tsupj =

n+ l
i= l

q i - qi - 1
) f j (t p +
( q
base

t - ti )

APVIII16

Prueba de Incremento
Tiempo efectivo de Agarwal
(Tiempo Equivalente)
Definicin:

t ef

tp t
t ef =
tp + t

Aplicacin:
Modificar la escala del tiempo para ajustar
datos de una prueba de incremento a
curvas tipo de pruebas de decremento
APVIII17

PRUEBA DE INCREMENTO
Tiempo equivalente
Log
pws

tef

dpws
dtef
t

pws
'
Log

tef
APVIII18

PRUEBA DE INCREMENTO
Efecto del tiempo equivalente
Lineal

Log

1/2
1

dpws
dtef
t

Bilineal

1/4
1

pws
'
-1/2

Esfrico
1

Log

tef

APVIII19

PRUEBAS DE INCREMENTO
Tiempo de Superposicin
tsup =
i=1

(qi - qi-1)

Log (

qn

tp +

t - ti-1

tp +

t - ti

Derivada con respecto a tsup


t pequeo

dpws
dtsup t

t grande

dpws
dtsup

p'

t|
2

p''|

Funcin de
1a. derivada
Funcin de
2a. derivada
APVIII20

PRUEBA DE INCREMENTO
Derivada con respecto a tsup
Log
t

p'

Lineal

Bilineal
1/2

t2

p''

1/4
1

dpws
dtsup

Radial
1
-1/2

Esfrico

Log

APVIII21

PRUEBA DE INCREMENTO
Deconvolucin
Estimar la funcin influencia a partir de
datos de presin y de gasto
Respuesta de presin :
tp
pws (

q' (

t)=

p 1 ( tp +

t-

) d

o
APVIII23

PRUEBAS DE INCREMENTO
Deconvolucin
- Analizar datos del periodo de decremento
pws ( t) , q (t)

p1 (t)

o t tp

- Estimar la funcin influencia a partir de


datos de incremento
n

p1 (tp+ t-ti-1/2) =

pws ( t) -

i=2

p1 (tp+

t-ti-1/2) x

q(t1)
( q(ti) - q (ti-1) )
APVIII24

PRUEBA DE INCREMENTO
Respuesta de presin q -constante
pi

pw

pws

q p1 ( tp )

pwf

q p1 (tp + t)

q p 1 ( tp )

t
0

tp

APVIII25

PRUEBAS DE INCREMENTO
Si tp >> t

p1 (tp+ t)

p1 (tp)

pws ( t) = q p1 (tp) - q p1 ( t)
pi - pws ( t)
p1 ( t) =
pq ( t) =

pi - pwf (tp)
pws ( t) - pwf (tp)
q
pws ( t) - pwf (tp)
pws

APVIII26

PRUEBA DE INCREMENTO
Si tp >>

y q - constante.

- Usar pws para el anlisis de diagnostico


de flujo y ajuste de curva tipo
- Este anlisis es aplicable a cualquier
tpo de flujo
- A medida que t crece el error en el
anlisis aumenta .

APVIII27

PRUEBAS DE INCREMENTO
Si tp >> t

, q - constante
1.151 q

pws ( t) = pwf (tp) +

kh

+ log

ctrw

Flujo Radial

0.80907
2.303

log
+ (

t
p)dao

Grfica de MDH
pws
1

1.151 q
m=
kh

log

APVIII28

PRUEBA DE INCREMENTO
Deconvolucin (Desuperposicin)
Si tp << t

y q - constante

La funcin influencia puede estimarse de:


pi - pwf (tp)
p1 (tp) =
q
pi - pws (tp) + q p1 (tp)
p1 (2tp) =
q

p1 (ntp) =

pi - pws ((n-1) tp) + q


q

p1 ((n-1) tp)
APVIII29

PRUEBA DE INCREMENTO

DECONVOLUCION
Comentarios :
El proceso de deconvolucin nicamente se
puede aplicar si se tienen datos de gasto
y presin del perodo de flujo antes de
cierre.
El proceso de deconvolucin se simplifica
s el gasto antes del cierre es constante.
La funcin influencia puede estimarse para
tiempos de cierre pequeo ( t 0.1tp )
La funcn influencia puede calcularse para
tiempos mltiples de tp

APVIII30

PRUEBA DE INCREMENTO
Mtodo de Impulso
Soliman
Antecedentes
Lineal
1/2
1

Log

-1/2

Bilineal

1/4
-3/4

Radial
-1
1

Log tp+

APVIII31

PRUEBA DE INCREMENTO
Mtodo de Impulso
Antecedentes

Ayoub et al.
pD

Log pi-pwf
(tp+t)(pi-pws)

tDpD'

Log tp+

APVIII32

PRUEBA DE INCREMENTO
Mtodo de Impulso
Respuesta de Presin de un Impulso
( pws(t))imp= pi-(pws)imp

= 24 Q

p1' (t)

pi
( pws)imp

pws
(pws)imp

Q (Volumen del Impulso)

APVIII33

PRUEBA DE INCREMENTO
Mtodo de Impulso Cinco-Ley et al.
pws
(pws)imp
qimp
Q

q
tp

2tp

3tp

Lim pws( t) = (pws)imp


t
8

t
t

APVIII34

Esfrico
Radial
Bilineal

10

Lineal
100

( pws-( pws)imp)/( pws)imp x 100

PRUEBA DE INCREMENTO
Mtodo de Impulso
Diferencia entre la respuesta
de incremento y de impulso

t/tp

APVIII35

PRUEBA DE INCREMENTO
Mtodo de Impulso
Derivada de la Funcin Influencia
p1'

=
tp/2 + t

(pi - pws(

t))

24 Q

Funcin de Primera Derivada


t

p1'
tp/2 + t

(pi - pws( t))


(t
/2
+
t)
=
p
24 Q
APVIII36

PRUEBA DE INCREMENTO
Mtodo de Impulso
Ventajas
* Vlido para cualquier tipo de flujo
* La derivada se estima sin un proceso
de diferenciacin.

Desventajas
* Requiere la presin inicial
* Vlido para

t 2 tp
APVIII37

PRUEBA DE INCREMENTO
Mtodo de Impulso
pi
k 1

pw
q
q q1 2

qN

3N

t
tp

2N

2tp

t
t

t1 t2 t3 tN tN+1

tp

2tp

3tp
APVIII38

PRUEBA DE INCREMENTO
Mtodo de Impulso
k2N

Tiempos largos
p1( tk+tp/2)

t
Tiempos cortos
p1( tk+tp-tj-1/2)
t
( pi - pws(

( pi - pws(
24 Q
k<2N
=
N-1

tk ))

tk ))

qj

j=1

p1( tk+tp-tj-1/2)
t

qN
APVIII39

PRUEBA DE INCREMENTO
Mtodo de Impulso
Derivada de la Funcin Influencia
p1'

=
tp/2 + t

(pi - pws(

t))

24 Q

pi es desconocida
p1''

=
tp/2 + t

- pws(

t)

24 Q

Segunda derivada
APVIII40

PRUEBA DE INCREMENTO
Mtodo de Impulso
pi desconocida
Tiempos largos
2
p1( tk+tp/2)
t

Tiempos cortos
2

p1( tk+tp-tj-1/2)
t
- pws
'(

k2N
- pws' (

24 Q
k<2N
=

N-1

tk ) /

tk )

qj

t j=1

qN

p1( tk+tp-tj-1/2)
t2
APVIII41

PRUEBA DE INCREMENTO
Mtodo de Impulso
pi es desconocida
t2

Funcin de Segunda derivada

(t

)2
/2
+
t
p

p1''

tp/2 + t )2 pws' (

p1''
t)

24 Q

APVIII42

Curva tipo de segunda derivada para


Flujo radial con almacenamiento y dao

tD2 |p''D|

1.0E+02

1.0E+01

1.0E+00

1.0E-01
1.0E+00

1.0E+01

1.0E+02

1.0E+03

1.0E+04

tD/CD
APVIII43

PRUEBA DE INCREMENTO
Mtodo de Impulso
Estimacin de la presin inicial
pws vs

t2| p1''|

Log
t2| p1''|

n
1
Log t

del ltimo perodo de flujo detectado


pi = pws( t) - (

t+tp/2) pws'(

t) / (n-1)
APVIII44

Ejemplo

Prueba de decaimiento
Fall Off
qi = -2124 STB/D

pw

tp = 1.13 horas

pi

pi = 4203 lb/plg2
0

tp

t(h) pws(psi) t(h) pws(psi) t(h) pws(psi)


.25
.50
.75
1.0
1.5

4287.8
4262.8
4249.6
4241.7
4232.1

2.0
3.0
4.0
6.0
8.0

4224.9
4218.0
4214.6
4211.3
4208.3

12 4206.8
16 4206.2
20 4205.6
24 4206.8
APVIII45

Ejemplo
Anlisis

Prueba de decaimiento
Fall Off

Diagntico de flujo:
* Derivada con respecto al tiempo de
cierre
* Derivada con respecto al tiempo de
superposicin
* Mtodo de impulso
APVIII46

Prueba de decaimiento
Fall Off
Ejemplo
Derivada con respecto
Diagnstico de flujo:
al tiempo de cierre
1E+02
t tp'
.375
.625
.875
1.25
1.75
2.50
3.50
5.00
7.00
10.0
14.0
18.0

37.5
33.0
27.6
24.0
25.2
17.25
11.90
8.25
10.50
3.75
2.1
2.7

tp'

1E+01

-1

1E+00
1E-01

1E+00

1E+01

1E+02
APVIII47

Prueba de decaimiento
Ejemplo
Derivada con respecto
al tiempo de
Diagntico de flujo:
superposicin
1E+03
t dt/dt
dp/dtsup
sup
.375
.625
.875
1.25
1.75
2.50
3.50
5.00
7.00
10.0
14.0
18.0

109
116
111
114
144
123
113
100
166
83
66
86

dp/dtsup
1E+02

Flujo radial

1E+01
1E-01

1E+00

1E+01

1E+02
APVIII48

Prueba de decaimiento
Ejemplo

Mtodo de Impulso

Usar de funcin de primera derivada puesto que


la presin inicial es conocida
(t p1')tp/2+t = ( tp/2 +t ) (pi - pws (t)) / 24 Q
pi = 4203 Lb/plg2

t p = 1.13 horas

Q = q tp = 2124 x 1.13 / 24 = 100 Barriles


24 Q = 24 x 100 = 2400

APVIII49

Ejemplo

Prueba de decaimiento

Diagnstico de flujo:

tp/2+t
0.815
1.065
1.315
1.565
2.065
2.565
3.565
4.565
6.565
8.565
12.56
16.56
20.56

tp1'

2.87x10-2
2.65x10-2
2.55x10-2
2.50x10-2
2.50x10-2
2.33x10-2
2.21x10-2
2.20x10-2
2.27x10-2
1.89x10-2
1.98x10-2
2.20x10-2
2.22x10-2

1E-01

Mtodo de impulso
Funcin de derivada de
la funcin influencia

t p1'
1E-02

1E-03
1E-01

Flujo radial

1E+00

1E+01

tp/2+t

1E+02
APVIII50

Prueba de decaimiento
Flujo radial
Grfica de Horner

Ejemplo
4300

pws

4250

4200

m = 114 Lb/plg 2/ciclo

pi = 4203 Lb/plg2

4150
1

(tp+t) / t

10
APVIII51

Prueba de decaimiento
Ejemplo

Mtodo de Impulso

Si no se conoce la presin inicial se debe


estimar la funcin de segunda derivada
t2 Ip'' I = (tp/2 +t)2 pws '(t) / 24 Q
en donde pws' es la derivada de presin de
cierre con respecto al tiempo de cierre

APVIII52

Prueba de decaimiento
Mtodo de Impulso
Ejemplo
Presin inicial desconocida
Grfica de la funcin de
2
tp/2+t t I p'' I
segunda derivada
0.94
3.68x10 -2
1.16
1.44
1.81
2.31
3.05
4.06
5.56
7.56
10.56
14.56
18.56

0.1

-2

2.98x10
2.73x10 -2
2.91x10 -2
2.77x10 -2
2.70x10 -2
2.34x10 -2
2.13x10 -2
3.57x10 -2
1.74x10 -2
1.32x10 -2
2.15x10 -2

t Ip''I
0.01

0.001
0.1

n=0

Flujo radial

tp/2+t

10

100
APVIII53

Ejemplo

Prueba de decaimiento
Mtodo de Impulso

Estimacin de la presin inicial

pi = pws(t) - (tp/2+t) pws'(t) / ( n-1 )


del diagnstico de flujo :

n=0

pi = pws(t) + (tp/2+t) pws'(t)

APVIII54

Prueba de decaimiento
Mtodo de Impulso

Ejemplo
Estimacin de la presin inicial
t
.375
.600
.875
1.25
1.75
2.5
3.5
5.
7.
10.
14.

pws'

(tp/2+t) pws'

-100
-52.8
-31.6
-21.2
-12.4
-6.9
-3.4
-1.65
-1.5
-0.375
-0.15

-94.
-61.5
-45.5
-38.5
-28.7
-19.9
-13.82
-9.18
-11.34
-3.96
-2.18

pi
4181
4194
4199
4198.4
4199.3
4201.5
4202.4
4203.8
4203.5
4203.3
4203.12
APVIII55

Prueba de Incremento
Diseo
* Objetivo
* Duracin
. t ews
. t bh
. t eh
. t 0.1 tp o t > 3 tp
* Herramienta
. Alta resolucin
. Medidor de flujo
. Cierre en el fondo
APVIII56

Presin Esttica p
Definicin
Presin promedio en el rea de drene
de un pozo en el momento del cierre

APVIII57

Presin Esttica p
Estimacin
* Mtodo de Matthews-Brons-Hazebroek
(MBH)
* Mtodo de Miiler-Dyes-Hutchinson
(MDH)
* Mtodo de Dietz
* Mtodo de Ramey-Cobb
pws(t)

p
APVIII58

Presin Esttica p
Mtodo de MBH
pws vs t
Horner

p*

1
pws

p*
Log (tp+t)/t

MBH
p

p = p* - m pDMBH / 2.303
APVIII59

Mtodo de MBH
pDMBH = f (tpDA, forma del Area,
posicin del pozo)
tpDA = k tp / ct A
MBH
pDMBH

tpDA
APVIII60

Grfica de MBH
6
5
4

pDMBH 3
2
1
0
-1
1E-02

1E-01

tpDA

1E+00

1E+01

APVIII61

Grfica de MBH
6
5
4

pDMBH

3
2
1
0
-1
1E-02

1E-01

tpDA

1E+00

1E+01

APVIII62

Grfica de MBH
6
5
4

pDMBH 3
2
1
0
-1
1E-02

1E-01

tpDA

1E+00

1E+01

APVIII63

Grfica de MBH
5
4
4

1
4

pDMBH 2

4
1

0
-1

4
1

-2
1E-02

1E-01

tpDA

1E+00

1E+01

APVIII64

Presin esttica
Mtodo de MDH

pws

m
1

Log t
p = pws(t) + m pDMDH(tDA) / 1.1513
tDA = k t /ct A
Supone flujo pseudo estacionario
antes del cierre

APVIII65

Grfica de MDH
5
4

pDMDH

3
2
1
0
1E-05

1E-04

1E-03

1E-02

1E-01

1E+00

tDA
APVIII66

Presin esttica
Mtodo de Dietz

MDH

p
pws

Log t

(t)p

(t)p = tp / CA tpDA = ct A / k CA
Supone flujo pseudo estacionario
antes del cierre

APVIII67

Presin esttica
Mtodo de Ramey-Cobb
Horner
p
pws

Log tp+t /t

( tp+t /t ) p

(tp+t/t)p = CA tpDA = k tp CA / ct A
Supone flujo pseudo estacionario
antes del cierre

APVIII68

Presin Dinmica de Bloque


Malla de Simulacin

* Presin esttica
* Presin de bloque o celda
APVIII69

Presin Dinmica de Bloque pd

pd
Area
de
drene

pwf
x
Bloque
APVIII70

Presin Dinmica de Bloque pd


Estimacin a partir de la grfica MDH
pd
pws

(t)pd
Log t
( t)pd = 200 ct ( x)2 / k
APVIII71

Presin media del yacimiento

p1

p3

p2
4

p4

py = pi Vpi / Vpt
i=1

APVIII72

IX.
IX. PRUEBA
PRUEBA DE
DE INTERFERENCIA
INTERFERENCIA

Objetivo:
Objetivo:
Presentar,
Presentar, discutir
discutir yy analizar
analizar los
los
mtodos
mtodos de
de interpretacin
interpretacin para
para
pruebas
pruebas de
de interferencia
interferencia yy de
de un
un
solo
solo pulso.
pulso.

APPIX01

IX.
IX. PRUEBA
PRUEBA DE
DE INTERFERENCIA
INTERFERENCIA

Definicin:
Medicin contnua en un pozo de observacin
de la respuesta de presin causada por un
cambio del gasto (caudal) en otro pozo (activo).

Pozo
Activo

Pozo de
Observacin

p(t)
APPIX02

PRUEBA
PRUEBA DE
DE INTERFERENCIA
INTERFERENCIA

q
Pozo
Activo

p(t)
Pozo de
Observacin

t
APPIX03

PRUEBA DE INTERFERENCIA

Ventajas
Caracterizacin de la zona localizada entre
el pozo activo y el pozo de observacin

Desventajas
Produccin diferida en varios pozos
Respuesta de presin pequea
APPIX04

Interpretacin
Datos

Filtrado

Anlisis
APPIX05

PRUEBA DE INTERFERENCIA

Interpretacin
Mtodo :
Ajuste de curva Tipo
Modelos de Flujo:
* Flujo Lineal
* Flujo Radial
* Flujo Esfrico

APPIX06

PRUEBA DE INTERFERENCIA
Curva Tipo

1E+02

Lineal
1E+01

F1(pD)

Radial

1E+00

Esfrico
1E-01

1E-02
1E-01

1E+00

1E+01

F2(tD)

1E+02

1E+03
APPIX07

Definicin de Variables Adimensionales

Flujo

F 1(pD)

F 2(tD)

Lineal

kbhp / lqBx

kt / ctx2

Radial

khp / qB

kt / ctr2

Esfrico

krp / sphqB

kt / ctr2

APPIX08

PRUEBA DE INTERFERENCIA
Curva Tipo

1E+02

Lineal
1E+01

F1(pD)

Radial

1E+00

Esfrico
1E-01

1E-02
1E-01

Punto de Ajuste

t (h)
1E+00

1E+01

F2(tD)

1E+02

1E+03
APPIX09

Anlisis de Prueba de Interferencia


Resultados del ajuste:
(p)M
( t )M

(F 1)M
(F 2)M

Geometra
de flujo

Estimacin de parmetros
Lineal:
Radial:

kbh = lqB(F1)M / (p)M


ctbh =kbh(t)M / x2(F2)M
kh = qB(F1)M / (p)M

cth =kh(t)M /r2(F2)M


Esfrico: k = sphqB(F1)M / (p)M
ct = k(t)M / r2(F2)M
APPIX10

Pruebas de Interferencia
Zona en estudio

Activo

Observacin
r

La zona que afecta una prueba


de interferenciaes una elipse
(Vela)

APPIX11

Pruebas de Interferencia
Factores que complican el anlisis de una
prueba de interferencia:
* Rudo en la informacin
. Respuesta de presin pequea
. Efectos de pozo (temperatura,
segregacin de fluidos, almacenamiento)
* Tendencia de presin
* Corta duracin de la prueba
APPIX12

Pruebas de Interferencia
q
q

Activo

to
p

t
pinterf

Observacin
t
APPIX13

Ejemplo. Analizar prueba de Interferencia


Datos:

t (hrs)
20
30
40
50
60
70
80

q = 1200 BPD
= 1.2 cp
h = 150 pies
Dp (lb/plg2)
1.2
3.6
6.5
9.5
11.5
14.
17.

Bo = 1.3 bls/stb
= 0.08
r = 900 pies
t(hrs)
90
100
110
120
140
160
180

Dp(lb/plg2)
19.5
21.5
23.
24.5
28.
32.
36.
APPIX14

PRUEBA DE INTERFERENCIA
Mtodo de El Khatib

Flujo Radial

pD(rD, tD) = E1( rD2/ 4tD ) / 2


tDpD' = exp(-rD2/ 4tD) / 2tD
t p' =qB exp(-ctr2/ 4kt) /2 kh
Log t p' = Log (qB/ 2kh)
- ( ctr2/ 9.212 k) t-1
APPIX15

PRUEBA DE INTERFERENCIA
Mtodo de El Khatib

Flujo radial

b =qB/ 2kh
1
m = -ctr2/ 9.212 k

Log tp'

1/t
APPIX16

PRUEBA DE INTERFERENCIA
Mtodo de El Khatib

Flujo radial

Estimacin de parmetros
k h = q B / 2 b
ct h = - 9.212 kh m / r2

APPIX17

PRUEBA DE INTERFERENCIA
Funcin Integral I(p)
Definicin:
t

p( ) d

I(p) = ( 1 / t )
0

Anlisis: Mtodo de ajuste de curva tipo


Ventajas: Datos suavizados

APPIX18

PRUEBA DE INTERFERENCIA
Curva Tipo para la Funcin Integral I(p)
tD/rD2

I(pD) = (1/tD/rD2)

pD( ) d
0

pD es la respuesta de presin adimensional


para el modelo de flujo presente durante la
prueba de interferencia.

APPIX19

PRUEBA DE INTERFERENCIA
Curva Tipo para la Funcin Integral I(p)
1E+01

Lineal

Radial

1E+00

Esfrico

1E-01

1E-02

1E-03
1E-02

1E-01

1E+00

1E+01

F2(TD/RD2)

1E+02

1E+03
APPIX20

PRUEBA DE INTERFERENCIA
Diseo de una prueba
* Duracin
t D / r D2 > 1
t > ctr2 /k
* Condiciones
- Mximo gasto posible (cte)
- Herramienta de alta resolucin
- Medir tendencia de presin
APPIX20A

PRUEBA DE UN PULSO
Activo

t0

tp

t
Observacin

p
ppulso

APPIX21

PRUEBA DE UN PULSO
Activo
q

t0

tp

t
Observacin

p
ppulso

APPIX22

PRUEBA DE UN PULSO
Objetivo:
Facilitar la deteccin de la seal de
presin en el pozo de observacin.
Interpretacin:
Ajuste de curva tipo
Log
pD

tpD/rD2

p
t

Log tD/rD2

APPIX23

PRUEBA DE UN PULSO
Flujo Lineal
1E+01

Curva Tipo
tpD/xD2
4

1E+00

pDL/tDL

1E-01

.1
.04
.02

1E-02
1E-01

1E+00

.2

.4

1E+01

tDL/xD2

1E+02

1E+03
APPIX24

PRUEBA DE UN PULSO
Flujo Radial
1E+01

Curva Tipo

tpD/rD2
1E+00

pD
1E-01

1E-02
1E-01

.4
.2
.1
.04
.02
1E+00

10
4
1 2

1E+01

tD/rD2

100
40
20

1E+02

1E+03
APPIX25

PRUEBA DE UN PULSO
Flujo Esfrico
1E+00

1E-01

pDsph rD
1E-02

1E-03
1E-01

.1
.04
.02

.2

.4

1E+00

Curva Tipo

4
1 2

10 20 40 100
tpD/rD2

1E+01

tD/rD2

1E+02

1E+03
APPIX26

PRUEBA DE FORMACION (DST)

PRUEBA DE FORMACION (DST)

3. Deteccin y Evaluacin de los


Lmites de un Yacimiento

Objetivo:
Presentar, discutir y analizar los
mtodos disponibles para
detectar y evaluar mediante
pruebas de presin los diversos
elementos que constituyen los
lmites de un yacimiento

CY3-001

Lmites de un yacimiento

Pozo

Aceite

Gas

Falla

Agua
CY3-002

Lmites de un yacimiento
Elementos:
Fallas
Acuamientos
Discontinuidades
Contacto agua-aceite
Contacto gas-aceite
Estratos semipermeables
Volumen poroso

CY3-003

Fallas Geolgicas
Clasificacin de acuerdo a su
comportamiento hidrodinmico:
Impermeables

Semipermeables

Conductivas

CY3-004

Fallas Geolgicas
Caracterizacin Hidrodinmica
Deteccin
Localizacin
. Posicin
. Orientacin
Caractersticas hidrulicas
. Conductividad
. Dao
. Longitud
CY3-005

Fallas Geolgicas
Fallas impermeables

df

CY3-006

Fallas Impermeables
Comportamiento de flujo
Simulacin
Pozo

Flujo radial

Imagen

Real
df
2 df

p (t) = ( p)real + ( p)imagen


Efecto de la falla

CY3-007

Fallas Impermeables
Comportamiento de flujo
Radial

Transicin

Semiradial

Zona de expansin
CY3-008

Fallas Impermeables
Evaluacin
Prueba de decremento
Prueba de incremento
Prueba de inyeccin
Prueba de abatimiento
Prueba de interferencia
CY3-009

Fallas Impermeables
Prueba de decremento
Comportamiento de presin
pw(t) = m(Log(t) + Log( /rw2) - 3.2275 + 0.87 S)
+ 2.303 m E1( df2/
donde

m = 162.6 q B

t)

/kh

E1 = Integral exponencial
CY3-010

Fallas Impermeables
Prueba de decremento
Comportamiento a tiempos cortos
t 0.4 df2/

Perodo de flujo radial

pw(t) = m(Log(t) + Log( /rw2) - 3.2275 + 0.87 S)

Funcin de derivada

p = m / 2.303
CY3-011

Fallas Impermeables
Prueba de decremento
Comportamiento a tiempos largos
t 20 df2/

Perodo de flujo semiradial

pw(t) = 2 m Log(t)
+ m (Log( /rw2) +Log(
- 3.2275 + 0.87 S)
Funcin de derivada

/4df2)

p = m / 1.151
CY3-012

Fallas Impermeables
Prueba de decremento
Grfica semilogartmica

2m

pw
1
m
1
Log t

CY3-013

Fallas Impermeables
Prueba de decremento
Grfica de diagnstico
Log
t pw'

2 ciclos

Log 2

Log t

CY3-014

Fallas Impermeables
Prueba de decremento
Estimacin de df

* Interseccin de rectas
semilogartmicas
* Ajuste de curva tipo
* Desuperposicin

CY3-015

Fallas Impermeables
Prueba de decremento
Interseccin de rectas
1

pw
2m
1
m

tint
Log t

df = 0.01217 (

tint)1/2
CY3-016

Fallas Impermeables
Prueba de decremento

Respuesta de presin
pWD = 1/2 ( Ln tD + 0.80907 ) + S
+ 1/2 E1( dfD2/tD)
Funcin de derivada
tD pWD' = 1/2 + 1/2 e
dfD = df / rw

- 1/(tD/dfD2)

CY3-017

Fallas Impermeables
Prueba de decremento
Curva Tipo

Log
tD pWD'

Log tD/dfD2

CY3-018

Fallas Impermeables
Prueba de decremento
Ajuste de Curva Tipo

100

t p'

Log
tD pWD'

Punto de
ajuste

10
.1

t (hrs)

Log tD/dfD2

10

CY3-019

Fallas Impermeables
Prueba de decremento
Ajuste de curva tipo
Estimacin de
parmetros

Ajuste
(tD pWD')M
(t

pW')M

kh=

(tD/dfD2)M
(t)M

df =

(tD pWD')M

141.2 q B
(t

pW')M

2.637x10-4

(t)M

(tD/dfD2)M
CY3-020

Fallas Impermeables
Prueba de decremento
Desuperposicin

pw
2m
1

pfalla

m
Log t

pfalla satisface la solucin de lnea fuente


CY3-021

Fallas Impermeables

Log pD

pfalla

Prueba de decremento
Desuperposicin

M
Punto de
ajuste
t
Log tD/rD2

CY3-022

Fallas Impermeables
Prueba de decremento
Desuperposicin
Estimacin de
parmetros

Ajuste
(pD)M
( pfalla)M

kh=

(tD/rD2)M
(t)M

df =

(pD)M

141.2 q B
(

pfalla)M

2.637x10-4

(t)M

4 (tD/rD2)M
CY3-023

Fallas Impermeables
Prueba de decremento
Posicin de la falla:
La falla es tangente a
un crculo de radio df

df
* No es posible determinar la
orientacin de una falla con
una sola prueba de presin
CY3-024

Fallas Impermeables
Prueba de decremento
La posicin de una falla se puede hallar
con un mnimo de tres pruebas en pozos
no alineados
2
Pozo 1

df 1

Pozo 2

df 2

Pozo 3

df 3

df 1

df 2
df 3
3

CY3-025

Fallas Impermeables
Prueba de incremento
Respuesta de presin

pws( t) = ( p(tp+ t))real


+ ( p(tp+ t))imagen
- ( p( t))real
- ( p( t))imagen
CY3-026

Fallas Impermeables
Prueba de incremento

pws vs

Normalizacin
p1 vs t

Aplicar la metodologa para


pruebas de decremento
CY3-027

Fallas Impermeables
Prueba de incremento

Grfica de Horner
pi

tp5 > tp4 > tp3 > tp2 >tp1

2m

pws

tp1

tp2
tp4

tp3
1

tp5
.001

.01

Log

t / (tp +

0.1

t)

1
CY3-028

Fallas Impermeables
Prueba de interferencia
2

Activo

Observacin

1
Observacin

1. La prueba proporciona resultados


cuantitativos
2. La prueba proporciona resultados
cualitativos

CY3-029

Fallas Impermeables
Prueba de interferencia

Activo

Real
Simulacin

Imagen

ri

Observacin
CY3-030

Fallas Impermeables
Prueba de interferencia
Respuesta de presin
Observacin
pDO(tD)

Real
=

Imagen

pDA(rD,tD) + pDI(rDi,tD)

pDO(tD) = 1/2 E1(1/4tD/rD2) + 1/2 E1(1/4tD/rDi2)


donde

rDi = ri / rw
CY3-031

Fallas Impermeables
Prueba de interferencia
Curva Tipo

Log
pD

ri/r

4
8

Log tD/rD2

CY3-032

Fallas Impermeables
Prueba de interferencia
Ajuste de curva tipo
4
Log
pD

ri/r

p
M
t
Log tD/rD2

CY3-033

Fallas Impermeables
Prueba de interferencia
Ajuste de curva tipo
Estimacin de
parmetros

Ajuste
(pD)M
( p)M
(tD/rD2)M
(t)M
(ri /r)M

kh=
=

141.2 q B

(pD)M

( p)M
r2 (tD/rD2)M
2.637x10-4 (t)M

ri = r (ri /r)M
CY3-034

Fallas Impermeables
Prueba de interferencia
Posicin de la falla

ri
O

La falla es tangente
a la elipse
Posibles localizaciones
del pozo imagen

CY3-035

Fallas Impermeables
Prueba de interferencia
Posicin de la falla

La falla es tangente a la elipse


y

2
2 1/2
((r
/2)
-(r/2)
)
b= i

r
O

a = ri /2
Ecuacin

x2/a2 + y2/b2 = 1
CY3-036

Fallas Impermeables
Determinacin de la posicin de una falla
Se requiere informacin de un mnimo
de 3 pruebas de un solo pozo o de
interferencia en pozos no alineados

A
O
O
CY3-037

Fallas Impermeables
Fallas mltiples

Perodos de flujo
* Radial
* Transicin
* Radial-sectorial

Log
t p'
Log (360/ )
Log t

CY3-038

Fallas Impermeables
Fallas paralelas

Perodos de flujo
* Radial
* Transicin
* Lineal

Log

1/2
1

t p'

Log t

CY3-039

Fallas Impermeables
Fallas en U

Perodos de flujo
* Radial
* Transicin
* Lineal

Log

1/2
1

t p'

Log t

CY3-040

Fallas Impermeables
Flujo Lineal

Falla
impermeable
df
Tipos de flujo:
* Radial
* Lineal
* Semilineal

CY3-041

Fallas Impermeables
Flujo Lineal
Anlisis

Diagnstico
p'
Log t

2mlf

1
1/2

p
1/2

mlf

Log t

t1/2
CY3-042

Fallas Impermeables
Flujo Esfrico

df

Falla
impermeable

Tipos de flujo:
* Esfrico
* Transicin
* Semiesfrico
CY3-043

Fallas impermeables
Flujo Esfrico
Anlisis

Log t

p'

Diagnstico
1
1

2msph

-1/2

1
msph

-1/2
1

Log t

t-1/2
CY3-044

Fallas Semipermeables
Flujo Radial

Yaxley

La falla exhibe resistencia al flujo debido a


que su permeabilidad (kf) es menor que la
permeabilidad de la formacin (k)

df
k > kf

k
kf

CY3-045

Fallas Semipermeables
Flujo Radial
La resistencia de la falla se caracteriza por
medio del parmetro adimensional:

= kf df / k bf

Transmisibilidad
adimensional de
la falla

o tambin:
Sf =

/2

Dao de la falla

bf es la amplitud de la falla

CY3-046

Fallas Semipermeables
Flujo Radial
Comportamiento de la presin
m

1
p(m,

pw
1

Log t
CY3-047

Fallas Semipermeables
Flujo Radial

1.151m p

A
CY3-048

Fallas Semipermeables

tD pD'

Curva Tipo

Falla
impermeable

0
.01
8

Flujo Radial

.1 .05

tDdf
CY3-049

Fallas Conductivas
Caracterstica
La falla permite flujo dentro y a travs
de su plano de acuerdo a su conductividad y dao

CY3-050

Fallas Conductivas
Parmetros de caracterizacin

Falla
kf , bf
df

Conductividad
adimensional
Factor de dao

Dao
kd , bd

FCD = kfbf / k df
Sf = (k bd / kd df

) ( /2)
CY3-051

Fallas Conductivas
Falla daada de conductividad
infinita

Falla
impermeable

tDf pwD'

100

-1

Sf

10

50

tDf
CY3-052

Fallas Conductivas
Falla daada de conductividad
infinita
Comportamiento de presin
pw = A - (1 + S / 2
2
f

q B ct

df2

3.73x10-6 k2 h

1
t

A es una constante que respresenta la


mxima caida de presin que se observa
en el pozo.
CY3-053

Fallas Conductivas
Falla daada de conductividad
infinita
3.73x10-6 k2 h mcp

df =

pw

q Bct

(1+Sf2/2

mcp
1

1/t

CY3-054

Fallas Conductivas

Falla
impermeable

tDf pwD'

102

-1

103
104
105
FCD
1

tDf

Falla conductiva daada

Sf

1/4

50

10

CY3-055

Fallas Conductivas

Falla
impermeable

tDf pwD'

102

-1

103
104
105
FCD
1

tDf

Falla conductiva daada

Sf

1/4

50

10

CY3-056

SAL.
111

114

115
105

103

6000

408

429
125 5500

121
6400

109

107

101 B

120

127
129

428
448

123
145

5500

301 A

119

117

149

5200

426
422

446

147

447

468

6400

169

62

6000

488
6400

58

56

14 D

NOMENCLATURA.

5500

8
26

5200

43

44

67
65

63
62

6400

83

49
69

6000

45

22 A

24

6000

6400

PENETRACION PARCIAL.

47

42

RADIAL HOMOGENEO.

FLUJO RADIAL COMPUESTO.

27

23 A

DOBLE POROSIDAD.

FLUJO LINEAL O BILINEAL .

29

3
2A

18

438

12

25

6000
6400

7
5

14

16

15

5800

13 A

38
439

32

34

36 A

6000

33

54

189

459

53

466

5500

167

444

5500

89

348

Acuamientos

También podría gustarte