Logo EBM-9000 EBM-9000 Ecuación de Balance de Materiales · simulador didáctico
Listo para simular
Modelo de tanque: se representan volúmenes a condiciones de yacimiento, no la posición real de los fluidos. Pasa el puntero sobre cada elemento.
Año 0.0 de 25 Arrastra para recorrer la simulación

Teoría: la Ecuación de Balance de Materiales

Unidades:

1. ¿Qué es el balance de materiales?

La ley de conservación de la masa es la base de los cálculos de balance de materiales. Un balance de materiales es una contabilidad del material que entra y sale de un sistema. El cálculo trata al yacimiento como un gran tanque y usa cantidades que sí se pueden medir para determinar la cantidad de un material que no se puede medir directamente.

Las cantidades medibles son: los volúmenes acumulados producidos de petróleo, agua y gas; presiones de yacimiento precisas; y las propiedades de los fluidos (PVT) obtenidas de muestras de los fluidos producidos. Con ellas se determinan, por ejemplo, el petróleo original en sitio N, el gas original G o la intrusión de agua We.

💡 En la pestaña Simulación el tanque de la izquierda es exactamente esa idea: un recipiente cuyo volumen poroso está ocupado por gas, petróleo y agua. Cuando la presión cae, los fluidos y la roca se expanden y esa expansión es la que "empuja" la producción.

2. ¿Para qué sirve?

  • Estimar volúmenes de forma independiente. Da una estimación del volumen de petróleo, agua y gas del yacimiento que se puede comparar con la estimación volumétrica (mapas, registros) o con la de un simulador numérico. Comparar la volumetría calculada por un simulador con la de un estudio de balance de materiales es una forma de validar ambas con fuentes independientes.
  • Identificar los mecanismos de empuje. La magnitud de cada término de la ecuación indica la contribución relativa de los distintos mecanismos de producción que actúan en el yacimiento (expansión del petróleo y su gas disuelto, capa de gas, acuífero, compactación, inyección).
  • Predecir el comportamiento futuro del yacimiento y ayudar a estimar la eficiencia de recuperación acumulada.

Referencias clásicas sobre estos usos: Dake [2001], Craft et al. [1991], Ahmed [2000], Towler [2002] y Pletcher [2002] (ver Bibliografía).

3. Ecuación general para un yacimiento de petróleo

La forma de la ecuación depende de si el yacimiento es predominantemente de petróleo o de gas. Para un yacimiento de petróleo, la ecuación general es la ecuación de balance de materiales de Schilthuis (1936), expresada en la forma dada por Guerrero (1966):

(2.6)\begin{aligned} & N\left(B_t - B_{ti}\right) + N m B_{ti}\left(\frac{B_{gc} - B_{gi}}{B_{gi}}\right) \\ & + N \frac{B_{ti} S_{wio}}{1 - S_{wio}}\left(\frac{B_{tw} - B_{twi}}{B_{twi}}\right) + N \frac{m B_{ti} S_{wig}}{1 - S_{wig}}\left(\frac{B_{tw} - B_{twi}}{B_{twi}}\right) \\ & + N\left(\frac{1}{1 - S_{wio}} + \frac{m}{1 - S_{wig}}\right) B_{ti}\, c_f\, \Delta P \\ & = N_p B_o - N_p R_{so} B_g + \left[G_{ps} B_g + G_{pc} B_{gc} - G_i B_g'\right] - \left(W_e + W_i - W_p\right) B_w \end{aligned}

Esta ecuación se deduce conservando el volumen (no la masa directamente) y por eso Dake [2001] la llama balance de materiales volumétrico. Se considera una ecuación general porque se puede aplicar a un yacimiento de petróleo con capa de gas y con acuífero; si no hay capa de gas basta con hacer m = 0, y si no hay acuífero ni inyección, We = Wi = Gi = 0.

Cada término tiene un significado físico: el lado izquierdo reúne los cambios de volumen (expansiones) y el lado derecho la producción y la inyección de fluidos. Antes de mirar la tabla de significados, conviene nombrar cada expansión.

4. Los términos de expansión

Para mostrar el significado físico de la ecuación (2.6) se definen primero los términos D, que son cambios de volumen por cada barril de petróleo original (RB/STB):

(2.7)\begin{aligned} D_o &= B_t - B_{ti}, \\[2pt] D_{go} &= m B_{ti}\left(\frac{B_{gc} - B_{gi}}{B_{gi}}\right), \\[2pt] D_w &= \frac{B_{ti} S_{wio}}{1 - S_{wio}}\left(\frac{B_{tw} - B_{twi}}{B_{twi}}\right), \\[2pt] D_{gw} &= \frac{m B_{ti} S_{wig}}{1 - S_{wig}}\left(\frac{B_{tw} - B_{twi}}{B_{twi}}\right), \\[2pt] D_r &= \left(\frac{1}{1 - S_{wio}} + \frac{m}{1 - S_{wig}}\right) B_{ti}\, c_f\, \Delta P \end{aligned}
Do · petróleo y su gas disuelto
¿Cuánto cambió el volumen del petróleo inicial y su gas asociado?

Usa el factor volumétrico total Bt = Bo + (Rsoi − Rso)·Bg: el petróleo que queda más el gas que se liberó de él. Por encima del punto de burbuja Bt = Bo y Do es solo la expansión del líquido; por debajo, el gas liberado domina.

Dgo · gas libre (capa de gas)
¿Cuánto cambió el volumen del gas libre?

La capa de gas inicial ocupa m·N·Bti barriles de yacimiento; al caer la presión cada pie cúbico ocupa más (Bgc > Bgi) y la capa se expande hacia la zona de petróleo.

Dw, Dgw · agua connata
¿Cuánto cambió el volumen del agua connata inicial?

El agua que ya estaba en los poros de la zona de petróleo (Swio) y de la capa de gas (Swig) también se expande. Con Btw ≈ Bw, el paréntesis vale cw·ΔP.

Dr · roca (volumen poroso)
¿Cuánto cambió el volumen poroso de la formación?

Al bajar la presión del fluido, la roca soporta más carga y el poro se reduce (compresibilidad cf). Ese volumen "perdido" expulsa fluido hacia los pozos.

5. Forma compacta y significado físico de cada término

Sustituyendo (2.7) en (2.6) se obtiene la ecuación general en la forma:

(2.8)N\left[D_o + D_{go} + D_w + D_{gw} + D_r\right] = N_p B_o - N_p R_{so} B_g + \left[G_{ps} B_g + G_{pc} B_{gc} - G_i B_g'\right] - \left(W_e + W_i - W_p\right) B_w

Los términos del lado derecho representan la producción e inyección de fluidos; los del lado izquierdo representan cambios de volumen. Su significado físico (Tabla 2-1 de Fanchi) convertido en la pregunta que responde cada término:

TérminoSignificado físicoPregunta a la que responde
N·DoCambio de volumen del petróleo inicial y su gas asociado¿Cuánto cambió el volumen del petróleo inicial y su gas asociado?
N·DgoCambio de volumen del gas libre¿Cuánto cambió el volumen del gas libre?
N·(Dw + Dgw)Cambio de volumen del agua connata inicial¿Cuánto cambió el volumen del agua connata inicial?
N·DrCambio del volumen poroso de la formación¿Cuánto cambió el volumen poroso de la formación?
Np·BoProducción acumulada de petróleo¿Cuánto volumen de yacimiento ocupa el petróleo producido acumulado?
Np·Rso·BgGas acumulado producido en solución con el petróleo¿Cuánto gas se produjo disuelto en el petróleo? (ya está dentro de Np·Bo, por eso se resta)
Gps·BgGas de solución acumulado producido como gas liberado¿Cuánto gas de solución se produjo como gas liberado?
Gpc·BgcProducción acumulada de gas de la capa de gas¿Cuánto gas de la capa de gas se ha producido?
Gi·Bg′Inyección acumulada de gas¿Cuánto gas se ha inyectado?
We·BwIntrusión acumulada de agua¿Cuánta agua ha entrado del acuífero?
Wi·BwInyección acumulada de agua¿Cuánta agua se ha inyectado?
Wp·BwProducción acumulada de agua¿Cuánta agua se ha producido?
📌 Lectura rápida. Todo lo que está a la izquierda es energía del yacimiento (volumen que "aparece" al expandirse fluidos y roca). Todo lo que está a la derecha es volumen retirado por los pozos, corregido por lo que se inyecta o entra del acuífero. En la simulación, la intrusión y la inyección se dibujan junto a las expansiones (lado izquierdo) como fuentes de vaciamiento: es la misma ecuación con los términos WeBw, WiBw y GiBg′ pasados de lado.

6. Suposiciones de la deducción

  1. El sistema está en equilibrio de presión (una sola presión media representa todo el yacimiento).
  2. El sistema es isotérmico.
  3. Los datos PVT disponibles son representativos de los fluidos del yacimiento.
  4. Los datos de producción son confiables.
  5. La segregación gravitacional de las fases se puede despreciar (el modelo de tanque no tiene dimensiones).

La importancia relativa de los mecanismos de empuje que se obtiene de la ecuación (2.6) se discute en el capítulo 7 de Fanchi y, con los índices de empuje que deben sumar uno, en Pletcher [2002].

7. Balance de materiales para un yacimiento de gas

La ecuación general para un yacimiento de gas se deduce de la (2.6) reconociendo primero que la relación

(2.9)G\, B_{gi} = N\, m\, B_{ti}

define el gas original en sitio G (el volumen de la capa de gas a condiciones de yacimiento es el mismo escrito con gas o con petróleo). Sustituyendo (2.9) en (2.6):

(2.10)\begin{aligned} & N\left(B_t - B_{ti}\right) + G B_{gi}\left(\frac{B_{gc} - B_{gi}}{B_{gi}}\right) + N \frac{B_{ti} S_{wio}}{1 - S_{wio}}\left(\frac{B_{tw} - B_{twi}}{B_{twi}}\right) \\ & + G \frac{B_{gi} S_{wig}}{1 - S_{wig}}\left(\frac{B_{tw} - B_{twi}}{B_{twi}}\right) + \left(\frac{N B_{ti}}{1 - S_{wio}} + \frac{G B_{gi}}{1 - S_{wig}}\right) c_f\, \Delta P \\ & = N_p B_o + \left[G_{ps} B_g + G_{pc} B_{gc} - G_i B_g'\right] - N_p R_{so} B_g - \left(W_e + W_i - W_p\right) B_w \end{aligned}

Como el balance de un yacimiento de gas no incluye petróleo en sitio, N = 0 y Np = 0, y la ecuación se reduce a:

(2.11)G B_{gi}\left(\frac{B_{gc} - B_{gi}}{B_{gi}}\right) + G \frac{B_{gi} S_{wig}}{1 - S_{wig}}\left(\frac{B_{tw} - B_{twi}}{B_{twi}}\right) + \left(\frac{G B_{gi}}{1 - S_{wig}}\right) c_f\, \Delta P = \left[G_{pc} B_{gc} - G_i B_g'\right] - \left(W_e + W_i - W_p\right) B_w

La compresibilidad del agua y la de la formación son pequeñas comparadas con la del gas; por eso la (2.11) suele escribirse en la forma simplificada:

(2.12)G B_{gi}\left(\frac{B_{gc} - B_{gi}}{B_{gi}}\right) = \left[G_{pc} B_{gc} - G_i B_g'\right] - \left(W_e + W_i - W_p\right) B_w

Si además no hay acuífero ni inyección, (2.12) se convierte en la conocida recta P/z vs Gp de los yacimientos volumétricos de gas.

8. Simbología y unidades

Unidades de campo entre corchetes según la nomenclatura de Fanchi; la columna SI se obtiene con las conversiones exactas 1 bbl = 0.158987 m³, 1 ft³ = 0.0283168 m³ y 1 psi = 6.894757 kPa. En la práctica hispanoamericana: STB = BN/BF (barriles normales o fiscales), SCF = PCN/PCS (pies cúbicos normales o estándar), RB = BY (barriles de yacimiento).

SímboloDefiniciónUnidad

9. La ecuación con los datos de la simulación

Primero la ecuación y luego su reemplazo con los datos del instante seleccionado de la simulación (el mismo que muestra la pestaña Simulación). Si cambias el sistema de unidades arriba, todos los valores se convierten automáticamente; los volúmenes de yacimiento resultantes se expresan en MMRB.

Año 0.0 Datos: ejemplo

10. ¿Cómo genera los datos esta simulación?

EBM-9000 no resuelve un simulador numérico: usa el propio balance de materiales como modelo de tanque, con estas reglas didácticas:

  • La presión declina desde Pi hasta la presión final con una curva prescrita (rápida al inicio, más lenta al final).
  • Las propiedades PVT (Bo, Rso, Bg, Bw) se interpolan entre los valores inicial y final que das; por encima del punto de burbuja Rso es constante y el petróleo solo se expande con co. Se toma Bgc = Bg′ = Bg y Btw = Bw (sin gas disuelto en el agua).
  • La intrusión de agua sigue el modelo de estado estable de Schilthuis (proporcional a ∫(Pi − P)dt), escalada para alcanzar el We final indicado. La inyección crece linealmente desde el instante de inicio.
  • La RGP de producción y el corte de agua siguen curvas suaves hasta los valores finales; el gas libre producido se toma de la capa de gas y del gas liberado en proporción a los volúmenes disponibles.
  • En cada paso se calcula el vaciamiento total (expansiones + intrusión + inyección) y se reparte en petróleo, gas y agua producidos, de modo que la ecuación (2.8) se cumple exactamente en todo instante: por eso las dos barras del panel siempre miden lo mismo.
  • La simulación termina en el punto muerto: producción de agua (corte de agua ≥ 90 %) o depletación (presión de abandono).

11. Autores de la teoría

Ralph J. SchilthuisHumble Oil & Refining Co. · 1936

Publicó "Active Oil and Reservoir Energy" (Trans. AIME 118), donde formuló por primera vez la ecuación general de balance de materiales para yacimientos de petróleo y el modelo de intrusión de agua en estado estable que lleva su nombre.

E. T. GuerreroUniversidad de Tulsa · 1966

En Practical Reservoir Engineering presentó la forma de la ecuación de Schilthuis que adopta Fanchi: expansiones por barril de petróleo original y producción neta a condiciones de yacimiento.

John R. FanchiColorado School of Mines · Texas Christian University · 2006

Autor de Principles of Applied Reservoir Simulation (3.ª ed.), la fuente de esta pestaña. Reúne la ecuación general, sus términos de expansión, la Tabla 2-1 de significado físico y las formas para yacimientos de gas.

Laurence P. DakeShell · 1978, 2001

En Fundamentals of Reservoir Engineering y The Practice of Reservoir Engineering desarrolló la deducción volumétrica ("balance de materiales volumétrico") y su interpretación por mecanismos de empuje.

D. Havlena y A. S. Odeh1963–1964

Reordenaron la ecuación como la ecuación de una línea recta (F vs. Eo, Eg, Ef,w), el método gráfico estándar para estimar N, m y We a partir de datos de producción.

Craft, Hawkins y Terry · Ahmed · Towler · Pletcher1991 · 2000 · 2002 · 2002

Textos y artículos citados por Fanchi para profundizar: derivaciones paso a paso, ejemplos numéricos, índices de empuje y mejoras a los métodos de balance de materiales (efecto de acuíferos débiles y de la compresibilidad de la formación).

Bibliografía y fuentes consultadas

  1. Fanchi, J. R. (2006). Principles of Applied Reservoir Simulation (3.ª ed.), cap. 2 "Basic Reservoir Analysis", sección 2.3 Material Balance, pp. 16–20. Gulf Professional Publishing / Elsevier, Burlington, MA. ISBN 978-0-7506-7933-6. shop.elsevier.com
  2. Schilthuis, R. J. (1936). "Active Oil and Reservoir Energy". Transactions of the AIME, 118(1), 33–52. doi:10.2118/936033-G
  3. Guerrero, E. T. (1966). Practical Reservoir Engineering. The Petroleum Publishing Company, Tulsa, OK.
  4. Dake, L. P. (2001). The Practice of Reservoir Engineering (ed. revisada). Developments in Petroleum Science 36. Elsevier, Ámsterdam. ISBN 978-0-444-50671-9. shop.elsevier.com
  5. Dake, L. P. (1978). Fundamentals of Reservoir Engineering. Developments in Petroleum Science 8. Elsevier, Ámsterdam.
  6. Craft, B. C., Hawkins, M. F. y Terry, R. E. (1991). Applied Petroleum Reservoir Engineering (2.ª ed.). Prentice Hall, Englewood Cliffs, NJ.
  7. Ahmed, T. (2000). Reservoir Engineering Handbook. Gulf Publishing Company, Houston, TX.
  8. Towler, B. F. (2002). Fundamental Principles of Reservoir Engineering. SPE Textbook Series Vol. 8. Society of Petroleum Engineers, Richardson, TX. ISBN 978-1-55563-092-8. store.spe.org
  9. Pletcher, J. L. (2002). "Improvements to Reservoir Material-Balance Methods". SPE Reservoir Evaluation & Engineering, 5(1), 49–59. SPE-75354-PA. doi:10.2118/75354-PA
  10. Havlena, D. y Odeh, A. S. (1963). "The Material Balance as an Equation of a Straight Line". Journal of Petroleum Technology, 15(8), 896–900. SPE-559-PA. doi:10.2118/559-PA
  11. Havlena, D. y Odeh, A. S. (1964). "The Material Balance as an Equation of a Straight Line — Part II, Field Cases". Journal of Petroleum Technology, 16(7), 815–822. SPE-869-PA. doi:10.2118/869-PA
  12. Society of Petroleum Engineers. The SI Metric System of Units and SPE Metric Standard. Factores de conversión usados en esta app. spe.org
Logo EBM-9000

EBM-9000

Simulador didáctico de la Ecuación de Balance de Materiales

Objetivo

Explicar de forma didáctica e interactiva la Ecuación de Balance de Materiales (EBM), mostrando qué pregunta responde cada término y cómo la expansión de los fluidos y la roca, la intrusión de agua y la inyección compensan el vaciamiento que produce la extracción de petróleo, gas y agua.

Resumen

EBM-9000 es una aplicación web educativa que anima, en cinco minutos, la vida productiva de un yacimiento tipo tanque: la presión declina, el petróleo y el gas se expanden, el acuífero avanza y los pozos producen hasta un punto muerto por depletación o por producción de agua. El usuario puede reproducir, pausar y recorrer la simulación, usar el ejemplo incluido o ingresar sus propios datos, y consultar en la pestaña de teoría la formulación general de Schilthuis presentada por Fanchi, la simbología, la sustitución numérica con los datos simulados y la conversión automática entre unidades de campo y Sistema Internacional. Está construida con HTML, CSS y JavaScript puros para desplegarse directamente desde GitHub Pages o Netlify.

Fotografía de Edgar Izurieta

Ing. Edgar Izurieta

Ingeniero en Petróleos | Especialista en Datos e Inteligencia Artificial | Aspirante a Ingeniero de Reservorios / EOR | Creador de Contenido

Fuentes y herramientas

  • Teoría: Fanchi, J. R. (2006). Principles of Applied Reservoir Simulation, 3.ª ed., sección 2.3 (Elsevier). Ecuación original de Schilthuis (1936); ver la bibliografía completa en la pestaña de teoría.
  • Tecnología: HTML, CSS y JavaScript sin dependencias de compilación; ecuaciones renderizadas con KaTeX; tipografías Exo 2 e Inter (Google Fonts).
  • Despliegue: GitHub Pages (flujo de trabajo incluido en el repositorio) o Netlify conectado al repositorio edgardoomer/EBM-9000.