2.
BALANCE DE MATERIA EN YACIMIENTOS DE PETRLEO SUBSATURADO
2.1
CONSIDERACIONES GENERALES.
La ecuacin de balance de materia, fue presentada por primera vez por Schilthuis, y
consista en un balance volumtrico, donde se igualaba el vaciamiento del reservorio
debido a la produccin, con la expansin de fluidos dentro el yacimiento.
Este concepto contina vigente desde entonces; bajo ciertas restricciones, el
balance de materia es una aplicacin de la ley de conservacin de masas. Las
principales restricciones son las siguientes:
a) No existe transferencia de masa entre fases
b) El reservorio es homogneo e isotrpico,
c) El reservorio es volumtrico y el sistema opera a volumen constante, aunque con
ciertas suposiciones se puede incluir la intrusin de agua.
Aunque el planteamiento original data de la dcada de los 40, su aplicacin es an
importante, podra suponerse que la tcnica de la simulacin matemtica de
reservorios confinara a la obsolescencia a este mtodo pero en realidad lo utiliza
como auxiliar muy confiable. Por ejemplo, se usa para verificar la consistencia de los
resultados de la corrida inicial, por medio de efectuar comparaciones de los
volmenes originales obtenidos por el simulador y los obtenidos por BM. Una
importante empresa que opera en Sudamrica la utiliza en la deteccin de la
comunicacin del yacimiento con otros niveles petrolferos que no fueron
considerados al disear el modelo matemtico, el BM puede detectar si hay un
intercambio de fluidos con estos niveles intercomunicados.
Los resultados obtenidos por el mtodo de balance de materia no deberan coincidir
forzosamente con los obtenidos por el mtodo volumtrico, aunque es razonable
esperar cierta similitud. La razn es que los clculos volumtricos trabajan con datos
por lo general aproximados y en base a suposiciones ms drsticas; el mtodo
volumtrico se utiliza con preferencia en las primeras evaluaciones de reservas,
muchas veces cuando solo hay un pozo descubridor y los datos que ste aporta son
generalizados a todo el yacimiento. Por otra parte la informacin geolgica y
petrofsica no es de gran exactitud.
A diferencia del mtodo volumtrico, el balance de materia requiere mas datos y de
mayor precisin. Al trabajar con la historia de produccin y presin, implcitamente
considera los cambios importantes que acontecen en el yacimiento y los efectos de
la heterogeneidad.
Cuando es extrada una cierta cantidad de fluidos y se produce una cada de
presin, acontecen cambios volumtricos cuyo anlisis permite deducir ecuaciones
tiles.
El caso ms general en un reservorio de petrleo es aquel donde inicialmente una
parte del gas est disuelto en el lquido y otra parte forma una capa de gas libre.
Ingeniera de Reservorios II
Ecuacin de Balance de Materia
La figura 8.1 muestra los cambios volumtricos que ocurren en el yacimiento donde:
A = Expansin del volumen de petrleo + gas disuelto,
B = Expansin del volumen del casquete de gas libre,
C = Contraccin del volumen poral de hidrocarburos debido a la expansin del
agua connata y la roca.
Casquete
de gas libre
Petrleo +
Gas disuelto
Fig. 8.1 Expansiones y contracciones de fluidos en el reservorio
Una roca reservorio est formada por granos de arena de sedimentacin unidos
entre s; entre ellos existen espacios que forman los poros de la roca. Al disminuir la
presin el material rocoso se descomprime, o en otras palabras se dilata, con lo que
el espacio de poros disminuye. Por otra parte, los granos de roca estn cubiertos por
la fase humectante que por lo general es el agua que al bajar la presin, tambin se
dilata, contribuyendo an ms al decremento del volumen poral.
El proceso de produccin de hidrocarburos crea un vaciamiento dentro el reservorio,
el cual es inmediatamente llenado por la expansin de los mismos. Esta expansin
se debe a que por efecto de la extraccin la presin en el yacimiento va
disminuyendo en forma continua.
Por lo tanto, la ecuacin de balance de materia se expresar a condiciones de
yacimiento.
Entonces, el balance de materia expresado en forma literal es:
Vaciamiento
del yacimiento
por la produccin
2.2
Expansin
del petrleo y
gas disuelto
Expansin
del gas libre
del casquete
Expansin del
agua connata y
reduccin del
volumen poral
8.1
ANLISIS DE EXPANSIONES VOLUMTRICAS.
Se puede expresar la ecuacin literal, como ecuacin algebraica, utilizaremos la
nomenclatura siguiente:
N=
Np =
Petrleo inicial in situ, bls
Produccin acumulada de petrleo, bls
Ingeniera de Reservorios II
Rp =
m=
Ecuacin de Balance de Materia
Relacin gas petrleo acumulada, pcs/bls
Relacin de volmenes originales de gas petrleo.
Fig. 8.2 Relacin de volmenes entre gas libre de casquete y petrleo
2.2.1
Expansin del lquido y gas disuelto.
Si se considera que el espacio que almacena al petrleo no cambia de volumen,
entonces la expansin del petrleo es la diferencia entre el volumen actual y el
volumen inicial, sea:
NBo NBoi N ( Bo Boi )
8.2
La expansin del gas liberado, es la diferencia de volumen entre el gas que
inicialmente estaba disuelto en el lquido y el gas que actualmente contina disuelto,
esto es:
NR si B g NR s B g N ( Rsi Rs ) B g
8.3
2.2.2 Expansin del gas libre del casquete.
El vaciamiento que resulta de extraer petrleo del yacimiento tambin es llenado en
parte con la expansin o dilatacin del gas libre contenido en el casquete original, al
disminuir la presin del reservorio.
El volumen inicial del gas libre en reservorio es:
mNBoi
El cual expresado en condiciones de superficie es:
mNBoi
B gi
Esta cantidad a una presin menor p, tendr el siguiente volumen dentro el
reservorio:
mNBoi B g
B gi
Entonces, la expansin del casquete de gas ser la diferencia entre el volumen
actual y el inicial, a condiciones de yacimiento, sea:
mNB oi (
Bg
- 1)
B gi
8.4
Ingeniera de Reservorios II
Ecuacin de Balance de Materia
2.2.3 Contraccin del volumen poral total.
Esta contraccin se debe a que el volumen del agua connata se expande por efecto
de la descompresin y tambin los granos de la roca, ocasionando la expulsin de
los fluidos que consiguientemente van a ocupar el espacio o vaciamiento dejado por
la produccin.
En consecuencia, la reduccin del volumen poral destinado a almacenar
hidrocarburos (VPHC) tiene la siguiente expresin:
d (VPHC ) dVw dV f
8.5
donde:
dVw = Expansin del agua connata,
dVf = Reduccin del volumen poral de la formacin,
los signos (-) se emplearon para indicar contraccin.
Fig. 8.3 Efecto de contraccin de los poros por dilatacin de la roca al bajar la presin
Aplicando la definicin de compresibilidad c
infinitesimal se tiene:
1 dV
, para una cada de presin
V dP
d (VPHC ) (c wV w c f V f ) dp
8.6
En esta ecuacin Vw representa el volumen total de agua connata, en tanto que V f es
el volumen de poros total de la formacin.
El volumen poral disponible para almacenar hidrocarburos VPHC es igual a
Ah (1 S wc ) . Por lo tanto, el volumen de poros total es:
Vf
VPHC
1 S wc
8.7
Similarmente, el volumen total de agua connata en el yacimiento es:
S wc
Vw VPHC
1 S wc
Reemplazando las ecuaciones (8.7) y (8.8) en (8.6), resulta:
d (VPHC ) VPHC
c f c w S wc
1 S wc
dp
8.8
8.9
Ingeniera de Reservorios II
Ecuacin de Balance de Materia
Por otra parte, considerando que el volumen poral total es el reservorio donde se
almacenan los hidrocarburos incluyendo el gas libre de casquete y el petrleo, se
puede establecer que:
VPHC mNBoi NBoi (1 m) NBoi
8.10
Reemplazando este valor en la ecuacin (8.9), se tiene que la contraccin del
volumen poral del reservorio es:
c w S wc + c f
d(VPHC) = (1 + m) NB oi (
)dp
8.11
1 - S wc
2.2.4 Vaciamiento del reservorio.
El proceso de extraccin de hidrocarburos causa un vaciamiento en las zonas que
rodean a los pozos productores; stos espacios son inmediatamente llenados por
las expansiones anteriormente analizadas.
El vaco que deja el petrleo es volumtricamente igual a la produccin acumulada,
que expresada a condiciones de yacimiento es:
NpBo
8.12
El vaco que deja el gas disuelto extrado, es la diferencia entre el gas producido y el
gas remanente que contina disuelto en el petrleo, o sea: N p(Rp-Rs); expresando
este volumen a condiciones de reservorio, resulta:
Np(Rp-Rs)Bg
8.13
2.2.5 La Ecuacin General de Balance de Materia.
Combinando las anteriores ecuaciones de acuerdo a la ecuacin literal expresada
como (8.1), se tiene:
N p Bo N p ( R p Rs ) B g N ( Bo Boi ) N ( Rsi Rs ) B g mNBoi (
(1 + m) NBoi (
Bg
B gi
1)
c w S wc + c f
)dp
1 - S wc
Efectuando algunas transformaciones algebraicas, se tiene:
N p [ B o + ( R p - R s ) B g ] = NB oi [(
+ (1 + m)(
Bo - B oi + ( R si - R s ) B g
B
) + m( g - 1)
Boi
B gi
c w S wc + c f
)p + ( W e - W p ) B w
1 - S wc
8.14
Donde dp ha sido reemplazado por p para variaciones finitas de la presin. El
trmino (W e-W p)Bw es el influjo neto de agua en el yacimiento.
2.3
LINEARIZACIN DE LA EBM.
Ingeniera de Reservorios II
Ecuacin de Balance de Materia
Fue planteada inicialmente en la dcada de los 60 por Odeh y Havlena. La
linearizacin permite expresar la ecuacin de balance de materia en una forma ms
sencilla y de fcil aplicacin a los distintos tipos de empuje.
Para este efecto se hacen los siguientes cambios de variables:
Vaciamiento del reservorio:
F = NpBo+Np(Rp-Rs)Bg + W pBw
Expansin del petrleo y su gas disuelto:
Eo = (Bo-Boi)+(Rsi-Rs)Bg
8.15
8.16
Expansin del gas libre del casquete:
E g = B oi (
Bg
- 1)
B gi
Reduccin en el volumen poral:
c w S wc + c f
)p
E f,w = (1 + m) NB oi (
1 - S wc
8.17
8.18
Reemplazando las relaciones (8.15), (8.16), (8.17) y (8.18) en la ecuacin general
de balance de materia (8.14), sta adopta la siguiente forma:
F = N(Eo+mEg+Ef,w)+W eBw
8.19
que es la forma lineal de la EBM
El objetivo principal del balance de materia es la determinacin de los volmenes
originales in situ, el factor de recuperacin y otros parmetros como ser, la relacin
de volmenes m y el volumen de intrusin de agua.
La forma linearizada es cmoda para estudiar los diferentes tipos de empuje. A
continuacin se harn las consideraciones correspondientes a los siguientes
mecanismos:
Empuje por gas disuelto
Empuje por la expansin del gas libre
Empuje por intrusin del acufero
.
2.4
EMPUJE POR GAS DISUELTO.
Tambin llamado mecanismo por agotamiento, donde el principal elemento es el gas
disuelto en el petrleo.
Se distinguen dos etapas:
a) Cuando el yacimiento es sub-saturado.
b) Cuando la presin cae debajo del punto de burbuja.
a) Petrleo sub-saturado.
Ingeniera de Reservorios II
Ecuacin de Balance de Materia
El trmino sub-saturado indica que el petrleo no est completamente saturado de
gas lo que significa que todava puede disolverse cierta cantidad en l. En estas
condiciones el yacimiento se encuentra sobre el punto de burbuja y por lo tanto los
hidrocarburos se encuentran formando una sola fase. Obviamente no existe
casquete de gas original.
En esta primera etapa los lquidos, petrleo y agua, se encuentran comprimidos al
momento de descubrir el reservorio. A medida que se desarrolla la produccin, los
fluidos salen a superficie principalmente por la descompresin; este perodo tiene
una corta duracin y se puede notar que el yacimiento declina muy rpidamente.
Si un yacimiento es clasificado como productor con mecanismo de gas disuelto,
obviamente es porque no existe el empuje de un acufero activo porque si lo hubiera,
ste sera el predominante y el que caracterice el mecanismo productor. Por esta
razn, el parmetro W e es cero.
Puesto que tampoco existe casquete de gas inicial, el valor de la relacin de
volmenes m es cero. Entonces, la ecuacin de balance de materia (8.19) para un
yacimiento con empuje de gas disuelto sub-saturado, en su forma lineal es:
F = N(Eo+Ef,w)
8.20
Por otra parte, puesto que no hay liberacin de gas en reservorio, el factor de
solubilidad Rs permanece constante, y adems es igual a la relacin gas/petrleo
obtenida de la produccin.
Esto significa que: Rp = Rsi = Rs
Por lo tanto, expandiendo la ecuacin lineal (8.20) y aplicando las consideraciones
anteriores, la EBM se simplifica notablemente del siguiente modo:
Bo - Boi c w S wc + c f
+
)p
8.21
N p B o = NB oi (
1 - S wc
B oi
El factor de compresibilidad del petrleo es por definicin:
1 dV o
co = V o dp
El mismo que puede ser expresado en funcin de los factores volumtricos de la
siguiente forma:
1 dBo
co =
Bo dp
Si se trata de un intervalo finito de presin, se tiene:
1 Bo - Boi
co =
Boi p
Reemplazando correspondientemente en la ecuacin (8.21):
Ingeniera de Reservorios II
N p Bo = NB oi ( co p +
Ecuacin de Balance de Materia
c w S wc + c f
)p
1 - S wc
8.22
En la zona petrolfera de un yacimiento subsaturado, evidentemente slo existe agua
connata y la fase petrleo. Si se expresa esto en funcin de sus saturaciones, se ve
que:
So+Sw = 1
por lo tanto, haciendo operaciones algebraicas, resulta:
co S o + c w S wc + c f
)p
8.23
N p B o = NB oi (
1 - S wc
Factor de Recuperacin
La relacin Np/N que representa la fraccin de petrleo recuperado constituye el
factor de recuperacin FR, que tambin puede expresarse en porcentaje.
Al inicio de la vida productiva no se conoce el valor de la produccin acumulada, y
por lo tanto no se puede determinar el factor FR.
Una aproximacin muy utilizada es la correlacin de Arps que tiene la siguiente
expresin:
pb 0.1741
(1 - S w ) 0.1611
k
0.0979 0.3722
FR = 41.8150[
][
]
]
]
Sw [
Bo
1.00 o
pa
donde:
pb =
punto de burbuja, psi
pa =
presin de abandono, psi
2.6
EMPUJE POR AGUA
El empuje natural del acufero, diferente de la inyeccin de agua fue considerado al
tratar la ecuacin general de balance de materia. Una cada de presin en el
yacimiento debido a la produccin de petrleo ocasiona que el acufero se expanda
y se introduzca dentro del yacimiento desplazando al petrleo.
Esta expansin de agua que se introduce al yacimiento se conoce con el nombre de
intrusin o influjo de agua, que en su forma ms simple se puede representar por la
expresin:
W e = cpW i
Donde c, que es el factor compresibilidad total es la suma de la compresibilidad del
agua y la de los poros debido a que el espacio poral de la zona del acufero est
ntegramente lleno de agua.
A menos que el volumen W i sea muy grande, la influencia del acufero es pequea y
generalmente despreciable. Si W i es muy grande, el influjo de agua es considerable
y merece un anlisis especial.
Ingeniera de Reservorios II
Ecuacin de Balance de Materia
Fig. 8.4 Empuje de agua, dos incgnitas
Usando la tcnica de Odeh-Havlena y asumiendo que B w = 1 la ecuacin de balance
de materia es:
F = N(Eo+mEg+Ef,w) + W e
8.26
Cuando hay accin del acufero, la tendencia es que la presin del reservorio se
mantiene casi constante; por esto, siendo p pequeo, el trmino E f,w es tambin
pequeo y puede despreciarse. Si adems no existe gas libre, la ecuacin se reduce
a:
F= NEo + W e
8.27
El empuje de agua es un mecanismo cuya energa se origina en la intrusin de
agua, ya sea por la expansin del volumen o por la energa potencial. Al tratar el
problema de intrusin de agua es necesario definir el modelo de acufero como
principal parmetro.
Para reproducir con sta ecuacin el historial de produccin existe cierta
incertidumbre en cuanto a la determinacin del influjo de agua. Para el clculo de W e
se requiere conocer la geometra y el modelo del acufero que puede obtenerse
mediante la informacin geolgica, aunque sta no siempre se encuentra completa
puesto que la perforacin llega a atravesar la zona de agua nicamente en los
pozos exploratorios y algunos pozos de avanzada y no as en los pozos de
desarrollo.
Por ejemplo, si se aplica el modelo radial debe definirse la relacin entre radios
acufero/reservorio re/rw. El tiempo adimensional y el factor de intrusin pueden
estimarse o asumirse. Entonces es posible aplicar la relacin planteada por Hurst y
van Everdingen cp We(tD).
2.6.1 Empuje de agua sin Casquete de Gas Original
Ingeniera de Reservorios II
Ecuacin de Balance de Materia
La ecuacin anterior puede entonces expresarse de la siguiente manera:
c pWe ( t D )
F
=N+
Eo
Eo
8.28
Se pueden seguir los siguientes pasos:
1) asumir un re/rw y un tiempo adimensional tD.
2) Calcular p W e(tD),
3) dibujar F/Eo versus cp W e(tD)/Eo.
Si los valores asumidos son correctos, la grfica resultante ser una lnea recta,
siendo N la ordenada al origen y C el valor de la pendiente. Pero antes de hallar la
solucin correcta habr muchos intentos con valores errados, las grficas no sern
lneas rectas sino curvas. Si la curva est situada por encima de la recta que
representa la solucin, el valor de p We(tD) es muy pequeo. Por otra parte, si la
curva est por debajo, quiere decir que el valor p We(tD) es demasiado alto. Una
curva semejante a una S alargada indica que el ajuste sera mejor si se asumiera si
se asumiera influjo lineal de agua.
Havlena y Odeh plantean una prueba llamada de consistencia que se describe del
siguiente modo: Se escogen varios valores de tD prximos al punto que da la
mnima desviacin estndar y se lee su valor.
Para cada valor seleccionado de tD, se calculan N y C como funciones de tiempo
real. Entonces se construyen grficas de N versus tiempo real y C versus tiempo real
y por medio de mnimos cuadrados, se hace el ajusta a una lnea recta para cada
grfica.
Se determinan las pendientes de las lneas rectas de N y C y se dibujan versus su
correspondiente valor de delta tD en una misma grfica. La interseccin de estas
dos curvas da el valor ms probable de tD y tericamente esta interseccin debera
coincidir con la pendiente cero.
En efecto, si re/rw y tD son correctos, y si tambin los dems datos, entonces N y C
no deberan variar con el tiempo.
BIBLIOGRAFIA
1. D Havlena, A.S. Odeh, "The Material Balance as an Equation of a Straight Line"
Journal of Petroleum Technology, August 1963, pag 896-900.
2. [Link], "Fundamentals of Reservoir Engineering", Elsevier Science Publisher,
1978.
3. [Link], A.H. Harvey, M. Honarpour, "Introduction to Petroleum Reservoir
Analysis", Gulf Publishing Company, Houston, Texas, 1989.
4. F. Kuppeand, S. Chugh and P. Connel, "Material Balance for Multi-layered.
