1.3.1 Simulación de Procesos Normales con Simuladores de Petróleo Negro.

Los simuladores de Petróleo Negro fueron los primeros desarrollados, son aún
los más usados. Estos simuladores pueden modelar el flujo de agua, petróleo y gas y
pueden contabilizar la solubilidad del gas en el petróleo dependiente de la presión,
pero ellos no pueden modelar los cambios en la composición de petróleo o gas.
Cuando los simuladores modernos de petróleo negro se usan apropiadamente, el
usuario puede confiar en que se tendrán respuestas realistas si los datos de entrada
son completos y razonablemente exactos. Las primeras versiones de estos
simuladores no garantizaban esta confianza, pero los simuladores de petróleo negro
actuales son tan confiables que se usan para verificar nuevos métodos analíticos.

1.3.2 Aplicaciones de Modelos de Campo Completo

Los simuladores de petróleo negro que tienen los mejores incentivos económicos
modelan reservorios completos o los segmentos mayores del reservorio. Tales
modelos pueden ser aplicados rentablemente a muchos tipos de estudios, tal como
modelo de inyección de agua de costado, inyección de gas, operación competitiva y
unitaria, ubicación y medición de pozos de desarrollo. Pueden también usarse para
estudiar modelos de espaciado de pozos, intervalos de terminación, caudales de
producción, medición gas-lift, y otras opciones operativas.

1.3.3 Modelos de Corte Transversal y Cilíndricos

Los modelos de corte transversal y cilíndrico (a veces llamados de campo simple,

radial, o r-z) son usados para simular segmentos de un reservorio. Los modelos son
comúnmente usados como entrada para los modelos areales.
Modelos radiales fueron usados en estudios del campo donde la conificación a
pozos individuales tiene un efecto significante en el comportamiento general del
reservorio. Un simulador radial capaz de modelar fracturas producidas manualmente
fue usado para estudiar la producción de gas de baja permeabilidad.

1.3.4 Simuladores de Propósito Especial

Después de la aparición de los simuladores de petróleo negro, se desarrollaron

otros programas para resolver problemas que eran matemáticamente más complejos
que los problemas del petróleo negro. Hay disponibles simuladores para modelo
composicional, termal, y procesos químicos; que se aplican normalmente en estudios
de proyectos EOR.
Los simuladores composicionales pueden modelar el rendimiento de reservorios
de petróleo volátil y gas condensado en los cuales las composiciones de las fases
varían ampliamente con la presión. Un uso típico es la optimización del reciclaje de
gas en la capa de gas o un reservorio de gas condensado.
Los simuladores de procesos termales pueden modelar reciclaje y flujo de vapor

Ing. Hermas Herrera Callejas Página: 1 de 21

Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

Simuladores de procesos químicos se han desarrollado para inyección de

polímeros, flujo de surfactante y soluciones alcalinas. Estos métodos de
desplazamiento de petróleo involucran cambios en propiedades de los fluidos
complejos y reacciones químicas

1.4 Planificación de un Estudio de Simulación.

Un estudio comprensivo puede tomar un año o más para completarse y, algunas

veces, puede requerir intensas demandas del hardware de computadoras y personal
calificado. Estudios menos comprensivos requieren menos recursos pero usualmente
deben ser conducidos bajo severas restricciones de tiempo. Ambos tipos de estudios
deben seguir planes claros y prácticos para asegurar que ellos proporcionen la
información correcta al equipo de administración del reservorio con el detalle apropiado
y, sobre todo, a tiempo para ser usada efectivamente.
La mayoría de los estudios involucran esencialmente los mismos tipos de
actividades, aunque la distribución del esfuerzo entre las actividades variará de un
proyecto a otro. La Fig. 1.1 grafica las actividades más significativas.
Tiempo Actividades
Definición del problema
Adquisición y Revisión de datos
Selección del enfoque
Descripción del reservorio y diseño del modelo
Soporte de programación
Ajuste Histórico
Predicción del Comportamiento y Análisis
Emisión de Reportes
Fig. 1.1 – Actividades en un típico estudio de simulación de reservorio

1.4.1 Definición del Problema

El primer paso en la conducción del estudio es definir el problema de

rendimiento del reservorio y el problema operativo asociado. Para hacer esto, uno
debe reunir suficiente información del reservorio y su ambiente operativo para
identificar qué proyecciones de rendimiento son necesarias, cuándo serán necesarias
y como pueden contribuir a la administración del reservorio.
Tan pronto como los datos experimentales adecuados han sido acumulados,
deben definirse clara y concisamente objetivos prácticos para el estudio. Muchas
decisiones que deben tomarse a medida que progresa el trabajo, seguirán
naturalmente si se tiene cuidado en desarrollar objetivos realistas. Aunque esta fase
del proyecto es relativamente corta, puede tener un mayor impacto en la eficiencia con
la cual es conducido el proyecto.
Una vez que los objetivos y el alcance del estudio están claros, debería
conducirse un análisis preliminar de la mecánica de fluidos del reservorio. El análisis
puede obtenerse de evaluaciones anteriores o puede involucrar cálculos manuales,
pequeños modelos de simulación o simplemente una revisión del rendimiento de otros
reservorios que parecen ser similares en los aspectos apropiados al reservorio en
estudio. El propósito del análisis es identificar el mecanismo principal de agotamiento
Ing. Hermas Herrera Callejas Página: 2 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

y reconocer los factores que afectarán el rendimiento del reservorio tales como
conificación, gravedad, heterogeneidad. El análisis debería ser lo suficientemente
completo como para permitir que los objetivos del estudio sean restablecidos en
términos de la mecánica del reservorio. Entonces debería ser posible determinar el
nivel de complejidad en el modelo del reservorio, para comenzar el diseño del modelo
e identificar los datos necesarios para construirlo.

1.4.2 Adquisición y Revisión de los Datos

Los datos normalmente deben ser revisados y reorganizados una vez hayan
sido recolectados porque habrán sido obtenidos por un número de razones
relacionadas a su resguardo y normalmente no habrán sido depurados y organizados
tan bien como para ser usados de inmediato. Debido a que una revisión detallada de
los datos puede ser morosa y tediosa, el esfuerzo debe ser cuidadosamente enfocado.
La revisión puede ser tan extensiva como sea necesario pero diseñada para evitar
trabajo innecesario.
La revisión de los datos disponibles debería revelar siempre brechas e
inconsistencias que necesitan ser resueltos. Algunas veces será necesario decidir si
hay los suficientes datos y si su calidad es la adecuada para construir un modelo de
reservorio con suficiente exactitud a fin de lograr los objetivos del estudio. Si los datos
disponibles, modificados y complementados por la experiencia y el buen juicio, son
inadecuados, los objetivos deben ser disminuidos o bien deben recolectarse datos
adicionales. Debemos tener presente, sin embargo, que muchos objetivos pueden ser
conseguidos aún con datos pobres para la evaluación de la sensitividad del
rendimiento del reservorio y la descripción del reservorio en un rango de valores
pensados para alcanzar los valores reales.

1.4.3 Selección del Enfoque de Estudio

Habiendo definido los problemas de mecánica de fluidos, uno debe decidir qué
modelos de simulación son los más adecuados para resolverlos. No siempre es
necesario ni deseable intentar modelar el reservorio completo. En estudios de
conificación, cuspidificación o severo decrecimiento en corrida de pozo individual, por
ejemplo, deben usarse modelos seccionales o segmentados.
Los factores que influyen en el enfoque de estudio son:
1) Disponibilidad de simuladores que pueden resolver adecuadamente los problemas
mecánicos del reservorio.
2) Cambios de programación que deben hacerse al simulador para modelar pozos y
facilidades.
3) El tipo y número de corridas del simulador necesarias para alcanzar los objetivos
del estudio.
4) Tiempo, mano de obra, computación y recursos financieros disponibles para el
estudio.
5) La necesidad de capacidades especiales de edición
6) Disponibilidad de recursos periféricos necesarios para concluir el estudio a tiempo.

1.4.4 Descripción del Reservorio y Diseño del Modelo

Ing. Hermas Herrera Callejas Página: 3 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

El diseño de un modelo de simulación es influenciado por el tipo de proceso a

ser modelado, la dificultad del problema de la mecánica de fluido, los objetivos del
estudio, la calidad de los datos descriptivos del reservorio, restricciones de tiempo y
presupuesto y el nivel de credibilidad necesario para asegurar aceptación de los
resultados del estudio.
Como se mencionó previamente, los simuladores más comunes son los de
petróleo negro que usan datos tabulares de propiedades de los fluidos y funciones de
permeabilidad relativa convencional. La simulación de algunos procesos requiere el
uso de simulador de propósito especial, a menudo soportado por programas periféricos
para predecir tales factores como pérdida de calor en boca de pozo o comportamiento
de la fase de hidrocarburos.
Las restricciones de tiempo o costo frecuentemente imponen compromisos en
el tipo de simulador a usar y en el diseño del modelo de reservorio. El número de
bloques-malla y el detalle incluido en el tratamiento de los pozos individuales son tal
vez las dos áreas más comunes para el compromiso. El ingeniero tendrá que
determinar, dentro del esqueleto de los objetivos del estudio, el nivel de sofisticación
realmente necesario para resolver el problema.

1.4.5 Soporte de Programación

Después de que se ha escogido el simulador apropiado y diseñado el modelo,

es necesario confeccionar partes del programa específicamente para el problema a
mano. Las modificaciones más comunes están relacionadas con la administración de
los pozos y la edición de resultados.
Los programas de administración de pozos pueden ayudar a automatizar la
simulación usando lógica que traduce las condiciones operativas del campo y
restricciones a condiciones de límites matemáticos para el simulador. Una rutina
comprensiva de manejo de pozos puede incorporar suficiente lógica y guías de
operación para hacer que la mayoría de las decisiones necesarias para simular
realistamente la administración del reservorio no debería ser más complicada de lo
necesario. En algunos casos es necesario solo especificar caudales de pozos o
presiones de la capa de arena; la mayoría de los simuladores deberá tener estas
opciones disponibles. Otras veces, la rutina de cálculo debe tomar en cuenta las
restricciones impuestas por facilidades de superficie, boca de pozo e hidráulica de flujo
de la línea, elevación artificial, objetivos de producción, y obligaciones reguladas o
contractuales.
Ocasionalmente, puede necesitarse acoplar un simulador en la línea de
producción de flujo multifase del reservorio de modo que la presión, facilidades de
producción e inyección sea impuesta al modelo de reservorio durante todo el tiempo
de simulación. Es práctica común incluir lógica en la rutina de administración de pozos
para adición, cierres, o reacondicionamiento de pozos.
En general, los aspectos más complicados de una rutina de administración de
pozos se requieren durante la fase de predicción. Durante la fase de history matching,
los caudales de pozos proporcionados al simulador serán corregidos por la producción
histórica e inyección. Por tanto, el history matching puede comenzar y el modelo puede
ser sintonizado mientras se programa la rutina de gestión de pozos.
Ing. Hermas Herrera Callejas Página: 4 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

1.4.6 Ajuste Histórico

Después de que se ha construido el modelo del reservorio, debe ser probado

para ver si puede reproducir el comportamiento del campo. Generalmente la
descripción del reservorio usada en el modelo es validada corriendo el simulador con
datos históricos de producción e inyección y comparadas las presiones calculadas y
movimientos de fluidos con el comportamiento real del reservorio. Una prueba más
severa es hacer que el simulador calcule el rendimiento pasado de pozos individuales
así como presiones históricas y movimientos de fluidos. Los datos usados en el history
matching varian con el alcance del estudio, pero usualmente incluyen presión de
reservorio y datos de producción. El proceso de history matching puede ser tan simple
que involucre solo parámetros de acuífero, o puede ser un procedimiento complejo que
involucre tanto el reservorio como el acuífero.
Los parámetros de entrada al modelo deben ser ajustados hasta obtener una
aceptable coincidencia entre resultados del modelo y datos del campo. Tener en
cuenta que cuando se ajustan parámetros del reservorio durante el history matching el
objetivo es describir el reservorio real tan exactamente como sea posible con los datos
disponibles. Típicamente uno ajusta:
1) Permeabilidad del reservorio para coincidir gradientes de presión del campo.
2) La permeabilidad y extensión del área de las arcillas y otras zonas de baja
permeabilidad para coincidir el movimiento vertical de fluidos.
3) Relaciones de Permeabilidad/Saturación relativas para coincidir distribución de
saturación dinámica y gradientes de presión.
4) Tamaño del acuífero, porosidad, espesor y permeabilidad para coincidir la cantidad
y distribución de influjo natural de agua.
Los reservorios son tan complejos que ninguna cantidad práctica de datos
descriptivos del reservorio pueden proporcionar el detalle necesario para predecir el
comportamiento futuro con seguridad completa y exactitud. La validez de los datos
descriptivos del reservorio, usados en un modelo, pueden evaluarse calculando el
rendimiento pasado del reservorio y comparando el rendimiento calculado con el real.
Este proceso de History Matching o coincidencia con el histórico descubrirá las
debilidades en la base de datos del reservorio que pueden ser fortificados alterando el
modelo del reservorio hasta que pueda predecir exactamente el rendimiento pasado.
Usado de esta manera, un simulador puede ser una herramienta poderosa de
descripción del reservorio, siempre que los cambios de datos y adiciones sean
consistentes con la base de datos conocida del reservorio.
Los datos que son emparejados son presión, WOR, GOR, relación gas/agua,
movimiento del fluido de contacto, tiempos de arribo de agua y gas, y saturaciones
medidas en muestras, registros de pozos, y pruebas químicas.
El History Matching puede ser tedioso, caro y frustrante, primeramente porque
el rendimiento del reservorio puede ser complejo, respondiendo a numerosas
iteraciones que, como un todo, puede ser dificultoso entender. Como consecuencia ha
habido investigación en programas para simulación inversa o coincidencia histórica
automática. Los programas son raramente usados, sin embargo, porque para la
mayoría de los estudios la coincidencia histórica automática es menos eficiente que la
coincidencia manual.
Ing. Hermas Herrera Callejas Página: 5 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

Aunque el objetivo primario del ajuste histórico es sintonizar el modelo, los

resultados del trabajo pueden contribuir a uno o más objetivos secundarios. Por
ejemplo, la coincidencia histórica contribuirá a un entendimiento de la distribución real
y movimiento de fluidos del reservorio y proporcionará una base racional para la
descripción del reservorio, incluyendo el petróleo y gas in situ, en partes del reservorio
donde no hay datos. Un modelo bien emparejado que se mantiene actualizado puede
también ser una excelente ayuda en la vigilancia del reservorio.

1.4.7 Predicción del Comportamiento y Análisis de Resultados

Una vez que se ha obtenido una coincidencia aceptable con los datos históricos, el
modelo puede ser usado para predecir el comportamiento futuro del campo y alcanzar
los objetivos establecidos para el estudio. Los tipos de predicción de rendimiento que
pueden generarse por una corrida incluyen:
1. Caudales de producción de petróleo
2. Rendimiento WOR (Relación Agua/Petróleo) y GOR (Relación Gas/Petróleo)
3. Requerimientos de pozos y de pozos reacondicionados
4. Rendimiento de la presión del reservorio
5. Posición de los frentes de los fluidos
6. Eficiencia de recuperación por área
7. Información general concerniente a requerimientos de facilidades (Ej. Medición de
gas liberado, manejo del agua, y separación a baja presión)
8. Estimación de recuperación final.
Uno de los aspectos más dificultosos para hacer predicciones es la evaluación de
los resultados de las corridas del computador. Debido a que los simuladores de
reservorios pueden generar cientos de miles de líneas de salida, debe tenerse cuidado
en concentrarse en los resultados que son necesarios para alcanzar los objetivos del
estudio de simulación.
La exactitud de las predicciones de rendimiento obviamente dependerá de las
características del modelo y de cuan exacta y completa sea la descripción del
reservorio. Es importante, por tanto, dedicar cierto tiempo estimando la calidad de la
simulación para determinar si es el adecuado para su uso.

1.4.8 Emisión de Reportes

El paso final en un estudio de simulación es la compaginación de los resultados

y conclusiones en un reporte claro y conciso. El formato usado para los reportes puede
variar desde un breve memorando entre oficinas para un estudio pequeño hasta un
reporte completo multivolumen con despliegues a color para un estudio completo del
campo. Independiente del estilo y tamaño del reporte, este debe definir los objetivos
del estudio, describir el modelo usado, y presentar los resultados y conclusiones en un
contexto apropiado al estudio.

1.5 Conceptos de Modelaje

Es necesario entender algunos conceptos básicos del modelaje de reservorios

para seleccionar el simulador apropiado, escogiendo entre las opciones disponibles y
Ing. Hermas Herrera Callejas Página: 6 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

controlando la calidad de los resultados de la simulación. Los conceptos involucrados

son simples y fáciles de entender aunque su representación matemática pueda ser
compleja.

1.6 Concepto de Bloques-Malla e Intervalos de Tiempo

En general, las ecuaciones diferenciales parciales que describen el flujo de fluidos

en el reservorio no pueden resolverse analíticamente; sin embargo pueden ser
resueltas numéricamente, reemplazando las ecuaciones diferenciales por ecuaciones
de diferencia finita. Implícito en una ecuación de diferencia finita es la discretización,
la subdivisión de la distancia y el tiempo en incrementos definidos y específicos. En
otras palabras, para usar ecuaciones de diferencia finita es necesario tratar el
reservorio como si estuviera compuesto de elementos de volumen discretos y calcular
los cambios dentro de cada elemento de volumen sobre cada uno de los muchos
intervalos de tiempo discretos. Los elementos conceptuales de volumen de un
reservorio son frecuentemente referidos como bloques-malla; los intervalos de tiempo
son llamados time steps.

1 2 3 4 5 6
Fig. 1.2 Un tanque de pozo agitado análogo de un simulador con sistema de malla

Aunque la subdivisión de un reservorio es una abstracción, es cualitativamente

correcto visualizar los bloques-malla como pozos agitados en tanques con lados
permeables. Para extender esta analogía, visualicemos los contenidos del bloque-
malla como uniformemente distribuidos dentro del bloque y los caudales a los que
fluyen los fluidos hacia dentro o fuera determinados por las permeabilidades de los
lados del bloque y las diferencias de presión entre los bloques-malla adyacentes. En
esencia el problema matemático se reduce a un cálculo del flujo entre los bloques-
malla adyacentes.
La Fig. 1.2 ilustra la analogía del tanque agitado para un modelo en el que el flujo
del fluido está restringido a una dimensión (1D). Observe que en esta analogía el
caudal y/o la composición del fluido que fluye de un tanque a otro está determinado
por el contenido del tanque superior.
Como implica la analogía del pozo con tanque agitado, las propiedades dentro
del bloque-malla no varían con la ubicación dentro del bloque, las ubicaciones dentro
de un bloque no están definidas. Por ejemplo, en cualquier momento, un bloque solo
tiene un valor de saturación de cada fase y cualquier propiedad que es dependiente
Ing. Hermas Herrera Callejas Página: 7 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

de la saturación (presión capilar y permeabilidad relativa). Para representar

variaciones en las propiedades del reservorio, las propiedades de los bloques-malla
difieren de un bloque a otro. Hay cambios abruptos o saltos en las propiedades de un
bloque-malla a otro. El modo abrupto con el que cada propiedad cambia entre los
bloques vecinos es una función del tamaño del bloque-malla.

(a) (b) (c)

(d) (e)

Δx
o Δy

(f)
Fig 1.3 - Modelos típicos usados en
Simulación de reservorios
(a) De Tanque, (b) 1D, (c) 1D radial (g)
(d) Corte transversal, (e) Areal
(f) Corte transversal radial, y (g) 3D
(g)

La precisión con la que un reservorio puede ser descrito en un modelo y la

exactitud con la que el flujo del fluido del reservorio puede calcularse dependerá del
número de bloques-malla usados en el modelo. En la práctica, el número de bloques
estará limitado principalmente por el costo de los cálculos y el tiempo disponible para
preparar los datos de entrada y para interpretar los resultados. Como una
consecuencia, el tamaño y la complejidad del modelo de reservorio deben ser
considerados cuidadosamente. El modelo debería contener suficientes bloques (y
dimensiones) para simular el reservorio y su comportamiento adecuadamente, pero
debería ser tan pequeño y simple como sea posible.
La Fig. 1.3 ilustra algunos de los muchos tipos de modelos que son necesarios
en la simulación del reservorio.
Ing. Hermas Herrera Callejas Página: 8 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

1.7 Consecuencias de la Discretización.

Para ilustrar las consecuencias de usar bloques-malla para modelar un

reservorio, veamos un modelo 1D de cuatro bloques-malla de un reservorio hipotético
de dos pozos que está siendo inundado. La Fig. 1.4 muestra la vista de plano del
reservorio, un boceto del sistema de malla para el modelo y diagramas de la
distribución de saturación de agua en el reservorio y en el modelo relativamente
temprano en la inundación. Note que en el reservorio la saturación del agua es una
función continua de la distancia, pero en el modelo hay un salto abrupto en la
saturación del agua entre los bloques 1 y 2 y de nuevo entre los bloques 2 y 3. Un
diagrama similar de presión en el modelo demostraría una distribución escalonada,
mientras un diagrama de la presión real del reservorio sería una curva continua.

Saturac agua Reservorio Modelo de simulador

75 75

50 50

25 25

0 0
Distancia 1 2 3 4
(a) (b) Nro. de bloque
Fig. 1.4 - Modelo de empuje hidráulico de 4 bloques-malla: (a) reservorio hipotético y su distribución de saturación
de agua a un cierto tiempo en el influjo de agua y (b) modelo de 4 bloques-malla y distribución simulada de
saturación de agua.

1.7.1 Ponderación de la Movilidad

Para calcular el flujo de agua y petróleo de bloque a bloque en nuestro ejemplo

de inundación, debemos asignar movilidad de agua y petróleo aplicables al flujo a
través de los límites hipotéticos entre bloques. Pero la movilidad de una fase es una
función de la saturación, y hemos mostrado que la saturación en dos bloques
adyacentes pueden diferir significativamente. Esto lleva a preguntarnos qué saturación
usar al definir movilidades para el flujo entre los dos bloques. Ningún valor único de
saturación de agua, por tanto ningún valor único de movilidad de agua o petróleo,
puede describir adecuadamente el comportamiento del flujo sobre el amplio rango de
saturaciones que puede existir en las regiones del reservorio representado por dos
bloques-malla adyacentes. Pero en un simulador, debemos usar un valor único de
movilidad cada vez que calculamos el caudal de flujo de una fase entre dos bloques-
malla.

Ing. Hermas Herrera Callejas Página: 9 de 21

Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

Una decisión racional puede ser usar un promedio de las movilidades de los dos
bloques. Actualmente hay cuatro métodos a considerar.
(1) Ponderación aguas arriba. La movilidad a usarse en el cálculo del flujo entre dos
bloques adyacentes es determinada enteramente por la movilidad en el bloque
aguas arriba (bloque 1 para el flujo entre los bloques 1 y 2 del modelo mencionado)
(2) Ponderación aguas abajo. En la cual la movilidad es determinada enteramente por
la movilidad en el bloque aguas abajo.
(3) Ponderación mixta. En la cual alguna combinación de movilidades en bloques
aguas arriba y aguas abajo son usadas.
(4) Métodos de extrapolación o interpolación. Tal como ponderación aguas arriba en
dos puntos, las movilidades en uno o más bloques-malla vecinos se usan para
lograr un estimado de movilidades en la interfase de dos bloques-malla.

1.7.2 Dispersión Numérica

La dispersión numérica es un artificio de las actuales técnicas de análisis

numérico que pueden causar severas distorsiones en simulaciones de procesos en los
cuales ocurren relativamente rápidos cambios de saturación. Tales cambios son
comunes en muchas pérdidas de desplazamiento de reservorios. Para ilustrar el
problema nuevamente usaremos el modelo de flujo de agua de cuatro bloques y
consideraremos flujo entre los bloques 2 y 3 y entre bloques 3 y 4 bajo las condiciones
representadas por la Fig. 1.4. Los bloques 3 y 4 aun no han sido invadidos por agua,
y en esos bloques el agua tendrá poca o ninguna movilidad. La saturación de agua en
el bloque 2 está por encima de la baja saturación inicial y el agua tendrá considerable
movilidad en este bloque. Si se usa ponderación total de movilidad aguas arriba
(método 1), la movilidad en el bloque 2 definirá la movilidad a usarse en el cálculo del
flujo del bloque 2 al bloque 3. Por tanto, en el próximo time step, algo de agua fluirá al
bloque 3 porque el agua es móvil en el bloque 2. El flujo de agua al bloque 3
incrementará la saturación de agua y causará que el agua llegue a ser móvil en este
bloque. Como resultado, el agua fluirá del bloque 3 al bloque 4 en el time step
subsiguiente. En otras palabras, en esta simulación al menos algo de agua se moverá
hacia delante de un bloque-malla en cada time step, independiente del tamaño del time
step. En el modelo de 4 bloques ilustrado acá, cálculos con ponderación aguas arriba
moverá algo de agua hacia el bloque-malla del pozo productor dentro de 4 time steps.

1.7.3 Efectos de Orientación de la Malla

En modelos multidimensionales, la dispersión numérica tiende a un fenómeno

interesante y algunas veces incómodo en que el rendimiento calculado es influenciado
por la orientación de la malla relativa a las ubicaciones de los pozos de inyección y
producción. El efecto de la orientación de la malla usualmente no es importante
excepto en simulaciones en que la fase desplazante es mucho más movible que la
fase desplazada (como en flujos de vapor de petróleo pesado)
1

2 3 4 5
A
Ing. Hermas Herrera Callejas Página: 10 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

6 7

4 5

1 3
B 2

Fig. 1.5 – Vías de flujo para flujo paralelo y diagonal en un sistema de malla rectangular.

La Fig. 1.5 ilustra el problema. Es un bosquejo del sistema de malla de un modelo

para simular inyección de vapor en un reservorio de petróleo viscoso. Esta parte del
modelo contiene un pozo de producción y dos pozos de inyección de vapor. Los pozos
A y B.
En el simulador, el vapor del pozo A se moverá en vía directa al productor. Sin
embargo, el vapor del pozo B debe zigzaguear en su ruta al productor. No-solo la vía
de flujo del pozo B es más larga, sino que el vapor del pozo B verá más petróleo que
el vapor del pozo A. Por tanto, en adición a llegar más tarde, el vapor del pozo B
parecerá barrer más eficientemente que el vapor del pozo A. Si la malla se rotara 45
grados (0.8 radianes), los rendimientos calculados para los dos pozos serían
reversados.
Las presiones así como las saturaciones son distorsionadas por la orientación de
la malla. El uso de un gran número de bloques-malla, que normalmente controlaría la
dispersión, reducirá pero no eliminará el efecto de la orientación de la malla.
Afortunadamente, el efecto no es importante en la mayoría de los simuladores de
desplazamiento agua/petróleo. La orientación de la malla puede ser importante en
algunas simulaciones gas / petróleo, y a veces puede causar un problema serio en la
simulación del flujo de vapor o desplazamientos de gas soluble.

1.8 Funciones Explícita e Implícita.

Hemos discutido algunas consecuencias de las discontinuidades en saturación y

presión que existen en las interfases entre bloques-malla en un modelo de reservorio.
Similarmente, las saturaciones y presiones son discontinuas en el tiempo: las
condiciones del reservorio son definidas solo al final de cada time step. Para ilustrar,
la Fig. 1.6 compara la propagación de saturación de agua real con la calculada en el
tiempo en una región del reservorio hipotético representado por el segundo bloque-
malla del modelo de flujo de agua de 4 bloques.

Saturac de agua

75 75

50 50

25 25

Ing. Hermas Herrera Callejas Página: 11 de 21

Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

0 0
Tiempo 1 2 3 4 5
(a) (b) Time step
Fig. 1.6 - Comparación de saturación vs tiempo (a) en el reservorio y (b) en un simulador

Durante el time step 3, la saturación de agua en el segundo bloque-malla se

incrementó en cerca del 10 %. Las propiedades, que son función de la saturación, tales
como la movilidad y la presión capilar, habrán cambiado significativamente sobre este
time step. En el reservorio real, durante el tiempo modelado por el time step 3, estas
funciones habrán asumido muchos valores. Obviamente, ningún valor de estas
funciones es correcto por todo el tiempo representado por el time step 3. En el
simulador, sin embargo, un solo valor para cada variable dependiente de la saturación
debe usarse para describir condiciones por todo el time step.
Como con la ponderación de la movilidad, no hay un método mejor para
seleccionar valores a utilizarse durante un time step para las funciones dependientes
de la saturación que cambian con el tiempo durante una simulación. Pueden usarse
varios procedimientos, pero todos son variaciones de los tres siguientes:
1. Procedimientos Explícitos usan valores de saturación conocidos al comienzo de
un time step. Estas saturaciones se conocen de los cálculos del time step previo (Ej.
El fin del time step 2 es el comienzo del time step 3).
2. Procedimientos Implícitos usan la movilidad y la presión capilar calculada como
funciones de la saturación al final de time step. Los valores no son conocidos hasta
que los cálculos para el time step hayan sido completados.
3. Procedimientos Semi-Implícitos usan estimados de movilidad y presión capilar
desarrollados asumiendo que las funciones serán de línea recta de la saturación
durante un time step. Las saturaciones no son resueltas independientemente, por tanto
el procedimiento es aún parcialmente implícito.

1.8.1 Procedimientos Explícitos, Semi-Implícitos e Implícitos.

Para ilustrar el uso de los tres procedimientos, calculemos el cambio en

saturación de agua en el bloque-malla 1 (el bloque del pozo de inyección) que ocurre
en el modelo de flujo de agua de 4 bloques-malla durante los primeros time steps del
comienzo de la inyección de agua. En la discusión el superíndice n se refiere al time
step inmediatamente precedente y n+1 se refiere al time step que se está procesando.
Inicialmente el flujo hacia fuera del bloque será petróleo. Pronto, sin embargo,
tanto agua como petróleo fluirá fuera del bloque 1 con la WOR creciente a medida que
procede el desplazamiento. El único flujo en el bloque-malla 1 es agua del pozo de
inyección. Como consecuencia, el cambio de saturación en el bloque durante un time
step será completamente determinado por la fracción de petróleo que fluye del bloque-
malla 1 al bloque-malla 2 durante el time step. En forma de ecuación, el cambio de
saturación en el bloque 1 durante el time step n+1 es:
So  ( Son 1  Son )  ( S wn 1  S wn )   qtf o ,...... 1.1
Donde: S n , S n 1 = Saturación al comienzo y fin del time step, respectivamente
q = Caudal de flujo de los fluidos fuera del bloque 1
t = Duración del time step
Ing. Hermas Herrera Callejas Página: 12 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

fo = Fracción de petróleo en los fluidos que fluyen fuera del bloque 1

Note que q se expresa en volumen poral por unidad de tiempo y que el signo
negativo implica flujo de salida del bloque 1.
Los procedimientos explícitos, semi-implícitos e implícitos difieren en la manera
en que ellos evalúan f o . Para simplificar la ilustración y facilitar los cálculos, se hacen
las siguientes suposiciones:
1. La ponderación de la movilidad es completamente aguas arriba; por tanto, todas
las funciones de saturación son aquellas del bloque-malla 1
2. La fracción de flujo de petróleo puede expresarse analíticamente como función de
la saturación de agua.

1.8.1.1 Procedimiento Explícito

En un procedimiento explícito, f o para el cálculo del flujo durante un time step es

función de las saturaciones en el bloque al comienzo del time step
f o  f on ................................ 1.2
La ecuación 1.1 llega a ser
( S on 1  S on )  qtf on ..................................................... 1.3
Donde f on se conoce del time step previo. Por tanto todos los términos, excepto S on 1 son
conocidos, así que S on 1 puede calcularse directamente.

1.8.1.2 Procedimiento Semi-Implícito

En un procedimiento semi-implícito f o se estima por extrapolación de los valores

conocidos, f on , al final del time step:
f
f o  f on  o ( S wn 1  S wn ) ............................................ 1.4
S w
La ecuación 1.1 llega a ser:
 f 
( Son 1  Son )  qt  f on  o ( S wn 1  S wn ) ...................... 1.5
 S w 
La pregunta aquí es donde tomar la pendiente f o / S w . El enfoque usado en la
formulación más simple del procedimiento consiste en estimar el cambio de saturación
que ocurrirá durante el time step y usar el valor promedio de la pendiente sobre ese
rango de saturación. (Mencionaremos una manera de estimar el cambio de saturación
en la sección 1.8.1.4). De esta manera la pendiente puede estimarse al comienzo del
time step sin información del valor eventual de S on 1 . El procedimiento semi-implícito
introduce un término adicional: S wn 1 ; pero debido a que el coeficiente f on + f o / S w es
definido, no se añade mucho trabajo al cálculo.

1.8.1.3 Procedimiento Fully Implicit

Ing. Hermas Herrera Callejas Página: 13 de 21

Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

En un procedimiento fully implícit, la fracción de petróleo de los fluidos que salen

del bloque se evalúa al final del time step.
f o  f on 1 .......................................................... 1.6
La ecuación 1.1 se convierte en:
n 1
S w  S wn  qtf on 1 ............................................ 1.7
En esta instancia, f on 1 es una función de la saturación al final del time step n+1 y
no puede calcularse independientemente de S wn 1 .

1.8.1.4 Comparación de Soluciones

Se deben especificar algunas condiciones adicionales antes de que pueda

calcularse el rendimiento del bloque-malla 1.
1. El Volumen Poral del bloque es 10.000 bbl (1.590 m 3 )
2. El caudal de inyección de agua es 100 BD (16 m 3 /día). Basado en el volumen poral,
este es un caudal de entrada alto (0.01 VP/D) y lleva ventaja al método de solución.
1
3. f o 
S  0.25 2.5
1 ( w )
0.75  S w
4. La saturación de agua inicial es 0.25 VP.
La saturación de agua en el bloque-malla 1 después de 50 días de inyección se
calcula por cada uno de los tres procedimientos con time steps de 0.01, 0.1, 1, 10, 25,
y 50 días.
Los procedimientos de cálculo son como sigue:
1. Procedimiento explícito. f on se calcula con el valor conocido de S wn y la ecuación 1.3
se resuelve directamente para S wn 1
2. Procedimiento semi-implícito. Al comienzo de cada time step el cambio de
saturación esperado durante el time step se estima con el procedimiento explícito. Este
estimado se usa en el cálculo del promedio de la pendiente de f o / S w . Entonces la
ecuación 1.5 se resuelve para S wn 1 .
3. Procedimiento implícito. Un método iterativo se usa para resolver la ecuación 1.7.
Primero se estima S wn 1 con el procedimiento explícito. Luego se calcula f on 1 a S wn 1 . Se
usa f on 1 para calcular un nuevo valor estimado de S wn 1 . El procedimiento se repite
hasta que S wn 1 no cambie entre iteraciones, o sea hasta que la solución converja.
TABLA 1.1 – SATURACIÓN DE AGUA EN EL BLOQUE-MALLA 1
DESPUES DE 50 DIAS DE INYECCIÓN
Tamaño del Saturación de agua
Time step (%VP) .
(Dias) Explícito Semi-implícito Implícito
0.01 56.57180 56.56648 56.56647
0.1 56.58438 56.55182 56.55182
1 56.73216 56.40623 56.40656

Ing. Hermas Herrera Callejas Página: 14 de 21

Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

10 56.33037 54.98752 55.03765

25 62.50000 52.63948 53.00287
50 75.00000 50.00000 50.00000
Los resultados de los cálculos se resumen en la Tabla 1.1. Cuando se usaron
time steps cortos, los tres procedimientos dieron respuestas casi idénticas. Con time
steps mayores, el procedimiento explícito sobreestimó el cambio de saturación en el
bloque-malla. Los procedimientos semi-implícito e implícito dieron respuestas
comparables para todos los tamaños de time steps, pero aún estos métodos dieron
resultados que variaron con el tamaño del time step usado en el cálculo. (Note que la
mejor medida de la longitud del time step es VP’s de flujo en vez de días. En este
problema, 10 días es un time step grande). Normalmente se impone un número de
restricciones durante un cálculo para controlar el tamaño de un time step.
Los datos presentados en la Tabla 1.1 no hacen el punto pero en un cálculo semi-
implícito la saturación cambia, el lugar donde se toman las pendientes de las funciones
de saturación pueden afectar la calidad de la respuesta.

1.9 Tratamiento de la Saturación Vertical y Distribuciones de Presión.

El bloque-malla tiene espesor específico en adición a sus dimensiones areales;

pero consistente con la analogía de tanques-batidos. Las variaciones verticales en la
distribución del fluido, presión, y propiedades de la roca no pueden ser definidos
directamente dentro de un bloque-malla único. La Fig. 1.7 ilustra el problema para un
bloque-malla que representa parte de la transición agua petróleo efectuada en un
modelo areal bidimensional (2D). En el reservorio, la saturación de petróleo es alta en
la parte superior de la zona y la movilidad del petróleo será alta. Similarmente la
movilidad de agua será alta en la parte baja de la zona. Por contraste, la saturación de
agua en el bloque-malla modelo se define por un solo valor promedio que no distingue
entre las partes alta y baja de la zona. Para modelar el comportamiento de flujo con
razonable precisión, el simulador debe asignar al bloque malla funciones de movilidad
que sucedan propiamente para las movilidades reales de los fluidos asociados con
esta distribución específica de fluidos a esta saturación específica.
Segmento del reservorio Bloque-malla del simulador
Profundidad

Valores simples de saturación

de agua igual al promedio
volumétrico de la curva de la
izquierda
(a) Saturación de Agua (b) Saturación de Agua
Fig. 1.7 – Distribución de saturación vertical en un segmento del reservorio es modelado como un solo
valor en un bloque-malla

Ing. Hermas Herrera Callejas Página: 15 de 21

Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

Por supuesto, la zona podría ser modelada exactamente con un modelo

tridimensional (3D) teniendo bloques-malla de suficiente espesor para definir la
distribución de saturación vertical. Hay una alternativa, sin embargo, en muchos
modelos de reservorio, un bloque grueso puede hacerse comportar como si fuera
construido de muchos bloques delgados usando un juego revisado de funciones de
saturación (permeabilidad relativa y presión capilar). Como se discutió anteriormente,
las funciones son derivadas en una manera que cuenta para el flujo a través de un
número de capas hipotéticas, cada cual con su propia saturación de fluido, espesor y
permeabilidad. Se calcula la contribución de cada capa hipotética a la movilidad de
petróleo y agua. La movilidad total es la suma de las movilidades de todas las capas.
La movilidad total resultante para cada fase es una movilidad ponderada promedio
permeabilidad y espesor para una distribución específica de saturación y un valor de
saturación promedio que definen condiciones en una parte del reservorio a un tiempo
específico. Repitiendo el cálculo para el rango completo de distribución de saturación
esperado en una parte específica del reservorio producirá un juego de funciones de
movilidad. El juego de valores así derivados se llama una ‘seudo función’.
Para que las seudo funciones sean satisfactorias es esencial usar la distribución
de saturaciones correcta al momento de derivarlas. Hay dos juegos de condiciones a
considerar:
1. En el caso más general, las condiciones son dinámicas y las fuerzas de viscosidad,
gravedad y capilaridad todas influyen en la distribución de la saturación vertical.
2. En muchas instancias las seudo funciones satisfactorias pueden ser desarrolladas
asumiendo condiciones de equilibrio, en el cual la gravedad y capilaridad controla
la distribución de saturación vertical
Para desarrollar seudo funciones para condiciones dinámicas, la distribución
vertical de saturación y permeabilidad deben ser modelados. El tipo de modelo más a
menudo usado con estos propósitos es un modelo de corte transversal 2D con una
adecuada malla vertical definida. El modelo de corte transversal debe ser diseñado
para simular las principales direcciones del flujo, por tanto pueden necesitarse varios
modelos diferentes. La validez de un juego particular de seudo funciones dinámicas
pueden establecerse usando las funciones en un modelo 1D para simular el
rendimiento del modelo de corte transversal.
Si el reservorio se comporta como si las fuerzas capilar y de gravedad estuvieran
en equilibrio, las seudo funciones de equilibrio vertical pueden calcularse
analíticamente por integración vertical y las curvas de permeabilidad relativa y presión
capilar para diferentes posiciones de la zona de transición como si se moviera a través
del bloque-malla.

1.9.1 Funciones de Pozos

Es claro que la falta de definición vertical dentro de un solo bloque-malla hace

imposible simular directamente el rendimiento de un pozo en un modelo areal 2D en
el cual el pozo es representado por un bloque-malla. Un solo bloque-malla obviamente
no puede modelar una terminación parcial y una sola capa de bloques no puede
modelar conificación de agua o gas. Estas limitaciones pueden resolverse por el uso
de ‘funciones de pozos’ que son seudo funciones derivadas de modelos de pozos
individuales finamente enmallados o tipos de pozos. Las funciones de pozos
Ing. Hermas Herrera Callejas Página: 16 de 21
Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

relacionan movilidades de fluidos cerca de la boca de pozo a las saturaciones de

bloques-malla. El modelo usado para derivar funciones de pozos debe representar
apropiadamente el intervalo de terminación, propiedades del reservorio cerca de boca
de pozo, y la mecánica de flujo multifase en el pozo (ya sea por cuspidificación,
conificación o flujo estratificado). En otras palabras, las seudo funciones asignadas a
un pozo deben representar el mecanismo de producción y la distribución geométrica
de saturaciones que pueden esperarse en la vecindad del pozo. Esto significa que
puede necesitarse un juego separado de seudo funciones para cada terminación.

1.9.2 Re-Solución de Gas

Dentro de un bloque-malla, la relación gas/petróleo GOR se expresa

normalmente como función de la presión. Esto vale para presión constante o
declinante, porque en una escala de tiempo del reservorio el establecimiento de un
equilibrio en la distribución de saturación durante la evolución de gas es efectivamente
instantáneo. Si se establece una fase de gas libre, sin embargo, algún gas segregará
y no permanecerá en contacto directo en el petróleo en el cual estaba disuelto. Si la
presión entonces sube, el equilibrio no necesariamente será restablecido
‘instantáneamente’. Porque la distribución areal y vertical del gas libre dentro de un
solo bloque-malla no está definido, un modelo de reservorio no puede rigurosamente
calcular la cantidad de gas que se re-disolvería. Si la re-solución se piensa que es
importante, deben usarse ya sea bloques malla más delgados o algunas veces
arealmente más pequeños o bien el proceso debe primeramente ser estudiado en un
modelo externo de malla fino y los resultados impuestos al simulador mediante
relaciones modificadas entre la relación GOR y la presión.

1.10 Administración de Pozos.

En adición a las funciones de pozos mencionadas anteriormente, un simulador

de reservorio debe contener un programa que traduce datos de producción / inyección
a las condiciones operativas deseadas del campo y restricciones en controles para el
modelo de reservorio. Tales programas se llaman rutinas de administración de pozos.
Una rutina de administración de pozos establece relaciones de pozos o presiones,
implementa políticas de operación, y satisface restricciones operativas a los niveles de
intervalo de producción, de pozo, grupo de pozos, reservorio, y campo.
En su forma más simple, una rutina de administración de pozos asigna caudales
específicos o presiones a bloques de producción individuales en el simulador a tiempos
especificados. Esta forma simplificada es a menudo referida como una ‘rutina de
caudal’. Rutinas más sofisticadas pueden también hacer y ejecutar decisiones que
simulan muchas de las decisiones operativas requeridas para administrar el reservorio
real. Por ejemplo, una rutina de administración de pozos puede perforar,
reacondicionar, o re-terminar pozos, calcular la hidráulica de pozos y líneas de flujo,
instalar ascenso artificial, controlar caudales de gas o agua, mantener objetivos de
producción, y dirigir re-inyección de gas o agua. La rutina puede también imponer
restricciones de producción que en la práctica puede resultar de tales controles
externos como contratos de producción, guias de operación, o reglas estatutarias.

Ing. Hermas Herrera Callejas Página: 17 de 21

Modelos de Simulación y Laboratorio Capítulo 1 – Fundamento Conceptual

1.11 Métodos de Solución.

La solución de ecuaciones de flujo consume gran parte del esfuerzo de

computación en una simulación típica. Como consecuencia, la selección de un método
efectivo de solución de la ecuación para un problema particular puede ser un paso
importante en el manejo del costo y control de la dificultad del estudio. La mayoría de
los simuladores proporcionan dos o más métodos para resolver las ecuaciones y
permiten al usuario seleccionar el método más apropiado para el problema.
Los métodos de solución son directos e iterativos. En un método directo, las
incógnitas son eliminadas una a una hasta que las ecuaciones son resueltas. Excepto
por posibles errores de redondeo, las soluciones son exactas. La base de un método
iterativo es el desarrollo de una solución ‘aproximada’ al sistema de ecuaciones. La
aproximación es reemplazada sistemáticamente hasta que las respuestas ‘converjan’
dentro de una tolerancia determinada de las respuestas correctas.
Alguna forma de solución directa – el método que normalmente utilizamos para
resolver pequeños sistemas de ecuaciones simultáneas a mano- está incluido en la
mayoría de los simuladores. El método es usualmente una opción apropiada para
modelos de tamaño pequeño a intermedio. Las respuestas son exactas y no se
requiere intervención del usuario. Con modelos que tienen un gran número de bloques-
malla, sin embargo, la solución directa puede no ser práctica porque el número de
cálculos y los requerimientos de almacenamiento durante los cálculos puede llegar a
ser excesivos. El significado de ‘grande’ aquí es influenciado por el tipo de problema
que se está resolviendo y por las características del software y hardware de
computación.
La mayoría de los simuladores también contienen uno o más métodos iterativos
para resolver las ecuaciones. El esfuerzo de cómputo para una iteración simple de un
método iterativo puede ser solo una pequeña fracción de aquellos requeridos para una
solución directa. Pero un factor adicional debe considerarse: el número de iteraciones
requeridas antes de que se alcance el criterio de convergencia. Sin usar el método
iterativo en el problema específico que se está resolviendo, puede no ser posible
estimar el número de iteraciones necesarias para una solución. Una complicación
adicional de algunos métodos iterativos es la necesidad de seleccionar parámetros de
iteración que son usados para acelerar la convergencia. Los valores de los parámetros
usualmente son influenciados por las características del problema particular que se
está resolviendo.
Para cualquier estudio particular la selección del método de solución puede no
ser clara. Puede requerirse alguna experimentación con el modelo de reservorio antes
que el método más apropiado sea seleccionado.

Ing. Hermas Herrera Callejas Página: 18 de 21