Balance de Masa en Petróleo Saturado
Balance de Masa en Petróleo Saturado
Petróleo Saturado
Laura Sánchez M.
[Link]@g
[Link]
Balance de Masa – Petróleo Saturado
VP VP . ( 1 – Swirr )
Soi VP . Soi = POIS . Boi =
POIS . Boi + GOIS . Bgi
• m = GOIS . Bgi .
POIS . Boi
Con esto resuelto armemos las expansiones como en los casos anteriores.
Balance de Masa – Petróleo Saturado
Expansiones actuantes:
donde Bt = Bo + Rl . Bg ,
Rl = Rsi - Rs
• ΔVP = Cp . VP . ΔP
• ΔVw = Cp . Vw . ΔP = Cp . ( VP . Swirr ) . ΔP
• POIS . (Bt – Bti) + GOIS . (Bg – Bgi) + [Link] .(1 + m).(Cp + Cw . Swirr).ΔP
1 - Swirr
Expansión del petróleo Expansión del gas Expansión de los poros y el agua intersticial
original y su gas disuelto libre original en las zonas del casquete y de petróleo
Balance de Masa – Petróleo Saturado
PEROOOOOOOO…
Balance de Masa – Petróleo Saturado
Petróleo en superficie
Petróleo en fondo
Balance de Masa – Petróleo Saturado
Gp Gp . Bg
Np . Bo
Np
Petróleo en fondo
Es decir que al multiplicar por Bo y Bg, para llevar los volúmenes a fondo, se está
considerando dos veces al gas liberado por el petróleo.
Para corregir esta duplicación, debe descontarse una de ellas. El gas liberado, en
condiciones de fondo puede expresarse como:
• Np . Bo + Gp . Bg - Np . Rs . Bg
Expansión del Expansión del gas Expansión de los poros y el agua Entrada
petróleo original libre original intersticial en las zonas del de agua
y su gas disuelto casquete y de petróleo
Np . Bo + Gp . Bg - Np . Rs . Bg + Wp
Versión Amyx
Balance de Masa – Petróleo Saturado
Versión Dake
Balance de Masa – Petróleo Subsaturado
Por encima de No hay gas libre Por no haber casquete, no hay VP o Entrada
Pburb: Bt = Bo originalmente Vw en zona de gas de agua
Np . Bo + Gp . Bg - Np . Rs . Bg + Wp
=0
Balance de Masa – Petróleo Subsaturado
Ce
Esta ecuación puede utilizarse entre la presión original del reservorio y la presión
de burbuja. En el momento en que se traspone ese límite y se libera gas en fondo,
esta ecuación deja de ser válida.
No hay petróleo Expansión del gas Escrita de esta forma, la fórmula Entrada
originalmente libre original presenta una indeterminación de agua
=0
Np . Bo + Gpeq . Bg - Np . Rs . Bg + Wp
Si seguimos operando:
Gpeq . Bg + Wp
Si seguimos operando:
=0
GOIS.(Bg – Bgi) + ([Link] + [Link]).(Cp + Cw . Swirr).ΔP + We =
1 – Swirr
Gptotal . Bg + Wp
Que es la ecuación de balance de masa para yacimientos de gas, tanto de gas seco
como de gas húmedo y gas con condensación retrógrada.
Esta ecuación incluye la expansión de poros y agua intersticial en las zona de gas (no
existe zona de petróleo) que en yacimientos de gas podría ser despreciable. De
querer hacerlo la ecuación se reduce aún más:
• Eo = Bt – Bti
• Eg = Bg – Bgi
• F = Np . Bo + Gp . Bg – Np . Rs . Bg + Wp
Balance de Masa – Petróleo Saturado - Histórico
• POIS . Eo + GOIS . Eg + We = F
• POIS + GOIS . Eg / Eo + We / Eo = F / Eo
35
F / Eo
30 Acuífero Fuerte
25
20 Acuífero Débil
15 We = 0
η tg η = GOIS
10
POIS5
0
0 2 4 6 8 10
Eg / Eo
Balance de Masa – Petróleo Saturado - Predictivo
La diferencia es que en modo histórico, de cada momento del que se tienen datos
de presión y acumuladas, las incógnitas son las mismas. Para cualquier instante
de la vida del reservorio, POIS y GOIS son invariantes.
Es decir que dos puntos en el tiempo son suficientes para formar un sistema
determinado (aunque sería conveniente tener más de dos).
Balance de Masa – Petróleo Saturado - Predictivo
Pero cuando se usa en modo histórico, en que POIS y GOIS son conocidos, y las
acumuladas de petróleo y gas son incógnitas, estas varían para cada momento:
Por lo que el sistema nunca queda determinado. Cada punto futuro de presión
propuesto agrega una nueva ecuación con dos nuevas incógnitas.
• qg = ( 2 . π . h . kg / µg ) . ∂Pg / ∂r
• qo = ( 2 . π . h . ko / µo ) . ∂Po / ∂r
Integrando entre el radio de drenaje y el radio del pozo despreciando los efectos
capilares, y dividiendo ambas ecuaciones:
• [Link] = qg = kg . µo .
qo µg . ko
• GOR = Qg = kg . µo . Bo + Rs
Qo µg . ko Bg
Balance de Masa – Petróleo Saturado - Predictivo
• Yacimiento sin entrada de agua, con casquete gasífero actuando por difusión
• Yacimiento sin entrada de agua, con casquete gasífero actuando por segregación
Balance de Masa – Petróleo Saturado - Predictivo
Sg
Pi P<Pi
Swirr
Swirr
ti=0 Soi ΔP t>ti
So
• So = (POIS – Np) . Bo / VP
• So = (POIS – Np) . Bo . (1 – Swirr) = 1 – Np . Bo . (1 – Swirr)
[Link] POIS Boi
Balance de Masa – Petróleo Saturado - Predictivo
• Yacimiento sin entrada de agua, con casquete gasífero actuando por difusión
Casquete Casquete
Zona de Gas
Zona de Petróleo
Swirr
Swirr
ΔP
Pi Soi P<Pi
ti=0 So t>ti
• So = (POIS – Np) . Bo / VP
• So = (POIS – Np) . Bo . (1 – Swirr) = 1 – Np . Bo . (1 – Swirr)
[Link] POIS Boi
Balance de Masa – Petróleo Saturado - Predictivo
• Yacimiento sin entrada de agua, con casquete gasífero actuando por difusión
Casquete Casquete
Zona de Gas
Zona de Petróleo
Swirr
Swirr
ΔP
Pi Soi P<Pi
ti=0 So t>ti
Como puede verse, se llega a la misma ecuación que en caso anterior, con la
diferencia de que la entrada de masa del casquete a la zona de petróleo ayuda a
mantener la presión, aumentando la saturación de petróleo.
• Yacimiento sin entrada de agua, con casquete gasífero actuando por segregación
Casquete
Zona de Gas Casquete
Zona de Petróleo GOIS . (Bg – Bgi)
Swirr
Swirr
Pi Soi ΔP P<Pi
ti=0 So t>ti
• So = (POIS – Np) . Bo .
VP – GOIS . (Bg–Bgi) / (1 – Swirr)
Balance de Masa – Petróleo Saturado - Predictivo
• Yacimiento sin entrada de agua, con casquete gasífero actuando por segregación
Reemplazando VP:
• Yacimiento sin entrada de agua, con casquete gasífero actuando por segregación
• F= 1 .
1 – m . (Bg / Bgi - 1)
• So = 1 – Np . Bo . (1 - Swirr) . F.
POIS Boi
Donde F sólo tiene sentido físico cuando es positiva y mayor que uno.
Esta ecuación engloba los 3 casos, a condición de que se asigne a m valor 0 cuando
el casquete actúe por difusión.
Balance de Masa – Petróleo Saturado - Predictivo
Pero la relación entre todas las variables no es lineal, por lo tanto su solución no
siempre es trivial.
• Babson (1944)
• Tarner
• Tracy
• Muskat
• Pirson o diferencias finitas
• Bugari (2009)
Balance de Masa – Petróleo Saturado - Predictivo
Método de Tracy:
• Np’ . Фn + Gp’ . Фg = 1
Donde la única incógnita es GOR2 ya que todas las demás variables son conocidas
o fácilmente calculables con información de un PVT.
Recordemos que la primera ecuación complementaria justamente dejaba
planteado el incremental de petróleo en función del GOR de un momento futuro.
Balance de Masa – Petróleo Saturado - Predictivo
Si se está trabajando con Excel, puede plantearse todo el cálculo en modo iterativa
de forma tal de no tener que ajustar manualmente ningún valor, sino que el
programa encuentre el punto de equilibrio de cada escalón y corroborar su
coherencia.
ΔGp1-2
Gp’2
Gp2
Siendo (1) un momento conocido y (2) el momento desconocido que queremos hallar
Balance de Masa – Petróleo Saturado - Predictivo
Método de Tarner:
ΔNp1-2
Siendo (1) un momento conocido y (2) el momento desconocido que queremos hallar
Balance de Masa – Petróleo Saturado - Predictivo
Método de Bugari:
Se utiliza una única ecuación complementaria, pero antes hay que modificar la
ecuación de balance de masa, y definir nuevas variables:
• Gp = FGp + RGp
• 1 = FGp / Gp + RGp/Gp
• RGp = Np . Rs
Balance de Masa – Petróleo Saturado - Predictivo
Método de Bugari:
Método de Bugari:
RGp / Gp vs Np
1,00
0,90
0,80
0,70
0,60
0,50
0,40
0,30
0,20
0,10
0,00
0 2000 4000 6000 8000 10000 12000 14000
• Cálculo volumétrico
• Balance de materia
Cálculo volumétrico
RECURSO
NO DESCUBIERTO
COMERCIAL NO COMERCIAL
En Producción
Desarrolladas Cerradas
No En Producción
Reserva
Detrás de la Cañería
No desarrolladas
RESERVAS Y PECES
PROBADA DESARROLLADA: El pez está en el bote. Usted lo pesó puede olerlo y puede comerlo.
PROBADA NO DESARROLLADA: El pez está en su anzuelo dentro del agua junto al bote y usted
está listo para agarrarlo con una red. Usted puede decir lo grande que parece (siempre se ven
más grandes dentro del agua).
PROBABLE: Hay peces en el lago. Usted puede aún ser capaz de verlos, pero usted no ha
pescado ni uno.
POSIBLE: Hay agua en el lago. Alguien le dijo a usted que hay peces en ese lago. Usted tiene su
bote sobre el tráiler, pero en vez de ir a pescar usted puede ir a jugar al golf.
• Cálculo volumétrico
𝑉𝑟 ∗𝜃∗(1−𝑆𝑤𝑖)
𝑃𝑂𝐼𝑆 = 𝐵𝑜𝑖
• Balance de materia
𝑝
•¿Qué es el Recurso?
•¿Dónde se encuentra?
•¿Cómo se cuantifica?
• Descubiertas
• Recuperables
• Comerciales
• Remanentes
(Para un proyecto de
desarrollo determinado)
Descubrimiento
Descubrimiento de una acumulación para la cual un pozo exploratorio ha
establecido la existencia “significativa” de hidrocarburos
▫ Locaciones futuras
▫ Desarrollo de un proyecto de secundaria nuevo
▫ Cambio de presión de red en yacimientos de gas
▫ Work Over’s o pescas de alto costo estimado
▫ Reperforaciones
A menos que
exista una
P1
P2 justificación P2
técnica (ambas P2
P1 P3
direcciones)
P3
P3
1000
100
10
1
1938 1944 1950 1956 1962 1968 1974 1980 1986 1992 1998 2004 2010
qoP[m³/DC] NpP[Mm³]
DATOS HABITUALES:
• Inicio de Prod. = 11/1939 • Qo @ 1/2010 = 3.7 m3/d
• Qoi = 394 m3/d • Np @ 1/2010 = 1351 km3
• Qmax = 499 m3/d
Desarrolladas en Producción
Np = 2750 m3
EUR 1P = 4978 m3
Comportamiento Reserva P1 = 2228 m3
estable
10
0.1
2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019
qoCP[m³/DC] qoP[m³/DC] qlP[m³/DC]
Np: producción acumulada de petróleo // EUR: estimated ultimate recovery // P1 = 1P: reservas probadas
Desarrolladas en Producción
Reservas Probadas + Probables Desarrolladas en
Producción (ó Reservas 2P en Producción)
R.e-1005
100
Np = 2750 m3
EUR 2P = 6357 m3
Reserva 2P = 3607 m3
Reserva 1P = 2228 m3
10 Reserva P2 = 1379 m3
Límite: Fin de
la Concesión
1
0.1
2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019
qoIP[m³/DC] qoCP[m³/DC] qoP[m³/DC] qlP[m³/DC]
Np = 2750 m3
EUR 3P = 7768 m3
Reserva 3P = 5019 m3
Reserva 2P = 3607 m3
10 Reserva P3 = 1412 m3
Límite: Fin de
la Concesión
1
0.1
2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019
qoRP[m³/DC] qoIP[m³/DC] qoCP[m³/DC] qoP[m³/DC] qlP[m³/DC]
P3: reservas posibles // 3P = 2P + P3 = P1 + P2 + P3: reservas probadas + reservas probables + reservas posibles
Desarrolladas en Producción
10
0.1
2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019
qoRP[m³/DC] qoIP[m³/DC] qoCP[m³/DC] qoP[m³/DC] qlP[m³/DC]
Desarrolladas en Producción
De P1 a P2
P3
P2
• Reservas 2P considerando el
naranja + el verde
• Reservas Probadas
▫ LKG = - 3340 m
• Reservas Probables
▫ LKG + 1esp = - 3350 m
• Reservas Posibles
▫ Spill Point = - 3380 m
No Desarrolladas
LKH
Práctica Común: Reservas
Posibles hasta el Spill Point
1,0%
1,1% 1,2%
6,4%
1,4%
21,7%
Noroeste
Cuyana
Neuquina
Golfo San Jorge
Austral
Pcia. De Tierra del Fuego
Estado Nacional
67,3%
Volumen % de Reservas de Gas por Cuenca – Argentina
2012
10,1%
16,4% 0,2%
10,0%
Noroeste
Cuyana
5,5%
Neuquina
Golfo San Jorge
Austral
Pcia. De Tierra del Fuego
42,4% Estado Nacional
15,4%
Reservas Probadas de PETRÓLEO - Argentina
Noroeste
Cuyana
Neuquina
Austral
Noroeste
Cuyana
Neuquina
Austral
• Analógicos
• Analíticos
• Simulación Numérica
Tipos de pronósticos
Historia de producción
Análisis Declinatorio
• Generalmente se emplea en forma empírica (aunque hay
aportes teóricos interesantes –Fetkovich)
• Constancia en el régimen de explotación
• Puede hacerse a nivel:
capa, pozo, formación, yacimiento
• Principalmente a partir de curvas:
log Q Q
t Np
Análisis Declinatorio
(Log Scale)
Q
t Np
Análisis Declinatorio
• En 1924, Cutler, después de una intensa investigación, postula que la mejor
aproximación no era una línea recta en escala semilogarítmica, sino una
línea recta en escala log-log (es decir, ambos ejes representados
logarítmicamente), pero haciendo algunas ajustes del eje X.
• Esta sería la primera vez en aplicarse declinaciones hiperbólicas para
pronósticos en la industria del petróleo.
(Log Scale)
t (Log Scale)
Análisis Declinatorio
• En 1925, Larkey propone utilizar el método de mínimos cuadrados
para hallar la hiperbólica que mejor se ajusta a los datos históricos.
(ln q ) 1 q
Dn .
t q t
• Siendo ∂q/∂t la pendiente en un eje lineal.
Análisis Declinatorio
• Es decir, por ejemplo, que para una declinación exponencial:
Q Q
∂logq1
∂q/∂t1
(Log Scale)
∂t1
∂logq2 q1 ∂q/∂t2
∂t2 q2
t [años] t [años]
Q Q
t [años] t [meses]
Análisis Declinatorio
• Por lo tanto la declinación nominal respeta la relación:
qi q f
De
qi
1 q
Dn . k .q n
q t
n<1
• Hiperbólica (n≠0; n≠1) Exponencial
t
Análisis Declinatorio
Constanza Menghi
Análisis Declinatorio
Fariña y Lavia
n ni .ea.t
log Q
Exponencial
No es tan pesimista como
t
la exponencial.
Cálculo de “De” Método con Escuadra
Punto 1
Punto 2
Franco Molini
Cálculo de “De” Método con Escuadra
Punto 1
Punto 2
Franco Molini
Cálculo de “De” Método con Escuadra
Punto 1
Punto 2
Franco Molini
Cálculo de “De” Método con Escuadra
Punto 1
Punto 2
Punto 1
trasladado
Punto 2
trasladado
Franco Molini
Cálculo de “De” Método con Escuadra
Punto 1
Punto 2
Punto 1
trasladado
Punto 2
trasladado
qi=10
1,5 Δq
qf=8,5
1 año
𝒒𝒊 − 𝒒𝒇 𝟏𝟎 − 𝟖, 𝟓
𝑫𝒆𝒂𝒏 = = = 𝟏𝟓%
𝒒𝒊 𝟏𝟎
Cálculo de “De” Método con Escuadra
qi=10
1,5 Δq
qf=8,5
1 año
𝒒𝒊 − 𝒒𝒇 𝟏𝟎 − 𝟖, 𝟓
𝑫𝒆𝒂𝒏 = = = 𝟏𝟓%
𝒒𝒊 𝟏𝟎