Existen diferentes tipos de reservas:

Reservas Probadas: aquellas que han sido cuantificadas por medio de

pozos, equipos y métodos técnicos específicos que garantizan su existencia.

Reservas Probables: aquellas que han sido certificadas por medio de

pruebas de producción, pero que por encontrarse dentro de los limites
geológicos conocidos de un yacimiento, son susceptibles de ser probadas
perforando pozos adicionales y haciendo pruebas de producción.

Reservas Posibles: aquellas de posible existencia pero que por falta de

información fehaciente no pueden garantizarse con exactitud.

Reservas por Descubrir: son cantidades de hidrocarburos que en una

fecha determinada, se estiman. Podrían existir en acumulaciones en
yacimiento aun no descubierto, pero que se presume su existencia sobre la
base de la información de geología en la superficie, sensores remotos,
gravimetría, sísmicas y otros métodos.

Los métodos para la estimación de reservas son:

 Método Volumétrico.
 Balance de Materiales.
 Curva de Declinación de Producción.
 Simulación Matemática.
 Método probabilístico.

MÉTODO VOLUMÉTRICO

El método volumétrico permite la estimación del petróleo y del gas

originalmente en sitio (POES) y (GOES) a partir de la determinación del
volumen de la roca que conforma el yacimiento, la capacidad de
almacenamiento de la roca y de la fracción de hidrocarburo presente en los
poros de dicha roca. Este requiere de la definición de la parte sólida del
yacimiento y el conocimiento de las propiedades físicas de los fluidos y área
de drenaje.

1
La ecuación para determinar el método volumétrico es la siguiente:

Donde:

V  Volumen de la roca (Acre-pie)

Porosidad promedio (fracción)

Swc  La saturación promedio del agua connata.

Boi  Factor volumétrico del petróleo a presión inicial By/Bn

N Petróleo original en el sitio POES.

G 43560 V  (1- Swc)

Bgi

G= Gas original en sitio.

V= volumen de la roca.

= Porosidad promedio (fracción).

Swc= La saturación promedio del agua connota.

Bgi= factor volumétrico del gas a presiones iniciales( PCY/PCN).

2
 Estimación Del Volumen De La Roca.

El cálculo del volumen de la roca puede hacerse mediante mapas

estructurales o isópacos, que son esenciales para determinar los hidrocarburos
en el sitio.

Los mapas estructurales identifican el tope y la base de las formaciones y/o

arenas del yacimiento que son utilizadas para determinar el volumen total de
la roca ; de la misma forma en el mapa isópodo se representa la cantidad de
arena neta permeable o zona efectiva de hidrocarburos dentro de ese volumen
total identificado en el mapa estructural.

Su solución puede obtenerse a través de tres métodos: por el método

gráfico, el método tabular y el método practico.

 Método Gráfico

El método gráfico para el cálculo del área bajo la curva, consiste en

dividir la curva en pequeñas porciones que se asemejan a figuras geométricas
conocidas.

Una vez calculada el área de cada una de estas porciones, el volumen de roca
será la suma de todas ellas.

Figura 3. Porción del área bajo la curva semejante a un trapecio

3
  Método Tabular:

El método tabular para el cálculo de área bajo la curva, consiste en leer en

dicha curva los valores de área correspondientes a una serie de espesores
igualmente espaciados.

El volumen encerrado entre dos isópacas se calcula a través de la

aproximación piramidal y por la aproximación trapezoidal.

*Aproximación Piramidal: en este caso se aproxima el volumen entre dos

curvas de igual espesor, por el de una pirámide de la siguiente manera.

Donde:
d: diferencia de los valores de dos líneas isópacas consecutivas (constante)
An: Área en acres encerrada por línea isópaca inferior.
An+1: Área en acres encerrada por la línea isópaca superior.

*Aproximación Trapezoidal: este método aproxima el volumen entre dos

curvas de igual espesor por una relación similar a un trapecio.

V= d (An + An+1)

2
En general:

Donde:
d: Intervalos de la líneas isópacas en pies.
Ao: Área encerrada por línea isópaca 0.
A1, A2,..., An = Áreas encerradas por las líneas isópacas consecutivas.
Tave = Espesor promedio de la base de la estructura.

4
Fig. 4 Sección esquemática en láminas horizontales para el cálculo de método
trapezoidal.

Las A0, A1 y A2…An, corresponden a las áreas dentro de un límite isópaca. El

volumen total bruto es igual a la sumatoria del volumen de todos los
trapezoides equivalentes.

Vb = ∑ [ (AI + Ai+1)/2] * (Ahi, i+1)

Donde:

Vb = volumen total bruto.

Ai = área de la cara inferior de cada lámina horizontal
Ai+1 = área de la cara superior de cada lámina horizontal
Ah, i+1 = epesor de la lámina (diferencia de espesor entre líneas isópacas que
corresponden a AI y Ai+1).

La aproximación trapezoidal introduce un error del orden del 2 % cuando la

razón entre sus áreas es de 0,5. Por ello, si dicha razón es inferior al valor de
0,5, se recomienda utilizar la aproximación piramidal.

 Método Práctico

Este método es fundamentalmente una variación del método trapezoidal.

La diferencia fundamental es que para fines del cálculo el yacimiento se
divide en láminas verticales.

5
 Cada lámina vertical tiene un área base que se mide con planímetro y que
equivale al área entre dos líneas de isoespesor. La altura promedio
correspondiente a cada área así medida es el espesor algebraico promedio
entre las dos líneas de isoespesor que limitan cada área medida. Así, se ha
aproximado el volumen de cada lámina vertical por el volumen del trapezoide
equivalente.

Fig 5. Sección esquemática para el método práctico.

El volumen bruto de cada lámina se computa como el del trapezoide

equivalente a(i, i + 1) [ (hi + hi+1) / 2.

El volumen bruto total es igual a la sumatoria del volumen de todos los

trapezoides equivalentes tomando las láminas forma vertical:

Vb =  (Ai, Ai+1) * ( hi + hi+1) /2

Donde:
Vb = volumen bruto total
Ai, Ai+1 = área entre dos líneas isópacas de espesor hi y hi+1.
hi = espesor correspondiente a una línea isópaca.
hi+l = espesor correspondiente a la siguiente línea isópaca.

6
 El método práctico tiene una precisión del mismo orden de magnitud que
el método trapezoidal, la cual es menor que la del método piramidal.

Mapas Estructurales

La estimación del volumen de roca mediante mapas estructurales

requiere el cálculo del área encerrada por líneas de igual profundidad, tanto
para el mapa del tope, como para el mapa de la base de la estructura.

Así, un gráfico de profundidad versus estas áreas produciría dos curvas

entre las cuales estará comprendido el volumen de roca del yacimiento.

De tenerse identificados los contactos gas-petróleo y agua-petróleo,

estos serían los límites para el cálculo del volumen de roca ocupada por
petróleo, que es el requerido para el cálculo POES.

Curvas de Profundidad vs Áreas obtenidas de los mapas del tope y la

base del yacimiento. En dicho gráfico se señala un contacto gas-petróleo y un
contacto agua-petróleo, como límites de la zona ocupada por petróleo.

Propiedades Promedio Ponderadas De La Roca Y Los Fluidos.

Obtenido el volumen bruto del yacimiento, el próximo paso lo constituye

la estimación del volumen poroso interconectado, el cual representa a
condición de yacimiento el espacio susceptible de contener fluidos.

7
 Así:
Vp = Vb *

Donde:

Vp = Volumen poroso interconectado total.

 = porosidad interconectada promedio ponderada volumétricamente.

Vb = Volumen total de la roca.

a) Porosidad promedio ponderada volumétricamente

Con el fin de ponderar volumétricamente la porosidad (o cualquier otra

propiedad del yacimiento), se requiere un mapa isópaco del mismo, así como
un mapa de isovalores de la propiedad que se desea ponderar. Superponiendo
ambos y diseñando una retícula, se podrán asignar valores de la propiedad
dada a cada bloque de la retícula. Finalmente se planimetría el área de cada
retícula, se asigna el espesor promedio correspondiente, y se computa el valor
promedio ponderado volumétrico de la propiedad.

Por ejemplo la porosidad se tiene:

 =  ( hi * i*Ai)
 (hi * Ai)
Donde:
 = Porosidad promedio ponderada volumétricamente.
Ai = Área de cada celda o bloque de la retícula.
hi = Espesor promedio asignado a cada celda o bloque de la retícula.
i = Porosidad asignada a cada celda o bloque de la retícula.

b) Saturación promedio ponderada volumétricamente

Utilizando el procedimiento descrito anteriormente, se puede pensar en

ponderar la saturación de cualquier fluido (gas, agua y/o petróleo) utilizando

8
el mismo concepto de retícula pero con dos posibilidades:

 Ponderar utilizando el volumen bruto de cada celda.

 Ponderar utilizando el volumen poroso de cada celda.

Aunque lo último es lo más deseable, no siempre se dispone de un mapa

de isopropiedades. El cálculo se efectúa mediante la aplicación del siguiente
formato. En este caso se muestra para la determinación de la saturación de
petróleo.

So = ( ( Vbi * i*Soi) / ( Vbi * i*Soi )

Donde:
So = Saturación promedio ponderada volumetricamente.
Vbi = Volumen bruto de cada celda.
i = porosidad de cada celda.
Soi = saturación de petróleo en cada celda.

Cálculos de volúmenes originalmente en sitio para cada uno de los

yacimientos.

Yacimiento volumétrico de gas seco.

El método volumétrico es el más generalizado, ya que presenta la única

forma de estimar un gas original en sitio (GOES) en yacimientos nuevos,
recién descubiertos sin historia de producción. Consiste en obtener el volumen
de arena gasífera, mediante el planímetro de los mapas isópacos, que
multiplicados por los datos Petrofísicos de porosidad y saturación inicial de
hidrocarburo se obtiene el GOES.

9
 El método consiste en determinar el volumen del espacio poroso en el
yacimiento que contiene gas, a condiciones estándar. El volumen estándar de
GOES es simplemente el volumen de tres factores:

1. El volumen poroso del yacimiento

2. Saturación del gas
3. Factor volumétrico del gas

Volumen poroso: el volumen poroso ocupado por el gas también puede

cambiar debido a la intrusión de agua en el yacimiento.

Saturación del gas: es la fracción del espacio en el yacimiento ocupado por el

gas libre.

Factor volumétrico: es el factor que representa el volumen del gas libre a

presión y temperatura del yacimiento, por unidad volumétrica del gas libre a
condiciones normales.

 PCY   PCN 
GOES  43560   *  * 1  Swi  *Vb ac. p  * Bgi  
 ac. p   PCY 

Vb  A * h Sgi  Swi  1 Sgi  1  Swi

Pi  PCN 
Bgi  35.35 *   Vb: Volumen bruto; es decir, el volumen total del
Zi * T  PCY 
yacimiento.

Gr43560 * * 1  Swi *Vb * Bga Gr: Gas Remanente

Gp  GOES  Gr Gas Producido

Yacimiento no-volumétrico de gas seco.

Gr 43560 *  * Sgr * Vb * Bga

Gp  GOES  Gr

10
Yacimiento de petróleo sub-saturado (P >Pb) volumétrico.

 BY 
7758  *  * 1  Swi  * Vb
 ac. p 
POES 
Boi
7758 *  * 1  Swi  * Vb
Nr  (P>Pb) So  Swi  1
Bo

7758 *  * 1  Swi  Sg  * Vb
Nr  (P<Pb) Sg  So  Swi  1
Bo

Yacimiento de petróleo sub-saturado (P > Pb) no volumétrico.

7758 *  * 1  Swr  * Vb
Nr  (P>Pb)
Bo

7758 *  * 1  Swr  Sg  * Vb
Nr 
Bo
(P<Pb) Si la presión se estabiliza por debajo de Pb.

Yacimiento de petróleo saturado (P < Pb).

Las ecuaciones generales son:

7758 *  * 1  Swi  Sgi  * Vb

N
Boi

7758 *  * 1  Swi  Sg  * Vb
Nr 
Boi

43560 *  * 1  Swi  Soi  * Vb

G
Bgi

Pero como se dificulta el calculo para la Sg, a través de los registros podemos
conocer el volumen bruto de la zona de gas y el volumen bruto de la zona de

11
petróleo, estos datos nos permiten modificar las formulas de la siguiente
manera:

43560 *  * 1  Swi  * VpZG

G
 PCY 
Bgi  
 PCN 
7758 *  * 1  Swi  Sg  * Vb
Nr 
Boi

Yacimiento de petróleo saturado (P < Pb) no volumétrico.

7758 *  * Sor * Vb
Nr  Sor  1  Swr  Sg
Bo

43560 *  * Sgr * Vb
Gr 
Bg

Cálculos de reservas de cada uno de los yacimientos.

Yacimiento de gas seco.

Grab  43560 * * 1  Swi *Vb * Bgab

Re serva  GOES  Grab

Re serva  Bgi ( PCN / PCY )  Bgab 

FR  *100    *100
Goes  Bgi 

 Bgab ( PCY / PCN )  Bgi 

FR    *100
 Bgi 

12
 BALANCE DE MATERIALES

El método de balance de materiales constituye una de las herramientas

más usadas en la interpretación y análisis de los yacimientos.

Su aplicación permite el cálculo del petróleo originalmente en sitio

(POES), y la identificación y grado de importancia relativo de los mecanismos
de empuje.

A partir de un balance de fluidos, se desarrollara la ecuación general de

balance de materiales así como también un análisis de los parámetros de dicha
ecuación.
La E.B.M relaciona a los fluidos producidos con los fluidos
originalmente en sitio, también considera que los fluidos que se producen en
la superficie son los que se expanden en el yacimiento, producto de la
disminución de la presión y a la compresibilidad tanto de los fluidos como de
la formación.

La ecuación de balance de materiales fue presentada por primera vez

por Schilthuis en 1936.

Este método se basa en la premisa que el volumen poroso de un

yacimiento permanece constante o cambia en una manera predecible con la
presión del yacimiento, cuando el petróleo, agua y/o gas es producido.

Ley De La Conservación De La Masa

Reducción de
Fluidos Expansión de Expansión
volumen
Producidos petróleo + gas de gas libre
poroso del
en solución
= + + hidrocarburo
Fundamentos del método

El E.B.M se fundamenta en el principio de conservación de la energía.

El volumen de control sobre el cual será aplicado este principio es el
yacimiento.

13
 También se fundamenta en que el volumen poroso de un yacimiento
(volumen central), permanece constante o puede ser determinado cada vez que
se produce una reducción de la presión del yacimiento, como consecuencia de
la producción de los fluidos.

En este sentido, un balance de los fluidos del yacimiento podría ser

expresado de la siguiente manera:

 El volumen de los fluidos presentes en el yacimiento en un momento

determinado será igual al volumen de los fluidos iniciales menos el
volumen de los fluidos producidos.

 En este balance el volumen de los fluidos debe calcularse a una

misma condición de presión y temperatura para que tenga validez.

Suposiciones.

 El yacimiento es considerado como un tanque, y por esto es visto

como un modelo de dimensión cero.

 Las propiedades de los fluidos y las rocas se consideran

uniforme.

 Las presiones y las saturaciones se distribuyen en forma

continua.

La ecuación de balance de materiales se utiliza para:

 Determinar los volúmenes de fluidos originales en el

yacimiento.

 Establecer la ocurrencia de una intrusión de agua hacia el

yacimiento.

14
  Predecir el comportamiento de presión-producción.

Usos y Aplicaciones De La Ecuación De Balance De Materiales

Uso:

 El balance de materiales se puede aplicar a ciertas áreas

del yacimiento.
 Permite realizar estudios de agotamiento.
 Establecer los mecanismos de recobro.
 Determinar N y G para confirmar los cálculos volumétricos y así
verificar posibles extensiones del yacimiento.
 Establecer el empuje hidráulico y la tasa de intrusión del agua.
 Verifica la existencia de una zona de gas libre.

Aplicaciones:

 Presión de yacimiento representada por un valor ponderado.

 Petróleo y Gas libre en equilibrio en el yacimiento.
 Datos PVT representativos.

Limitaciones De La Ecuación De Balance De Materiales.

Considera el comportamiento del petróleo y el gas en el laboratorio

similar al del yacimiento.

La EBM no toma en cuenta el factor geométrico del yacimiento, como

tampoco la migración de los fluidos.

15
 El uso de un único PVT introduce errores en los cálculos del Balance
de Materiales.

La presión y temperatura promedio del yacimiento genera una nueva

fuente de error.

En presencia de un acuífero activo (Pi  P) no se puede aplicar el

balance de materiales.

Deducción de la ecuación general

Para la deducción de la ecuación es necesario considerar un yacimiento

con una capa de gas y un acuífero activo. Se debe suponer que después de
cierto tiempo de haber sido puesto en producción, la presión del yacimiento
ha decaído desde la presión inicial (Pi) hasta una presión promedio P.

En este intervalo de presión, el balance volumétrico del yacimiento puede

expresarse de la siguiente manera:

Volúmenes
═ = Expansión
+* + Intrusión
producidos de fluidos de agua

Ilustración:
Gas

Petróleo Contacto gas- petróleo

Agua Agua Contacto agua- petróleo

Volúmenes producidos:

16
 Debido a que el yacimiento que se está considerando posee una capa de
gas, la cual se encuentra en estado de saturación por lo que cualquier
disminución de presión producirá una liberación de gas. Los volúmenes
producidos se pueden encontrar de la siguiente manera:

A condiciones de superficie

Para la caída de presión ∆P =Pi- P

Se obtiene en la superficie:

Volumen acumulado de petróleo, Np (Bls)

Volumen acumulado de gas, Gp (Pies3)
Volumen acumulado de agua, Wp (Bls)

Se obtendrá también un valor de la relación gas- petróleo acumulado Rp,

definido:

A condiciones de yacimiento

Estos volúmenes deben ser puestos en yacimiento mediante factores

volumétricos y la razón gas disuelto- petróleo y evaluados a la presión
promedio P.

Npo= Petróleo + gas en solución.

NpRsg = Gas disuelto producido.
NpRpg = Gas total producido
WpBw = Agua producida

La suma de estos factores representa los volúmenes acumulados de fluidos

que han salido del yacimiento durante la caída de presión P; sumando y
realizando artificios matemáticos, la expresión nos queda:

Np  o + (Rp  Rs) g  Wpw

17
 Expansión de los fluidos y rocas

Como consecuencia de la caída de presión producto de la presión de

los fluidos, el yacimiento y todos los factores que lo componen experimentan
cambios de volumen. Así, habrá cambio de volúmenes en el petróleo y su gas
disuelto, en la capa de gas, en el agua connata y el volumen poroso del
yacimiento. A continuación, se describen los cambios de cada uno de estos
factores:

Petróleo + gas disuelto

 El volumen inicial de petróleo en la zona de petróleo es Nboi. La

reducción de petróleo debido a la caída de presión es:

N(Bo –Boi)

 El volumen de gas inicialmente disuelto en el petróleo es NRsi.

 El volumen de gas disuelto en el petróleo cuando la presión inicial del

yacimiento se reduce hasta la presión promedio P es NRs.

El volumen de gas liberado como consecuencia de la caída de presión

∆P es:
N(Rsi – Rs ).Bg

 El cambio total de volumen en la zona de petróleo es la suma de la

reducción del volumen de petróleo, mas la expansión del volumen del
gas liberado.

N(Bo –Boi) + (Rsi – Rs).Bg.

Capa de gas (m)

Se define como la relación entre los volúmenes originales de gas y

petróleo contenido en la capa de gas y en la zona de petróleo,
respectivamente:

18
  El volumen de petróleo original en la zona de petróleo es NBoi.

 El volumen del gas original en la capa de gas obtenido de la ecuación

anterior es m. NBoi.

 La suma de estos volúmenes será el volumen total original de

hidrocarburo en el yacimiento es:

Nboi (1 + m)= Nboi +mNBoi

Finalmente, para obtener la expansión de la capa de gas, el volumen

original en la capa de gas (mNBio), puede expresarse en condiciones de
superficie como mNBoi/Bgi.

Al producirse la caída de presión P, este volumen se expandirá a

mNBoi/Bgi –Bg. entonces la diferencia entre el volumen anterior menos en
volumen original de gas en la capa de gas (mNBoi) será la expansión de la
capa de gas:

Agua connata

La expansión de la saturación de agua connata (Swc) promedio en la

zona de hidrocarburos pueda determinarse a partir de la definición de
compresibilidad del agua:

dVw =Cw.Vw.P, en donde:

Cw =Compresibilidad de agua, en Lpc-1

Vw = volumen de agua connata en Bls, expresado de la siguiente manera:

.P = Pi – P, caída de presión en lpc.

19
 Si sustituimos el valor de Vw en la primera expresión tenemos:

Volumen poroso

Similar a la expansión de agua connata, el cambio de volumen poroso (Vp) se

puede expresar como:

Donde:

Cf =Compresibilidad de la formación Lpc-1

Intrusión de agua

La expansión de agua del acuífero es comúnmente referida como

intrusión de agua acumulada y se denota We, a condiciones de superficie o
WeBw a condiciones de yacimiento.

La figura muestra en forma esquemática los cambios de volumen que tienen

lugar en el yacimiento al producir una reducción de presión P:

3
Pi P
1
P

4
W Acuífero

1: Expansión de petróleo + gas disuelto.

2: Expansión de la capa de gas.
3: Expansión del agua y reducción del volumen poroso
4: Intrusión de agua.

20
 Ecuación General De Balance De Materiales

Como ya los factores han sido descritos de los mismos permitirá escribir la
ecuación general de balance de materiales:

Aplicación

La ecuación de B.D.M desarrollada anteriormente debe evaluarse

siempre entre la presión inicial del yacimiento y cualquier otra presión P
donde se entregan valores de la producción acumulada de petróleo, gas y agua.
A pesar de que la presión aparece solo explícitamente en los términos de
expansión de agua y reducción del volumen poroso (P = Pi – P), está
implícita en los datos PVT, Bo, Bg, Rg que deben ser calculados a esa presión.

Los parámetros Bo, Bgi y Rsi se obtendrán de los datos PVT a la

presión inicial del yacimiento.

Nota importante:

Si hubiese inyección de gas y/o agua en el yacimiento, los volúmenes

acumulados de estos fluidos a condiciones de yacimiento, deben ser sumados
al lado derecho de la ecuación de balance de materiales.

21
Parámetros de la ecuación general

Estos parámetros contenidos en la ecuación general de balance de

materiales pueden ser separados en dos categorías: aquella que ocupa los
parámetros relativamente conocidos y las que agrupan las potenciales
incógnitas.

Parámetros relativamente conocidos

En la E.B.M podemos encontrar:

Historia de producción

De los parámetros que tiene que ver con la historia de producción, la

producción acumulada de petróleo (Np) es la que presenta el menor grado de
incertidumbre.

Los valores medidos de gas producido acumulado (Gp), presentan un

alto grado de incertidumbre, ya que en yacimientos alejados a los centros de
recopilación, el gas es venteado o usado como generador de electricidad, sin
medirse, o en otros casos su medición es realizada sin mucha precisión.

En relación con el agua producida acumulada (Wp), por razones

económicas su grado es mayor que en el caso anterior.

Datos PVT

La aplicación de la E.B.M requiere que los datos PVT, como son los
factores de los fluidos (Bo, Bg, Bw), la razón gas-petróleo disuelto (Rs) y las
compresibilidades del agua y de la formación (Cw y Cf), hayan sido válidas y
corregidas para ser usadas en las condiciones de operación del yacimiento.

Datos Petrofísicos

En estos datos, la saturación de agua connata promedio del yacimiento

(Swc), obtenida de registros de pozos, referida a cierto dato.

22
Historia de presiones

Esta tendencia es obtenida a partir de la medición de las presiones

individuales de los pozos, referida a cierto dato.
En el de yacimientos con un alto valor de la constante de difusividad
hidráulica (K / D Mc), la presiones de los pozos graficada como función del
tiempo, permitirán definir claramente la tendencia de la historia de presiones.
En yacimientos donde el valor de la constante de difusividad hidráulica
sea baja, la tendencia no podrá definirse claramente. En este caso será
necesario encontrar maneras para promediar las presiones individuales de los
pozos, obtenidas dentro de un período de tiempo, para calcular el valor
representativo del yacimiento, para ese período. Para explicar mejor este punto
observemos la siguiente ilustración:

La figura presenta el gráfico de comportamientos de presión de un

yacimiento basándose en las presiones individuales de los pozos. La figura (a)
presenta el caso de yacimientos con un alto valor de la constante de
difusividad hidráulica y la figura (b), el caso donde dicha constante tiene un
valor bajo:

(b) (a)

Potenciales incógnitas

En la ecuación de balance de materiales las potenciales incógnitas son

POES (N) el tamaño de la capa de gas (m) y la intrusión de agua acumulada
(We).

23
 El POES (N) contabiliza el petróleo que contribuye a la declinación de
presión del yacimiento.

La intrusión de agua se considera cuando hay un contacto agua-petróleo

en un yacimiento. Otro elemento que lleva a pensar en la existencia de una
intrusión de agua al yacimiento es la observación de una suave o inexistente
declinación de la presión del yacimiento y/o el incremento de la producción de
agua en los pozos buzamiento abajo.

Aplicación De La Ecuación De Balance De Materiales De Acuerdo Al

Tipo De Yacimiento

BALANCE DE MATERIALES EN YACIMIENTOS DE GAS

Al igual que en los otros casos la aplicación de la ecuación de balance

de materiales se rige por la ley de la conservación de la materia y se expresa
de la siguiente manera:
Masa Masa de gas
inicialmente remanente
Masa de gas
= de gas en el - en el
producido
yacimiento yacimiento

Para hacer el análisis y deducción de la fórmula de balance de materiales

para el gas deben tenerse en cuenta las siguientes consideraciones:

 Coexisten dos fases en el yacimiento (agua y gas)

 Los efectos de expansión de la formación y agua connota son
despreciables.
 En el yacimiento el petróleo producido es cero.

La ecuación general de balance de materiales se reduce a:

GpBg + WpBw = G (Bg-Bgi) + We (We ‡ 0)

24
BALANCE DE MATERIALES EN YACIMIENTOS SUBSATURADOS
(P > PB)

Suposiciones:

 Hay producción de agua.

 Hay intrusión de agua.
 No hay capa de gas original.
 Compresibilidad cero.

En este caso la Presión de trabajo se encuentra por encima de la

presión de burbujeo. En el yacimiento sólo se tiene petróleo
(líquido).

Si se asumen compresibilidades cero, entonces, el volumen

inicial disponible para hidrocarburos al inicio es igual al final.

25
26
 Cuando hay producción e intrusión de agua y asumiendo las
compresibilidades cero, el balance de materiales puede expresarse de la
siguiente manera:

Suposiciones:

 Hay expansión de la roca + agua connata.

 Hay producción de agua.
 No hay intrusión de agua.
 No hay capa de gas original.

La expansión de la roca y el agua siempre ocurre a medida que la

presión del yacimiento disminuye.

Anexando a la ecuación anterior la expansión de la roca y el agua

connata y partiendo del principio de que los fluidos producidos son iguales a
la expansión de los fluidos del yacimiento (Balance de materiales)

Donde:

Fluidos Producidos = A

Expansión del Petróleo + gas disuelto = B

27
Expansión de la Roca + agua connata = C
Expansión del acuífero = D

Pero:

A = NpBo + WpBw

B = N(Bo – Boi)

C = ∆Vw + ∆Vf

D = We

Determinación de ∆Vw y ∆Vf

Partiendo de la definición de compresibilidad, se tiene:

Determinación de Vp

Si conocemos la saturación inicial de agua (Swi), entonces, se dice que:

28
 Donde:

SHci : Saturación inicial de hidrocarburo.

Vp: Volumen Poroso.

Soi: Saturación inicial del petróleo.

Determinación de Vw

Sustituyendo las ecuaciones obtenidas de ∆Vf y ∆Vw , se tiene que la

expansión de la roca + agua connata, va a estar expresada de la siguiente
manera:

29
sustituyendo A, B, C y D en la ecuación general se tiene que:

También se puede expresar esta ecuación en base a la compresibilidad de

petróleo (Co) o compresibilidad efectiva (Ce), para realizar esta deducción hay
que considerar que:

30
Pero como 1 – Swi = So, entonces la ecuación quedaría definida:

SATURACIÓN DE GAS LIBRE EN EL YACIMIENTO

La cantidad de gas en solución del tiempo t a condiciones normales es

N, Rsi (gas original en el yacimiento). Esta cantidad de gas debe ser en
condiciones normales, al tiempo t igual a la cantidad de gas en el yacimiento
mas el Gas producido a ese tiempo t (el gas producido incluye todo el gas
medido en el separador y en el tanque).

El gas al tiempo t remanente en el yacimiento es el que se encuentra aun

en solución mas el gas libre que queda en el yacimiento, bien sea formando
una capa de gas secundaria, o bien difundido dentro de la masa del petróleo.
Efectuando un balance de gas a condiciones de yacimiento, se obtiene lo
siguiente:

Gas Gas original Gas original Gas Aun

liberado en =
disuelto en -
disuelto en -
en la
el el el superficie.
yacimiento yacimiento yacimiento

GL= [NRsi - (N-Np)Rs – RpNp]Bg

31
 Expresando el gas libre en el yacimiento G, como la saturación del gas,
Sg, se tiene considerando:

Sgi GOES capa Sg

Producción
Capa de Gas GRem Libre

Soi So
N GOES Sg=0 NRem
Swi oil solución Sw

Pi , Rsi (PCN/BN) Pab , Rsab

Sg= Vg = GRem solución a CY (BY)

Vp V Poroso a CY (BY)

GRemtotal = GRemlibre + GRemsolución

GOEStotal = GOEScapa + GOESsolución … También viene dado por:
GOEStotal = Gptotal + GRemtotal
Rp = Gp Gp = RpNp
Np
Sustituyendo en la Ecuación se tiene:
GOEScapa + GOESsolución = RpNp + GRemtotal
GOEScapa + GOESsolución = RpNp + GRemlibre + GRemsolución
GRemlibre = GOEScapa + GOESsolución – RpNp – Gremsolución
m = GOEScapa(PCN)Bgi (BY/PCN) Despejando GOES se tiene:
N (BN) Boi (BY/BN)
GOEScapa (PCN)= mNBoi
Bgi

32
GRemlibre(BY) = [ (mNBoi/Bgi) + NRsi – RpNp - NRemRsab]Bg
(N – Np)
Vporoso = Vt*Φ
N = 7758(BY/A-P)Vt(A-P) Φ(1-Swi)/Boi(BY/BN)

Vporoso(BY) Deapejando Vporoso se tiene:

Vporoso = NBoi/(1-Swi) ; Sustituyendo en Sg = Vg/Vp se tiene:
Sg = [ (mNBoi/Bgi) + NRsi – (N – Np)Rs - RpNp ] Bg
NBoi/ (1 – Swi)
Para yacimientos subsaturados m = 0
Calculando Sg de otra manera se tiene:
So = Vo/Vp = NRem
(NBoi/(1 – Swi) …. En By
So = (N-Np)(BN)Bo(BY/BN) So = (N-Np)Boi(1 – Swi)
Nboi/ (1 – Swi) NBoi
Sg + So + Sw = 1 Sg = 1 – So – Sw
Sg = [ 1 – (N-Np)Bo(1 – Swi) - Swi
NBoi

Esta ecuación demuestra la relación directa que hay entre la saturación

del gas (Sg) y el recobro del petróleo (Np/N), es decir a mayor saturación de
gas libre en el yacimiento, mayor es el recobro del petróleo, razón por la cual
se debe producir con la menor relación gas-petróleo que se pueda
operacionalmente, y menor aún, si las condiciones del yacimiento son muy
atractivas, iniciar a una temprana edad un proceso de inyección de agua y/o
gas según sea el caso.

BALANCE DE MATERIALES EN YACIMIENTOS SATURADOS

(P < Pb)

33
 En este caso la presión ha disminuido por debajo de la presión de
burbujeo y el gas disuelto en el petróleo comienza a liberarse.
Asumiendo las compresibilidades cero:

Como:

VHCi = VHCf, entonces:

Nboi = (N – Np)Bo + GI.Bg

Determinación de GI.

El gas libre proviene del petróleo, entonces si al gas originalmente

disuelto en el petróleo se le resta el gas disuelto a las nuevas condiciones y el
gas producido, se obtiene el gas libre en el yacimiento.

34
 En otro caso, el yacimiento se encuentra inicialmente saturado, con una
capa de gas libre, que debe ser considerada en la ecuación de balance de
materiales.

Donde:

Fluidos Producidos: A

Expansión del petróleo: B

Expansión del gas libre: C

Expansión de la Roca + agua connata: D

Expansión del acuífero: E

Entonces:

35
Gp: gas producido.

Se define “m” como el tamaño de la capa de gas, y se tiene que:

despejando G de “m” se establece que:

sustituyendo G en C:

36
Utilizando la fórmula de compresibilidad, despejando ∆v y adaptándola a los
parámetros Cf y Cw, se tiene:

Sacando factor común ∆p (Nboi + GBgi)

Finalmente la ecuación de balance de materiales nos queda:

37
 FACTOR DE RECOBRO DEL PETRÓLEO

Factor de recobro de petróleo (FRP): De todo el volumen de

hidrocarburos presente en un yacimiento (POES), solo una fracción puede
ser conducido hasta la superficie, debido a la existencia de las fuerzas
retentivas de fluidos en el yacimiento. Además, solo una fracción de este
último puede hacerse en condiciones de rentabilidad económica. El factor
que agrupa estas dos fracciones se denomina Factor de Recobro de Petróleo
(FRP) que es característico de cada yacimiento.
El factor de recobro es función de las propiedades de las rocas y los
fluidos del yacimiento, su nivel de presión y del mecanismo de producción
predominante en el yacimiento, y se obtiene por medio de correlaciones
establecidas provenientes de yacimientos con características similares.
El producto de este factor por el POES de un yacimiento constituye las
reservas recuperables (Npr), lo cual en términos matemáticos puede
representarse de la siguiente manera:

Partiendo de la ecuación de Balance De Materiales se tiene:

Np[Bo + (Rp-Rs)Bg] + WpBn = N[(Bo-Boi) + (Rsi-Rs)Bg + mNBoi((Bg/Bgi) – 1) +

(NBoi + GBgi) (CwSwi+Cf/(1-Swi))∆p+We

Yacimiento Subsaturado m = 0
Yacimiento Volumétrico (We = 0)

38
Se desprecia Cw, Cf Quedando la ecuación como:
Np[Bo + (Rp-Rs)Bg] = N[(Bo-Boi) + (Rsi-Rs)Bg]
Por lo tanto el factor de recobro del petróleo queda expresado de la siguiente
manera:

FRP = N/Np =(Bo – Boi) + (Rsi – Rs) Bg/ (Bo + (Rp – Rs)Bg)

Esta relación demuestra que hay una relación inversa entre el recobro del
petróleo (FRP) y la relación gas petróleo acumulada (Rp). De aquí se concluye
que para obtener un alto recobro primario, tanto gas como sea posible debería
permanecer en el yacimiento, lo cual requiere que la razón gas-petróleo debe
mantenerse tan baja como sea posible.

Técnica de Havlena y Odeh (Usos)

La aplicación de la ecuación general de balance de materiales se

simplifica mediante la técnica de Havlena y Odeh.

Está técnica se basa en ver la mencionada ecuación como la ecuación de

una línea recta, donde después de definir esta tendencia, el cálculo de la
pendiente y el intercepto, permite obtener algunos parámetros hasta
desconocidos (incógnitas).

Ecuaciones

A partir de la ecuación general de balance de materiales, utilizando la

técnica de Havlena y Odeh, se expresa de la siguiente manera:

39
40
Aplicaciones de la Técnica de Havlena y Odeh

 Yacimientos Subsaturados

 No hay capa de gas inicial.

 Pi >Pb

Caso 1: sin intrusión de agua (We = 0)

Todos los datos de F, Eo y Efw se pueden obtener de una historia de

producción. Si se toman estos como variables y se construye una gráfica, se
obtendrá una línea recta y se podrá calcular N.

41
 Yacimientos Subsaturado sin intrusión de agua

Caso 2: Con intrusión de agua (We ╪ 0) We es conocida

Se divide entre (Eo + Boi.Efw):

b=N

42
 Ecuación de una recta que parte del punto b y tiene una pendiente igual a 1.

Yacimiento Subsaturado con intrusión de agua

 Yacimientos Saturados (Pi < Pb)

Para yacimientos Saturados Efw es despreciable en comparación a la

expansión del petróleo y el gas, es por ello que se hace cero.

43
Caso 1

 Sin intrusión de agua (We = 0)

 Sin capa de gas inicial (m = 0  G = 0)
 Pi = Pb

Yacimiento Saturado sin intrusión de agua y sin capa inicial de gas.

Caso 2

Sin intrusión de agua

Con capa inicial de gas (m ╪ 0 → G ╪ 0)
Pi = Psat

Si se conoce “m” (se desconoce G)

44
45
Yacimiento saturado sin intrusión de agua y con capa inicial de gas (conocida)

Si “m” es desconocida (G se conoce)

Se divide toda la ecuación entre Eo:

n=G

46
 Yacimiento saturado sin intrusión de agua y con capa inicial de gas
(desconocida)

Procedimiento:

 Asumir un valor de “m”

 Calcular F/Eo y Eg/Eo

 Graficar F/Eo Vs Eg/Eo

 Ajustar “m” hasta que se obtenga una línea recta

 Calcular el valor de la pendiente (N).

Caso 3

Con intrusión de agua (We ╪ 0, conocida)

 Con capa de gas inicial (m ╪ 0)
 Pi = Psat

47
Dividiendo la ecuación entre Eo:

Yacimiento saturado con intrusión de agua y con capa de gas inicial

CURVAS DE DECLINACIÓN DE PRODUCCIÓN

48
 La curva de declinación de producción representa un método dinámico
para la estimación de las reservas recuperables de un yacimiento.

Su característica dinámica proviene del hecho de que utiliza la historia

de producción de los fluidos, concretamente de petróleo, por pozos o por
yacimientos, para la estimación de sus reservas recuperables.

La aplicación de este método parte de que existe suficiente historia de

producción como para establecer una tendencia de comportamiento, por lo que
la predicción del yacimiento puede hacerse a partir de la extrapolación de
dicha tendencia.

En general se busca un tipo de gráfico donde la tendencia se presenta de

forma lineal para facilitar su extrapolación. Este procedimiento lleva implícito
una suposición básica. "Todos los factores que han afectado al yacimiento
en el pasado lo seguirán afectando en el futuro".

En Venezuela este método es ideal porque:

1.- Los Yacimientos son grandes y tienen extensos periodos de vida extensa.

2.- No existen restricciones severas de las tasas de producción por parte del
M.E.M.

El estudio de la declinación de producción de un pozo o yacimiento en

particular puede hacerse según dos procedimientos:

49
 1. Método grafico: El método grafico consiste en lo siguiente:

a) La vida del yacimiento se representa gráficamente en distintos tipos de

papel (normal, semilog, etc.), con el objeto de obtener la curva mas
sencilla a usar, en la mayoría de los casos la línea recta, que permita
hacer extrapolaciones y así poder efectuar predicciones.

b) Una vez determinada la representación grafica as conveniente, efectuar

las correspondientes extrapolaciones (preedición) hasta las condiciones
de abandono.

Entre las variables a considerar en este método es:

Tasa limite económico (T.L.E.): es la base para llevar a cabo el

abandono de un yacimiento, siendo un factor muy variable y particular de cada
yacimiento, su expresión matemática viene dada por:

T.L.E: = C (Bs / Mes*Pozo) / P (Bs / BN ) * 30.4 (días / mes)

Donde:

C = Costo de Producción por pozo en 1 mes

P = Entrada de dinero por BN de petróleo producido antes de ISLR. En

Venezuela es el precio de venta.

Factores que afectan las curvas de declinación: entre los muchos factores
que afectan las curvas de declinación se tienen los siguientes:

50
  Períodos desiguales de tiempo en las medidas.

Algunas mediciones hechas en los pozos, a intervalos desiguales de

tiempo, dificultan la definición de la tendencia de su comportamiento. Habrá
entonces, periodos de tiempo donde se tendrá que interpolar para completar la
tendencia, con los riesgos que esto implica.

 Cambio en la producción de los pozos.

Cuando el cambio de productividad de un pozo es significativo o más

de un pozo del yacimiento experimenta cambios similares, ya sea por la
aplicación de métodos de recuperación secundaria o por los mecanismos de
empuje propios del yacimiento, también se afectara la tendencia del
comportamiento de producción del yacimiento.

 Completación de nuevos pozos.

Cada vez que se complete un nuevo pozo, la producción del yacimiento

se incrementara, lo cual altera la tendencia del comportamiento anterior a
dicha completación. En este caso, habrá que esperar hasta que se observe una
nueva tendencia para proceder a la aplicación de este método.

Sin embargo, se necesitara una tendencia, podría tratarse una paralela a

la tendencia anterior por el nuevo valor de la tasa de producción del
yacimiento.

 Interrupción de los programas de producción.

Esporádicamente ocurren cierres de producción en los pozos, estos

cierres pueden ocurrir debido a carencia de mercados, caída del fluido

51
eléctrico, etc. Una vez que se produce la apertura del yacimiento se obtendrán
muy probablemente nuevos niveles de producción lo que limita la aplicación
del método.

 Veracidad de información disponible.

Este factor tendrá un efecto importante cuando no se conozca la

metodología usada para asignar los valores de producción a los pozos y por
ende, al yacimiento en estudio.

 Prorrateo.

En países donde existe restricción en las tasas de producción, los

yacimientos no producen a su potencial y, por tanto, el método no podrá
aplicarse.

2. Método matemático: por medio de este método se trata de

encontrar una expresión matemática en base a la información disponible hasta
la fecha y luego utilizar dicha expresión para predecir (reservas existentes,
tiempo de abandono, etc.). Realmente es el mismo método grafico pero
obteniendo la mejor línea recta (u otra curva) que pase a través de los puntos
datos. Para hallar esta mejor línea recta se utilizan diferentes métodos
estadísticos entre los que tenemos:

52
 a. Método de la pendiente – intersección: consiste en trazar una
recta cualquiera a través de los puntos graficados y entonces determinar su
pendiente y su intersección con el eje Y. Estos valores definen la ecuación
de esa recta.

b. Método de los puntos notables: de la información disponible ya

graficada se eliminan aquellos puntos que posiblemente y en relación con la
mayoría presenten error (muy distante de la posible solución) y con los
puntos se estima no tengan error, se determina la ecuación de la línea recta.

c. Métodos de los promedios: considera a los puntos que están

sobre la línea recta como una serie de puntos y a los que están abajo como
otra serie de puntos.

d. Métodos de los mínimos cuadrados: considera a todo los puntos

involucrados en la declinación de producción y con la ayu8da de
procedimientos estadísticos (aplicación de mínimos cuadrados) se determina
la ecuación de la línea recta.

53
 Existen tres tipos de curvas de declinación de producción: La Curva de
Declinación Exponencial, La Curva de Declinación Hiperbólica y la Curva de
Declinación Armónica.

Tipos De Curva De Declinación De Producción

1) Curvas de Declinación Exponencial

Se dice que ocurre una declinación exponencial de las tasas de producción

cuando la variación de la tasa de producción con el tiempo expresada como
una fracción de ella misma es una constante.

Matemáticamente, este tipo de declinación se expresa de la siguiente

manera:

 dq 
D 
dt 
q

Donde:
q = tasa de producción, Bls/ Días

t = tiempo de producción (días).

D = constante de declinación exponencial (días).

Dada la historia de producción de un pozo o yacimiento, el

reconocimiento de declinación del tipo exponencial puede hacerse mediante su
definición.

54
 Este cálculo se presenta en la tabla siguiente, lo que permite obtener
como resultado un conjunto de valores de la constante de declinación D.

t q dq Dt dq/dt q D
día Bls
dia
Bls
dia Día Bls
dia
Bls
dia día-1
dia

t0 q0
t1 q1 q1  q0 t1  t 0 q1  q0 q1  q0 D1
t1  t 0 2

t2 qz q2  q1 t 2  t1 q2  q1 q 2  q1 D2
t 2  t1 2

t3 q3 q3  q2 t3  t 2 q3  q 2 q3  q 2 D3
t3  t 2 2

tn qn qn  qn1 t n  t n 1 qn  qn1 q n  q n 1 Dn
t n  t n1 2

Si estos valores se representan alrededor de un cierto valor (D),

estaremos en presencia de una declinación exponencial y podría determinar el
valor promedio de la constate D de la siguiente forma:
n

D i
D i 1
n
Si los valores no se representan alrededor de un valor, tendríamos que
probar otro tipo de declinación.

Desarrollo de la ecuación de la Curva de Declinación Exponencial

Al desarrollar la ecuación que define la declinación exponencial,

podremos encontrar expresiones para la tasa de producción y la producción
acumulada de petróleo.

Tasas de Producción

55
 La ecuación de la tasa de producción bajo este tipo de declinación se
obtendrá a partir de su definición:

 dq 
 dq 
D 
dt 
  Ddt  
q  q

Esta ecuación se integra desde los valores iniciales (t i, qi), hasta los
últimos valores (t, q) de la producción:

q dq t
qi q
  D   dt
ti

Al resolver la integral tenemos que:

ln  q   D  t  ti 
 qi 

Al despejar q, la ecuación queda:

q  qi  e  Dt ti  Asumiendo el ti = 0

q  qi  e  Dt

En forma logarítmica se obtiene la recta:

Log q    D  t  Log qi  Donde: m = - D/ 2,303

2,303

Con la cual se puede graficar y predecir el comportamiento de la tasa de

producción en función del tiempo para el pozo o yacimiento en estudio.

56
Procedimiento para obtener la grafica de la tasa de producción en función del
tiempo

Pasos Acciones
1 Identificar si la declinación es de tipo exponencial utilizando la
tabla de historia de producción del pozo o yacimiento.
2 Graficar la tasa de producción (Bls/días) del pozo o yacimiento
en función del tiempo (días) en papel semi - logarítmico.
3 Observar el comportamiento de D:
- Si el lineal continuar con el siguiente paso.
- Si no es lineal, la curva no es de tipo exponencial y se
cancela este procedimiento.
4 Se calcula la pendiente de la recta trazada, de la recta
logarítmica se sabe que la m = - D / 2,303 de la cual se despeja
y se calcula D.
5 Se extrapola la línea recta hasta el tiempo cuando se quiera
predecir la producción de la tasa de petróleo del pozo o
yacimiento, o hasta la tasa de producción limite económico o de
abandono para predecir el tiempo cuando se alcanzara esta tasa.

Función de la ecuación de la tasa de producción

 Identificar el tipo declinación exponencial.

 Predecir el comportamiento de la tasa de producción de petróleo del

pozo o yacimiento.

Producción acumulada de petróleo (Np)

Por definición la producción acumulada de petróleo en forma

diferencial es:

57
dNp  q  dt

Al despejar q la ecuación queda:

q  dNp
dt

Sustituyendo la expresión de la tasa de producción para la declinación

exponencial, y despejando dNp:

qi  e  Dt  dNp dNp = qi e-Dt dt

dt

Integrando esta ecuación desde los valores iniciales del tiempo o de un

tiempo dado hasta los valores finales de la producción acumulada y el tiempo:

Np t

 dNp  q  e
 Dt
i  dt
u u

Al resolver esta integral se obtiene:

 e  Dt 1  qi  e  Dt qi
Np  qi     ó Np  
 D D D D

Se sustituye la ecuación de la producción acumulada de petróleo,

obteniéndose la recta:

q qi
Np    Donde m = - 1/D
D D

58
 Con lo cual se puede graficar y predecir el comportamiento de la
producción acumulada de petróleo en función de la tasa de producción del
pozo o yacimiento.

Procedimiento para obtener la grafica de la producción acumulada de

petróleo en función de la tasa de producción

Pasos Acciones
1 Identificar si la declinación es de tipo exponencial utilizando la
tabla de la historia de producción del pozo o yacimiento.
2 Graficar la producción acumulada de petróleo (Bls) en función de
la tasa de producción de petróleo (Bls/días) del pozo o yacimiento
en papel normal.
3 Observar el comportamiento de D:
- Si es lineal continuar con el siguiente paso.
- Si no es lineal, la curva no es de tipo exponencial y se
cancela este procedimiento.
4 Se calcula la pendiente de la recta trazada, de la ecuación de la
producción acumulada de petróleo se sabe que la m = - 1/D de la
cual se despeja y se calcula D.
5 Se extrapola la línea recta hasta la tasa de producción del pozo o
yacimiento o hasta la tasa de petróleo límite económico o de
abandono para predecir el volumen de reservas que podrían
recuperarse cuando se alcance dicha tasa.

2) Curva de declinación hiperbólica

Es aquella donde la variación del inverso de la constante de declinación

con el tiempo es una constante.

La definición matemática de la declinación hiperbólica es la siguiente:

59
 q
a 1 
D dq
dt
Donde:

D = Tasa de declinación (días -1)

q= Tasa de producción (Bls/días)
t = Tiempo de producción (días)
a = Inverso de declinación (días)

Dada la historia de producción de un pozo o yacimiento, el tipo de

declinación hiperbólica se puede reconocer usando su definición. La tabla
siguiente muestra la secuencia de cálculos para obtener los valores de "a".

t q dq dt dq/dt q a
q a b  a
dq t
dt
Bls Bls Bls Bls
día dia dia día dia dia día
dia

t0 q0
t1 q1 q1  q0 t1  t 0 q1  q0 q1  q0 a1 ∆a1 b1
t1  t 0 2

t2 qz q2  q1 t 2  t1 q2  q1 q 2  q1 a2 ∆a2 b2
t 2  t1 2

t3 q3 q3  q2 t3  t 2 q3  q 2 q3  q 2 a3 ∆a3 b3
t3  t 2 2

tn qn qn  qn1 t n  t n 1 qn  qn1 q n  q n 1 an ∆an bn

t n  t n1 2

Si los valores de a se encuentran ordenados en forma creciente diríamos

que la declinación es de tipo hiperbólica y procedemos a calcular los valores
de b. De lo contrario tendríamos que probar otro tipo de declinación.

60
 La variación de los valores de b con el tiempo, se representarán
alrededor de un cierto valor cuyo promedio puede calcularse como:
n

b i
b i 1
n 1

Desarrollo de la ecuación de la curva de declinación hiperbólica

Al desarrollar la ecuación de la declinación hiperbólica podremos

encontrar expresiones para la tasa de producción y la producción acumulada
de petróleo.

Tasa de producción

La ecuación de la tasa de producción bajo este tipo de declinación se

obtendrá derivando la declinación hiperbólica con respecto al tiempo:

 
 q 
d
 dq 
 
  dt   b
da
dt dt

Si separamos las variables e integramos:

 
 q 
 d
 dq  
 b  dt
 dt 

Se obtiene:
q
 bt  c
dq
dt
Donde "c" es la constante de integración que al evaluarla para t = 0, será
igual a a1.
El intervalo de a1 es la constante de declinación inicial de Di.
Entonces, la ecuación anterior quedaría así:

61
 q  Di  bt  1

dq Di
dt

Se separan variables y se integra desde los valores iniciales (t 0 = 0, q0)

hasta los últimos valores (t, q) de la producción:

D
t q
dq
0 Di  bt i 1  dt  q q
i

Resolviendo la integral se obtiene:

Ln (Di b t + 1)-1/b = Ln (q/qi)

Simplicando:

q = qi (Di b t + 1)-1/b

q-b = qi-b (Di b t + 1)

Con la cual se puede predecir el comportamiento de la tasa de

producción en función del tiempo para el pozo o yacimiento en estudio.

Dado que la curva de declinación hiperbólica no se puede ajustar a una

línea recta no se hace ningún procedimiento para graficar la tasa de
producción en función del tiempo ya que en la curva es muy difícil hacer una
extrapolación.

Función de la ecuación de la tasa de producción

62
 La ecuación de la tasa de producción de petróleo es útil para predecir
cual será el nivel productivo de un pozo o yacimiento a través del tiempo.

La predicción de la tasa de producción de petróleo se realiza al sustituir

el tiempo al cual se quiere estimar dicha tasa en la ecuación anterior.

Con la ecuación de la tasa de producción y conociendo la tasa de

abandono se puede despejar el tiempo de abandono:

q-b = qi-b (Di b t + 1) donde q será qa

((qa/qi)-b - 1) / Di b = ta

Producción acumulada de petróleo

La ecuación de producción acumulada de petróleo para la declinación

hiperbólica puede deducirse al sustituir la ecuación de la tasa de producción
de petróleo desarrollada anteriormente, en la producción acumulada de
petróleo.

dNp  q  dt

Se integra desde un tiempo inicial hasta un tiempo final:

t
pero q  Di  bt  1
1
Np   q  dt , b

Por lo tanto:
t
Np  qi  Di  bt  1
1
b  dt
0

63
 Resolviendo esta integral se tiene:

b 1
 Di  bt  1 b  1
qi
Np 
Di b  1  

Pero como:
b
Di  bt  1   q 
 qi 

La expresión anterior se reduce a:

qi  q 1b 
Np      1
Di b  1  qi  
Dado que la curva de declinación hiperbólica no se puede ajustar a una
línea recta no se hace ningún procedimiento para graficar la producción
acumulada de petróleo en función de la tasa de producción de petróleo ya que
en la curva es muy difícil hacer una extrapolación.

3) Curva de declinación armónica

La declinación armónica es un caso particular de la declinación

hiperbólica, la cual ocurre cuando b=1.

La definición matemática de este tipo de declinación es la misma que la

declinación hiperbólica.

La declinación armónica puede ser reconocida al construir una tabla

idéntica a la que se escribió cuando se describió declinación hiperbólica.

Tasa de producción para una curva de declinación Armónica

64
 La ecuación de la taza de producción de petróleo para esta declinación se
obtiene a partir de la ecuación para la declinación hiperbólica, al asignarle al
parámetro “b” el valor de uno (1). Así tendremos que:

qi
q
Di  t  1

La tasa de producción en función del tiempo debe graficarse en papel Log

– Log para poder ajustar los resultados a una línea recta para poder hacer las
estimaciones necesarias.

Función de la ecuación de la tasa de producción

Predecir la taza de producción de petróleo a un tiempo determinado.

Una vez definida la taza de producción de petróleo al abandono (q a), la

ecuación anterior permite calcular el tiempo en que se calcularía dicha tasa.

1  qi 
ta    1
Di  qa 

Producción acumulada de petróleo

La ecuación de la producción acumulada de petróleo se obtendrá al

sustituir la ecuación de la taza de producción de petróleo y se tiene:

dNp= q.dt = qi. (DI.t +1)-1 dt

Al integrar ambos miembros:

t
ln Di  t  1
1 qi
Np  qi  dt 
0
Di  t  1 Di

65
Pero el factor Ln (Di t + 1 ) es igual a Ln (qi / q)

Entonces:

qi  qi 
Np  ln  
Di  q 

Luego la ecuación lineal para la producción acumulada de petróleo será:

ln qi   i ln q  log qi   23,03 i log q 

qi q qi q
Np  ó Np  23,03
Di Di Di Di

La producción acumulada de petróleo debe graficarse en papel semi –

Log para poder ajustar los resultados a una línea recta para poder hacer las
estimaciones necesarias.

SIMULACIÓN DE YACIMIENTOS

La simulación: es la representación del funcionamiento de un

determinado proceso por medio de la computadora

Simulación de yacimiento

La simulación de yacimientos ha sido practicada desde el comienzo de

la ingeniería de petróleo en la década de los 30’s. Antes de 1960, estos
cálculos consistieron básicamente en métodos analíticos, balance de materiales
0-dimensiónales y cálculos de modelos unidimensionales del tipo Buckey y

66
Leverett. Como tal el término simulación se hizo común en los años 60’s,
como un método predictivo involucrado con sofisticados programas de
computación. Estos programas presentan amplias ventajas ya que permiten la
solución de un conjunto de ecuaciones diferenciales que describen flujo en
varias dimensiones (2-D y 3-D), transiente y multifásico en un medio poroso
heterogéneo. Esta ventaja fue hecha posible por la rápida evolución a gran
escala, alta velocidad y desarrollo de métodos numéricos para resolver grandes
sistemas de ecuaciones diferenciales.

La simulación de yacimiento es el proceso de inferir el comportamiento

real a partir del comportamiento de un modelo el cual puede ser físico,
conceptual o numérico, y consiste en describir el yacimiento y predecir el
desempeño futuro del yacimiento bajo una variedad de escenarios. Los
simuladores son ampliamente usados en el análisis del yacimiento y el estudio
de métodos para el mejoramiento de la recuperación de hidrocarburos. La
observación del comportamiento del modelo en diferentes condiciones de
producción permitirá seleccionar el esquema óptimo de explotación del
yacimiento.

Propósito de la Simulación de un Yacimiento

 Es utilizada para estimar la recuperación de un esquema de

producción existente para un futuro.
 Evalúa los efectos sobre la recuperación de condiciones de
operador activa.

 Comparar la economía de diferentes métodos de recuperación.

 La simulación numérica ha sido utilizada para estimar el
rendimiento de inyección térmica en un ambiente de yacimiento donde los
procesos son muy complejos y muchos parámetros del yacimiento afectan
los resultados.

Aplicaciones más comunes del Simulador de yacimiento.

 Determinar el comportamiento de un yacimiento bajo un proceso

de intención particular o agotamiento natural.
 Evalúa las ventajas de un proceso de inyección de agua de flanco
contra un proceso de agua por arreglos.

67
  Determina el efecto de ubicación de los pozos y el
espaciamiento.
 Investiga el efecto sobre el recobro de variaciones en las tasas de
inyección y/o producción. Investiga el efecto sobre el recobro de la
perforación interespaciada.

Tipos de Simuladores de Yacimientos.

Existen varios tipos de simuladores de yacimiento. Estos son

generalmente clasificados en Simulador de petróleo negro, composicional,
térmico y químico, dependiendo de las características de flujo de fluidos y el
comportamiento de la transferencia de masa y calor.

A.- Simulador de Petróleo Negro. Son frecuentemente utilizados para

simular procesos isotérmicos, flujo simultáneo de petróleo, gas y agua debido
a fuerzas gravitacionales, viscosas y capilares. El término petróleo negro es
usado cuando se señala que la fase de hidrocarburo se considera como un
líquido único y gas, y no hay cambio en la composición química, por lo tanto
la composición de las fases permanece constante aunque la solubilidad del gas
en petróleo y agua se toman en cuenta.

B.-Simulador Composicional:
Estos simuladores consideran la variación de la composición de las
fases con la presión, estos son utilizados para desarrollar estudios en
yacimientos de gas condensado y petróleo volátil.

C.-Simuladores Térmicos:
Este tipo de simulador toma en cuenta tanto el flujo de fluidos como la
transferencia de calor y reacciones químicas. Los usos prácticos de estos
simuladores están ubicados dentro de los procesos de simulación de inyección
de vapor de agua, agua caliente y procesos de combustión en sitio.

D.-Simuladores Químicos:
Estos simuladores consideran el flujo de fluidos, el transporte de masa,
debido a dispersión, adsorción, filtración, cinética de reacción y cambios del
comportamiento de fases. Son usados en procesos de inyección de sulctantes,
polímeros, emulsiones, sistemas gelificantes y flujos de compuestos alcalinos.

68
 Cada programa de simulación van a depender de la empresa que lo
realice. Cabe destacar:
Sclumberger: a través de geoquest ofrece su suit de :

 Eclipse 100 – Petróleo negro.

 Eclipse 300 – Composicional.
 Eclipse 500 – Químico y Térmico.

CMG: (Computer modellin group) tiene:

 IMEX – Petróleo nogro.

 GEM – Composicional.
 STARS – Químico y Térmico.

METODO DE SIMULACION NUMERICA.

Consiste en simular numéricamente el yacimiento y esta basado en la

disgregación del yacimiento en un número de bloques, lo cual permite
considerar sus heterogeneidades y predecir su comportamiento. La validez
de este método requiere de una buena definición geológica del yacimiento y
de las características petrofisicas y de los fluidos

Modelo del Yacimiento.

El modelo integrado de yacimiento es aquel que caracteriza el
yacimiento incorporando de manera consistente aspectos sísmicos, geológicos
propiedades de rocas y fluidos, historia de producción y/o inyección e historia
de presiones.
En un modelo de simulación se debe representar las propiedades de roca
y fluido, las características de los pozos y su relación con el yacimiento. Debe
tener la descripción geológica-estructural del mismo, en general debe tener
reflejado todas las características del yacimiento que se va a estudiar.

Tipos de Modelos según la Geometría.

A. Modelo de Cero Dimensiones.

69
 El modelo más simple que se puede construir es el cero dimensional o
una celda, mejor conocido como balance de materiales. El balance de
materiales se usa normalmente para estimar fluidos inicialmente en sitio o la
presión del yacimiento. Este tipo de modelo asume que todas las condiciones
del yacimiento existen para una presión promedio y que el sistema completo
está en condiciones de equilibrio. Este modelo es muy útil al comienzo del
estudio para realizar revisiones rápidas de consistencias de datos como PVT,
restricciones de pozos, etc.

Figura n° 2.11; Modelo de Simulación cero dimensional

B. Modelos de Una Dimensión.

La orientación de los bloques puede ser horizontal, vertical o con cierto

ángulo de inclinación. Estos modelos dan una buena representación del
movimiento de fluidos globalmente, así como distribución promedio de
presiones. Los efectos de caída de presión del pozo sobre el comportamiento
global del yacimiento no pueden ser, generalmente, simulados con éstos
modelos, debido a que la menor unidad del yacimiento (un bloque) es muy
grande comparado con el volumen del yacimiento que está afectado por la
presión en el pozo. Sin embargo, estos efectos pueden ser simulados por un
modelo 1-D radial. Los modelos 1-D son útiles cuando el espesor del
yacimiento, h, es pequeño en comparación con su longitud; el petróleo se
drena por un sistema de pozos casi equidistantemente espaciados; es decir,
formando filas paralelas al contacto agua - petróleo; los efectos de
conificación se desprecian. En muchos otros casos, modelos 1-D son
representaciones pobres del yacimiento.

70
Figura n° 2.12; Modelos de una dimensión (1-D).

C. Modelos de Dos Dimensiones.

Para modelar la eficiencia de barrido de un fluido desplazante es necesario

utilizar modelos 2-D. Este puede ser un modelo radial, un modelo transversal
para simular conificación y segregación gravitacional, o un modelo área para
simular efectos de barrido. El modelo 2-D radial es útil para determinar la tasa
crítica de producción a la cual ocurrirá conificación, para predecir el
comportamiento futuro de un pozo conificado y para evaluar los efectos de
barreras de lutitas o permeabilidad vertical baja. Otro uso de los modelos 2-D
radial es en el análisis de pruebas de presiones. Probablemente, el uso más
extensivo de los modelos 2-D areales es para determinar los patrones óptimos
de inyección de agua o gas. También son útiles para determinar la posición de
pozos.

Figura n° 2.13; Modelos de dos

dimensiones (2-D)

D. Modelo de Tres Dimensiones.

Estos modelos pueden tomar en cuenta casi todas las fuerzas presentes
en el yacimiento. Considera, no solamente los efectos de barrido areal, sino
también los efectos gravitacionales. Sin embargo, pueden ser muy difíciles
para modelar fenómenos locales (tales como conificación) donde se requieren
bloques muy pequeños para una representación adecuada. Los modelos 3-D
radiales son una generalización del 2-D radial, en el cual se pueden tomar en

71
cuanta penetración parcial del pozo en la arena productora y cañoneo parcial,
así como cualquier otro parámetro que dependa de la profundidad. En el uso
de simuladores sofisticados se deberá siempre pensar cuidadosamente los pro
y contra de cada tipo de modelo. Usando 2-D se puede ahorrar tiempo pero se
pueden obtener resultados irreales debido a que la situación es mucho más
compleja para ser representada por una aproximación simplificada. Todo
depende de los datos disponibles, de la complejidad del yacimiento, del patrón
de pozos, de la distribución de producción entre pozos y otros elementos como
las completaciones de los mismos.

Figura 2.14 Modelos de tres dimensiones

Información requerida en Simulación.

 Datos del Yacimiento: geometría, tamaño de los bloques,

porosidad permeabilidad, profundidad de cada bloque, espesor de
la arena, curva de permeabilidad relativa y presiones capilares.

 Propiedades PVT: factor volumétrico, relación gas-petróleo y

viscosidad, ubicación de pozos, completación índice de
productividad e historia de producción.

Simulación matemática:

 Se basa en principios de balance de materiales.

72
  Toma en cuenta la heterogeneidad de yacimiento directo de
flujo de los fluidos.
 Toma en cuenta la localización de los pozos productores e
inyectores y sus condiciones operacionales.
 Los pozos pueden operarse y cerrarse de acuerdo a condiciones
específicas.
 Se pueden prefijar las tasas o las presiones de fondo o ambas.
 El yacimiento se divide en en multiples celdas o bloques (
tanques)
 Los cálculos se efectúan para las fases de petroleo, gas y agua e
intervalos discretos.

Etapas de un estudio integrado (esquema)

Valor de la Simulación.

 Permite conocer el comportamiento del yacimiento a diferentes

escenarios de explotación.
 Permite validar todos esos escenarios antes de seleccionar.
 Permite evaluar tecnologías nuevas de alto riesgo y planificar el
desarrollo de yacimientos complejos.

Costo de la simulación.

El costo de la simulación es muy pequeño comparado al número de

barriles recuperados, sin embargo, el proceso de organizar los datos y
construir el modelo pudiera ser todavía más valioso que su resultado, ya que
permite entender la naturaleza de los yacimientos.
El costo de la simulación incluye la mano de obra , la cual es menos de
un medio céntimo de dólar por barril.

Abuso de la Simulación

73
  Expectativas irreales.

 Tendencia a creer que “la respuesta es infalible, ya que salió del

computador”.

 Justificación suficiente para la Simulación.

 Se pide un trabajo que no es justificado con la data disponible.

 El objetivo de estudio no está claro.
 Estudios convencionales pueden dar resultados similares.

 íDescripción irreal del yacimiento.

Para hacer el cotejo de la historia se manipulan parámetros

(permeabilidad relativa y propiedades de los hidrocarburos.) hasta hacerlos
irreales.

Diez reglas de oro para ingenieros de simulación.

1. Entienda su problema y defina sus objetivos.

2. busque la sencillez
3. Entienda la interacción entre las diferentes partes.
4. No crea que siempre lo más complicado es mejor.
5. Conozca sus limitaciones y confíe en su juicio.
6. Sea razonable en sus expectativas.
7. Cuestiónese el ajuste de los datos en “cotejo histórico”
8. No suavice los extremos.
9. Ponga atención a las mediciones y use escalímetro
10.No deseche el necesario trabajo de laboratorio.

ENUNCIADO DE LOS PROBLEMAS.

Problema nº 1

74
Se tiene las siguientes características de un yacimiento:

POZO Φ(%) AREA DE AREA(ACRES)

PRODUCCION
UDO 1 13.3 A0 169.73
UDO 2 14.0 A1 156.62
UDO 3 14.5 A2 135.02
UDO 4 12.0 A3 117.92
UDO 5 12.0 A4 95.72
UDO 6 13.2 A5 69.55
A6 30.05

El crudo presente en el yacimiento tiene 29 ºAPI y se considera constante. El

yacimiento formo una pequeña capa de gas al iniciarse la producción.
Identifique y determine el tipo de yacimientos:

a) Volumen poroso del yacimiento.

b) Reservas originales de petróleo.

Otros datos adicionales:

Pi= 3600 Lpca; T= 191 ºF; Rs= 1200 PCN/BN; Swi= 25%

Problema nº2

Las propiedades PVT del fluido de yacimiento volumétrico de petróleo de

la arena “R” se presenta en la tabla nº 1. Cuando la presión del yacimiento
disminuye desde la presion inicial 2500 Lpc, a una presión promedia de
1600 Lpc, la producción correspondiente es de 26,0 MMBF. La RGP

75
acumulativa a 1600 Lpca es 954 PCS/BF y la RGP instantánea es 2250
PCS/BF. La porosidad promedia es 18 % y la saturación promedia de agua
connota es también 18%. La cantidad de agua producida es insignificante y
las condiciones normales so 14.7 Lpca y 60 ºF.

P (Lpca) Bo (by/Bf) Rs ( PCN/BF) Z

2500 1.270 575 0.870
2300 1.290 575 0.845
2200 1.300 575 0.835
2100 1.275 544 0.827
2000 1.272 513 0.820
1800 1.243 456 0.815
1600 1.215 381 0.803

76
 BIBLIOGRAFÍA.

 AMIX, James W: "Petroleum Reservoir Engineering. Physical

Properties". McGraw Hill. Classic Textbook Reissue. 1960.

 BRADLEY, Howard: "Petroleum Engineering Handbook". Socíety of

Petroleum Engineers, Richardson TX. 1992.

 CENTRO INTERNACIONAL DE EDUCACIÓN Y

DESARROLLO. "Caracterización física de los Yacimientos." CIED.
1997.

 B.C. Craft y M.F. Hawkins: “Ingeniería aplicada de yacimientos

petrolíferos”. Editorial Tecnos. Madrid.

 CIED — PDVSA: “Tópicos para ingenieros”. Caracas, enero de 1999.

 MANUCCIV, Jesús E: “Car

 acterización Física de yacimientos”. Mannyna consultores.

 T. E. W, Nind: “Fundamentos de producción y mantenimiento de

pozos petrolíferos.”. Editorial LIMUSA. México, 1987.

 CEIP. CONACEP. Maturín, octubre de 2000.

77
 SIMBOLOGÍA

°API: Grados API

Bt: Factor volumétrico total.
Bo: Factor volumétrico del petróleo.
Bg: Factor volumétrico del gas.
Bw : Factor volumétrico del agua.
By: Barriles de yacimiento.
BN: Barriles normales.
BF: Barriles fiscales.
Cf: Compresibilidad de la formación o de la roca.
Co: Compresibilidad del petróleo.
Cw: Compresibilidad del agua.
EBM: Ecuación de balance de materiales.
Gp: Gas producido.
Lpca: Libras por pulgadas cuadradas absolutas.
N: Petróleo original en el sitio.
Np: Petróleo producido.
Pb: Presión de burbujeo.
Py: Presión en el yacimiento.
PCN: Pies cúbicos normales.
Rs: Relación gas disuelto-petróleo.
Rsi: Relación gas disuelto-petróleo inicial.
RGP: Relación gas petróleo.
Rp: Relación gas petróleo acumulativa.
So: Saturación del Petróleo.
Sg: Saturación de gas.

78
Sor: Saturación residual de petróleo.
Vt: Volumen total
Vp: Volumen poroso.
We: Intrusión de agua.
P: Variación de presión.

79