DETERMINACION DE LA ECUACION QUE REPRESENTA LA ENTRADA DE
AGUA EN EL YACIMIENTO
Para aplicar la EBM se debe determinar la entrada de agua, We, como funcin del
tiempo si el Balance de Materia se usa para determinar los volmenes originales
de aceite y gas. We, depende del tamao del acufero y de la cada de presin
entre el acufero y el yacimiento. En este trabajo se presentarn dos mtodos para
determinar We: el mtodo de flujo en rgimen permanente de Schilthuis y el
mtodo de flujo en rgimen pseudopermanente de Fetkovich
Determinacin de los parmetros del acufero mediante Balance de Materia
La constante del acufero y los volmenes originales de aceite o gas tambin se
pueden determinar de grficas de Balance de Materia. La EBM se puede escribir
de la forma:
( )
fw gi fgi OI fgi e g fgi O foi P W
E B G B N W E G E N W B F + + + + = +
en donde, para determinar la entrada de agua se tiene:
( )
fw gi fgi OI fgi E g fgi O foi P W
E B G B N Pi W B E G E N W B F + + A + + = +
0
para yacimientos de aceite bajosaturado, la ec. se reduce a:
( )
A + + = + Pi W B E B E N W B F
E fw gi O foi P W 0
la ec. 4.9 se puede rearreglar para que quede:
fw oi
n
j
e
foi
fw gi O
P W
E B E
Pi W B
N
E B E
W B F
+ +
A
=
+
+
0
1
0
}
Mtodo de Schilthuis
El modelo ms simple para la entrada de agua es el mtodo de Schilthuis el cual
considera que el ritmo de entrada de agua es proporcional al decremento de
presin (pi-p).
( ) P Pi K
dt
dWe
=
Integrando la ecuacin anterior.
( ) P Pi K We
t
=
}
0
Donde k es una constante de proporcionalidad en unidades de RB/(D-psi). En el
SI, k puede tener unidades de res m3/(s Pa). Se usa k en lugar de k para no
confundirlo con la permeabilidad. Es difcil predecir un valor a priori para k, pero
se puede medir a partir de los datos de presin y produccin del yacimiento. Una
manera de hacer esto es calculando primero We de la EBM.
( )
fw gi fgi OI fgi g fgi O foi P W
E B G B N E G E N F W B We + + =
La ecuacin requiere que Nfoi y Gfgi sean conocidos de volumetra. Despus, por
medio de una grfica We/t vs. (pi-p), se determina k a partir de la pendiente de
la grfica, m. Si los volmenes originales de aceite y gas son conocidos, k se
puede determinar si la produccin del yacimiento se puede limitar por un periodo
en donde la presin se mantenga relativamente constante. Esto es, dp/dt=0, lo
cual significa que el volumen de fluidos producidos es igual al volumen de entrada
de agua.
t
Wp B F
dt
W B F d
dt
dWe
W
p w
A
A + A
=
+
=
) (
para yacimientos de aceite negro (Rv=0), esto se puede simplificar a:
( )
dt
dGp
BG
dt
dNp
B R B
dt
dWp
Bw
dt
dWe
g S
+ + =
0
para presiones mayores a la de burbujeo, se puede simplificar a:
dt
dNp
B
dt
dWp
B
dt
dWe
W 0
+ =
entonces el valor de dWe/dt se divide por (pi-p) para determinar k.
Mtodo de Fetkovich
Este mtodo asume que el tamao del acufero es conocido y que el agua que
fluye del acufero de presiona al yacimiento de acuerdo con la EBM. Debido a que
se asume un estado pseudopermanente en este mtodo, se debe especificar un
tamao del acufero finito.
Este mtodo se aplica de una manera simple en clculos de Balance de Materia. A
continuacin se describen los pasos para aplicar en el mtodo:
1. Se calcula el volumen de agua acumulable mximo de entrada, Wei,
( )
360
0
61458 . 5
2 2
| t h r r c p
w
R aq i i
ei
=
donde Wei est dado por RB si h, raq, y rR, estn en pies.
2. Se calcula el ndice de productividad, J, para el flujo del acufero al yacimiento
con:
360
0
75 . 0
007082 . 0
(
=
R
aq
r
r
In
kh
j
para acuferos finitos o con una presin constante en las fronteras del acufero:
360
0 007082 . 0
(
=
R
aq
r
r
In
kh
j
las dos ecuaciones tienen unidades de campo, J est dado por RB/(D-psi) en SI o
en res m3/(d Pa).
3. Inicialmente se toma la presin del yacimiento al final del tiempo n, pRn; la
presin media durante el tiempo n, pRn ; y la presin media del acufero al final del
tiempo n, paqn , para la presin inicial del yacimiento, pi. Se miden o se
especifican pRn y n t .
4. Se calcula la presin media del yacimiento para cada intervalo de tiempo:
2
1 Rn Rn
Rn
P p
p
+
=
5. Calcular la entrada de agua para cada intervalo de tiempo:
( )
(
= A
A
Wei
t Jpi
Rn n
i
ei
en
n
P P
P
W
W 0 1
1
6. Se calcula la entrada total acumulada de agua al tiempo actual:
A =
n
i j
ej e
W W
7. Calcular la presin media en el acufero al trmino del tiempo actual:
(
=
ei
e
i n
W
W
p P 1
8. Regresar al paso 4 y repetir para el siguiente intervalo de tiempo.
Para los ejemplos de este trabajo se supusieron tanto Jw como Wei, empleando
un proceso iterativo para el mtodo de Fetkovich.
Funcin Influencia del Acufero
Para un acufero de geometra arbitraria, que produce a gasto constante, la
relacin entre la cada de presin y el gasto es proporcional, a esa constante de
proporcionalidad se le conoce como Funcin Influencia del Acufero:
( ) ( ) ( ) t F t e t p P
w o
=
para un acufero de geometra arbitrara pero que no produce a gasto constante la
relacin entre la cada de presin, el gasto y la Funcin Influencia del Acufero se
puede obtener aplicando el principio de superposicin, el cual establece que un
sistema es lineal la respuesta a la suma de dos estmulos es la suma de las
respuestas a cada uno de ellos, el cual se representa mediante la integral de
convolucin
( ) ( )
( )
( )
}
= =
t
w w o
d
d
dF
t e
dt
t dF
t e t p P
0
t
t
t
t
Aplicando la transformada de Laplace a la ecuacin anterior
( )
( )
( ) { }
( )
)
`
=
)
`
dt
t dF
L t e L
dt
t dF
t e L
w w
desarrollando:
( )
( )
( ) ( ) ( ) ( ) 0 F s sF s e
dt
t dF
t e L
w w
=
)
`
Considerando que F(0) = 0 y (0) = 0 w e la ecuacin anterior se puede escribir
como:
( )
( )
( ) ( ) ( ) ( ) ( ) s F e s se
dt
t dF
t e L
w w w
0 =
)
`
Aplicando la transformada inversa de Laplace:
( )
( )
( )
}
=
t
w
w
t F
dt
t de
dt
t dF
t e
0
donde:
p 0 = presin inicial en el contacto original acufero-yacimiento, psia.
p = presin en el contacto original acufero-yacimiento, psia.
e w = gasto de agua en el contacto original acufero-yacimiento, PCD para
yacimientos de gas y BPD para yacimientos de aceite.
F = Funcin Influencia del Acufero, psia/PCF, para yacimientos de gas y psia/BPD
para yacimientos de aceite.
Hutchinson y Sikora, 1959.
El primer mtodo para describir el comportamiento de un acufero sin la necesidad
de suponer sus propiedades fue propuesto por Hutchinson y Sikora en 1959,
adems de que estableci que la Funcin de Influencia del Acufero es una
solucin de la ecuacin de difusin aplicando el principio de superposicin para
fluidos ligeramente compresibles (agua o aceite bajosaturado) para gasto
constante.
Establecen que el valor de la FIA debe siempre ser positivo, adems de cncavo;
para calcular el valor de la FIA proponen un mtodo iterativo, basado en la
discretizacin:
( )
=
A = A
+
n
l j
j w n
F e P
j n 1
el cual toma como primera aproximacin para el valor de la FIA a:
n n
F
=
Donde:
n = nmero de mes en que se est calculando la Funcin Influencia del Acufero,
FIA.
El proceso iterativo mejora a partir de un factor de correccin:
( )
n
cal n
n
p
p
s
A
A
=
donde:
( )
cal n
p A es el valor de la cada de presin calculada con la Ecuacion.
Con los valores calculados de DF y los valores obtenidos para el factor de
correccin, se estiman los valores de la funcin de influencia para la siguiente
iteracin:
( )
( )
,..... 3 , 2 , 1 . = = n
S
F
F
n
cal n
mejorada
n
El proceso se repite hasta que los valores de n S tiendan a 1.
Katz y Cols., 1962
Katz y Cols, 1962, present una extensin y modificacin del trabajo presentado
por Hutchinson y Sikora para determinar los valores de la Funcin Influencia del
Acufero para yacimientos de gas con acufero asociado.
Su mtodo se basa en la solucin de la Ecuacin de Difusin para fluidos
ligeramente compresibles que producen a gasto constante:
( ) ( )
D D w l
t P e K t p p =
0
Donde:
h k
k
w
w
l
15 . 25
=
2
00633 . 0
b wt w
w
D
r c
t k
t
|
=
wf w wt
c c c + =
donde:
m w = viscosidad del agua del acufero, cp.
k w= permeabilidad del acufero, md.
h = espesor del acufero, pie.
t = tiempo, das.
c wt= compresibilidad de acufero, psi-1.
c w =compresibilidad del agua, psi-1.
c wf= compresibilidad de la formacin en la zona del acufero, psi-1.
El valor de la Funcin Influencia es:
( ) ( )
D D l
t P K t F =
Katz y cols establecieron que para el clculo de la Funcin Influencia se tiene que
considerar que la porcin inicial de la curva es ms sensible a errores provocados
por malas mediciones de la presin, por lo que propuso que los primeros datos
estn sujetos a la desigualdad siguiente:
| |
c
c
A
s A s
(
1 1
2
1
1
1
n
n n
Z
Z
n
n
Z
Los autores observaron que el valor de que funciona mejor es el de 0.02.
Katz y cols establecen adems de que la FIA debe ser siempre positiva y cncava,
esdecir, debe de cumplir con las condiciones siguientes:
0
) (
0
) (
0 ) (
2
2
s
>
>
dt
t F d
dt
t dF
t F
NOMENCLATURA
A= pendiente
b= punto de interseccin
b= ordenada al origen de una recta
B= constante de correlacin
Bg= factor de volumen de formacin del gas al tiempo t, L3/L3, RB/scf
Bgi= factor de volumen de formacin del gas inicial, L3/L3, RB/scf
Bo= factor de volumen de formacin del aceite al tiempo t, L3/L3, RB/STB
Bob= factor de volumen del aceite en el punto de burbujeo, L3/L3, RB/STB
Bofb= factor de volumen del aceite en el punto de burbujeo a partir una prueba de
separacin flash, L3/L3, RB/STB
Bod= factor de volumen del aceite a partir una prueba DL, L3/L3, RB/STB
Bodb= factor de volumen del aceite en el punto de burbujeo a partir una prueba
DL, L3/L3, RB/STB
Boi= factor de volumen de formacin del aceite inicial, L3/L3, RB/STB
Btg= factor de volumen de formacin del gas total al tiempo t, L3/L3, RB/scf
Bto= factor de volumen de formacin del aceite total al tiempo t, L3/L3, RB/STB
Bw= factor de volumen de formacin del agua, L3/L3, RB/STB
c= compresibilidad isotermal, m/Lt2, psi
C= constante de correlacin
cf= compresibilidad isotrmica de la formacin, Lt2/m, psi-1
ct= compresibilidad isotrmica total definida en la, Lt2/m, psi-1
cw= compresibilidad isotrmica del agua, Lt2/m, psi-1
D= constante de correlacin
Efw= factor de expansin del agua de formacin, L3/L3, RB/STB
Eg= factor de expansin del gas, L3/L3, RB/scf
Eo= factor de expansin del aceite, L3/L3, RB/STB
F= factor definido por la ecuacin 2.34
F= fluidos producidos del yacimiento, L3, res bbl
G= gas inicial total en el yacimiento, L3, scf
GI= gas de inyeccin acumulado al tiempo t, L3, scf
Gfg= gas en fase libre al tiempo t, L3, scf
Gfgi= gas inicial en fase libre, L3, scf
Gp= produccin acumulada de gas al tiempo t, L3, scf
Gpg= produccin acumulada de gas de la fase vapor al tiempo t, L3, scf
IDD= ndice de empuje por agotamiento
ISD= ndice de empuje por segregacin
IFD= ndice de empuje de formacin
IWD= ndice de empuje por entrada de agua
J= funcin J de Leverett
k= permeabilidad promedio, L2, md
krg= permeabilidad relativa del gas
kro= permeabilidad relativa del aceite
krw= permeabilidad relativa del agua
K= constante definida por la ecuacin 2.24
m= pendiente de una recta
m= masa, m, lbm
M= peso molecular, m
n= nmero de moles, n
N= aceite inicial total en el yacimiento, L3, STB
Nfo= aceite en fase libre al tiempo t, L3, STB
Nfoi= aceite inicial en fase libre, L3, STB
Np= produccin acumulada de aceite al tiempo t, L3, STB
p= presin actual del yacimiento, m/Lt2, psia
pi= presin inicial del yacimiento, m/Lt2, psia
BIBLIOGRAFIA
- FACULTAD DE INGENIERA Dr. Fernando Samaniego Verduzco TESIS
ECUACIN GENERAL DE BALANCE DE MATERIA Y UNA
HERRAMIENTA COMPUTACIONAL PARA SUS CLCULOS
- FACULTAD DE INGENIERIA TESIS LA FUNCIN INFLUENCIA DEL
ACUFERO Y SU UTILIZACIN EN LA PREDICCIN DEL
COMPORTAMIENTO DE YACIMIENTOS DE GAS Y ACEITE
BAJOSATURADO
