UNIVERSIDAD AUTONOMA GABRIEL RENE

MORENO

SIMULACION DE RESERVORIO
Principios Básicos y Conceptos

Autor: Msc. Ing. José Luis Rivero Sandoval

Santa Cruz-Bolivia
 Introducción

Simulación de Yacimientos es un texto universitario dirigido a estudiantes de pregrado y

postgrado en Ingeniería de Petróleos los cuales reúne las experiencias y los conocimientos
adquirido a lo largo de los años en la carrera de docente y de sus estudios de postgrado. El texto,
dirigido a estudiantes de pregrado y posgrado en Ingeniería de Petróleos, consta de quince
unidades que tratan de resaltar los aspectos más importantes de la Simulación de Yacimientos,
una ciencia que se encuentra en su mayor auge, dada la necesidad de describir lo más
exactamente posible el comportamiento de un yacimiento de hidrocarburos.

En el presente libro se presentan los conceptos fundamentales de matemáticas, física e

ingeniería de yacimientos requeridos para desarrollar un simulador de un yacimiento al igual
que la metodología para llevar a cabo un estudio de simulación de yacimientos numérica y
analíticamente. El texto menciona algunos conceptos modernos considerados en la actualidad
en el área.

Las primeras cinco unidades dan la clasificación de los simuladores, tipo de modelos,
enmallado, errores en la información requerida para un estudio de simulación, ajuste histórico
y efectos de orientación de la malla. Las unidades seis, siete y ocho dan los conceptos
fundamentales y las bases matemáticas de la ecuación de difusividad, clasificación de las
ecuaciones en diferencias y problemas de valores de frontera. De la unidad nueve a la quince,
se estudia el tratamiento numérico de las ecuaciones de flujo al igual que los esquemas que se
han introducido en la literatura para resolver más eficientemente los problemas de flujo en el
medio poroso. Allí se consideran los métodos para resolver las ecuaciones algebraicas
resultantes de la aplicación de las diferencias finitas.

En este libro se presentaran las formas de

Simulación de Reservorio
2
 Contenido

CAPITULO I Introducción al Método Numérico

1.1.- Introducción al Método Numérico

1.2. Definición
1.3. Breve Historia de la Simulación
1.4.- Objetivos de la Simulación
1.4.1. Utilidad de la Simulación
1.4.2. Ajuste de Simulador con la Historia del Yacimiento
1.4.3. Resultados de una Simulación
1.5. Etapas para Desarrollar un Modelo
1.6. Como Trabaja un Modelo
1.7.- Consenso en el modelaje.
1.8.- Manejo de Estudios de Simulación.

CAPITULO II Proceso del Modelaje

2.1.- Introducción
2.2.- Análisis básicos del Reservorio.
2.3.- Pre-requisitos.
2.4.- Modelaje computacional.
2.5.- Elementos principales en el estudio de Simulación de Reservorios
2.6.- Escala y muestreo del Reservorio
2.4.- Modelaje computacional.
2.5.- Elementos principales en el estudio de Simulación de Reservorios
2.6.- Escala y muestreo del Reservorio
2.7.- Escalas integradas – Unidades de Flujo
2.8.- Estructura del Reservorio
2.8.1.- Giga Escala.
2.8.2.- Mega Escala.

CAPITULO III Información Requerida del Proceso de Modelaje

3.1 Introducción
3.2. Descripción Física Del Yacimiento
3.3. Mecanismos De Desplazamiento
3.4. Propiedades Petrofísicas
3.4.1.- Porosidad
3.4.1.1. La Porosidad primaria
3.4.1.2. La Porosidad secundaria
3.4.2. Permeabilidades
3.4.3 Saturación
3.3.4. Permeabilidad relativa
3.3.5. Presión Capilar

Simulación de Reservorio
3
3.3.6.-Propiedades de los fluidos.
3.3.7.- Tipos de fluidos.
3.3.8.- Modelaje del fluido
3.3.9.- Muestreo de fluidos
3.3.9.1. Toma de Muestras del Reservorio
3.3.9.2.- Recomendaciones para el muestreo de Pozo
3.3.9.3-Tipos de Muestreo
3.3.9.4.-Muestra de fondo
3.3.9.5.- Muestra de Superficie
3.3.9.6- Ventaja y Desventaja en diferente tipos de muestreo
3.3.9.7.-Muestreo en Superficie
3.3.9.8.- Muestreo en Fondo
3.3.10.- Determinación de los Volúmenes In –situ
3.3.10.1.- Ecuación de Balance de Materia
3.3.10.2.-.-Determinación del Mecanismo de Producción
3.3.10.3..-Linealizacion de la Ecuación de Balance de Materiales
3.3.10.4.-.-Metodo Propuesto
3.3.10.5..- Método Presentado por Roberto Aguilera SPE100451
3.3.10.6.- Determinación del Volumen In-Situ para Reservorios de Petróleo

CAPITULO IV Clasificación de los Simuladores

4.1.- Introducción
4.2.- Tipos de Simuladores
4.2.1.- Modelo tanque o dimensión cero
4.2.2.- Modelo 2D
4.2.2.- Modelo 3D
4.2.3.- Enmallado
4.3.- Tipo de Yacimientos
4.3.1.- Nivel de Simulación
4.3.2.- Tipos de Yacimientos
4.3.2.1.- Simuladores de Gas
4.3.2.2.- Simuladores de Gas Condensado
4.3.2.3.- Simuladores de Petróleo Negro
4.3.2.4.- Simuladores Geotérmicos
4.3.2.5.- Simuladores de recuperación química
4.3.2.6.- Simuladores de recuperación con miscibles
4.3.2.7.- Simuladores de recuperación Térmica
4.4.- Comparaciones de Modelamiento
4.4.1.- Métodos Analógicos
4.4.2.- Métodos Experimentales
4.4.3.- Métodos Matemáticos
4.5.- Exactitud de las Soluciones
4.5.1.- Error de Redondeo
4.5.2.- Error de balance de materia

Simulación de Reservorio
4
 4.5.3.- Errores no lineales
4.5.4.- Error de Inestabilidad
4.5.5.- Errores de Truncamiento
4.5.6.- Errores de Dispersión numérica

Simulación de Reservorio
5
CAPITULO I Introducción al Método Numérico

1.1. Antecedentes

La explosión electrónica en las últimas dos décadas ha transformado la simulación de

yacimientos de algo inaccesible y oculto en una herramienta muy importante que, entre muchas
cosas, le permite al ingeniero tener un mejor entendimiento de la dinámica del flujo de fluidos
en yacimientos muy complejos y las características de la dinámica del flujo de fluidos en
cercanías al pozo, la interacción del pozo (horizontal, vertical o desviado) con el yacimiento, el
modelamiento adecuado de las estructuras geológicas, fallas y acuñamiento y la complejidad de
la caracterización del yacimiento.

Los simuladores numéricos de yacimientos se usan muy ampliamente ya que permiten resolver
problemas que no se pueden resolver por otros medios. Aunque recientemente la aplicación de
métodos semianalíticos, como el de las líneas de flujo (streamlines) han tomado auge,
especialmente para simular yacimientos estratigráficamente complejos. Su versatilidad se debe
a que utilizan menor esfuerzo de cómputo, fundamentalmente radica en el desacoplamiento del
problema de flujo de fluidos de 3D a 1D.

El potencial de la simulación arrancó a finales de los 40’s. Con el compromiso de aunar

esfuerzos en el análisis numérico y el desarrollo de métodos prácticos de cálculo.

Inicialmente, los primeros simuladores fueron lo suficientemente grandes para justificar

estudios costosos y la necesidad de contar con un comportamiento detallado a mediano o largo
plazo. Esto, por supuesto, convirtió la simulación en una herramienta indispensable en el manejo
de los yacimientos.

Anteriormente para calcular la recuperación de hidrocarburos, por ejemplo, se utilizaban

métodos de balance de materia como los de Schilthuis, Tarner, Muskat, Pirson y Tracy en los
cuales se considera al yacimiento como un tanque con propiedades promedio, tanto de presión
como de propiedades petrofísicas y PVT de los fluidos. Sin embargo, esta suposición de
homogeneidad a lo largo de todo el yacimiento, aunque se ha demostrado que puede ser valida,
muchas veces no existe, por lo cual se pensó en dividir el yacimiento en una serie de bloques o
celdas, asignándole a cada una de ellas propiedades promedio y aplicar la ecuación de balance
de materia para cada bloque, acoplado a la ecuación de Darcy que es una ecuación de flujo que
permite determinar la interacción entre los bloques.

A esto, es decir, el dividir el yacimiento en una serie de bloques para su estudio se le conoce en
forma general como simulación y los aspectos nuevos que presenta es que, como puede
fácilmente suponerse, se requiere de una gran cantidad de cálculos (hay que utilizar balance de
materia en cada bloque) por lo que se hace indispensable el uso de una computadora para
llevarlos a cabo.

Simulación de Reservorio
6
Lógicamente el primer problema que surge es obtener la información necesaria para cada
bloque, lo cual se desconoce y una vez estimado continúa tendiendo cierto grado de
incertidumbre. Sin embargo, suponiendo que se pueda conseguir dicha información, sin lugar a
dudas, esta es la mejor manera de llevar a cabo el estudio de un yacimiento cuando éste no es
homogéneo. Con ello no se quiere decir que esto sea lo mejor o deba de aplicarse indistintamente
a cualquier problema, pues la experiencia ha demostrado que el método de balance de materia
simplificado, bien aplicado, en determinados casos puede proporcionar resultados acertados y
económicos.

Actualmente, la disponibilidad de computadores y el enorme progreso que ellos han venido

teniendo recientemente, hace de la simulación un instrumento práctico en la toma de decisiones
y planeación durante la vida de un yacimiento.

1.2. Definición

La simulación de yacimientos es una ciencia que combina la física, la matemática, lageología,

la ingeniería de yacimientos y la programación de computadores para desarrollar herramientas
que pronostiquen el comportamiento de los yacimientos de hidrocarburos bajo diferentes
condiciones de operación.

Simular quiere decir “dar la apariencia de”. Luego, esta ciencia es indispensable en virtud a
que se requiere obtener predicciones exactas del desarrollo de un yacimiento. Dicha necesidad
nace del hecho que un proyecto de recuperación de un campo de hidrocarburos involucra una
inversión de cientos de millones de dólares y presenta varios riesgos que están asociados con el
desarrollo seleccionado y por tanto se precisa la evaluación y minimización de dichos riesgos.
Los factores que contribuyen al riego incluyen:

• Complejidad del yacimiento debido a las propiedades de heterogeneidad y anisotropía de las

rocas,

• Variaciones regionales del flujo de fluidos y características de las curvas de permeabilidades

relativas,

• Complejidad del mecanismo de recuperación de hidrocarburos,

• Aplicabilidad de otros métodos predictivos limitados e inapropiados.

El último factor es el único controlable por el ingeniero pero requiere experticia y práctica
adecuada.

1.3. Breve Historia de la Simulación

Prácticamente la simulación de yacimientos se ha venido empleando desde los inicios de la

ingeniería petrolera, es decir, hacia 1940. Antes de 1960, los cálculos usados para predecir del
comportamiento del yacimiento pronosticar la recuperación o comparar alternativas económicas

Simulación de Reservorio
7
entre diversos métodos de recuperación consistían en su mayoría de métodos analíticos tales
como: el método de balance de materia o simulador de cero dimensiones y el método de
Buckley-Leverett o modelo de una dimensión.

El termino “simulación” se hace común a los principios de 1960, refiriéndose con él a métodos
de predicción desarrollados en programas de computadoras relativamente sofisticados. Dichos
programas representaban un mayor adelanto debido a que permitían la solución de un conjunto
de ecuaciones expresadas en diferencias finitas que describían flujo multifásico a través de un
medio poroso heterogéneo en dos y tres dimensiones. Este adelanto se hizo posible gracias a la
evolución tan rápida que tuvieron las computadoras y sobretodo el desarrollo de métodos
numéricos capaces de resolver grandes sistemas de ecuaciones en diferencias finitas.

Durante los años 60´s, los esfuerzos de la simulación fueron dedicados en gran medida a los
problemas de dos fases (gas y agua) y, en tres fases, así como modelos de petróleo.

La simulación de métodos de recuperación se limitaba esencialmente a los problemas de

agotamiento natural y de mantenimiento de presión. Con esto era posible el desarrollo de un
modelo de simulación único, capaz de dirigirse a la mayoría de los problemas de yacimientos
que se tenían. Este concepto de un modelo general siempre ha sido atractivo debido a que
significa reducción en el costo de su preparación y de su uso y, potencialmente, en el costo del
desarrollo del modelo y de su mantenimiento.

Sin embargo, durante los años 1970´s el panorama cambió radicalmente. El aspecto económico
motivó a que se buscara la forma de obtener una mayor recuperación, llevándose a efecto
proyectos de pruebas de campo o pruebas piloto encaminadas al estudio de procesos de
recuperación mejorada. Esto condujo a la simulación de nuevos procesos que iban más allá del
depresionamiento convencional y del mantenimiento de presión, tales como la inyección de
miscibles, la inyección de vapor, la inyección de productos químicos y la combustión in-situ.
Con esto, al manejo relativamente cómodo de dos componentes hidrocarburos (gas y petróleo)
en flujo simple inmiscibles, había que agregarle entonces la influencia de la temperatura, agentes
químicos y los efectos del comportamiento complejo del equilibrio entre fases. La proliferación
que tuvieron estos métodos de recuperación en los años 1970´s motivo la orientación del
concepto de modelo único o general hacia modelos individuales desarrollados para representar
cada una de estas nuevas técnicas.

Las investigaciones realizadas durante este tiempo, dieron como resultado un avance
significativo en la formulación de modelos de simulación y de métodos numéricos para la
solución de sistemas de ecuaciones. Estos avances permitieron simular procesos de recuperación
de lo más complejo y/o reducir el costo de tiempo del computador.

En la actualidad el enfoque de la simulación es el de afinar los avances que se han obtenido y

volver a tender hacia un simulador general aplicable a todos o a la mayoría de los procesos de
recuperación que interesen. El éxito depende en gran parte, de la obtención de ecuaciones de

Simulación de Reservorio
8
estado que representen el comportamiento PVT de los componentes de un sistema de fluidos en
tres fases bajo un rango de presiones y temperaturas bastante amplio.

Aunado a esto, se vienen desarrollando métodos de enmallado más eficaces que permitan
capturar con más exactitud los detalles locales del medio poroso y de heterogeneidades
normalmente halladas en un yacimiento. Tal es el caso, de refinamiento de mallas
convencionales para describir mejor el pozo dentro del yacimiento, o las nuevas mallas PEBI –
Bisección Perpendicular- que permiten modelar más realísticamente los procesos intrincados de
flujo dentro de una fractura, pozo horizontal, o la caracterización adecuada de fallas, lentes, y
discontinuidades.

1.4.- Objetivos de la Simulación

La simulación de yacimientos es un proceso mediante el cual el ingeniero con la ayuda de un

modelo matemático, integra un conjunto de factores para describir con cierta precisión el
comportamiento de procesos físicos que ocurren en un yacimiento.

Básicamente, un modelo matemático de simulación de yacimientos, consiste en un número

determinado de ecuaciones que expresan el principio de conservación de masa y/o energía,
acoplada con ecuaciones representativas de flujo de fluidos, temperatura y/o la concentración
de estos fluidos a través de medios porosos. Las Ecuaciones aplicadas son ecuaciones
diferenciales en derivadas parciales no lineales, su solución es posible únicamente en forma
numérica y de manera discreta, es decir, en un número de puntos preseleccionados en tiempo y
en espacio y no de una manera continua.

La no linealidad delas ecuaciones se debe principalmente a:

 La heterogeneidad en el yacimiento.
 La relación no lineal entre la saturación con la presión capilar.
 Las propiedades PVT de los fluidos son funciones no lineales de la presión, composición
y temperatura.

Los modelos matemáticos requieren el uso de un programa de computo debido a la cantidad de

cálculos tan grande que se realizan al efectuar una simulación. El objetivo primordial al hacer
uso de la simulación es predecir el comportamiento de indeterminado yacimiento y con base a
los resultados obtenidos, optimizar ciertas condiciones para aumentar la recuperación.

La selección del modelo a utilizar, además del aspecto económico, está en función de lo que
sede sea simular y de la información con que se cuente para realizar la simulación, pero una
regla general es utilizar el modelo más simple capaz de resolver el problema planteado.

1.4.1. Utilidad de la Simulación

La simulación de yacimientos constituye la herramienta más poderosa con que cuenta el

ingeniero siempre y cuando la geología y las propiedades de los fluidos están propiamente

Simulación de Reservorio
9
caracterizadas y el modelo matemático de simulación ha sido probado y calibrado
adecuadamente. Mientras que físicamente el yacimiento puede producirse una sola vez y lomas
probables es que no sea en la forma más adecuada, dado que un error cometido en el proceso
afectará cualquier cambio subsiguiente, el modelo permite “producir” un yacimiento varias
veces y en muy diferentes maneras, con lo cual se pueden analizar varias alternativas y
seleccionar el mejor escenario.

El observar el comportamiento del modelo bajo diferentes condiciones de operación, ayuda a

seleccionar un conjunto de condiciones de producción óptimas para el yacimiento. Más
específicamente, con la ayuda de la simulación, se puede hacer lo siguiente:

 Conocer el volumen original de hidrocarburo.

 Tener una buena idea del movimiento de los fluidos dentro del yacimiento.
 Determinar el comportamiento de un campo de hidrocarburo bajo diversos mecanismos
de desplazamiento, como puede ser: la inyección de agua, la inyección de gas,
depresionamiento natural o el uso de algún método de recuperación mejorada.
 Determinar la conveniencia de inyectar agua en un yacimiento de petróleo por los
flancos en lugar de utilizar un patrón determinado de pozos inyectores o viceversa.
 Optimizar los sistemas de recolección.
 Determinar los efectos de la localización de los pozos y su espaciamiento. De esta
manera desarrollar un campo con base en una información limitada, pudiéndose
determinar donde perforar nuevos pozos.
 Estimar los efectos que tiene el caudal de producción sobre la recuperación.
 Calcular la cantidad de gas que se obtiene de un número determinado de pozos
localizados en puntos específicos.
 Definir valores de parámetros en el yacimiento, para llevar a cabo estudios económicos.
 Obtener la sensibilidad de los resultados o variaciones en las propiedades petrofísicas
del yacimiento o las propiedades PVT de sus fluidos cuando no son bien conocidas.
 Realizar estudios individuales de pozos.
 Conocer la cantidad de gas almacenado.
 Hacer un programa de producción.
 Simular un proceso físico específico: resultados de una inyección de agua, una prueba
de presión, etc.

1.4.2. Ajuste de Simulador con la Historia del Yacimiento

Si la información con que se cuenta para llevar a cabo una simulación es amplia y decalidad,el
objetivo de la simulación tenderá a satisfacerse la predicción del comportamiento. Si por el
contrario, la información es incompleta o no muy confiable, los simuladores solo podrán
utilizarse para comparar semicuantitativamente los resultados al explotar el yacimiento de
diferentes maneras. De cualquier forma, la apropiación que proporciona el simulador puede
mejorarse mediante el ajuste de éste a medida de que se vaya obteniendo mayor información del
yacimiento.

Simulación de Reservorio
10
Lo primero que se hace para ajustar el simulador con la historia del yacimiento, es calcularel
comportamiento de éste usando la mejor información disponible. De esta manera los resultados
obtenidos de la simulación se comparan con aquellos obtenidos del campo, esto es, con los
valores reales. Si los resultados al compararlos no coinciden en una manera satisfactoria, se
hacen modificaciones en los datos utilizados y se efectúan otras corridas del simulador hasta
que se alcanza la aproximación deseada en los resultados. Cuando esto ocurre, el modelo ya
puede ser utilizado para predecir con cierto grado de precisión, del comportamiento del
yacimiento. Es importante notar que dicho comportamiento está influenciado por muchos
factores tales como: permeabilidades, distribución de saturaciones, espesores de las capas,
porosidades, permeabilidades relativas, etc. que nunca se conocen con exactitud a lo largo de
todo el yacimiento.

Por esto se debe tener en cuenta que al utilizar un simulador, después de haberlo ajustado a la
historia del yacimiento, no se puede asegurar la predicción que proporcione será exactamente el
comportamiento real de dicho yacimiento. Sin embargo, a medida de que el periodo ajustado
sea mayor, la predicción que se haga será mas confiable, lo que implica que el ingeniero deba
estar continuamente comparando la predicción hecha por el simulador con el comportamiento
presente y actualizar de ser necesario, las combinaciones de datos que maneja el modelo.

1.4.3. Resultados de una Simulación

Los resultados típicos que se obtienen de una simulación consisten de la distribución

depresiones y de saturaciones en cada una de las celdas en que ha sido dividido el yacimiento,
y de los volúmenes producidos y las relaciones agua-petróleo y gas-petróleo para los pozos
productores. Si hay inyección de fluidos se obtiene, el ritmo de inyección de los pozos o las
presiones necesarias para inyectar los volúmenes establecidos.

1.5. Etapas para Desarrollar un Modelo

El desarrollar un modelo es un proceso iterativo que consiste de las siguientes etapas:

 Descripción del yacimiento

 Determinar el tipo de mecanismo de desplazamiento
 Escribir el modelo matemático
 Desarrollar el modelo numérico
 Desarrollar el programa de cómputo
 Determinar la validez del modelo
 Ajustar el modelo con la historia del yacimiento
 Predecir el comportamiento futuro

El proceso iterativo mencionado se puede observar en la Fig. 2.1, ya que al ir avanzando ellas
diferentes etapas, es necesario regresar a modificar algo de las anteriores, como pueden ser las
suposiciones en las que se basó el modelo. Las razones de considerar varias suposiciones al
desarrollar un modelo, son las siguientes:

Simulación de Reservorio
11
  No obstante de haberse hecho todo lo posible por caracterizar al yacimiento de la mejor
manera, nunca podrá hacerse esta sino solo en una forma aproximada.
 Hacer el problema manejable.
 Reducir el costo de la simulación.

Obviamente la necesidad de utilizar suposiciones se hace cada vez menor, debido a los adelantos
e innovaciones que la ciencia va proporcionando día a día, especialmente en velocidad de
procesamiento de datos, la cual fue una limitante significativa antes de 1980.

1.6. Como Trabaja un Modelo

El procedimiento de cálculo en forma simplificada que utiliza un modelo, está dado por los
siguientes pasos:

 Se empieza con las celdas en las que se conoce su presión y su saturación inicial.
 Se selecciona un incremento de tiempo al cual el modelo vana efectuar los cálculos (los
incrementos de tiempo iniciales son generalmente, periodos de tiempo cortos alrededor
de uno o varios días, pero en los intervalos sucesivos los incrementos de tiempo pueden
ir aumentando hasta llegar a cubrir algunos meses).
 Calcular o asignar el volumen producido o inyectado, si es el caso, en cada pozo para el
intervalo de tiempo seleccionado.
 Calcular el flujo que hay entre las celdas durante el intervalo de tiempo utilizado y los
nuevos valores de saturación para cada celda.
 Seleccionar un nuevo intervalo de tiempo y repetir el proceso hasta que el modelo haya
hecho los cálculos para el tiempo total deseado.

De esta manera el simulador calculará la distribución de presiones y de saturaciones a lo largo

del yacimiento en función del tiempo. En la Fig. 2.2 se puede observar un diagramade flujo que
da idea de cómo trabaja un modelo.

En los capítulos siguientes se estudiará con más detalles cada una de las etapas de cálculo, de
manera que el lector se vaya familiarizándose con lo que realmente es la simulación de
yacimientos, para que tenga la capacidad de escribir modelos sencillos y pueda resolver
problemas en una y dos fases, así como en una y dos dimensiones.

Simulación de Reservorio
12
 Figura No 1.1 Etapa para desarrollar un Modelo de Simulación

Figura 1.2 Flujo de Modelo Matemático

Simulación de Reservorio
13
1.7.- Consenso en el modelaje. 1

El modelado de reservorios es la aplicación de un sistema de simulación en computadora para

la descripción del flujo en un reservorio.

El sistema de simulación es usualmente uno o más programas para computadora. Este sistema
será llamado de aquí en adelante el simulador del reservorio, y el conjunto de datos de entrada
será llamado “modelaje del reservorio”.

Muchas disciplinas contribuirán a la preparación del conjunto de datos de entrada.

La información es integrada durante el procesamiento del modelaje del reservorio, y la

concepción de los reservorios es cuantificada en el simulador del reservorio. La figura muestra
la contribución de las diferentes disciplinas

Interpretación
Petrofísica
Sísmica

Modelo Propiedades
Geológico de los fluidos

Curvas de
Pozos
Entubado

Modelo de
Simulación
Numérica
Efectos de
Medios Enrejado del
modelo

Figura 1.3 Contribución disciplinaria al modelaje de Reservorios.

El simulador es el punto de contacto entre las anteriores disciplinas. Sirve como un filtro que
selecciona entre todas las descripciones propuestas del reservorio. Este simulador no está
influenciado por argumentos o estilos de presentación. Provee una apreciación objetiva de cada
hipótesis, y reprime el poder de la influencia personal descrito por Milheim. Como un filtro de

1SCHLUMBERGER Publications, Martin Kaeser, Tobias Werther, “A Course in Reservoir

Simulation”, 2001, Pág.4-6 (1-3)

Simulación de Reservorio 14
hipótesis, el modelador de reservorios es a menudo el primero en saber cuando una hipótesis
propuesta sobre el reservorio es inadecuada.

Una de las más importantes tareas del modelador es lograr consenso en respaldo de la
representación de un reservorio. Esta tarea se hace más compleja cuando los datos disponibles
del desempeño del campo pueden asociarse a más de un modelo de reservorio. La no
singularidad del modelo será discutida en mayor detalle mas adelante. Esto significa que hay
más de una manera de percibir y representar los datos disponibles.

El modelador debe elegir a través de varias representaciones de reservorios y buscar consenso

entre todas las posibles.

Esto se hace a menudo rechazando una o más representaciones propuestas. Como consecuencia,
el elemento humano es un factor en el proceso, particularmente cuando los datos no respaldan
claramente la selección de una sola representación de reservorio de un conjunto de
representaciones competitivas.

El criterio dual de racionalidad y Ockham’sRazor son esenciales en este proceso, como también
lo es un entendimiento de como los individuos pueden contribuir mas efectivamente al
modelado.

1.8.- Manejo de Estudios de Simulación. 2

El manejo apropiado de estudios de la simulación reconoce estas etapas y permite tener un

tiempo para que el proceso de construcción del equipo madure.

Los estudios modernos de simulación de campos mayores son realizados por equipos que
funcionan como equipos proyectados en una matriz de manejo de organización.

Matriz de manejo es sinónimo de Manejo de proyecto, pero tiene dos características distintas:

 Organización funcional cruzada con miembros de diferentes áreas de trabajo que se

toman en un proyecto.
 Un empleado es responsable a dos o más superiores, lo cual puede causar dificultades
para superiores y empleados.

Para evitar potenciales problemas, el equipo debe constituirse tal como se muestra:

 Cada miembro del equipo es asignado a una tarea diferente.

 Todos los miembros trabajan en busca de la misma meta.

2
SPE 72160, “Integrated Reservoir Characterization and Simulation for Reservoir Management using a Web-Based
Collaborative Technical Workflow Manager”, HaiZuiMengetalii, 2001, Pág.3-4.

Simulación de Reservorio
15
Los miembros del equipo deberían tener roles únicos para evitar funciones redundantes. Si las
responsabilidades de dos o más miembros del equipo coinciden considerablemente, la confusión
puede suceder con respecto a las áreas de responsabilidad. Cada miembro del equipo debe ser
la clave de las tomas de decisión en cada disciplina particular, de otra manera las disputas
pueden no resolverse dentro del tiempo disponible para completar el estudio. A los equipos no
se les debería permitir derrochar tiempo debatiendo sobre utopías que nunca funcionaran.

Los equipos deben esforzarse eficazmente por llegar a un consenso, pero la presión del
acercamiento de las fechas topes requerirá a un miembro del equipo para que funcione como
líder. Hay muchas áreas en donde las decisiones no se tomaran por acuerdo general. Por esta
razón, los equipos deberán tener un líder con las siguientes características:

 Habilidades técnicas significativas

 Amplia experiencia
 Los líderes deberían tener autoridad técnica y monetaria sobre el proyecto.
 Si son percibidos como faltos de autoridad será imposible desempeñar su función.
 Por otro lado, los líderes deben evitar el control autoritario, de lo contrario debilitarán al
equipo y arrollaran a un grupo.
 De acuerdo con Maddox, los equipos y grupos difieren en la forma en que se comportan.

El comportamiento de un grupo exhibe las siguientes características:

 Los miembros piensan que están agrupadas únicamente para propósitos administrativos.
Los individuos trabajan independientemente, a veces en propósitos comunes.
 Los miembros tienden a concentrarse en ellos mismos porque no están lo
suficientemente inmersos en la planificación de los objetivos como una unidad.
Aproximan su trabajo a solamente una función de empleados contratados.

En contraste, las características del comportamiento de un equipo son las siguientes:

 los miembros reconocen su independencia y comprenden que las metas personales y las
metas del equipo son mejor logradas con apoyo cooperativista. El tiempo no es
desperdiciado esforzándose fuera de las labores propias o buscando ganancias
personales a expensas de los demás’.
 Los miembros sienten un sentido de pertenencia sobre sus funciones porque están
comprometidos con metas que ellos mismos ayudaron a establecer.

Observaciones similares fueron hechas por Haldorsen y Damsleth:

 los miembros de un equipo deben necesariamente entenderse entre si, respetarse, actuar
como el abogado del diablo con unos y otros, y mantenerse informados entre si.

Haldorsen y Damsleth defendieron que el objetivo de cada equipo debería ser:

 Innovación y creación de valor a través del acercamiento del equipo.

Simulación de Reservorio
16
  Orientación del cliente con enfoque en “mi es su comienzo”.

McIntosh apoyó la noción de que cada miembro del equipo debería ocupar un rol funcional, por
ejemplo, geólogo, ingeniero, etc. Un corolario es que los miembros del equipo pueden entender
sus roles porque estos han sido claramente definidos.

El manejo apropiado puede mejorar la probabilidad de que un equipo funcionara como debería.
Un sentido de pertenencia puede ser criado si los miembros de un equipo participan en el
planeamiento y la toma de decisiones.

Las visiones de un miembro del equipo deberían influenciar el trabajo y el horario de

actividades. Muchos problemas pueden ser evitados si expectaciones realistas son creadas en
los horarios del proyecto al comienzo, y después adheridos a lo largo del proyecto, expandiendo
el alcance del trabajo sin alterar la asignación de recursos o fechas tope, puede desmoralizar y
minar el concepto de equipo.

Finalmente, es de importante cuidado tener en mente cuando se realice estudios utilizando

equipos que: ‘‘menos ideas son generosas por grupos que por individuos trabajando solos’’, una
conclusión sostenida por evidencia empírica de la sociología.

Simulación de Reservorio
17
CAPITULO II Proceso del Modelaje

2.1.- Introducción

La mejor tecnología para hacer predicciones sobre el desempeño de los reservorios hoy en día
es modelar el flujo de fluidos en medios porosos utilizando programas de computadora
conocidos como simuladores.

2.2.- Análisis básicos del Reservorio. 3

La caracterización del reservorio y las evaluaciones ingenieriles son usualmente realizadas
como una parte de las prácticas estándares de las negociaciones, independientes del estudio del
reservorio en simulación. Las tareas asociadas con el análisis básico de un reservorio
proporcionan información necesaria para preparar los datos de entrada para un estudio en
simulación. Incluso proporcionan un concepto inicial del reservorio el cual puede ser utilizado
para diseñar un modelo de estudio.
El equipo de modelado necesita estar al tanto sobre estudios anteriores y debería relacionar el
desempeño del modelo con estudios previos siempre que sea posible.

2.3.- Pre-requisitos.
Varios prerrequisitos deberían ser satisfactorios antes de emprender un modelo de estudio. El
más importante, desde el punto de vista de los negocios, es la existencia de un problema de
importancia económica. Los objetivos del estudio de un modelo deberían rendir soluciones para
los problemas de importancia económica. Una vez que se especifican los objetivos de un estudio,
el modelador debería plegar toda la información disponible y los reportes relacionados con el
campo.
El termino modelador es usado en relación al texto “equipo de modelado” como sinónimo,
mientras no se haga ningún tipo de distinción al respecto. Si es que existen datos necesarios que
no están disponibles, el modelador debería determinar si estos datos pueden ser obtenidos, por
analogía con otros reservorios o por correlación. Los valores para todos los datos de entrada de
un modelo deben ser obtenidos, porque el simulador no correrá si es que no cuenta con un
conjunto completo de datos. En algunos casos, se debe asumir simplificaciones sobre el
reservorio ya que no se cuenta con suficientes datos disponibles para representar
cuantitativamente el sistema en gran detalle.
En adición a los objetivos claramente definidos, otro pre-requisito que deber ser satisfactorio
antes de comprometer un estudio en simulación es determinar que los objetivos de un estudio
no pueden ser logrados utilizando las técnicas más simples. Con técnicas de menor costo, como
el análisis de curvas de declinación o el algoritmo de desplazamiento de flujo de Buckley-
Leverett, no se proveen resultados adecuados, es entonces cuando se justifican métodos más
costosos y sofisticados.

3 SPE textbook:”BASIC APPLIED RESERVOIR SIMULATION”, Ertekin, Abou-Kassem,King,

vol.7, 2001, Pág.11.

Simulación de Reservorio 18
2.4.- Modelaje computacional. 4
Un modelo comprensivo de manejo de reservorios puede ser pensado como cuatro modelos
interactivos: el modelo de reservorio, el modelo de pozo, el modelo de perforación, y el modelo
de superficie. El modelo de reservorio representa el flujo en el reservorio.
El reservorio es modelado subdividiendo el volumen del reservorio en un arreglo, o enrejado,
de elementos de menor volumen. Muchos nombres son utilizados para denotar un volumen
individual, por ejemplo, bloque de grilla, celda, nodo. El conjunto de todos estos volúmenes es
conocido como enrejado o grilla.
Cada simulador de reservorio incluye un modelo de reservorio y un modelo de pozo.
El modelo de pozo es un término en las ecuaciones de flujo de fluido que representan la
extracción de fluidos del reservorio o la inyección de fluidos en el mismo. Los simuladores
comerciales completos incluso incluyen en su representación un modelo de perforación y un
modelo de superficie.
El modelo de superficie representa contención asociado a las estructuras superiores, tales como
plataforma y delimitadores.
Los algoritmos matemáticos asociados con cada modelo dependen de leyes de conservación de
la materia, momento y energía. Los más ampliamente utilizados asumen que el reservorio
permanece isotérmico. Si es que estamos modelando un reservorio en el cual los efectos
térmicos son considerados de importancia, en procesos tales como recuperación secundaria en
donde el calor es inyectado de alguna manera, entonces necesitamos usar un simulador que
cuente con variaciones de temperaturas y demás efectos termodinámicos asociados.
El conjunto de algoritmos es lo suficientemente complejo que, la alta velocidad de las
computadoras es la única manera práctica de resolver los procesos matemáticos asociados con
el estudio en la simulación de un reservorio.

2.5.- Elementos principales en el estudio de Simulación de Reservorios. 5

Los elementos esenciales de un estudio de simulación incluyen a la asociación de historiales
sobre el campo, hacer predicciones, incluir previsiones basadas en la estrategia operativa
existente, y evaluar escenarios alternativos de operación. Durante los ajustes de los historiales
sobre de campo, el modelador verificará y refinará la descripción del reservorio. Comenzará con
una descripción inicial del reservorio, el modelo es usado para asociar y predecir el desempeño
del mismo. Si es necesario el modelador modificará la descripción del reservorio hasta que se
llegue a una asociación aceptable.
La fase de asociación de historiales sobre el campo puede ser considerada un problema inverso
ya que la respuesta de por si ya existe. Sabemos cómo el reservorio se desempeñará, lo que
buscamos es entender porque. La tarea es encontrar el conjunto de parámetros del reservorio

4 SPE 89974, “Top-Down Reservoir modelling”, 2004, G.J.J Williams; M.Mansfield: D.G.
MacDonald; M.D.Busch, Pag 1-3.
5 SPE 89974, “Top-Down Reservoir modelling”, 2004, G.J.J Williams; M.Mansfield: D.G.

MacDonald; M.D.Busch,

19
Simulación de Reservorio
que minimizan la diferencia entre el desempeño del modelo y el desempeño histórico del campo.
Este no es simplemente un problema aislado, y hay más de una forma de asociar los datos
disponibles.
Una vez que se tiene una asociación histórica de datos, el próximo paso será hacer una
predicción, la cual es esencialmente una continuación de la ya existente práctica de
operaciones.Esta predicción nos da una línea de base para comparar las diferentes estrategias de
manejo de un reservorio. Se debe tener cuidado con la validez de las predicciones del modelo.
Una forma de obtener una idea de la exactitud de las predicciones es medir el éxito de las
previsiones hechas anteriormente. Aunque las previsiones del precio son esenciales para un
emprendimiento económico, está claro que existe una considerable incertidumbre asociada a las
previsiones del precio.
Las previsiones no toman en cuenta discontinuidades en los patrones históricos que surgen por
efectos inesperados. Esto es una verdad en el mundo físico, así como en el mundo social. En
adición a la incertidumbre en los parámetros económicos, existe incertidumbre en el desempeño
productivo del campo. Los simuladores no eliminan la incertidumbre, nos dan la posibilidad de
evaluar y mejorar el manejo del riesgo asociado con la predicción del desempeño productivo.
Un valioso pero intangible beneficio del proceso asociado con la simulación de reservorios es
la ayuda que nos proporciona para el manejo del reservorio. Una de las tareas críticas del manejo
de reservorios es la adquisición y mantenimiento de un banco de datos que permanezca al
corriente.
Un estudio en simulación puede ayudar a coordinar actividades y modelar recursos necesarios
para determinar el plan óptimo para operar un campo. Recolectar datos de entrada para un
modelo es una buena forma de asegurar que cada técnica importante es considerada, de las
diferentes disciplinas que contribuyen al manejo de un reservorio. Si el modelo es
particularmente sensible a un determinado parámetro, entonces se debe idear un plan para
determinar dicho parámetro con mayor exactitud, por ejemplo, las pruebas de laboratorio más
apropiadas.
Las Escalas Conceptuales del reservorio son muy importante considerarlas Una de las metas
más importantes del modelado es reducir el riesgo asociado con la toma de decisiones en un
ambiente donde se tiene un conocimiento algo limitado. El rango de aplicabilidad de datos
exactos y la integración de datos dependientes de una escala en un concepto cohesivo de
reservorio se discutirán en este capitulo.

2.6.- Escala y muestreo del Reservorio. 6 7

Un buen sentido de entendimiento sobre un reservorio puede ser obtenido considerando un área
de muestra del mismo, mediante diferentes técnicas. Como ejemplo supongamos que queremos
encontrar el tamaño del área de muestreo en el wellbore que tiene un radio de seis pulgadas. Si
asumimos un área circular, podemos calcular el área con πr2. El área resultante de la muestra es
de menos de un pie cuadrado. Para determinar la fracción del área de la muestra, normalizamos

6Alain Brie, “Carbonates Evaluation”, paper, 2001, Pág.2.

7SPE 77715, “On the Relationchip between Resistivity and Permeability Anisotropy, 2002,
Georgi, Bespalov, Tabarovsky, Baker Hughes & Schoen, Pág.2-3.

20
Simulación de Reservorio
el área de la muestra con respecto al área de drenaje del pozo, de unos muy modestos cinco
acres. ¿Qué fracción de área está directamente muestreada por el wellbore? El área de drenaje.
La fracción del área de muestreo es diminuta fracción del área de interés.
Una señal de registro de pozo se expandirá en el área que está siendo muestreada. Suponiendo
que con un registro de pozo se puede penetrar la formación hasta cinco pies desde el wellbore,
lo cual es una asunción razonablemente generosa. La fracción de área que está siendo
muestreada es ahora de aproximadamente cuatro parte en diez mil. El área de muestreo es un
área de drenaje de cinco acres, lo cual es un área de drenaje pequeña, es aún un fragmento de
uno por ciento.
Los núcleos y registros de pozo nos dan una visión muy limitada del reservorio. Una sección
sísmica expande la fracción de área de la muestra, pero la interpretación de datos sísmicos es
menos precisa. Los datos sísmicos son a menudo vistos como “datos de relleno” ya que
dependen de una determinada interpretación. La fiabilidad de las interpretaciones sísmicas
puede ser mejorada correlacionando datos fiables como son la información de núcleos y
registros. El rango de aplicabilidad de los datos de medición depende de la técnica de muestreo
utilizada. ¿Tomamos muestras de núcleos, medimos respuestas eléctricas del registro de pozo,
o detectamos energía acústica?. Los rangos son ilustrados en la figura 2.1.

Núcleo Registros Sísmica

0.0

20

12.7 40
cm

60

80

Figura 2.1 Rango de Datos en pruebas de muestreo.

Fayers y Hewett apuntaron que las definiciones de escalas no son universalmente aceptadas,
pero ilustran la escala relativa asociada con mediciones apropiadas en reservorios. El rango de
tamaño de las escalas va desde realmente grandes a microscópicas. Para reconocer las
variaciones en el rango deon aplicabilidad de los datos, han sido definidas cuatro escalas
conceptuales las cuales serán adoptadas para ser utilizadas en discusiones siguientes.
La escala Giga incluye información asociada con técnicas geofísicas tales como la arquitectura
del reservorio, teorías de carácter regional, tales como las placas tectónicas, datos sísmicos y
satelitales, que pueden ser interpretadas.
La escala Mega es la escala de caracterización del reservorio e incluye registros de pozo, las
pruebas de pozo y los análisis sísmicos en 3D. La escala Macro se enfoca en datos de muestra
a nivel del análisis de núcleo y el análisis de las propiedades de los fluidos.

Simulación de Reservorio
21
La escala Micro incluye datos de escala de los poros obtenidos mediante técnicas tales como el
análisis de secciones finas y mediciones de granulometría. Cada una de estas escalas contribuye
al modelo final del reservorio.

Micro

Macro

Giga
Mega
Figura 2.2. Integración de Escalas de Reservorio.
Toda información recolectada de las diferentes escalas debe ser integrada en una sola,
comprensiva y consistente representación del reservorio. La integración de los datos obtenidos
en diferentes escalas es un problema que a menudo es referido como “scale-up”o
“upscaling”.Los intentos de relacionar datos de dos diferentes escalas pueden resultar
dificultosos. Por ejemplo, la permeabilidad es a menudo obtenida de la prueba de presión
transiente y las rutinas de análisis del núcleo. Las respectivas permeabilidades, pueden parecer
no concordantes porque están representadas en dos escalas muy diferentes.
Una tarea importante del proceso de integrar escalas es desarrollar un buen entendimiento de
cómo los parámetros que han sido medidos varían con cada escala. El enfoque en detalle en uno
o más aspectos del proceso de modelado del reservorio puede opacar el concepto fundamental
del reservorio en un modelo de estudio.

2.7.- Escalas integradas – Unidades de Flujo. 8

Una manera de integrar datos disponibles en el contexto de “hacer una sola foto” es aplicar el
concepto de unidad de flujo. Una unidad de flujo es definida cono “un volumen de roca
subdivido de acuerdo con las propiedades geológicas y petrofísicas que influencian el flujo de

8SCHLUMBERGER Publications, “A Course in Reservoir Simulation”, Martin Kaeser, Tobias

Werther, 2001, Pág.7.

Simulación de Reservorio 22
fluidos a través de el”. Las propiedades geológicas y petrofísicas típicas son mostradas en la
tabla 2.1.

Propiedades geológicas Propiedades Petrofísicas

Textura Porosidad
Mineralogía Permeabilidad

Estructura sedimentaria Compresibilidad

Contactos Saturación de fluidos

Barreras de Permeabilidad

Tabla 2.1 Propiedades Geológicas y Petrofísicas

Un reservorio es modelado subdividiendo su volumen en un arreglo de volúmenes elementales
representativos (VER). El concepto de VER no es el mismo que el de unidad de fluido. Una
unidad de fluido es una parte contigua del reservorio que tiene similares propiedades de flujo
caracterizada por datos geológicos y petrofísicos. Las unidades de flujo contienen usualmente
uno o más VER. Por el contrario, los VER son el volumen de elemento del tamaño suficiente
para proporcionar un promedio de valores y parámetros estadístico y significativo describiendo
el flujo contenido en el volumen, pero lo suficientemente pequeño para proporcionar una
aproximación numérica significativa de las ecuaciones fundamentales de flujo.
El concepto de unidad de flujo es una forma efectiva de manejar la creciente base de datos
proporcionada por los geocientistas. Incrementando el refinamiento del análisis de los
geocientistas se puede obtener mayor información. Incluso en nuestros días, con 100000 bloques
de grillas de flujo, los modeladores no pueden usar toda la información proporcionada por
modelos geológicos en computadora, los cuales pueden estar basados en más de un millón de
puntos de grilla. Es aun necesario minimizar los modelos geológicos en unidades de flujo
representativas. Un entendimiento de “una foto grande”, incluso en un esquema simple, es un
recurso valioso para validar las ideas siendo cuantificadas en un modelo.
Richardson, esbozó varios tipos comunes de modelos de reservorios: ventilador de aguas
profundas, delta arenoso, un canal deltaico contrastado con una barra deltaica, etc.
Sus esbozos ilustraron la forma como el reservorio podría verse según un conjunto de
especificaciones. En muchos casos especialmente en campos relativamente pequeños, la mejor
imagen de reservorio puede ser solo una imagen cualitativa. Cuando se comience con estudios
detallados, la imagen cualitativa puede ser cuantificada. La confianza en el desempeño del
modelo se obtiene usando el modelo asociando al desempeño histórico del campo. Asociar el
historial y aprobar el modelo serán temas discutidos más adelantes. Como regla, los modelos
deberán ser actualizados y refinados cada vez que se obtengan y añadan más datos del campo.

Simulación de Reservorio
23
2.8.- Estructura del Reservorio
El tamaño físico y la forma del reservorio pueden ser inferidos por varios métodos que sirven
como recursos de información para definir la medida de la escala y la estructura del reservorio.
Estos recursos de información son brevemente revisados a continuación
2.8.1.- Giga Escala. 9
Las mediciones sísmicas discutidas por autores como Ausburn, McQuillin y Sherif, proveen
mucha de la información de la escala Giga que puede ser directamente utilizada para caracterizar
un reservorio. Históricamente los análisis sísmicos han sido de interés para establecer el tamaño
estructural del reservorio. La gente no creía que los datos sísmicos pudiesen proporcionar
suficiente detalle para tener información sobre toda la estructura del reservorio.
Pero esa forma de pensar ha cambiado con la aparición del monitoreo sísmico en 4D y la
geofísica de los reservorios. Ahora es importante introducir algunos conceptos geofísicos
básicos en el contexto de la función de manejo de reservorios. Las ondas sísmicas son
vibraciones que se propagan desde una fuente, tal como una explosión, a través de la tierra hasta
que encuentran una superficie reflectora y son reflejadas hacia un detector, tal como un geófono.

Respuesta Sísmica
Cuña de Arena

Cuña de Arcilla

Figura 2.3. Trazo sísmico para una cuña

Cada trazo representa la señal recibida por un detector, los cambios en la amplitud de ondas
sísmicas ocurren en los reflectores. Una reflexión sísmica ocurre en la interfase entre dos
regiones con diferentes impedancias acústicas. La impedancia acústica es fundamental como
parámetro sísmico. La impedancia acústica se define como Z   * V en donde:

ρ = es la densidad del medio.

V = es la velocidad de condensación de la onda en ese medio.
La figura 2.4 muestra la correlación entre la velocidad de la onda sísmica y la densidad de
diferentes tipos de roca.

9SCHLUMBERGER Publications, “Oilfield Review”, 2001, Pág.3-5.

Simulación de Reservorio 24
 Figura 2.4. Velocidad de Onda Sísmica a través de diferentes densidades de roca.
Un cambio en la impedancia acústica causara una reflexión de la onda sonora. La habilidad de
reflejar una onda sonora cambiando la impedancia acústica es cuantificada en términos del
coeficiente de reflexión.
El coeficiente de reflexión R en la interfase entre dos capas continuas es definida en términos
de la impedancia acústica como:

Z 2  Z1  2V2  1V1
R  Ec. 2.1
Z 2  Z1  2V2  1V1

Donde los coeficientes 1 y 2 están referidos a dichas capas contiguas. Las magnitudes de los
coeficientes de reflexión para interfases de capas típicas son ilustradas en la tabla 2.2. Los
valores de los coeficientes de reflexión entre materiales arenosos/calcáreos muestra que los
valores de los mismos pueden ser relativamente pequeños. En adicción al coeficiente de
reflexión se puede definir un coeficiente de transmisión.

Interfase Coeficiente de Reflexión

Arenosos a Calcáreos 0.04
Calcáreos a Arenosos -0.04
Fondo Marino 0.11 (suave) a 0.44 (duro)

Tabla 2.2. Coeficientes de Reflexión.

Un coeficiente es distinto de cero siempre que exista un cambio en las impedancias acústicas,
por haber un cambio en la velocidad de la onda que se propaga de un medio hacia otro, o porque

Simulación de Reservorio
25
las densidades del medio difieren. Si el cambio en la impedancia acústica es lo suficientemente
grande las reflexiones pueden ser medidas en la superficie. Es por esto que el gas tiende a
mostrarse como manchas luminosas en los datos sísmicos, hay un gran cambio en la densidad
del fluido.
Por el contrario, la presencia de un contacto petrolero/agua es más difícil de observar con
mediciones sísmicas porque las diferencias de las densidades entre ambos, agua y petróleo es
relativamente pequeña y como consecuencia da pequeñísimos cambios en la impedancia
acústica.
Los trazos sísmicos, amplitud sísmica versus desplazamiento en dos sentidos, o el tiempo que
toma que la onda sísmica se propague desde la fuente al receptor. Uno de los problemas
centrales en el procesamiento de datos sísmicos es determinar la conversión
tiempo/profundidad. La conversión de los datos de tiempo de viaje para formar profundidades
requiere que la velocidad asociada con cada zona geológica sea conocida o que se pueda saber
como la onda evoluciona con el tiempo.
Cuando la conversión tiempo/profundidad es aplicada a datos sísmicos, puede cambiar
profundidades de amplitudes sísmicas asociadas a trazos adyacentes. Los registros sónicos son
típicamente usados para calibrar los datos sísmicos cuando estos se usan para caracterizar un
reservorio. La respuesta del registro sónico delinea la cima y la base de una sección geológica
La longitud de onda de la onda sísmica es la velocidad de la onda dividida por su frecuencia.
Alternativamente, la longitud de onda es la velocidad en tiempos medios dados del periodo de
la onda. La frecuencia de la onda es una medida de la energía de la onda y es conservada
mientras la onda se propaga desde un medio hacia otro.
Cuando las ondas se superponen crean una perturbación, como se muestra en la figura 2.5. El
tiempo asociado a la duración de esta perturbación se denota como ∆t. esta perturbación tiene
una velocidad V en un medio, y el periodo T es la anchura de la misma cuando es dibujada en
un mapa de datos sísmicos. La longitud de la onda es igual a la velocidad V por el periodo T.Así,
la perturbación tiene un periodo de 10 milisegundos y la velocidad es 5000 pies por segundo en
un medio dado la longitud de onda es 50 pies.

Perturbacion

V = Velocidad en un medio
T = ∆t = Periodo de la perturbación

Simulación de Reservorio
26
 Figura 2.5. Perturbación Sísmica
Si los datos sísmicos no tienen el suficiente poder para mostrar los límites de reflexión de una
capa geológica, entonces las amplitudes de las ondas sísmicas pueden ser útiles par caracterizar
más las propiedades petrofísicas del reservorio. Por ejemplo supongamos una región del
reservorio caracterizada por una porosidad, permeabilidad K, espesor neto h, y una saturación
de petróleo So, la amplitud sísmica puede estar relacionada con la calidad de la roca (por ejemplo
Kh o Kh) o la capacidad productiva de petróleo. Cuando una correlación existe entre la
amplitud sísmica y una agrupación de parámetros petrofísicos, la correlación puede ser usada
de guía para la distribución de las propiedades del reservorio en áreas entre pozos.
Las figuras 2.6A y B muestran dos acercamientos a contornear un conjunto de valores a los
puntos de control. Las líneas de contorno lisas mostradas en la figura 2.6A son preferibles al
menos que las líneas de contorno onduladas estén acompañadas por datos adicionales.
Las correlaciones sísmicas pueden ser utilizadas para justificar el estilo de contorno mas
heterogéneo mostrado en la figura 2.6B.

Figuras 2.6 A y B Líneas de Contorno

2.8.2.- Mega Escala. 10

La escala Mega es la escala en la cual se comienza a integrar los registros y los datos de prueba
de pozo en un modelo funcional del reservorio. La Tabla 2.3 muestra el tipo de información que
puede ser obtenida a nivel de la escala Mega de los datos de registro de pozo. Las
interpretaciones más comunes de cada respuesta del registro están incluidas en la tabla. Por
ejemplo, una alta respuesta de rayo gamma implica la presencia de arcilla, mientras que una baja
respuesta de rayo gamma implica la presencia de arenas limpias o carbonatos. Una buena
combinación de herramientas de registro es usualmente necesaria para minimizar la ambigüedad
en la interpretación de registros.

Registro Variable Respuesta

Detecta Arcilla de radiación in situ

Rayo Gamma tipo de roca *Alto RG = arcilla

*Bajo RG = arenas limpias o carbonatos

10SCHLUMBERGER Publications, “Oilfield Review”, Volume 4, 2003, Pág.27-29.

Simulación de Reservorio 27
 Mide la Resistividad de Agua:
Tipo de
Resistividad *Alta Resistividad = Hidrocarburos
Fluido
*Baja Resistividad = Salmuera

Mide densidad del electrón detectando rayos gamma

Compton esparcidos. La densidad del electrón se
relaciona con la densidad de la formación. Buena para
detectar gas de baja densidad en comparación con roca
Densidad Porosidad o liquido.

*Baja Respuesta = Bajo contenido de Gas

*Gran Respuesta = Alto contenido de Gas

Mide la velocidad del sonido en un medio. La

velocidad del sonido es mas rápida en roca que en
fluidos.
Acústica
Porosidad *Mucho tiempo de viaje = Velocidad Baja = gran
(sónica)
cantidad de espacio poroso

*Poco tiempo de viaje = Velocidad Alta = poco

espacio poroso

Los electrones rápidos son detenidos por colisiones con

energías térmicas. Los neutrones térmicos son
capturados por los núcleos, y emiten rayos gamma
Contenido de detectables. Bueno para detectar gas.
Neutrón
Hidrogeno
*Gran Respuesta = Alto contenido de Hidrogeno

*Poca Respuesta = Bajo contenido de Hidrogeno

Mide el potencial eléctrico (voltaje) asociado con

movimientos de iones.
Potencial Asientos
Espontáneo permeables *Poca Respuesta = Arcillas impermeables

*Alta Respuesta = Asientos impermeables

Tabla 2.3 Respuesta de registros de pozo.

La Tabla 2.3 ilustra el tipo de información que puede ser obtenido a nivel de escala mega de los
datos de prueba del pozo. La tabla incluso muestra el tiempo más apropiado en la vida del
proyecto para realizar las pruebas. Usualmente es necesario realizar una serie de pruebas a
medida que el proyecto madura. Estas pruebas ayudan al operador a refinar su entendimiento
sobre el campo y a menudo motivar cambios en la forma como el pozo o el campo es operado.

Simulación de Reservorio
28
 Tipo de
Propiedades Etapa de Desarrollo
prueba

Comportamiento del Reservorio

Permeabilidad

Película
Prueba de Exploración y Apreciación
longitud de Fractura
Perforación del Pozo
Presión del Reservorio

Limite del Reservorio

Limites

Múltiples
Exploración y Apreciación
pruebas de Perfil de Presión
del Pozo
Formación

Comportamiento del Reservorio

Permeabilidad

Película Recuperación primaria,

Prueba de
secundaria y recuperación
Draw-down longitud de Fractura reforzada
Limite del Reservorio

Limites

Comportamiento del Reservorio

Permeabilidad
Película
Recuperación primaria,
Prueba de
longitud de Fractura secundaria y recuperación
Build up
reforzada
Presión del Reservorio
Limite del Reservorio
Limites

Separación de la formación Recuperación secundaria y

Presión recuperación reforzada

Simulación de Reservorio
29
 Pruebas de Permeabilidad
Proporción al
paso Película

Movilidad en varios bancos

Película

Pruebas de Presión del Reservorio Recuperación secundaria y

Fallo longitud de Fractura recuperación reforzada

Localización de Frente

Limites

Comunicación entre pozos

Pruebas de Tipo de comportamiento del Recuperación primaria,

Interferencia y Reservorio secundaria y recuperación
Pulso Permeabilidad entre pozos reforzada

Permeabilidad vertical

Propiedades de Capas individuales

Permeabilidad horizontal
Pruebas de Permeabilidad vertical Durante toda la vida del
Reservorios
Película Reservorio
capa
Presión promedio de Capas

Limites exteriores

Tabla 4.3. Información tenida a través de Escala Mega.

Las Tablas 2.3 y 2.4 ilustran unos pocos métodos usados para tratar información en escala Mega.
Los avances de la tecnología periódicamente añaden a la lista nuevas pruebas y herramientas.
En muchos casos, el constreñimiento presupuestario será el factor de control en la determinación
del número de pruebas a llevarse a cabo. El equipo de modelado debe trabajar con cualquier
información que tenga disponible.

Simulación de Reservorio
30
CAPITULO III Información Requerida del Proceso de Modelaje

3.1 Introducción

El éxito de un proyecto de ingeniería de yacimientos de basa en una buena información que

represente las condiciones que prevalecen en el yacimiento. Asípues la simulación sin ser la
excepción, requiere de una amplia descripción física delmismo y de los tipos de mecanismo por
medio de los cuales va a producir. Los resultadosque se obtengan de la simulación serán tan
buenos como los datos que se hayan empleadopara realizarla y, el tiempo que se pueda perder
en preparar esta información es un tiempobien empleado.

Hay que hacer notar que la información que debe tratarse de obtener con mayor precisión es
aquella que al variarla realizando diferentes corridas de simulación cause un cambiosignificativo
en los resultados obtenidos. Así por ejemplo, se sabe que una propiedad determinada varía en
un rango específico y al efectuar dos o tres corridas de simulación sevaría dicha propiedad
dentro de este rango y se obtienen resultados parecidos, se puedetomar como buena una de las
predicciones, o bien, relegar a segundo término esfuerzosadicionales para medir con precisión
dicha propiedad. Si por el contrario, variando esapropiedad se alteran los resultados
considerablemente, es necesario redoblar esfuerzos paraobtener con mayor aproximación dicha
propiedad.

La información que se requiere para efectuar una simulación es:

 Descripción física del yacimiento.

 Mecanismo o mecanismos de desplazamiento que operan en el yacimiento.
 Propiedades petrofísicas de las capas de interés.
 Propiedades PVT de los fluidos.
 Otros datos.

3.2. Descripción Física Del Yacimiento

Para obtener una descripción física del yacimiento es necesario llevar a cabo un
estudiogeológico de detalle que proporcione un conocimiento estratigráfico, estructural y
petrofísico, que permita de esta manera caracterizar al yacimiento perfectamente. Dichoestudio
geológico se completa con métodos geofísicos. La información de este tipo queinteresa a la
simulación es:

 Límites del yacimiento

 Características de la formación productora
 Características del acuífero
 Fallas
 Discontinuidad en las capas

Simulación de Reservorio
31
3.3. Mecanismos De Desplazamiento

Los cuatro mecanismos básicos que operan para recuperar los hidrocarburos del yacimientoson:

 Expansión del sistema roca-fluido

 Desplazamiento
 Segregación gravitacional
 Imbibición

La expansión del sistema roca-fluidos se provoca al haber un abatimiento de presión,

dandocomo resultado el movimiento de los fluidos a través del medio poroso del punto de
mayorpresión al punto de menor presión.

El desplazamiento se da con gas o con agua. Con gas puede ser empuje de gas disueltoliberado
o empuje de algún casquete de gas, ya sea natural o inyectado. Con agua puedeser agua de
inyección o bien entrada natural por la presencia de algún acuíferoconsiderable.

La segregación gravitacional se presenta en yacimientos de espesor considerable (o encapas

muy pronunciada) que tengan valores de permeabilidad altos en el sentidovertical y consiste en
el acomodo que tienen los fluidos de acuerdo con sus densidades.

La imbibición capilar se da generalmente en el sentido normal (perpendicular) al flujo ypuede

ser muy importante al inyectar agua en forma lateral en capas heterogéneas con variaciones
considerables en las permeabilidades verticales.

3.4. Propiedades Petrofísicas

Las propiedades petrofísicas se determinan en el laboratorio con pequeños núcleosobtenidos del

yacimiento y que se procuran sean representativos. Para asegurar una mayorprecisión en estos
datos se puede obtener información complementaria de estaspropiedades. Dicha información la
proporcionan los registros eléctricos y los análisis deprueba de presión. Además, existen
correlaciones publicadas para la obtención de estaspropiedades y pueden ser de utilidad en
determinado momento.

Los datos petrofísicos que se necesitan para efectuar una simulación son:

 Porosidades, φ
 Permeabilidades, k
 Saturaciones de agua, petróleo y gas, Sw, So, Sg
 Presión capilar entre diferentes interfaces, Pcw-o, Pcg-o, Pcg-w
 Permeabilidad relativa al agua, petróleo y al gas, krw, kro, krg
 Compresibilidad de la formación, cr

Simulación de Reservorio
32
3.4.1.- Porosidad

Es una de las propiedades fundamentales de la roca que se debe incluir en un modelo de

reservorio. La porosidad es la fracción de espacio nulo de un medio. Este espacio nulo está
conectado, esto es referido a su porosidad efectiva, de otra manera el espacio nulo es porosidad
inefectiva.

La porosidad puede también definirse como el radio del volumen del espacio poroso y el
volumen aparente de la roca reservorio. Es un parámetro a dimensional que puede ser expresado
en fracción o porcentaje.

La porosidad, es usualmente denotada por ϕ es la fracción de volumen vacío del medio, que es
0 1. La porosidad puede ser dependiente de la presión, pero en muchos casos se asume
que ϕ es constante.

Aun cuando llamamos medio poroso al objeto de nuestro estudio, debemos recordar que es roca.
Casi todas las rocas formadas naturalmente contienen poros, sin embargo, la distribución y la
fracción de volumen de estos poros determinan las propiedades de la roca, las cuales son los
parámetros que gobiernan el flujo de hidrocarburos en el reservorio.

Sin embargo, a pesar de tal simple definición, la porosidad puede ser un parámetro difícil de
cuantificar, siendo el volumen poroso del reservorio generalmente una red completa de espacios
de diferentes tamaños, dimensiones y orígenes. Como consecuencia de esta complejidad muchas
clasificaciones pueden ser consideradas.

Una clasificación general y simple del sistema poroso puede ser basada en los procesos
genéticos responsables de la formación de la porosidad. Desde este punto de vista, se pueden
distinguir 2 tipos fundamentales de porosidad, primaria y secundaria.

3.4.1.1. La Porosidad primaria: es la porosidad original preservada en los sedimentos luego

de la deposición y compactación inicial, es la porosidad original resultante de una deposición
de sedimentos. Es fuertemente dependiente de las características de textura del sedimento
(tamaño, forma y distribución de los granos) y tiende a decrecer con el tiempo y profundidad de
la sedimentación. Algunas veces es preservado por migraciones tempranas de hidrocarburos.

3.4.1.2.La Porosidad secundaria:está relacionada con esfuerzos tectónicos que afectan los
sedimentos después de la deposición y/o la circulación de aguas subterráneas. El primer proceso
tiende a formar fracturas, estilositas y pliegues, mientras que el último es responsable de
disolución, deposición, re cristalización y procesos de dolomitización que afectan a la roca
después de la deposición. La porosidad secundaria es normalmente más importante en rocas
carbonaticos que en sedimentos siliciclásticos, debido a la fragilidad de estos minerales y su

Simulación de Reservorio
33
relativa alta solubilidad. Es definida también como un incremento creciente en la porosidad
primaria debido a la disolución química y esfuerzos tectónicos.

Otra forma simple y popular de clasificar la porosidad de las rocas reservorios toma la imagen
de poros conectados o aislados, lo cual lleva a la distinción entre porosidad total y porosidad
efectiva o interconectada.

Figura 3.1. Muestras de Porosidad

La figura 3.1 ilustra el sistema poroso de una arenisca limpia con una simple historia
diagenética. Las rocas reservorios muestran generalmente cierto grado de porosidad aislada,
debido a la presencia de materiales cementantes que sellan ciertas porciones del espacio poroso.
La característica más importante de estos poros aislados es que estos no contribuyen al flujo de
fluidos.

Cuando se refiere a rocas carbonáticas, la descripción del sistema poroso puede ser aun más
complicada. En la mayoría de los casos, distintos tipos de porosidad coexisten y pueden ser
identificados en la misma roca carbonática.

Las figuras 3.2 y 3.3 muestran valores de porosidad diferente, de acuerdo al ordenamiento de
los granos, los efectos de cementación y sedimentación.

Figura 3.2 Empaque romboédrico de esferas porosidad 25%.

Simulación de Reservorio
34
 Figura 3.3 Efectos de cementación que afectan la porosidad.

Los valores de la porosidad también dependen del tipo de roca como se muestra en la tabla 3.1
hay dos técnicas básicas para medir la porosidad, análisis de núcleo en laboratorio, y registros.
Las mediciones de laboratorio tienden a ser más exactas, pero representan solo una fracción
pequeña del reservorio. Los cambios en las propiedades de las rocas pueden incluso ocurrir
cuando el centro es traído desde el reservorio hacia la superficie. Las mediciones del registro
del pozo representan una porción mucho mayor del reservorio que el análisis de centro, pero
típicamente da resultados menos exactos.

Idealmente una correlación puede ser establecida entre mediciones in situ tales como el registro
de pozo, y mediciones superficiales como en análisis de núcleo.

Rango de Porosidad Porosidad Típica

Tipo de Roca
(%) (%)

Arena Compacta 15-35 25

Arena no consolidada 20-35 30

Carbonatos, Roca Calcárea

intercristalina 5-20 15

Roca Calcárea oolítica 20-35 25

Dolomita 10-25 20

Tabla 3.1 Dependencia de la porosidad con el tipo de roca.

La descripción y la cuantificación de la porosidad es una etapa muy importante en el proceso de

la caracterización de Reservorio.

Simulación de Reservorio
35
3.4.2. Permeabilidades

La permeabilidad es definida como la habilidad de la formación para conducir fluidos. Sin lugar
a dudas, es la propiedad petrofísica más importante de un reservorio. La mayoría de los
parámetros utilizados en la definición económica de los proyectos de desarrollo están en cierta
forma relacionados con la permeabilidad del reservorio.

Desafortunadamente, la permeabilidad es también uno de los parámetros más difíciles de

describir en un estudio de reservorios. De hecho, la permeabilidad en muchos casos representa
una de las paradojas de los estudios convencionales.

La permeabilidad (absoluta), es denotada por K, es una propiedad dependiente de la roca,

definida por la habilidad de la roca de transmitir cierto flujo a ciertas condiciones. Se debe notar
que la permeabilidad no es necesariamente proporcional a la porosidad, puesto que la
orientación y la conexión de poros son esenciales.

La relación entre el tamaño de la abertura y la atracción molecular de la superficie rocosa juegan

un parte importante. La atracción molecular es la fuerza que una capa fina de fluida se adhiera
a la superficie de la roca a pesar de la fuerza gravitacional. La permeabilidad tiene significado
estadístico como una representación de un gran numero de poros. Una medición en escala micro
de la distribución de los granos muestra la diferente granulometría y las diferentes formas de
los granos que afectan la permeabilidad, como se muestra en las figuras 3.3 y 3.4. Usualmente
la permeabilidad disminuye a medida que el tamaño del grano disminuye. Esto puede ser visto
como una conveniencia matemática para describir el comportamiento estadístico de un
experimento de flujo dado.

Figura 3.3 Granos de mayor tamaño, mayor permeabilidad.

Simulación de Reservorio
36
 Figura 3.4 Granos pequeños, menor permeabilidad

Otra característica peculiar de la permeabilidad es que en contrastes con otras propiedades

petrofísicas, es que es una cantidad que tiene dirección y generalmente muestra anisotropía. De
hecho, la permeabilidad depende primeramente de las características de textura de la roca, que
son el resultado de distintos procesos de deposición. El arreglo de los granos que constituyen el
sedimento tiene un impacto importante en las características de flujo. El alineamiento y la
orientación de los granos, la presencia de horizontes de arcilla y de ciertas características
laminares, como ser estilolitas o fracturas son solamente características de las rocas reservorios,
responsables por la existencia de la anisotropía de permeabilidad.

Figura 3.5 Anisotropía de la Permeabilidad.

La unidad de la permeabilidad es el Darcy (D), o mili-Darcy (mD), y los valores típicos varían
en diferentes ordenes de magnitud, valores de 1mD - 10D son usuales en reservorios típicos. La

Simulación de Reservorio
37
definición precisa de un Darcy envuelve la transmisión de un fluido con 1cp de viscosidad en
una roca homogénea a una velocidad de 1cm/seg debido a un gradiente de presión de 1 atm/cm
Traduciendo esto a condiciones de reservorio, 1D es relativamente una permeabilidad alta.

Figura 3.6 Definición de la Permeabilidad

El aspecto matemático del comportamiento anisotrópico es que la permeabilidad es de hecho

descrita por un tensor, lo cual significa que la permeabilidad en las diferentes direcciones
depende de la permeabilidad de las otras direcciones. Sin embargo, para un cambio de base, K
puede ser diagonalizada por diferentes clases, por la estructura del reservorio, los ingenieros de
reservorios están usualmente interesados solo en la permeabilidad horizontal y vertical.

La permeabilidad es una constante física que describe el flujo en una muestra dada con un fluido
dado y un conjunto de condiciones experimentales. Si estas condiciones son cambiadas la
permeabilidad deja de ser aplicable. La distribución de permeabilidad y las curvas relativas de
permeabilidad son puestas en el modelo para reflejar el tipo de procesos que ocurren en el
reservorio

Los datos sísmicos pueden ayudar a definir la distribución de la permeabilidad entre pozos si es
que existe una buena correlación entre la amplitud sísmica y las mediciones cualitativas de la
roca que incluye la permeabilidad. No es usual encontrar que la permeabilidad se dirija solo en
un sentido. Cuando un modelo esta siendo diseñado, el equipo de modelado debería tomar en
cuenta la dirección de la permeabilidad. En principio los simuladores toman en cuenta todos
estos efectos. La creencia que la permeabilidad se dirige en una de tres direcciones ortogonales
conocidas como ejes principales del tensor. Este criterio tiene implicaciones que deben ser
consideradas en los modelos de estudio.

En muchos casos la permeabilidad vertical no es medida y debe ser deducida. Una regla para
esto es asumir que la permeabilidad vertical es aproximadamente un décimo de la permeabilidad
horizontal. Esto llega a ser razonable cuando no se llega a tener datos disponibles.

Simulación de Reservorio
38
3.4.3 Saturación

Es el espacio poroso de la roca reservorio esta lleno de fluidos, normalmente hidrocarburo y

agua. El espaciamiento relativo de estos fluidos depende en un gran número de factores que
están relacionados con las propiedades físicas tanto de la roca como los fluidos, así como
también la interacción roca-fluido.

La fracción de volumen ocupado por cada fase es denotada por la Saturación de la fase. De
modo que,

S i 1
todas_las_fases

La determinación de las condiciones de saturación de una formación del reservorio es una de

las tareas más importantes en un estudio de reservorios. De hecho, cada fase puede contener
más componentes. Que es una forma de distinguir entre diferentes hidrocarburos que están en
la fase gaseosa y en la oleica. La fracción de masa de cada componente i en la fase j es denotada
por Cij. En cada una de las fases, la fracción de masa debe alcanzar la unidad, así para los N
componentes tenemos:

N N N

C
i 1
i,g  C
i 1
i,o  C
i 1
i,w 1

Donde g, o y w son los subíndices comúnmente usados para las fases gaseosas, oleicas, y acuosa
(agua) respectivamente.

De hecho, no solo estas condiciones afectan los cálculos del hidrocarburo in situ, si no también
la mecánica de fluido y por lo tanto el desarrollo de producción esperado del campo.
Desafortunadamente, la saturación de fluidos es mas difícil de determinar que la porosidad,
como en la mayoría de los casos estas evaluaciones están sujetas a diferentes fuentes de
incertidumbre.

Un sin numero de técnicas están disponibles para establecer las condiciones de saturación de
una roca reservorio. Algunos de ellos están basados en mediciones directas de la cantidad de
fluidos presentados en el espacio poroso, otros están basados en mediciones indirectas realizadas
en núcleos o en ambiente de pozo. La saturación de los núcleos puede ser obtenida también en
laboratorio a través de mediciones de presiones capilares.

Simulación de Reservorio
39
3.3.4. Permeabilidad relativa

Los modelos de reservorio calculan la saturación como una función del tiempo. Consideremos
el caso de agua desplazando petróleo. Inicialmente el petróleo ocupa el interior del espacio
poroso y el agua se encuentra adyacente a la superficie de la roca. Cuando el flujo comienza, el
agua desplaza al petróleo a través de los espacios porosos interconectados. La medida de la
conectividad de estos espacios es llamada permeabilidad. El petróleo desplazado después del
flujo del agua es residual. Un comportamiento similar es visto por otras combinaciones, por
ejemplo, gas-petróleo, gas-agua y gas-petróleo-agua. El flujo en fases múltiples es modelado
incluyendo las curvas de permeabilidades relativas en el simulador. Un arreglo típico de curvas
de permeabilidades relativas se muestra en la figura 3.7 las curvas de permeabilidades relativas
representan mecanismos de flujo, tales como drenaje o procesos de imbibición.

Además, como la definición de permeabilidad involucra ciertos fluidos, los diferentes

fluidos experimentan diferentes permeabilidades en el mismo ejemplo de roca.

La permeabilidad de una de las fases depende de la saturación de las otras fases en el

mismo ambiente, así como también de la interacción de las fases con las paredes del pozo.

De modo que, nosotros introducimos una propiedad llamada permeabilidad relativa,

denotada por Kri, donde i = g, o, w, la cual es dependiente de la fase y describe como fluye una
fase en presencia de otras. De modo que, en general, nosotros asumimos que:

K ri  K ri ( S g , S o )

La permeabilidad experimental (efectiva) por fase i, es:

K ri  kri K
Es importante notar, que usualmente las curvas de permeabilidad relativa no son lineares, por
tanto la suma de permeabilidades relativas en un lugar específico (con una composición
específica) no es igual a la unidad.

Simulación de Reservorio
40
 Figura 3.7 Permeabilidad relativa Kro y Krw.

Una característica importante con la permeabilidad relativa, es que puede ser dependiente de la
historia. Este efecto, es llamado histéresis, lo que significa que las curvas de permeabilidad
relativa (dependientes de la Saturación) son diferentes para los procesos de drenaje (fase no-
mojante desplazando la fase mojante) e imbibición (fase mojante desplazando la fase no-
mojante, Ej. Inyección de Agua).

Además, para repeticiones de drenaje e imbibición, Ej. Inyección alternada de agua y gas
(WAG), las curvas cambian aun.

Para permeabilidades de dos fases (sistemas donde solo dos fases están presentes) se ha tenido
buenos datos experimentales en la dependencia funcional. Este, sin embargo, no es el caso de
los datos de permeabilidad relativa para tres fases. De modo que, usualmente se usan modelos
basados en las mediciones de dos fases.

La permeabilidad relativa debería obtenerse por experimentos que modelen mejor el tipo de
desplazamiento que se piensa dominara el desempeño del flujo del reservorio. Por ejemplo, las
curvas de imbibición agua-petróleo son representativas del flujo de agua, mientras que las curvas
de drenaje de agua-petróleo describen el movimiento de petróleo en una zona de agua. El equipo
de modelado necesita comprender que las curvas de permeabilidades relativas usadas en un
modelo de flujo son mas representativas del tipo de experimento que fue utilizado para medir
las mismas. Aplicar estas curvas a otro tipo de mecanismos de desplazamiento representa un
error significativo. En la práctica, la permeabilidad relativa es una de las cantidades físicas más

Simulación de Reservorio
41
útiles disponibles para realizar asociaciones históricas. Como consecuencia las curvas que son
ubicadas inicialmente en el modelo de reservorio son a menudo modificadas durante el proceso
de asociación histórica.

3.3.5. Presión Capilar

La presión capilar es la presión diferencial existente a través de la interfase curvada entre dos
fluidos inmiscibles considerados mojante y no-mojante con respecto a la superficie de
confinamiento. La presión capilar es directamente dependiente de la tensión interfacial y del
ángulo de contacto de la interfase con el cuerpo confinante e inversamente proporcional al radio
de curvatura en la interfase.

2Cos
Pc  (1)
r

Donde;

P = presión capilar, dinas/cm2

σ = tensión interfacial entre los fluidos, dinas/cm

ϕ = Angulo de contacto, grados

r = Radio del capilar (o de los poros), cm

La presión capilar ocurre cuando dos fluidos coexisten en el espacio poroso de una roca
reservorio y es definida como la diferencia de presión apreciable en las dos fases

 Presión capilar en las rocas reservorios

La presión capilar puede expresarse como una carga hidrostática

Pc  h(ρw  ρo )g (2)

Donde,

h = Altura por encima del nivel de agua libre, cm

∆ρ = (ρw – ρo) = Diferencia de densidades entre los dos fluidos, gr/cm3

g = Aceleración de la gravedad, cm/seg2.

Pc
h (3)
Δρg

Reemplazando (1) en (3)

Simulación de Reservorio
42
 2σσCos
h (4)
r(ρw  ρo )g

Por lo tanto, para una distribución de capilares la presión dará lugar a una distribución del
ingreso del fluido mojante en los capilares como se muestra en la figura 3.8 la posición relativa
de la elevación capilar esta dada con respecto al nivel de agua libre (FWL), es decir el punto
donde la presión capilar es cero.

Debe recordarse que aunque los conceptos de presión capilar fueron formulados en términos de
tubos finos, su aplicación en la práctica trata con una compleja red de poros interconectados en
el mátrix de la roca reservorio.

FWL

Figura 3.8 Presión Capilar

La tensión capilar o elevación es usada para establecer la zona inicial de transición en el

reservorio. La zona de transición petróleo-agua, por ejemplo, es la zona entre flujo puro de agua
y flujo puro de petróleo. Esto representa que parte del reservorio donde el cien por ciento pasa
a saturación de petróleo con saturación de agua irreductible. Pueden existir zonas de transición
en interfases entre cualquier par de fase inmiscibles.

La tensión capilar es usada primordialmente para determinar contactos iniciales de flujo y zonas
de transición. Incluso es usada en modelos de reservorios fracturados para controlar el flujo de
los fluidos entre la fractura y la roca madre. Si la tensión capilar es omitida las zonas de
transición no son incluidas en el modelo.

La figura 3.9 muestra como el contenido de fluido del bloque es determinado por la localización
del punto medio del bloque relativo al contacto entre dos fases. El punto medio del bloque es
mostrado como un punto en medio de los bloques de la figura 6.11 así, y el punto medio del

Simulación de Reservorio
43
bloque se encuentra por encima del contacto gas petróleo, el bloque entero es tratado como un
casquete (fase simple de gas con saturación de agua irreductible), incluso si en gran parte del
bloque se extiende sobre la columna del petróleo. Una representación mas exacta puede ser
obtenida disminuyendo el espesor de los bloques de la grilla, pero esto a menudo resulta en un
sustancial incremento en el costo de hacer correr la computadora. Los beneficios relativos de
aumentar la exactitud incrementando el costo deben ser considerados cuando se modelen las
zonas de transición.

Figura 3.9. Grillas en el contacto e Iniciación de un modelo.

La inclusión de una zona de transición en el modelo requiere una curva de tensión capilar
específica (Pc) como una función de saturación para cualquier zona de transición que sea
modelada: petróleo-agua, gas-petróleo o gas-petróleo-agua. La altura de la zona de transición
sobre un nivel de agua libre, es proporcional a la tensión capilar e inversamente proporcional a
la diferencia de densidades entre los dos fluidos. La altura de la zona de transición es una función
de saturación por que la tensión capilar depende de la saturación.

La zona de transición petróleo-agua es típicamente la de mayor espesor por que la diferencia

de densidades entre el gas y un liquido inmiscible. La figura 3.9 muestra también la iniciación
de un modelo conteniendo una curva de tensión capilar. Primero, la altura h sobre un contacto
especificado tal como el contacto agua-petróleo, es calculado con Pc y ∆P. la saturación de un
bloque con un punto medio a la altura h sobre el contacto es calculado de la relación tensión
capilar y saturación. Las zonas de transición complican la identificación de los contactos, porque
la definición de contactos no es universalmente aceptada.

Por ejemplo, el contacto agua-petróleo puede ser definido como la profundidad donde la tensión
capilar es 0. La profundidad del contacto agua-petróleo puede ser identificada usando una
repetición de la prueba de formación encontrando el punto de intersección entre la presión de la

Simulación de Reservorio
44
fase de petróleo y la presión de la fase de agua. Por el contrario el contacto agua-petróleo puede
ser definido como el punto mas profundo del reservorio en el cual aun se puede producir petróleo
libre de agua.

Las diferentes definiciones de contactos resultan en diferencias en las zonas de transición del
modelo, así que es importante saber cual definición es aplicable y quien tiene la autoridad de
juzgar la validez del modelo. En algunos casos, puede ser necesario preparar modelos con ambas
definiciones y entregar una definición como la base mientras la otra es vista como una
sensibilidad. La forma apropiada de incluir la tensión capilar en un modelo de estudio es corregir
las mediciones del laboratorio con las condiciones del reservorio. Esto se hace aplicando la
siguiente corrección:

Y cosο 
PC(res)  P(lab)η coor , η coot  res

Y cosο  lab

El problema con la corrección es que requiere datos que raramente son conocidos con certeza,
a saber la tensión interfacial y el ángulo de contacto a condiciones de reservorio.
Aproximaciones alternativas incluyen el ajuste de las curvas de la presión capilar par ser
consistentes con las estimaciones del espesor de las zonas de transición de los registros del pozo,
o asumiendo que el factor del ángulo de contacto queda fuera.

En este caso:

ηcoor  Yres Ylab

Si las mediciones en laboratorio de la tensión interfacial no están disponibles, podrán ser
estimadas de las correlaciones de Macleod-Sugden para componentes simples o de las
correlaciones de Weinaug-Katz para mezclas.

3.3.6.-Propiedades de los fluidos.

Las propiedades de los fluidos son tan importantes para el ingeniero de reservorios como las
propiedades de roca son importantes para el geólogo reservorista. Actualmente, el tipo de fluido
de reservorio es uno de los factores que más influencian el comportamiento de producción de
un reservorio y consecuentemente, la elección del plan de explotación más apropiado y la
infraestructura superficial de separación.

Las propiedades de fluidos de petróleo deben ser cuantificadas en un simulador de reservorio.

El rango de aplicabilidad del simulador del reservorio es definido, en parte, por los tipos de
fluidos que pueden ser modelados usando los algoritmos matemáticos codificados en el
simulador. Por esta razón, vale la pena considerar los tipos generales de fluidos que serán
encontrados en el ambiente del reservorio.

Simulación de Reservorio
45
En realidad, todas las aplicaciones de ingeniería de reservorios, desde las más simples hasta las
más complejas, requieren asumir algunos comportamientos de los fluidos de reservorio referente
a Presión-Volumen y temperatura (PVT). La importancia de estos parámetros, es que permiten
relacionar los volúmenes de producción de gas y liquido observados en superficie con su
correspondiente vaciamiento en el subsuelo, y desde este punto de vista, se puede considerar
una vinculación entre el reservorio y la producción de superficie.

Más que otro parámetro dinámico, las propiedades de los reservorios de hidrocarburos tienen
un impacto muy fuerte en los cálculos de ingeniería. Además, las características PVT cargan
comúnmente con un alto grado de incertidumbre, la cual puede ser relacionada con problemas
en el muestreo (o más general en la representatividad actual de los ejemplos de fluidos
disponibles), la existencia de variaciones verticales y/o laterales de la composición de fluidos
del reservorio, por la influencia de los procesos de migración debido a la alteración o
biodegradación entre otros.

Se asume que el vacío encontrado en el medio poroso es llenado con las diferentes fases
(componentes). Por propósitos de reservorios prácticos, usualmente tres fases son consideradas
(acuosa, oleica, y la fase gaseosa). Sin embargo, debido a condiciones extremas en un reservorio
hidrocarburífero, no es determinado exclusivamente cuales componentes de hidrocarburo
pertenecen a cual fase. 11

3.3.7.- Tipos de fluidos. 12

Una estimación de la composición elemental del petróleo está dada en la siguiente tabla:

Composición Elemental
Carbono 84%-87%
Hidrogeno 11%-14%
Azufre 0.6%-8%
Nitrógeno 0.02%-1.7%
Oxigeno 0.08%-1.8%
Metales 0-0.14%

Tabla 3.2 Composición elemental de un petróleo.

Puede verse de la tabla que los fluidos del petróleo son predominantemente hidrocarburos. Las
moléculas más comunes de hidrocarburos son las parafinas, naftalinas y aromáticas por la
relativa estabilidad de sus moléculas. Una parafina es un hidrocarburo saturado, que es, un
enlace simple entre dos átomos de carbono. Los ejemplos incluyen metano y etano. Las
parafinas tienen la formula general de CnH2n+2 para alcanos alifáticos (de cadena lineal) y CnH2n
para ciclo-alcanos. Los aromáticos son hidrocarburos no saturados con una estructura cíclica

11 Institut Français du Pétrole publications, COSENTINO LUCA, “Integrated Reservoir

Studies’’, 2001, Editions TECHNIP, Pág. 194-195.
12 IBIDEM, Pág. 195-196.

Simulación de Reservorio
46
que tiene múltiples enlaces entre átomos de carbono como el benceno. La estructura cíclica hace
a los aromáticos relativamente estables y no reactivos.

Un diagrama general PVT de una sustancia pura muestra el comportamiento de fases en función
a la presión, el volumen y la temperatura. Los tipos de propiedades de interés de un reservorio
desde el punto de vista de ingeniería pueden ser llevadas a un diagrama PT de comportamiento
de fase como el que se muestra en la figura 3.10. Muchos fluidos de un reservorio no exhiben
efectos térmicos significativos in situ, aunque reservorios condensados en arenas compactas
pueden presentar un gradiente composicional, que puede influenciar el rendimiento en función
de la profundidad de perforación del pozo.

Punto Crítico

Figura 3.10 Comportamiento de Fase.

El diagrama PT incluye la región de dos fases. La línea que separa la región de fase simple es
llamada fase envoltura o curva envolvente. La región de petróleo está a baja temperatura y en la
región de alta presión sobre la curva del punto burbuja separando las regiones de fase simple y
de dos fases. Si consideramos presiones en la región de la fase simple y nos movemos a la
derecha del diagrama dejando incrementar la temperatura hacia el punto crítico encontramos
petróleo volátil.

A temperaturas superiores al punto crítico pero menores que el cricondentérmico, los fluidos
del reservorio se comportan como condensados. El cricondentérmico es la máxima temperatura
a la cual un fluido puede existir en fases gaseosas y liquidas. Cuando la temperatura del
reservorio es superior al cricondentérmico encontramos reservorios gaseosos. Un sumario de
estos tipos de fluidos se da en la tabla 3.2.

Simulación de Reservorio
47
 Vaciamiento de Presión
Relación Gas-Oil en
Tipo de Fluidos Comportamiento en
Separador (MSCF/STB)
Reservorio

Gas Seco Sin líquidos en Restos de gas

Superficie
Gas Húmedo Restos de gas
> 100
Condensado Gas con liquido
3-100
Petróleo
Liquido con Gas
volátil 1.5-3
Petróleo
Liquido con poco Gas
Negro 0.1-1.5
Petróleo Despreciable formación de
Pesado ~0 Gas

Tabla 3.2 Clasificación de los tipos de fluidos.

Consideramos un reservorio conteniendo hidrocarburos que están a presión y temperatura

correspondientes a una fase simple de petróleo negro. Si la presión del reservorio declina a
temperatura constante, la presión eventualmente cruzara el punto de la curva de burbuja y pasara
a la región de dos fases gas/petróleo. Similarmente comenzando con un condensado en fase
simple y dejando que la presión del reservorio decline a temperatura constante, la presión del
reservorio cruzara la curva de rocío y entrara en la región de dos fases. En este caso, una fase
libre de líquido del gas condensado.

Una vez que el líquido cae es muy difícil recuperarlo. Un método de recuperación es el cíclico,
pero la eficacia en la recuperación será sustancialmente menor que 100%. Si dejamos caer la
presión aún más lejos es posible encontrar condensación retrograda para algunas composiciones
de hidrocarburos.

El diagrama PT también se aplica a los cambios de presión y temperatura. En el caso del flujo a
través de la perforación, el fluido se mueve de una temperatura y presión alta a una presión y
temperatura relativamente baja. Como resultado es común ver fluidos que están en fase simple
en el reservorio convertirse en fase doble al alcanzar la superficie. La figura 3.2 es un diagrama
Presión Temperatura que compara cuatro tipos de fluidos iniciales a distintas condiciones. El
fluido de un reservorio puede cambiar de un tipo de fluido a otro dependiendo en como el
reservorio es producido. Un buen ejemplo es una inyección de gas seco en un reservorio de
petróleo negro. La inyección de gas seco incrementa la cantidad relativa de los componentes de
bajo peso molecular en el petróleo negro.

Simulación de Reservorio
48
La tabla 3.2 muestra las diferentes composiciones para típicos tipos de fluidos. El gas seco
usualmente contiene solo componentes de menor peso molecular. El petróleo volátil continúa
añadiendo componentes de mayor peso molecular. La adición de componentes de mayor peso
molecular y la reducción de componentes de menor peso molecular eventualmente rinden un
petróleo negro. Si monitoreamos el contenido de metano (C1), vemos que tiende a decrecer a
medida que los fluidos cambian de gas seco a petróleo negro.

Componentes Petróleo Petróleo Gas y Gas seco

Volátil Condensado
C1 45.62 64.17 86.82 92.26
C2 3.17 8.03 4.07 3.67
C3 2.10 5.19 2.32 2.18
C4 1.50 3.86 1.67 1.15
C5 1.08 2.35 0.81 0.39
C6 1.45 1.21 0.57 0.14
C7+ 45.08 15.19 3.74 0.21
PM C7+ 231.0 178.00 110.00 145.00
Dens. Relativa 0.862 0.765 0.735 0.757
Color del Negro Anaranjado Café Ligero Acuoso
Líquido Verdoso Oscuro

Tabla 3.3 Diferentes composiciones para tipos de fluidos típicos.

3.3.8.- Modelaje del fluido. 13

El ingeniero debe además prestar gran atención al modelo PVT de los fluidos de reservorios y,
en lo que respecta, la integración de todas las fuentes de información puede proveer un medio
para reducir la incertidumbre relacionada con esta fase. En general, tres tipos de fuentes de datos
PVT están disponibles en un estudio de reservorios típico:

1.- Análisis experimental de laboratorio en fondo de pozo o muestras recombinadas de fluidos

de superficie.
2.- Datos de Producción del Campo.
3.- Correlaciones generalizadas.

En general el comportamiento de los fluidos es mejor modelado usando una ecuación cubicas
de estado. (EdE o EOS), utilizadas en simuladores comerciales, en adición a la presión (P),
volumen (V) y temperatura (T), las EdE contienen el gas como constante R y un conjunto de
parámetros ajustables {a, b} los cuales pueden ser funciones de la temperatura. Las EdE de la
tabla 3.4. Son llamadas cúbicas porque rinden una ecuación cúbica para el factor de
compresibilidad Z, Z = PV/RT. En el caso de un gas ideal Z = 1.

13Institut Français du Pétrole publications, COSENTINO LUCA, “Integrated Reservoir

Studies’’, 2001, Editions TECHNIP, Pág. 197-199.

Simulación de Reservorio 49
 Ecuaciones de Estado
a 1/2
RT T
Redlich-Kwong P  
V  b V V  b 

Soave-Redlich- RT aT 
P 
Kwong V  b V V  b 

aT 
Peng-Robinson P
RT
 V  b
V  b V V  b   b

Tabla 3.4 Ecuaciones de Estado

Las ecuaciones de estado son valiosas para representar propiedades de fluidos en muchas
situaciones. Por ejemplo, supongamos que se quiere modelar un sistema en el cual la producción
está dada por más de un reservorio y por más de un solo tipo de fluido. En este caso el simulador
más apropiado sería un simulador composicional porque un simulador de petróleo negro no
proveería la exacta representación del comportamiento del fluido.

Los dos tipos más comunes de modelo de fluidos de reservorio son los modelos de petróleo
negro y los modelos composicionales. Los modelos de petróleo negro están basados en la
asunción de que las propiedades de la fase saturada de dos fases de hidrocarburos (petróleo y
gas) dependen solamente de la presión. Los modelos composicionales incluso asumen dos fases
de hidrocarburos, y permiten la definición de muchos componentes de los hidrocarburos. De
forma contraria al simulador de petróleo negro el cual puede ser pensado como un simulador
composicional de dos componentes, un simulador composicional a menudo tiene de seis a diez
componentes. Los procesos de los simuladores utilizados para los modelos de superficie
requieren hasta 20 componentes o más.

El costo al utilizar un simulador composicional se incrementa dramáticamente porque se

incrementa el número de componentes modelados, pero los componentes adicionales hacen
posible un modelo más exacto para el comportamiento de los fluidos. Las ecuaciones de estado
deben ser usadas para calcular relaciones de equilibrio en un modelo composicional. Esto trae
consigo poner a punto parámetros tales como {a, b} de las EdE de la tabla 5.3. Varias técnicas
de regresión existen para poner a punto las EdE estas técnicas usualmente difieren en la elección
de parámetros de EdE que serán variados para emparejar los datos de laboratorio con las EdE.

Las figuras 3.11 y 3.12 muestran el comportamiento típico de las propiedades del gas y el
petróleo para un modelo de petróleo negro. Las propiedades de la fase gaseosa son el factor de
formación de volumen de gas (Bg), viscosidad del gas (µg), y solución GOR (Rso). Ambas
curvas, saturadas y no saturadas están incluidas como funciones de la presión solamente. Los
cambios de fase ocurren a presiones de saturación. El petróleo de fase simple se convierte en

Simulación de Reservorio
50
petróleo/gas de fase doble cuando la presión cae por debajo de la presión del punto de burbuja
(Pb), y el gas de fase simple se convierte en gas de fase doble condensado cuando la presión cae
por debajo del punto de rocío (Pd).

Figura 3.11. Propiedades del Petróleo.

Figura 3.12. Propiedades del Gas.

Los simuladores corren más efectivamente cuando los datos de las propiedades de los fluidos
son curvas parejas. Cualquier discontinuidad en una curva puede causar dificultades numéricas.
Ordinariamente, las propiedades reales de los fluidos son funciones parejas de presión a
excepción de los puntos donde ocurren transiciones de fases. Como asunto práctico, es
usualmente bueno graficar en PVT para verificar si las curvas son parejas. La mayoría de los
simuladores reducen la no linealidad del factor de formación de volumen de gas Bg usando la
inversa bg = 1/Bg para interpolar propiedades del gas.

Las propiedades del petróleo de un laboratorio deben usualmente estar relacionadas para el uso
en un simulador de petróleo negro. El flujo en el reservorio es un proceso lento correspondiente
a un proceso diferente del de laboratorio. Cuando el petróleo es producido, es enviado a la

Simulación de Reservorio
51
superficie a través de varios regímenes de presión. Las correcciones están diseñadas para
representar más adecuadamente los fluidos a medida que fluyen diferencialmente del reservorio
a la superficie. Las correcciones alteran la proporción de gas-petróleo de solución y el factor de
formación de volumen del petróleo.

Las propiedades del agua también deben ser ingresadas en un simulador. Idealmente las
propiedades del agua deberían ser medidas realizando análisis de laboratorio en muestras de
agua. Si las muestras no están disponibles, las correlaciones son lo suficientemente ordenadas
para describir el comportamiento del agua.

3.3.9.- Muestreo de fluidos.14

Todas las mediciones de laboratorio sobre las propiedades de los fluidos y los análisis
subsecuentes son inútiles si es que no representan adecuadamente los fluidos in situ. La meta
del muestreo de los fluidos es obtener una muestra que sea representativa del fluido original en
el reservorio. A menudo es necesario acondicionar el pozo antes que la muestra sea tomada. Un
pozo es condicionado produciendo cualquier fluido no representativo, tales como lodos de
perforación, o los que están alrededor de la perforación hasta que es reemplazado por el fluido
original del reservorio.

Los parámetros como la Presión, Volumen y Temperatura son los que gobiernan el
comportamiento de producción en un yacimiento de gas. Los análisis PVT de los fluidos del
reservorio consisten en una serie de procesos de laboratorio diseñados para proveer valores
de las propiedades físicas, requeridas en los métodos de cálculos conocidos como cálculos
de balance de materia y otros. Para obtener un buen estudio del comportamiento del fluido
en el reservorio es necesario que la muestra obtenida sea representativa y consistente, ya que
su validación resulta de vital importancia en una adecuada caracterización de reservorios.
Para realizar un buen chequeo de consistencia y validación de datos en una prueba PVT, se
requiere seguir un orden lógico de los procesos de análisis, los cuales son:
1. Balance de masa de la información reportada
2. Chequeo generales de parámetros
3. simulación termodinámica
La información de los estudios PVT es de gran importancia en la identificación de los
mecanismos de recuperación, en el comportamiento de flujo y la simulación composicional.En
este capítulo se revisan los procesos experimentales estándar realizados a muestras de fluido,
hay cinco principales procesos los mismos que son llevados a cabo con muestras de líquido del
reservorio.

3.3.9.1. Toma de Muestras del Reservorio

Las muestras de los fluidos para un análisis PVT deberán ser tomadas al inicio de su producción
del pozo, antes que ocurra una significativa caída de presión o por lo menos cuando la presión
14 Institut Français du Pétrole publications, COSENTINO LUCA, “Integrated Reservoir
Studies’’, 2001, Editions TECHNIP, Pág. 199-200.

Simulación de Reservorio 52
del reservorio esté por encima de la presión de roció y este en un sistema monofásico, ya que
existen reservorios que están cerca del punto crítico porque con una pequeña caída de presión
se da un aumento rápido de gas o liquido entrando a la región de dos fases, y cuando la presión
de reservorio haya declinado por debajo de la presión de roció, no es posible obtener muestras
representativas del fluido original. Por lo tanto para tomar una muestra representativa la presión
de reservorio deberá ser mayor a la presión de roció. En el caso contrario:
 Si el condensado retrógrado es inmóvil, la muestra presenta una composición más
pobre en componentes pesados que la original del reservorio.
 Si el condensado retrógrado es más móvil, la muestra resultante después de la
combinación daría una presión de roció mayor que la presión del reservorio. La
muestra presenta un contenido de líquido mayor que el original y no es
representativa.
De la misma manera ocurre para los reservorios petrolíferos. Salvando estos problemas,
podemos obtener una prueba representativa.
Los estudios PVT15,16 usualmente son basados en una o más muestras tomadas
durante una prueba de producción, las muestras de fondo del pozo pueden ser obtenidas por
wireline con un recipiente de alta presión durante la prueba de producción o un período de cierre,
alternativamente, las muestras del separador pueden ser tomadas durante una prueba de
producción. Las muestras de fondo del pozo es el método preferido para la mayoría de los
reservorios de petróleo, mientras las muestras recombinadas son tradicionalmente usadas para
reservorio de gas y condensado. Tomar muestras de fondo y de separador en pozo de petróleo
no es poco común. La ventaja de las muestras del separador es que ellas pueden ser
recombinadas en proporciones variadas para alcanzar una deseada presión de burbuja.
3.3.9.2.- Recomendaciones para el muestreo de Pozo
Para obtener un buen muestreo de pozo se deben considerar las siguientes recomendaciones:
1. La producción del pozo debe ser estable
2. Tener un alto índice de productividad y una baja caída de presión alrededor del pozo,
eliminando los daños de formación existente alrededor del pozo
3. Evitar la segregación de los fluidos más pesados al fondo de pozo, teniendo un caudal
adecuado para arrastrar estos fluidos.
4. La RGP, RGC, y la gravedad API deben ser representativas en el reservorio o campo.
5. Cuando el pozo produce agua es recomendable usar separadores trifásicos
6. El muestreo se deben realizar con las perforaciones que se encuentre alejadas de los
contacto ya sea gas/pet, gas/agua, pet/gas, debido a que las muestras no serian
representativas ya que ellas estarían contaminadas.
7. No se deben muestrear en pozos que ya tienen una conificacion de petróleo o agua.

15
API: ”recommended practice for sampling petroleum Reservoir fluids” , API (1966) 44.
16
Reudelhunber, F.O.:”Sampling Procedures for oil reservoir fluids”, J. Pet. Tech. (Dec. 1957) 9, 15-18.

Simulación de Reservorio 53
 8. Tiempo estable sin cambios brusco de temperatura en superficie, ya que ello conlleva el
incremento de condensado y gas. Por ejemplo cuando se tiene una baja temperatura
tenemos la mayor recuperación de líquido en superficie.
9. Cuando un pozo de gas se haya perforado con lodo base aceite es necesario darle un
tiempo de producción para evitar la contaminación del fluido.
3.3.9.3-Tipos de Muestreo
La obtención de la muestra de fluido para un análisis PVT es muy importante ya que de ello
dependen nuestros cálculos de reservas realizados, y la confiabilidad de los mismos,
dependiendo del lugar donde se tomen las muestras, estas se pueden clasificar en la forma
siguiente:
3.3.9.4.-Muestra de fondo
En este método, el pozo es cerrado se baja un equipo especial de muestreo dentro del pozo,
sujetado por un cable de acero (wireline) el fluido en el fondo es muestreado en un recipiente
de alta presión. A esto es llamado muestreo de fondo de pozo o muestra de subsuelo, (en la
figura 3.13, puede observarse un esquema de este método de muestreo).
Se recomienda tomar la muestra con el pozo fluyendo a un bajo caudal estabilizado, la presión
de fondo fluyente en el punto de muestreo debe ser mayor que la presión de saturación; de tal
manera que el fluido mostrado se encuentre en una sola fase. Se deben tomar 3 a cuatro muestras.
A todas estas muestras se le deben determinar la presión de saturación de campo a la temperatura
atmosférica la diferencia entre las medidas no deben ser superior a 30 psi.
El mejor lugar de muestreo es el punto más profundo en el pozo por donde pasa el fluido que
viene de la formación, y donde la presión sea mayor a la presión de saturación, y el nivel de
líquido acumulado debe estar por debajo del punto de control. Experimentalmente, la
composición de una muestra de fondo del pozo es determinada por los siguientes puntos (ver
figura 3.14):
 Flasheando la muestra a condiciones atmosféricas.
 Midiendo los volúmenes de gas, Vg , y del petróleo de superficie, Vo .
 Determinando las fracciones en peso normalizadas, Wgi , y Woi , de las muestras de
superficie por la cromatografía del gas.
 Midiendo el peso molecular del petróleo en superficie, Mo , y la gravedad específica,
GEo .
 Convirtiendo la fracción en peso normalizadas, Wgi a fracción molar normalizadas yi ,
xi .
 Recombinando matemáticamente la composición de la corriente del pozo.

Simulación de Reservorio
54
 Figura 3.13.-Diagrama del muestreo de fondo.

Figura 3.14.-Proceso para recombinar las muestras del separador de la primera etapa, para
obtener la
Composición de la corriente del pozo de una muestra del fondo de pozo; BHS = muestras de
fondo de pozo, GC = cromatografía del gas, FP = Depresión del punto de congelado, y
DM =densitométro.
3.3.9.5.- Muestra de Superficie

En este método la relación de producción como el líquido y gas son cuidadosamente controladas
en el separador. Las muestras son tomadas en el separador de alta, después de largos periodos
de flujo estabilizado. El gas y el líquido deberán ser recombinados en la misma proporción de
los caudales de flujo o de la relación gas petróleo o condensado medida en el separador,
obteniendo una muestra representativa del fluido del reservorio (ver figura 4-4).

Simulación de Reservorio
55
 Figura 3.15.-Diagrama del muestreo de superficie.

Las muestras deben obtenerse al comienzo de las operaciones de producción del reservorio,
preferiblemente en el primer pozo para que de esta forma, la muestra sea representativa del
fluido original que se encuentra en el reservorio. La composición del fluido obtenido en la
muestra depende del acondicionamiento adecuado antes de obtener la muestra. Por lo tanto es
indispensable cumplir con todas las recomendaciones para el muestreo enunciada en el punto
4.3. Kennerly y Reudelhumber2 recomiendan un procedimiento para acondicionar
debidamente el pozo. Un reservorio de petróleo deberá ser muestreado antes de que la presión
de reservorio caiga debajo de la presión del punto de burbuja del líquido. A presiones de
reservorio debajo del punto de burbuja, ningún método de muestreo dará una muestra
representativa de la mezcla del reservorio original. Muestras de fondo del pozo generalmente
contienen menos gas que líquido original, puesto que algo de gas tiene que ser desarrollado.
Hay diferentes compañías en el mundo que presentan distintas recomendaciones para realizar
muestreos en superficie para reservorio de petróleo y gas. Las principales recomendaciones son:
1. Las muestras deben ser tomadas bajo un modelo de flujo crítico para evitar
efectos de contrapresión.
2. Los caudales de muestreo deben ser mayores a los caudales críticos para arrastrar
el agua, lodo y condensado en el fondo de pozo.
3. Las muestras de gas y condensado deben ser tomadas simultáneamente.
4. El separador debe estar operando en condiciones estabilizadas de presión,
temperatura y flujo.
5. Determinar las condiciones del separador durante la toma de muestra como ser
(P, T, RGP, RGC, Qo, Qg).
6. Para el muestreo de gas se recomienda purgar el cilindro o recipiente.
7. Para el muestreo de líquido se recomienda usar técnica de desplazamiento para
evitar alguna contaminación.

Simulación de Reservorio
56
 8. Las muestras deben ser tomadas en el separador de alta o principal. La relación
Gas-Condensado se mide entre el gas de separador y el liquido del tanque.
Corrigiendo este valor por el factor de merma.
9. Considerar el estado del tiempo, para evitar las muchas variaciones en la
temperatura ambiente.
3.3.9.6- Ventaja y Desventaja en diferente tipos de muestreo
3.3.9.7.-Muestreo en Superficie
Ventajas
1. La muestra no se contamina con los fluidos acumulados en el fondo de pozo
2. Menor costo y menor riesgo comparado con el muestreo de fondo.
3. Se pude obtener un mayor volumen de muestra para estudios adicionales.
4. Las muestras son de fácil manejo en el campo
Desventajas
1. Los resultados obtenidos en la prueba dependen de la exactitud de la medición de los
caudales de flujo.
2. Resultados erróneos cuando se tiene separación gas-liquido deficiente
3. Pequeños errores de medición de caudal de flujo, donde la recombinación no generan
muestras representativas.
4. Cambios bruscos de temperatura ambiental en el caso que el fluido fuera gas
3.3.9.8.- Muestreo en Fondo
Ventajas
1. No requiere de medición de caudales de flujo
2. No se tiene problemas de separación gas-liquido
3. Es excelente para el muestreo de gases condensado y petróleo subsaturados, mientras no
se contamine la muestra en el fondo del pozo.
4. No se tiene problemas con los cambios climáticos en superficie en el caso de gas
Desventajas
1. No se toma muestras representativa cuando la presión fluyente es menor que la presión
de saturación.
2. no se recomienda el muestreo cuando el pozo tiene una columna de líquido en el fondo
del pozo.
3. Volumen de la muestra es pequeña
4. Fugas durante la sacada del muestrador hacia superficie.
5. Posibles pesca por ruptura de cable.

Simulación de Reservorio
57
3.3.10.- Determinación de los Volúmenes In –situ

El cálculo de reserva es muy importante ya que nos da una idea general de la reserva
disponible en el yacimiento, conociendo esta reserva podemos saber hasta cuando un
reservorio es económicamente rentable para explotarlo, de ahí la importancia de tener
valores precisos. Por lo tanto en cualquier proyecto a ejecutarse es necesario conocer los
volúmenes de hidrocarburos para ello existen métodos de cálculos dependiendo de la etapa
de maduración del proyecto los cuales mostraremos a continuación los cuales tenemos que
diferenciar las reservas de los recursos las cuales mostraremos de acuerdo a las definiciones
de la SPE-PRMS

Figura 3.16 Clasificacion por oportunidades de Comercializacion según SPE

Figura 3.17 Clasificacion por Criterio de descubrimiento y Comercial según SPE

Simulación de Reservorio
58
Los métodos de Estimación de Reservas o recursos están de acuerdo a la maduración del
proyecto, como se muestra en las Figuras 3.16, 3.17 los cuales mostraremos en la Figura
3.18. según podemos observar en la figura tenemos una etapa temprana o inicio del
proyecto, etapa de exploración o producción temprana cuyos volúmenes determinados
están basado por analogías de otros campos, seguidamente con el desarrollo de campo y
con la información proporcionada de pozos y sísmica podemos realizar, la determinación
volumétrica, cuyos valores tienen bastante incertidumbre principalmente cuando se tiene
un ambiente deltaico, fluvial, marino de poca profundidad donde no se puede definir la
continuidad de los niveles de arena.

El balance de Materia se aplica cuando ya se tiene una producción media, con datos de
levantamiento de presiones de reservorio, siempre y cuando ya se haya producido un
volumen acumulado de un 10 % del volumen In-Situ. El método de declinación es aplicada
en cualquier etapa de los reservorios petrolíferos no así en los reservorios de gas ya que
se tiene un caudal de entrega constante, solamente el método de declinación se aplica
cuando el caudal de gas está por debajo del caudal de entrega. En la parte final del
proyecto.

Para realizar la etapa de simulación de reservorio es necesario contar con suficiente

información de desarrollo como nuevos pozos datos de producción acumulados y datos de
presiones que sustente la valides de la información ya que los registro de presiones hayan
detectados los límites del campo.

Figura 3.18 Métodos de Estimación de Reservas o Recursos según maduración del proyecto

La determinación de las reservas o recursos es muy importante para definir la viabilidad de

cualquier proyecto ya que ello nos da la seguridad para las inversiones dando rangos de
incertidumbre.

Simulación de Reservorio
59
Las reservas son acumulaciones de hidrocarburo que se encuentra en el subsuelo y pueden ser
recuperadas comercialmente y su clasificación está de acuerdo a la Certidumbre de Ocurrencia,
Facilidades de Producción y Métodos de Recuperación lo cual podemos observar en tabla No
3.5

Tabla No 3.5 Clasificación de las Reservas de Hidrocarburo

Desafortunadamente las cifras de las reservas confiables son más necesarias durante las primeras
etapas de un proyecto, cuando se tiene disponible solamente una mínima información. Dado que
la información base es acumulativa durante el desarrollo de un campo, el ingeniero de
yacimientos tiene una creciente cantidad de datos con los cuales trabajar para un proyecto futuro
y maduro, este aumento de datos no solamente cambia los procedimientos para estimar las
reservas, sino correspondientemente mejora la confianza en los cálculos. Es decir:

1) Después del descubrimiento del campo, es necesario estimar las reservas para justificar
muchas decisiones que acarrearían gastos excesivos; de esta manera se podrían planificar las
perforaciones de pozos de desarrollo, plataformas marinas, líneas de tubería e instalaciones
superficiales.

2) Después del desarrollo, las necesidades de capital son altas, pero las reservas estimadas deben
aun soportar las inversiones, nuevas exploraciones, compra de propiedades y el financiamiento
general.

Así, mientras el proceso cambia desde los primeros hallazgos, que van desde la exploración,
continua la necesidad de la estimación de las reservas. Pero la calidad de dichas estimaciones
obviamente varía grandemente entre la etapa especulativa y la efectiva. Antes de que un
yacimiento sea explotado, los datos son muy limitados. Después del descubrimiento, los datos
se acumulan rápidamente; la calidad de la estimación de las reservas también mejora
rápidamente, porque los métodos que pueden usarse se determinan por el tipo y cantidad de
información disponible cuando se hace la estimación.
Las decisiones directivas son dictadas por resultados anticipados de inversión. En el caso del
gas, el ingeniero petrolero compara los costos estimados para algunas oportunidades de
inversión con el flujo de efectivo resultante de la producción.

Simulación de Reservorio
60
3.3.10.1.- Ecuación de Balance de Materia

El balance de materia en reservorios de gas, es un balance de moles de gas existente en el

reservorio. Matemáticamente el balance de materia está representado a través de una ecuación
denominada ecuación de balance de materiales (EBM). Los moles existente en un reservorio en
un determinado instante son la diferencia entre los moles original en el reservorio y los moles
producidos. Como los volúmenes de los fluidos producidos son generalmente medidos a una
determinada condición estándar de presión y temperatura. La ecuación de balance de materiales
es comúnmente escrita de manera que para cualquier instante, el volumen de los fluidos
existente en el reservorio sea la diferencia entre el volumen inicial y el volumen producido,
ambos medidos a esa condición de presión.
La ecuación de balance de materiales para reservorio de gas está sujeto a los tres principales
mecanismos de producción, tales como la expansión de fluido, expansión del gas y la entrada
de agua. Como principales utilizaciones prácticas de la ecuación de balance de materiales se
pueden citar:
 Determinación de volumen original de gas.
 Determinación de volumen original de condensado.
 Determinación de la entrada de agua proveniente de acuíferos.
 Prevención del comportamiento de reservorios.

Inicialmente la EBM será utilizada para la determinación de los volúmenes originales de gas
existente en un reservorio, dato que será la base para los estudios de previsión y comportamiento
futuro utilizando métodos analíticos de entrega de gas.
Conociéndose el volumen de reservorio, la porosidad de roca y la saturación de agua irreducible,
pueden ser calculados los volúmenes originales de gas a través del método volumétrico. El
volumen de gas, medido en condiciones estándar, está dado por:

Vr    Sgi
G Ec. 3.1
Bgi
V    (1  Swi)
G r Ec.3.2
Bgi
Donde Vr es el volumen total del reservorio, φ es la porosidad de la roca. Sg la saturación de
gas, Bg el factor volumétrico de gas, Sw la saturación de agua y el índice i se refiere a las
condiciones iniciales. El factor volumétrico de gas (Bg) , así mismo como otras propiedades de
los fluidos de la roca reservorio pueden ser determinados en laboratorio o estimados a través
de correlaciones empíricas. Como se demuestra en el capítulo 3 de propiedades de los fluidos.
En muchos casos la porosidad, la saturación de agua irresidual y el volumen del reservorio no
son conocidos con la precisión necesaria para el método volumétrico, por tanto no puede ser
aplicado, en esas situaciones la EBM puede ser empleada para calcular el Volumen In-Situ y la
entrada de agua, en el caso que el reservorio tenga el empuje del acuífero.

La ecuación de balance de material puede ser obtenida a partir del principio de conservación de
masa en el interior del reservorio, representado por la siguiente expresión:

Masa producida = masa inicial – masa actual.

Simulación de Reservorio
61
Si la composición del gas producido es constante, los volúmenes producidos y volúmenes
remanente en el reservorio son directamente proporcional a la masa, ya que para la ecuación de
estado de los gases se tiene:
Masa ZRT
V Ec. 3.3
Masa molecular p

Donde Z es el factor de compresibilidad del gas, R la constante universal de los gases, T la

temperatura y P la presión. Por lo tanto el balance de materiales es el balance de moles como se
puede visualizar en la ecuación 3.4
np = ni – n Ec. 3.4

Donde np, ni y n son los números de moles producidos, iniciales y remanente en el reservorio
respectivamente. El número de moles puede ser calculado a través de la ecuación de estado de
los gases reales.
p SC G p pV
PV  ZnRT np  ni  i i
Z SC RT SC Z i RT
Psc  G p Pi  Vi P V Psc P V P V
   Gp  i i 
Z sc RTsc Z i  RTi Z  R  T Tsc Z i  Ti Z  T

Si V=Vi y T=Ti
Psc V  P P Psc T  P P
 G p   i     G p   i  
Tsc T  Zi Z  Tsc V  Zi Z 
P P P T
 i  sc   G p
Z Z i Tsc V
Ec. 3.5
P
 a b  Gp Y = a + bX
Z

Reservorio Anormalmente Presurizado


V  Vi  Cw SW iVPi  C f VPi P  
V  Vi  Cw SW i  C f VPi P 
 
V  Vi  Cw SW i  C f 
Vi
P
 1  Swi   VW i  SW iVPi Vi  VPi S gi
 Cw SW i  C f   Vi
V  Vi 1  P VPi 
 1  Swi   1  S wi  Vi  VPi 1  S wi 
Psc 
P  V P  Vi 1  Cwf P
 Gp  i i 


V  Vi 1  Cwf P  Tsc Z i  Ti Z T
V P P 
Psc

 G p  i  i  1  Cwf P 
Ti  Z i Z

Tsc  Ec.3.6

Simulación de Reservorio
62
La hipótesis de que los efectos de la compresibilidad de la roca y el agua son despreciables no
válidos, debido a que el gas posee compresibilidades mucho mayor que las compresibilidades
de la roca o del agua. Pero pueden ocurrir casos de reservorios anormalmente presurizados
donde la roca reservorio generalmente es inconsolidada, y por lo tanto altamente presurizada,
aparte de eso en reservorios altamente presurizados la compresibilidad del gas es menor de lo
usual.

Los reservorios anormalmente presurizados pueden presentar gradientes de presión que están
en el orden de 0.86 psi/pie, siendo que los valores normales están en torno de 0.433 psi/pie lo
cual significa que la compresibilidad se reduce a la mitad de lo usual. En el primer caso, la
compresibilidad de la roca es de 6114 x10-6 psi-1, comparado con el segundo caso donde los
valores están en el orden de 995x10-6 psi-1 para las rocas consolidadas. Esto significa que sus
efectos no deben ser despreciados, principalmente en el inicio de la vida productiva del
reservorio. Cuando la compresibilidad del gas es menor.

La curva convencional de declinación de presión P/Z vs Gp para yacimiento de gas con

presiones anormales debe ser ajustada teniendo en cuenta las compresibilidades del agua
connota Cw y de la formación Cf con el fin de obtener un valor correcto del volumen In-Situ.
Típicamente un gráfico no corregido de P/Z vs Gp para un yacimiento de gas con presiones
anormales muestra dos pendientes como se muestra en la Figura No 3.19. Una pendiente suave
negativa al comienzo y una más inclinada a presiones menores.

Figura 3.19 Comportamiento de un Reservorio Anormalmente Presurizado

Reservorio Anormalmente Presurizado con entrada de Agua

VW  Cwf * P *Vi V  Vi  Cwf * P *V  We   wWP

Psc   
P  V P  Vi  Cwf * P *Vi  We   wWP
 Gp  i i 

Tsc Z i  Ti Z T

Psc   
P  V P  Vi 1  Cwf P  We   wWP
 Gp  i i 

Tsc Z i  Ti Z T {

Simulación de Reservorio
63
 P P  P
Psc
 
 T  G p  Vi  i  1  Cwf P   We   wWP 
Tsc  Zi Z  Z
 Psc  1 Pi P
  T  G p  We   wWP    1  Cwf P
P
 
 Tsc Z  Vi Z i Z

Psc T
  Gp 
Tsc Vi
P 1
Z Vi
We   wWP   Pi  P 1  Cwf P
Zi Z
 

P We   wWP  Pi Psc T

P
Z

1  Cwf P 
Z
 Vi
    Gp
Z i Tsc Vi

P W   wWP  Pi Psc T

 1  Cwf P  e
Z
      G p
Vi  Z i Tsc Vi
P W   wWP 

 1  Cwf P  e
Z
   a  bG p
Vi  Ec.3.7

P
 a b Gp
Z
Reservorio Volumétrico

P
Z

1  Cwf P  a  bG p 
Reservorio Anormalmente Presurizado

P W   wWP 
Z

 1  Cwf P  e 
  a  bG p
Vi 
Reservorio Anormalmente Presurizado Con entrada de Agua

Simulación de Reservorio
64
 Figura 3.20 Reservorio de Gas con entrada de Agua

3.3.10.2.-.-Determinación del Mecanismo de Producción

Existen tres mecanismos de producción los cuales son responsables de la recuperación de gas
en los yacimientos los cuales enunciaremos de acuerdo a su importancia:
 Expansión del Gas por declinación de la Presión
 Empuje del Agua proveniente del Acuífero activo adyacente
 Expansión del agua connota y reducción del volumen poroso por compactación al
ocurrir disminución de la presión de poros

La ecuación General (EBM) para el balance de materiales considerando los tres mecanismos
anteriormente mencionados se obtiene a partir del siguiente balance:

El volumen de fluido producidos = Gp*Bg+Wp*Bw

El espacio dejado por la producción es llenado por:

 Expansión del Gas
 Expansión del Agua Connata y reducción del volumen poroso
 Intrusión del Agua
Teniendo en cuentas los mecanismos de producción la ecuación de vaciamiento viene expresada
por:
Vaciamiento  G B g  B gi   V w  V p  We Ec.(3.8.)
La expansión del Agua connata, y la reducción del volumen poroso se pueden obtener de las
siguientes ecuaciones generales.
GBGI
Vw  C wVw Pi  P  V w  V p S wi V p  C f V p Pi  P  Vp 
1  SW I
Por lo tanto la suma de la expansión del agua y la reducción del volumen poroso esta expresada
a continuación
C  C S P  P 
GBgi
V w  V p  Ec. 3.9
1  S wi  f w wi i
Combinando las ecuaciones obtenemos el siguiente balance de materiales, considerando los tres
mecanismos de producción

Gp  B gi  B gi  C w S wi  C f  We  W p B w
 1     P  Ec. 3.10
G  B g  B
 g  1  S wi  GB g

 B  Expansión del Gas

1  gi 
 B g 

B gi  C w S wi  C f 
  P Expansión del Agua connata y reducción del volumen
B g  1  S wi 

Poroso

Simulación de Reservorio
65
 We  W p Bw Intrusión del Agua
GBg

3.3.10.3..-Linealizacion de la Ecuación de Balance de Materiales

La ecuación 3.5 del balance de materiales con entrada de agua no está definida como una línea
recta cuando se grafica los valores de p/Z en función de la producción de gas acumulada Gp, lo
que impide la extrapolación para obtener el volumen original de gas Gi. El método de Havlena-
Odeh de la línea recta nos permite la determinación simultánea de Gi y de la entrada de agua
acumulada We.
La ecuación 7.42 puede ser escrita como:
G P B g  W P BW W
G e Ec.3.11
ET ET

C Sw  C f 
E g  B g  Bi E fw  B gi
w
Pi  P  E t  E g  E fw
1  S wi
Cuya forma más general es: y=G+x :

G P B g  W P BW
y Ec.3.12
ET
We
x Ec.3.13
ET
EL gráfico de y vs. x forma una línea recta, con coeficiente angular unitario y coeficiente lineal
igual a Gi. La entrada de agua acumulada We debe ser calculado admitiéndose un modelo para
el acuífero en función al tiempo de producción, los distintos modelos de los acuíferos serán
presentados en el capítulo 8, también se muestra un estudio detallado de los modelos de entrada
de agua. Debido a que si tenemos valores incorrectos de entrada de agua We puede ser
visualizado y analizado.

Figura 3.21.- – Reservorio de gas con entrada de agua

Simulación de Reservorio
66
3.3.10.4.-.-Método Propuesto

Sabemos que el balance de materia del p/z vs producción acumulada en reservorio con entrada
de agua no es aplicable debido a que en la depleción se asume que no tiene entrada de agua. Sin
embargo si un reservorio es afectado por entrada de agua no se tiene una línea recta para
determinar volumen de gas. Por lo tanto debido a estos problemas presentados se presentó una
ecuación para determinar el volumen In-situ.
En cualquier instante del historial de producción, el cálculo del volumen original de gas (Gi)
debe producir el mismo resultado, a diferentes tiempos de la vida productiva de un reservorio
con entrada de agua que mostraremos a continuación:
La ecuación de balance de material p/z con entrada de agua puede ser descrita de la siguiente
manera:
n p  ni  nr
PscGp PiVi Pf (vi  We  ( BwWp )  V fw  Vin)
 
Tsc ZiTr ZfTr
PscGp PiVi PfVi Pf (We  BwWp  V fw  Viny )
  
Tsc ZiTr ZfTr ZfTr

PscGp Pf (We  BwWp  V fw  Viny )  Pi Pf 

    Vi
Tsc ZfTr  ZiTr ZfTr 
PscTrGp Pf (We  BwWp  V  Viny )  Pi Pf 
    Vi
fw

Tsc Zf  Zi Zf 

PscTrGp Pf (We  BwWp  V fw  Viny )  Pf 

 vs  
Tsc Zf  Zf 

De la gráfica determinamos  Pi  para obtener el volumen inicial para diferentes tiempos

 Zi 
de la vida productiva del reservorio obteniéndose el mismo valor original In-Situ Gi
 PscTrGP pf (We  BwWp  V fw  Viny   Pi Pf 
Vi     /   Ec. 3.14
 Tsc Zf   Zi Zf 

3.3.10.5..- Método Presentado por Roberto Aguilera SPE100451

La ecuación de balance de material presentada por Roberto Aguilera toma en cuenta los efectos
de los esfuerzos naturales para los reservorios fracturados ignorando los efectos de
compresibilidad de fractura la cual puede sobre estimar el volumen original in situ.

Simulación de Reservorio
67
Los efectos de la compresibilidad de fractura y matrix están incluido en este estudio usando
estudios de compresibilidad de laboratorios mediantes las muestras. Donde los cambios de
permeabilidad y porosidad pueden ser calculados con las siguientes formulas

Porosidad de la Fractura

Compresibilidad de la Fractura

m  1  2
i 1 i 1
* b Porosidad del Matrix

2 2
vi 1  i 1
 i 1
Coeficiente de Partición
2  m
i 1 i 1
t
i 1

Compresibilidad del Matrix Según Hall

Compresibilidad del Matrix Según Newman

Ec. 3.15

Simulación de Reservorio
68
 Determinacion del Volumen In-Situ Campo A
7000
P/Z Normal
P/Z (1-CwfDp)
6000
P/Z corr agua
Polinómica (P/Z corr agua)

5000
P/z , P/z(1-Cwf*Dp)

4000

3000

2000

1000

0
0 1000000 2000000 3000000 4000000 5000000 6000000 7000000 8000000 9000000
Gas Acumulado MMPC

Fig. No 3.22 Balance de Materiales Reservorio Naturalmente fracturado con entrada de agua

3.3.10.6.- Determinación del Volumen In-Situ para Reservorios de Petróleo

D. Havlena y A. S. han realizado un estudio del análisis de la ecuación de balance de materiales (E.B.M) en forma
de línea recta. En este método se arregla algebraicamente la E.B.M de tal manera que forme una línea recta. Los
grupos de variables se grafican en escala normal. La secuencia de puntos graficados, también como la forma general
del gráfico resultante, es de primordial importancia. En este método el análisis gráfico es esencial, para la
determinación adecuada del volumen original de hidrocarburos.
Inicialmente se arreglará la ecuación de balance de materiales a una línea recta considerando un yacimiento de
petróleo saturado que tenga capa de gas, entrada de agua, e inyección de gas y agua simultáneamente. Luego esta
ecuación se acondicionará según el mecanismo de empuje.
Partiendo de la ecuación general de yacimiento de petróleo tenemos

N
 
Np  Bo  Rp  Rs   Bg  Wp  Bw  We  Giny  Bginy  Winy  Bw iny

Bo - Boi  Rsi  Rs   Bg  m  Boi   Bg  1  1  m  Boi   SwCw  Cf   P
 Bgi   (1 - S w ) 
Ec..(3.16)
Si agrupamos en términos las variables de la Ec..(3.17) tenemos:
El primer término denominado F que agrupa las variables de producción

F  Np  Bo  Rp  Rs   Bg  Wp  Bw  Giny  Bginy  Winy  Bw iny

El Término Eo que expresa la expansión del petróleo y la disolución del gas

Eo  Bo - Boi  Rsi  Rs   Bg

El Término Eg que expresa la expansión de la capa de gas

Simulación de Reservorio
69
  Bg 
Eg  Boi    1
 Bgi 
El término Efw que expresa la expansión del agua inicial y la reducción del volumen del poro.

 S C  Cf 
Ef w  Boi   w w   P
 (1 - Sw ) 
Para efecto de linealización agrupamos y factorizamos E o, Eg y Efw en E

E  Eo  mE g  1  mEf w

En consecuencia tenemos un arreglo de términos en la expresión

F  NE  We Ec..(3.18)
Luego tomando la ecuación de una línea recta y acomodándola en la Ec..(3.18) tenemos:

Y  ax  b
La línea recta pasa por el origen por lo tanto b=0

Y a x
    
F  We N E
En el caso de no existir entrada de agua (We=0), correspondería:

Y a x
    
F N E
La presión del yacimiento declina, lo que hace que los factores volumétricos cambien. También el volumen de
hidrocarburos producidos aumentan y en caso de acuíferos, la entrada de agua al yacimiento aumenta. Estos sucesos
demuestran que desde que se empieza a calcular el volumen original en el yacimiento, en función a los factores
mencionados, los puntos de X, Y graficados si bien muestran una tendencia de línea recta, este conjunto de puntos
graficados no se alinean perfectamente, como se muestra en la fig. 6.1, razón por la cual tiene que tomarse el criterio
de una adecuada linealización que tome la mayor cantidad de puntos en una línea recta.
Para una correcta linealización es imprescindible ver la secuencia de puntos graficados para luego realizar el trazado
de una línea recta. Cuando se trace una recta manualmente y la única incógnita por despejar sea N se tiene que
considerar que la línea pase por el origen, pues la pendiente será el volumen de petróleo original (N). Cuando se
tiene dos incógnitas la intercepción en la ordenada (b) será igual a N, y la pendiente (a) será utilizada para encontrar
el valor de la otra incógnita por despejar, como se verá más adelante.

Simulación de Reservorio
70
 Y

Eo
F

Pendiente = Nm

N
0.0
Eg /Eo X

Fig.11.1 Determinación gráfica del volumen original de petróleo. Y We

Simulación de Reservorio
71
Capitulo IV Clasificación de los Simuladores

4.1.- Introducción

La Simulación de Yacimientos ha sido practicada desde el inicio de la Ingeniería de Petróleos.

En la década de los 40, el potencial de la simulación de yacimientos fue reconocido y muchas
compañías iniciaron el desarrollo de modelos analógicos y numéricos con la finalidad de
mejorar las soluciones analíticas existentes (cálculo de balance de materiales y desplazamiento
1-D de Buckley-Leverett). En la década de los 50, se llevaron a cabo investigaciones en lo que
respecta a solución numérica de ecuaciones de flujo. Como resultado, se obtuvieron programas
de computador para simulación de yacimientos, aunque sencillos pero útiles. Estos programas
representaron el mayor avance y usaron la solución de un conjunto de ecuaciones de diferencias
finitas para describir el flujo multifásico 2-D y 3-D en medios porosos heterogéneos. Fue la
primera vez que los Ingenieros de Yacimientos lograron resolver problemas complejos.

En la década de los 60, el desarrollo de la Simulación de Yacimientos estuvo dirigida a resolver

problemas de yacimientos de petróleo en tres fases. Los métodos de recuperación que fueron
simulados incluían agotamiento de presión y varias formas de mantenimiento de presión. Los
programas desarrollados operaban en grandes computadores (Mainframe) y usaban tarjetas para
el ingreso de datos. Durante la década de los 70, la tendencia cambió bruscamente, debido al
creciente número de investigaciones en procesos EOR, avances en técnicas de simulación
numérica y la disminución del tamaño e incremento de velocidad de los computadores. Los
simuladores matemáticos fueron desarrollados de tal manera que incluían procesos de inyección
química, inyección de vapor y combustión in situ. La investigación durante este período resultó
en avances significativos en lo que respecta a la caracterización de la física del hidrocarburo, en
el desplazamiento bajo la influencia de la temperatura, agentes químicos y comportamiento de
fase multicomponente.

Durante la década de los 80, el rango de las aplicaciones de la simulación de yacimientos

continuó expandiéndose. La descripción de yacimientos avanzó hacia el uso de la
GeoEstadística para describir heterogeneidades y proporcionar una mejor definición del
yacimiento. Se desarrolló la tecnología para modelar yacimientos naturalmente fracturados,
incluyendo efectos composicionales. Así mismo, el fracturamiento hidráulico y pozos
horizontales y su aplicación al monitoreo del yacimiento. Al inicio de la década de los 90, las
aplicaciones fueron hechas en grandes computadores, al final de la década se empezaron a usar
microcomputadores. Actualmente, computadores personales y una gran cantidad de sistemas de
simulación de yacimientos, proporcionan al Ingeniero, un medio económico y eficiente para
resolver complejos problemas de Ingeniería de Yacimientos.

Los avances recientes se han centrado principalmente en los siguientes temas:

1.- Descripción del yacimiento.

2.- Yacimientos Naturalmente Fracturados.
3.- Fracturamiento Hidráulico.
4.- Pozos Horizontales.
Referente a descripción del yacimiento, se están aplicando técnicas estocásticas sustentadas en
lo siguiente:

Simulación de Reservorio
72
a).- Información incompleta del yacimiento en todas sus escalas.
b).- Propiedades de roca variables.
c).- Relación desconocida entre propiedades.
d).- Abundancia relativa de muestras con información proveniente de los pozos.

Referente a yacimientos naturalmente fracturados, la simulación se ha extendido a aplicaciones

composicionales e inyección cíclica de vapor. Respecto a fracturamiento hidráulico, se ha
enfatizado en la predicción de la geometría de la fractura. Se dispone de varias técnicas para
predecir la distribución de los esfuerzos in situ, mejorando de esta forma la simulación del
crecimiento de la fractura en el sentido vertical y lateral.

El objetivo de la simulación de pozos horizontales es estudiar los efectos de longitud del pozo,
ángulo de inclinación, heterogeneidades locales, permeabilidad direccional, barreras y caída de
presión en el pozo. La simulación exacta de los fenómenos cerca al pozo, ha permitido estudiar
los efectos que tienen los pozos horizontales sobre la productividad, intersección de fracturas,
conificación y recuperación de hidrocarburos.

4.2.-Tipos de Simuladores

Los simuladores de yacimientos se pueden clasificar de acuerdo al número de dimensiones, al

tipo de yacimiento que se quiere simular o al proceso particular que se quiere estudiar. Ver
figura 4.1. De acuerdo al número de dimensiones El modelo más simple es el de dimensión cero
o modelo tanque ya que considera el yacimiento como un tanque. En este modelo las
propiedades del fluido y de la roca no varían de punto a punto y considera que el disturbio de
presión inducido en el yacimiento debido a la producción o inyección de fluidos se transmite en
forma instantánea a través de todo el sistema, razón por la cual los cálculos siempre se realizan
a condiciones estáticas.

Simulación de Reservorio
73
 Figura 4.1.- 17Clasificación General de los Simuladores

17
Simulación de Reservorio por Gildardo Osorio

Simulación de Reservorio 74
4.2.1.- Modelo tanque o dimensión cero

Con miras a tener un enfoque más realístico, se pueden unir dos o más de estos tanques,
asignando un valor único a las propiedades de la roca y del fluido a cada tanque y permitiendo
flujo de uno al otro a través de las caras adyacentes. Este constituye un modelo en una
dimensión. Aunque los modelos en una dimensión consideran al yacimiento más detalladamente
que el modelo de dimensión cero, solamente dan una idea general del movimiento de fluidos y
de la distribución de presiones en el yacimiento en función del tiempo. Estos modelos 1D son
usados para estudiar la sensibilidad del comportamiento del yacimiento a las variaciones de los
parámetros del mismo. Por ejemplo, la sensibilidad del petróleo recuperable a la relación de
movilidad, permeabilidad absoluta o la forma de las curvas de permeabilidad relativa.

Estos modelos son raramente usados en estudios de yacimientos para un campo entero, debido
a que no se puede modelar el barrido areal y vertical. Por ejemplo, no se pueden efectuar
cálculos confiables de la eficiencia del desplazamiento en regiones invadidas debido a que no
se puede representar los efectos gravitacionales que actúan perpendicularmente a la dirección
del flujo.

Simulación de Reservorio
75
4.2.2.- Modelo 2D

Los modelos areales cartesianos 2-D (x,y) son los más usados en los estudios de yacimientos.
Se usan principalmente para estudiar el yacimiento entero en casos donde el espesor de la
formación es relativamente pequeño o donde no hay una gran variación vertical en las
propiedades de los fluidos y la formación.

Los modelos areales usan normalmente sistema de coordenadas cartesianas (x,y), sin embargo
existen algunas aplicaciones que requieren sistemas de coordenadas radiales (r,θ ) o cilindricas.
Estos dos últimos sistemas proporcionan una mejor definición cerca a los pozos.

El modelo de geometría horizontal es particularmente útil en la simulación de eficiencias de

barrido, efectos de barreras y evaluación de arreglos geométricos de pozos de inyección. Los
modelos verticales permiten simular la variación vertical de permeabilidad, efectos de
estratificación y efectos de segregación de fluidos.

Simulación de Reservorio
76
El modelo de Secciones Transversales se usan para desarrollar:

a.- Para simular inyección de agua periférica o inyección de gas en la cresta con la finalidad de
proporcionar información sobre la uniformidad de la eficiencia de barrido.

b.- Se usan también para analizar el efecto de la gravedad, capilaridad y fuerzas viscosas sobre
la eficiencia de barrido vertical (conificación de agua).

Si la eficiencia de barrido areal, es un aspecto importante a ser tomado en cuenta, no se debe

usar este tipo de modelo para estimar el comportamiento total del campo.

El modelo Radial es usado para desarrollar funciones de pozo que permitan predecir el
comportamiento cuando se usen en modelos 2-D areales y 3-D y permiten evaluar el
comportamiento de los pozos cuando los efectos verticales dominan el comportamiento como
en el caso de la conificación de agua o gas. Los modelos 2-D radiales son muy usados para
simular la convergencia o divergencia del flujo en una región radialmente simétrica del
yacimiento. Además se usan estos modelos para estudiar el comportamiento de pozos en
yacimientos con empuje de agua de fondo, con capa de gas y yacimientos que tienen una delgada
columna de petróleo y se encuentran rodeados por agua o gas.

4.2.2.- Modelo 3D

Los modelos en tres dimensiones pueden ser cartesianos o radiales. Son los más versátiles ya
que permiten simular la variación de las propiedades de la roca y de las condiciones del fluido
en forma areal y vertical. Sin embargo su aplicación es limitada debido a su alto costo y a que
requiere una caracterización precisa del yacimiento.

Los modelos 3-D son usados donde la geometría del yacimiento es muy compleja como para
ser modelado por un 2-D. Los yacimientos en etapa de depleción avanzada tienen una
dinámica de fluido muy compleja y requieren ser modelados con un simulador 3-D. También
se usan modelos 3-D para simular el desplazamiento de fluidos donde los regímenes de flujo
son dominados por el flujo vertical. Un problema que está asociado a los modelos 3-D es el
tamaño. Un adecuado modelo puede tener tantos gridblock que consumiría mucho tiempo en
proporcionar resultados y retardaría la toma de decisiones

Simulación de Reservorio
77
4.2.3.- Enmallado
La exactitud y eficiencia de un simulador de yacimientos en sistemas complejos depende en
gran parte en la selección adecuada de la malla. Las mallas cartesianas son las más fáciles de
utilizar aunque presentan varias desventajas entre las cuales se destacan:
(a) inflexibilidad para describir fallas, pinchamientos, fracturas hidráulicas, pozos horizontales
y discordancias generales que se presentan en los yacimientos (Fig. 3.9),
(b) inflexibilidad al representar la localización del pozo
(c) inexactitudes inevitables debido a los efectos de orientación de la malla.
Para mejorar las deficiencias de las mallas cartesianas se empleó refinamiento de la malla, y
esquemas de diferencias finitas de 9 y 13 puntos. Los efectos de orientación de la malla han sido
reducidos en cierto grado con le uso de mallas hexagonales (en forma de colmena) La Fig. 3.9
muestra los nuevos tipos de mallas. Actualmente, tienen gran aplicación las mallas con
refinamiento local ya que permiten una mejor descripción del comportamiento de los pozos.

4.3.- Tipo de Yacimientos

En forma general, dependiendo de características físicas producto de la mecánica de las rocas

de los yacimientos, estos pueden dividirse en dos grandes grupos: yacimientos no fracturados y
yacimientos fracturados, siendo los estudios de simulación en estos últimos, los que presentan
mayor grado de dificultad debido a que las fracturas representan verdaderos canales de flujo que
modifican el comportamiento de los fluidos a través del medio poroso. En la actualidad, el
estudio de este tipo de formaciones está recibiendo considerable atención.

4.3.1.- Nivel de Simulación

Los estudios de simulación pueden realizarse a los siguientes niveles:

• Pozos individuales

Simulación de Reservorio
78
• Sector del yacimiento
• Yacimiento
Los estudios de simulación en pozos individuales serían más sencillos que los estudios de
simulación en un determinado sector del yacimiento y más aún en todo el yacimiento; sin
embargo, se debe comentar que existen estudios de simulación para un solo pozo con un grado
de dificultad muy elevado. Estos se ven en los casos de conificación y canalización de agua.

4.3.2.- Tipos de Yacimientos

Para determinar realizar la simulación de reservorio debemos considerar el tipo de yacimiento

en el cual se va realizar el proyecto ya que existen diferentes tipos de simuladores los cuales se
definen de acuerdo al tipo de fluido que contiene el yacimiento y los que utilizan proceso de
recuperación mejorada:

1.- De acuerdo al tipo de fluido en el yacimiento tenemos

 Simuladores de Gas
 Simuladores de Gas Condensado
 Simuladores de Petróleo Negro
 Simuladores Geotérmicos

2.- Simuladores de Recuperación Mejorada

 Simuladores de Recuperación Química

 Simuladores de Recuperación con Miscibles
 Simuladores de Recuperación Térmica

4.3.2.1.- Simuladores de Gas

Este tipo de simuladores se utilizan para llevar a cabo las predicciones del comportamiento de
un yacimiento de gas. Los estudios para este tipo de yacimientos son los más sencillos, si se
considera la presencia de una sola fase (gas). Los parámetros que pueden definirse con este tipo
de simulador son entre otros son: volumen de gas inicial, Caudal de producción y distribución
de presiones.

4.3.2.2.- Simuladores de Gas Condensado

Los Simuladores de Gas condensado o simuladores composicionales se utilizan para simular los
procesos de recuperación para los cuales no sean válidas las suposiciones hechas en modelo de
aceite negro. En esta categoría se incluyen los yacimientos de gas y condensado con
condensación retrograda y los yacimientos de aceite volátil, cuya composición varía
continuamente al existir pequeños cambios de presión y/o temperatura. Este tipo de simuladores
supone en cambio, que los fluidos contenidos en el yacimiento son una mezcla formada por n-
componentes.

Las propiedades de las fases gas – petróleo y su equilibrio se calculan por medio correlaciones
que están en función de la presión y de la composición y más recientemente por medio de

Simulación de Reservorio
79
ecuaciones de estado. Algunos ejemplos de procesos en los cuales son utilizados estos modelos
son los siguientes:

• Agotamiento de un yacimiento de petróleo volátil o de gas y condensado donde la composición

de fase y sus propiedades varían en una manera significativa, con presiones por debajo de la
presión de burbujeo
• Inyección de gas (seco o enriquecido) a un yacimiento de petróleo negro para lograr su
miscibilidad, ya sea total o parcial
• Inyección de CO2 a un yacimiento de petróleo

4.3.2.3.- Simuladores de Petróleo Negro

Este es el modelo más simple que puede utilizarse para estudios de agotamiento primario o
recuperación secundaria por medio de inyección de gas o de agua. Cuenta con los cuatro
mecanismos de desplazamiento básicos para la recuperación de petróleo. Los modelos de este
tipo se han utilizado durante más de treinta años y se basan en la suposición de que los fluidos
del yacimiento pueden representarse de solo tres pseudo componentes (petróleo, gas agua). Esta
posición funciona bien siempre y cuando el sistema durante el proceso de recuperación, quede
lejos del punto crítico y de la región de condensación retrograda y además, si los fluidos que se
inyectan, consiste de los mismos componentes que los fluidos que se encuentren en el
yacimiento. Los modelos de petróleo negro frecuentemente se utilizan para estimar los
siguientes efectos durante la recuperación de petróleo:

• Espaciamiento y arreglo de pozos

• Intervalos disparados
• Conificación del gas y/o el agua como función de la rata de producción
• Caudal de producción
• Mejorar el mecanismo de entrada de agua mediante inyección de la misma y ver la
conveniencia de inyectar por los flancos del yacimiento o inyectar con un arreglo de

4.3.2.4.- Simuladores Geotérmicos

Existen yacimientos cuya energía calorífica se emplea para la generación de energía eléctrica.
Aunque esto no tiene al parecer ninguna conexión con la industria del petróleo, un modelo que
se utilice en este tipo de estudios no puede quedar al margen de una clasificación general de
simuladores, de ahí la razón por lo que se menciona.

4.3.2.5.- Simuladores de recuperación química

En los últimos años se han desarrollados nuevos procesos para recuperar una mayor cantidad de
petróleo, lo cual ha originado la necesidad de contar con simuladores capaces de reproducir el
comportamiento de los yacimientos cuando se someten a este tipo de procesos, tal es el caso de
los simuladores de recuperación química. Dentro de este tipo de métodos de recuperación
mejorada, se pueden citar como los más importantes los siguientes:

• Desplazamiento de petróleo con soluciones miscelares y microemulsiones

• Desplazamiento de petróleo con polímeros

Simulación de Reservorio
80
• Desplazamiento de petróleo con surfactantes
• Desplazamiento de petróleo por combinación de los tres anteriores.

Como es de suponerse, los modelos que se utilizan en este tipo de estudios, presentan un mayor
grado de complejidad pues deben de considerar tanto la interacción que existe entre los propios
fluidos químicos, como la que hay entre dichos fluidos y el medio poroso.

4.3.2.6.- Simuladores de recuperación con miscibles

La miscibilidad es el fenómeno físico que consiste en la mezcla de dos fluidos en todas

proporciones, sin que se forme entre ellos una interfase. Existen diferentes fluidos que se
inyectan al yacimiento bajo esta condición y al estudio del efecto que produce cada uno de ellos
en la recuperación del petroleo se hace con la ayuda de un simulador. Entre los fluidos que se
utilizan en este tipo de procesos se puede citar:

• El gas enriquecido
• El bióxido de carbono, CO2
• El nitrógeno, N2

4.3.2.7.- Simuladores de recuperación Térmica

Este tipo de modelos se utiliza para simular el comportamiento de los yacimientos sujetos a
algún proceso de recuperación mejorada, por medio de métodos térmicos cuyo objetivo
principal es el de proporcionar energía calorífica al petróleo con el fin de disminuir su viscosidad
y de esta forma, facilitar su flujo hacia los pozos productores. Este tipo de métodos puede
clasificarse en:

• Inyección de fluidos calientes, que pueden ser agua caliente o vapor

• Combustión in-situ
• Calentamiento electromagnético

4.4.- Comparaciones de Modelamiento

Los métodos tradicionales para predecir el comportamiento de un yacimiento, se pueden

dividir en tres grandes categorías;

 Métodos analógicos,
 Métodos experimentales,
 Métodos matemáticos.

4.4.1.- Métodos Analógicos

Este método compara las propiedades de yacimientos maduros que tengan la misma ubicación
geográfica o posean condiciones petrofísicas similares al yacimiento en estudio. Antes de
perforar, cuando la información es limitada o ausente, el único método que disponen los
ingenieros para efectuar análisis económicos es la analogía. Para ello se utilizan los yacimientos
con características petrofísicas más parecidas dentro de la misma cuenca o provincia para

Simulación de Reservorio
81
estudiar el comportamiento de otro yacimiento en particular. Este método permite determinar
factores de recuperación, caudales de producción inicial, caudales de declinación, espaciamiento
de pozos y mecanismos de recuperación. Este método proporciona resultados confiables en
yacimientos muy similares con las mismas estrategias de explotación.

Estos se emplean muy raramente en la actualidad, sin embargo, existen dos puntos de vista de
consideración:
1. Históricamente, fueron importantes en los primeros estudios, particularmente la
incorporación de eficiencias de barrido en procesos de inyección de agua,
2. La diferencia entre redes de resistencia-capacitancia y modelos potenciométricos ilustra
la diferencia entre modelos continuos y discretos. Para simular el comportamiento de un
yacimiento, los modelos análogos usan la similitud entre el fenómeno del flujo de fluidos
a través de medios porosos con otros fenómenos físicos (flujo de calor, de electricidad y
flujo de fluidos entre dos capas paralelas o celdas Hele-Shaw).

4.4.2.- Métodos Experimentales

En estos métodos se miden parámetros físicos (caudal, presiones o saturaciones) en modelos de

laboratorio a escala y luego se apropian para todo el yacimiento. Ellos se subdividen en análogos
y físicos:

4.4.3.- Métodos Matemáticos

Estos son el tipo más comúnmente usados por los ingenieros modernos e incluyen balance de
materia, curvas de declinación, estadísticos (correlaciones) y analíticos (pruebas de presión,
Buckley-Leverett, etc.).

En la actualidad, existen dos tipos de modelos físicos.

 El primero no toma en cuenta la geometría de flujo que existe en el yacimiento, como

por ejemplo pruebas de desplazamiento en corazones de roca y pruebas en
empaquetamientos de arena. La desventaja es que los experimentos se corren a escala y
muchas veces no son representativos del yacimiento.
 El segundo tipo utiliza similitudes de conceptos geométricos, mecánicos y térmicos. Esto
es, la geometría, espesor, porosidad y permeabilidad del modelo y las propiedades de los
fluidos se escalan de tal forma que la forma y dimensiones del modelo sean las mismas
del yacimiento.

4.5.- Exactitud de las Soluciones

En cualquier tipo de trabajo de simulación por computador, es importante determinar la

exactitud de la solución generada. En este capítulo, son discutidos varios tipos de errores que
pueden causar errores en la solución.

Simulación de Reservorio
82
4.5.1.- Error de Redondeo

El error por redondeo puede ocurrir cuando se usa una precisión sencilla de exactitud donde es
requerida una precisión doble o una mezcla de variables de precisión sencilla y doble. El error
de redondeo ocurre debido al número de cifras significativas de cada variable definidas en el
simulador. Una precisión simple o doble tiene diferentes significados para diferentes
computadores. Es menos ambiguo indicar la precisión con una declaración tal como REAL*4,
la cual indica que 4 bytes son asignados para el número. Para una máquina de 8 bit, 32 bits o
números binarios serán usados para representar el número. La notación científica es usada. El
primer bit es almacenar el signo positivo o negativo, los pocos bits que quedan son usados por
la mantisa. El número de bits usados por la mantisa determinan la precisión del número.

REAL*4 6-7 cifras significativas

REAL*8 14-15 cifras significativas
REAL*16 30-31 cifras significativas

4.5.2.- Error de balance de materia

El error de Balance de Materia (EBM) es una medida de la consistencia, no de la exactitud. Su

causa principal surge de errores no lineales. Con el método IMPES, la EBM es una indicación
de la forma en que se soluciona la ecuación Ap = b. Un método para resolver , los EBM consiste
en usar una iteración de Newton o una iteración cíclica que eventualmente converja en cero los
EBM. La exactitud es un indicativo de la precisión de la solución. Dos de las formas para
determinar el EBM son:

1. EBM local, solo aplica a pasos de tiempo de producción sencillos que deben igualarse a cero.
La ecuación de balance de materia para el petróleo es:

2.-EBM Acumulado, la ecuación de balance de material para el aceite, es definida para que un
error positivo indique mucho fluido en el yacimiento:

Las ecuaciones de balance de material se conservan tal como están escritas. Sin embargo, las
ecuaciones escritas sobre el papel no se solucionan con exactitud en el simulador. Cuyas causas
son:

1. Ecuaciones no conservativas

Simulación de Reservorio
83
2. Residuos en métodos de matriz iterativa
3. Coeficientes no lineales (coeficiente que cambian con la solución, por ejemplo presión)
4. Errores de redondeo
5. Errores de sintaxis

4.5.3.- Errores no lineales

Los errores no Lineales ocurre cuando se usa una aproximación lineal para encontrar un valor
a un nivel de tiempo n + 1 de una función no lineal tal como son los factores volumétricos de
la formación . La ecuación que se resuelve es la de presión, la cual causa dos errores críticos.
La ecuación de presión con Rs es la siguiente:

Donde

La mayor fuente de error cuando se soluciona la ecuación anterior es el uso de aproximaciones

de la pendiente para cálculos de ct, especialmente cerca del punto de burbuja. Otro problema
es el cálculo del factor volumétrico de formación a un nivel de tiempo n + 1. Esencialmente,
cualquier cálculo a un nivel de tiempo n + 1 introducirá un error. Desde que estos valores no
sean realmente conocidos, los valores anteriores generalmente serán usados.

4.5.4.- Error de Inestabilidad

Los Errores de Inestabilidad es causado por variables dependientes de saturación explicita

en soluciones IMPES. Una solución por inestabilidad puede tomar pequeños pasos
de tiempo o volverse un modelo implícito completamente.

4.5.5.- Errores de Truncamiento

Los Errores de Truncamiento es causado por truncamiento en la serie de Taylor. En una

solución con error de truncamiento, el tamaño del paso de tiempo, ∆t, y el tamaño de la celda,
∆x, varían por ensayo y error hasta que la solución converja. La exactitud de la solución depende
de la magnitud del término de truncamiento. Las series de Taylor

Simulación de Reservorio
84
4.5.6.- Errores de Dispersión numérica

Los Errores de Dispersión Numérica son causados por discontinuidad en la saturación dentro de
una celda. Las soluciones para dispersión numérica son:
1.- Reducir los tamaños de celdas ∆x para celdas grandes
2.- Modificar los cálculos de kr,
3.- Seleccionar apropiadamente ∆t,
4.- usar pseudo kr.
5.-. Orientación de la Malla puede cambiar la solución final. Estos se presentan por la trayectoria
perpendicular o paralela que siguen las partículas durante el proceso. La orientación de la malla
es generalmente importante en cálculos de distribución de saturación en un modelo de inyección.
Típicamente, un sistema de malla diagonal tendrá mejores resultados de recobro (más óptimos).

Simulación de Reservorio
85
CAPITULO V Ajuste Histórico

5.1.- Introducción

El principal objetivo de un estudio de simulación es predecir el comportamiento futuro del

yacimiento con mayor exactitud que alguna otra técnica simple de predicción. Es evidente que
el comportamiento del modelo numérico debe ser similar al del yacimiento para que los
resultados sean aceptables. Debido a la incertidumbre inherente a los datos requeridos para
construir el modelo, se debe probar el comportamiento del modelo antes de ser usado para
predecir el comportamiento futuro.

La única forma de probar el modelo es simular el comportamiento pasado del yacimiento y

comparar los resultados con los datos históricos. El proceso de probar el modelo a través de
comparar el comportamiento pasado es usado también para identificar las inconsistencias del
modelo y corregirlo. El ajuste de historia es, por lo tanto, el proceso de refinar el modelo a través
del ajuste de parámetros de geología, roca y fluido, para producir la mínima diferencia entre los
datos de campo y los resultados del simulador.

Uno de los usos más comunes de la simulación de yacimientos para problemas de campo es el
ajuste histórico. Este es un proceso de estimación de datos de yacimiento mediante datos
arrojados por un simulador, los cuales generan un comportamiento del yacimiento similar a los
datos reales en el campo. Esto, algunas veces es llamado problema inverso. En otras palabras,
se inicia con la solución (datos reales de campo) y se prueba para definir el problema
(descripción del yacimiento).

El ajuste del historial comienza con una definición clara de los objetivos. Dados los objetivos,
es necesario adquirir información de entrada del modelo, especialmente el historial de
producción del campo. Un ejemplo de ajuste del historial es la comparación de resultados
obtenidos en una simulación con los datos obtenidos con mediciones reales, el ajuste es
utilizado para actualizar la simulación de reservorios.

Para lograr un Ajuste del historial exitoso, son utilizados lineamientos generales para el
desarrollo del Ajuste:

1.-Ajustar la presión de Reservorio La cual debe ser ajustada global y localmente. El ajuste
promedio de la presión de campo establece la cantidad global del modelo como un balance de
materiales en conjunto.

2.-Ajustar volumétricamente con el balance de Materiales e identificar el aporte del acuífero.

3.- Ajustar las variables dependientes de la saturación. Estas variables incluyen relación agua-
petróleo (WOR) y la relación gas-petróleo (GOR). El WOR y el GOR son frecuentemente las
variables de producción más sensitivas, en términos de tiempos de avance y la forma de la curva
de WOR o GOR.

Simulación de Reservorio
86
4.- Ajustar las presiones fluyentes de pozo. Una de las tareas esenciales de la etapa de
adquisición de datos es determinar cuál dato debe ser ajustado en el ajuste del historial. La
adquisición de datos es una parte esencial del modelo de inicio.

El modelo inicial es considerado cuando tiene todos los datos que necesita para calcular los
fluidos in place. El reservorio debe ser caracterizado en un formato que se pueda introducir al
simulador y que es aceptado por los comisionarios del estudio. La caracterización del
reservorio incluye la selección de una grilla y datos asociados para el uso en el modelo. Puede
también requerir el estudio de realizaciones múltiples de reservorio en caso de un estudio de un
modelo geoestático. Toda la corrección de datos de fluidos, como también las correcciones
rápidas aplicadas a los datos del diferencial PVT en una simulación de petróleo negro, debe ser
completada durante el proceso de iniciación.

En muchos casos, modelos conceptuales simples pueden ser útiles en la selección de una grilla
final para el modelo de estudio, especialmente cuando se determina el número de capas.
Supongamos que queremos definir el curso del fluido en un movimiento frontal en un campo
muy largo. En este caso, queremos la mayor cantidad de definición de posible. Esto quizás
signifique una pérdida de definición vertical. Una manera de solucionar este problema es poner
uno o más modelos de sección cruzada que representen diferentes partes del campo. Los efectos
de la conformidad vertical en estas regiones son modeladas en detalle por el cálculo

Otro aspecto de la iniciación de un modelo es el equilibrio. Este es el punto donde todo contacto
de fluido es estabilizado y los volúmenes de fluido son calculados. El resultado de los volúmenes
deben ser comparadas con otros estimados de fluido in place, notablemente volumétrico y
estimaciones de balance de material estimado. Allí habrá acuerdo razonable entre los diferentes
métodos. Finalmente el History matching o Ajuste del historial (HM) puede comenzar. 18 19

5.2.- Procedimiento del Ajuste del Historial

El Ajuste del historial es más un arte que una ciencia, porque varía dependiendo de la
complejidad del problema. Sin embargo hay algunas guías generales que pueden ayudar a ajustar
el historial hacia una complementación satisfactoria. Estas guías fueron presentadas por varios
autores como Crichlow (1977), Mattax y Dalton (1990), Thomas (1982), etc. Un conjunto de
guías es presentada en la tabla 5.1 los dos primeros pasos en la tabla toman precedente sobre los
dos últimos.

Si es que los dos primeros pasos no pueden ser logrados, hay un gran chance de que el modelo
sea inadecuado y será necesaria una revisión. Un modelo inadecuado quizás se deba a una
variedad de problemas. Así por ejemplo, si un modelo equivocado ha sido elegido, puede que
el reservorio sea pobremente representado, o que los datos del campo petrolero sean incorrectos
o estén incompletos. Entre las variables de datos ajustados en un estudio típico de petróleo negro

18 SPE textbook, “Gas Reservoir Engineering”, John Lee & Robert A. Wattenbarger, 1996
Pág.267-268.
19 Institut Français du Pétrole publications, COSENTINO LUCA, “Integrated Reservoir

Studies’’, 2001, Editions TECHNIP, Pág.272-273.

87
Simulación de Reservorio
o gas están las presiones, caudales de producción, WOR y GOR, y datos de trazado si es que
estos están disponibles. Estudios más específicos tales como los estudios composicionales o
termales son también datos únicos ajustados al proceso, como la composición del flujo de un
pozo o la temperatura de los fluidos producidos.

El match de presión es usualmente el primer match que será buscado en el

proceso del history matching. Una comparación estimada de las presiones del
fluido obtenidas de exámenes de pozos de un solo pozo en días sucesivos
muestra errores en el reporte histórico de las presiones que pueden ser 10 %
arriba de las presiones de Drawdown. Como una primera aproximación es
suficiente comparar el historial de presiones incorrecto directamente con el
modelo presiones, particularmente si el interés inicial es en las tendencias de
presión y no en los valores actuales de presión. Las correcciones de presión
deben ser aplicadas cuando se desee una sintonización fina del history match.

Paso Recomendaciones

Ajuste volumétrico con balance de materiales

I
identificando el soporte del acuífero.

Ajuste de la presión del reservorio. La presión será

ajustada globalmente y parcialmente. El ajuste promedio
II
de la presión del campo establece la cualidad global del
modelo es como un balance de materiales generales

El ajuste de la saturación depende de variables, las cuales

son: relación agua petróleo (WOR), relación gas petróleo
III
(GOR). WOR y GOR son a menudo las variables más
sensibles

IV ajuste de las presiones fluyentes de los pozos

Tabla 5.1 Recomendaciones para el HM

En la Figura 5.1 mostramos un ajuste volumétrico mediante las presiones de reservorio en la

misma figura se puede observar las distintas variables, como por ejemplo la reducción del NTG
o volumen neto saturado con hidrocarburo y la activación de los flancos de los acuíferos del
reservorio. En la figura 5.1 mostramos el comportamiento de la presión de reservorio en función
del historial de tiempo de producción. Observándose las distintas sensibilidades analizando cuál
de las corridas de simulación es la que más se ajusta al comportamiento real del reservorio una

Simulación de Reservorio
88
ves definido el modelo se procede a comparar la producción de gas del simulador con los datos
históricos de producción.

Los caudales de producción son usualmente de producciones mensuales

grabadas. El modelado especifica unos caudales o presiones de pozo, y luego
verifica que los caudales sean adecuadamente introducidos comparando la
observación de producción acumulativa con el modelo acumulativo de
producción. Después de que los caudales de una fase son especificados, los
caudales de todas las demás fases deben ser ajustados por el desarrollo del
modelo. En muchos casos los caudales observados serán promediados en
bases mensuales o trimestrales con los caudales calculados del modelo. Si la
historia del desarrollo del reservorio es extensiva, por tanto a menudo es
prudente poner gran seguridad en la validez de los datos más recientes del
campo cuando se está desarrollando el history match.

Las relaciones de fases como GOR y WOR, son indicadores sensitivos del
desarrollo del modelo. Las relaciones ajustadas proveen información sobre la
depleción de presiones y movimientos frontales. Los trazados son también
útiles para modelar el frente de los fluidos. Los trazos no necesitan ser
químicos o caros, ellos pueden a menudo ser cambiados en la salinidad del
agua producida. Los cambios de salinidad pueden ocurrir como el resultado
de mezcla, cuando es inyectada, y la salmuera ‘in situ’ tiene diferentes
salinidades, la muestra del análisis del agua de una base periódica es útil para
el rastreo de variaciones de salinidad en función del tiempo.
8000

7000

6000

HM SAL 1
DATA PEB
5000 DATA YPFB
DATA ULTIMA+Sheet2!
DATA JLR

4000 NTG 0.80

NTG 0.75
Flank Aq Props
Flank Aq Pc Krw
3000
Flank Aq Pc Krw Modif
SAL2B0.60
CORRIDA NTG0.55
2000

1000

0
24/7/1998 19/4/2001 14/1/2004 10/10/2006 6/7/2009 1/4/2012 27/12/2014 22/9/2017 18/6/2020

Figura 5.1 Ajuste Volumétrico en base al comportamiento de la presión

Simulación de Reservorio
89
 Figura 5.2.- Ajuste de la Producción del simulador Vs la Producción Histórica

5.3.- Parámetros de Ajuste del Historial

Un concepto fundamental del history matching es el concepto de ‘jerarquía de

incertidumbre’. La jerarquía de incertidumbre es una clasificación del modelo
de calidad de datos para introducir, esto permite que el modelador determine
cual dato es más y cual menos fiable. Las opciones para introducir datos, son
luego reprimidas por el principal, tal que el dato de menor fiabilidad sea
corregido.

La confiabilidad de los datos es determinada cuando son recolectados y

evaluados para su integridad y validez. Esto es como un paso importante para
establecer un sentimiento en el modelador para poder involucrar un dato que
debe ser cuidadosamente revisado. Las curvas de permeabilidades relativas
son a menudo ubicadas en el tope de la jerarquía de incertidumbre por que
son modificadas más a menudo que los demás datos. Las curvas de
permeabilidad relativa son a menudo determinadas desde el centro del fluido.
Como consecuencia la aplicabilidad final del juego de curvas al resto de la
región modelada es siempre una duda.

El volumen inicial de fluido quizás modifique o cambie una variedad de

parámetros de entrada, incluyendo los puntos finales de la permeabilidad
relativa y los contactos del fluido. El modelado computa el volumen original
del fluido en lugares que son restringidos por técnicas independientes como
estudios volumétricos y balance de materia. El esfuerzo para ajustar los datos
del pozo, puede requerir cambiar el intervalo de producción o el índice de
productibilidad de un intervalo de perforación. Si es dificultoso el ajuste del
desarrollo del pozo en una zona o un conjunto de zonas, los modeladores

Simulación de Reservorio
90
necesitaran mirar la variedad de posibilidades, incluyendo la terminación
inesperada y problemas de wellbore.

En un estudio por ejemplo, una inesperada alza de GOR desde el intervalo de

perforación que era conocido por permanecer debajo de los contactos de gas y
petróleo, esta alza era debido al escape de gas en el anillo entre la tubería y la
cañería. Este resultado fue confirmado al hacer circular cemento de atadura y
encontrando una fuga en el intervalo adyacente a la capa de gas y el wellbore.
El gas de la capa de gas adyacente fue encontrado en el wellbore provocando
así una producción de GOR más alta de la esperada. Este efecto puede ser
modelado por una variedad de opciones dependiendo del grado de exactitud
que se desee: por ejemplo puede ser modelado alterando el índice de
productividad (IP) en el modelo del pozo o diseñando un modelo conceptual
cercano al wellbore y preparando curvas de permeabilidad pseudo-criticas.

Esta opción de método influenciara la capacidad productiva del modelo. Así el

modelo de permeabilidad pseudo-relativa permitirá un GOR alto y que el pozo
sea recompletado, donde el IP podrá ser leído correctamente en el tiempo de la
reterminación del pozo para reflejar el mejoramiento en la integridad del
wellbore. Ajustes de mapa quizás sean necesarios. Estos usados para
considerar el último cambio recurrido. Por que los cambios del mapa requieren
un esfuerzo sustancial para re-digitalizar las modificaciones y preparar una
grilla revisada. El Pre-procesamiento de paquete y la ayuda de la computadora
en el modelado geológico están haciendo que los cambios de mapas sean un
método más aceptado en el history matching. En el caso de la geoestática, el
proceso de History matching puede actualizarse involucrando el uso de
diferentes modelos geológicos, llamados imágenes esthoestáticas.

Esto es notable porque provee una guía de cómo los cambios en algunos
parámetros de la historia de pozo (history matching) afectan el ajuste de la
saturación y gradientes de presiones. Un sumario es presentado en la tabla
5.1. Esta muestra por ejemplo que un cambio en el volumen del poro debido
al tiempo puede afectar la presión, como otro ejemplo los cambios en la
permeabilidad relativa son útiles para el ajuste de variaciones de saturación
en tiempo y espacio. Esto es porque los datos de la propiedad del fluido son
raramente cambiados para el ajuste de la historia del campo. Esto es porque
los datos de la propiedad del fluido tienden a ser medidos con más precisión
que los datos del modelo introducido.

El history matching no debe ser realizado haciendo modificaciones incorrectas

en los parámetros. Por ejemplo el ajuste de la presión debe ser realizado por el
ajuste de compresibilidad de roca, y todavía el valor final del ajuste debe ser
sin el conjunto de valores típicos asociados con el tipo de roca en la formación.
En general las modificaciones de los valores de parámetros deben ser
físicamente significativas.

Simulación de Reservorio
91
5.4.- Evaluación del Ajuste

Una manera para evaluar la historia del pozo es comparar, observar y calcular
parámetros. Típicamente, los parámetros observados y calculados son
comparados haciendo graficas de presión vs. Tiempo, presión acumulativa (o
inyección) vs. Tiempo, caudales de producción vs. Tiempo o GOR, WOR, o corte
de agua vs. Tiempo. Otras comparaciones pueden y deben ser hechas si los
datos son válidos y disponibles. Deben incluir por ejemplo modelos de
saturación vs. Saturación de well log y trazado de concentraciones (como la
salinidad) vs. Tiempo. En caso de simulación composicional, los componentes
dominantes (típicamente metano) deben ser graficados como función del
tiempo.

En muchos estudios, el indicador mas sensitivo en el desarrollo del modelo

son las graficas de GOR y WOR o cortes de agua vs. Tiempo. Estos gráficos
pueden ser utilizados para identificar áreas de problema. Por ejemplo
supongamos que graficamos todo alto/bajo WOR y GOR de pozos o graficamos
toda alta/baja presión de pozo, la revisión de estos gráficos pueden revelar una
agrupación de pozos con el mismo problema. Esto puede identificar la
presencia de un error sistemático o una falla en el modelo que necesita ser
corregida. Si la distribución es al azar, puede que las variaciones en el
desarrollo del pozo sean debidas a heterogeneidad que debe ser considerada.

5.5..- Pronostico del Comportamiento

El principal objetivo de un proyecto de simulación es pronosticar el comportamiento. Durante

el ajuste histórico, los caudales son especificados para cada pozo durante todo el periodo
histórico. Los caudales usualmente son desconocidos para el periodo de pronóstico, también
otras condiciones usualmente son especificadas.

La condición más común es especificar la presión de fondo fluyendo, Pwf, y permitir al

simulador calcular los caudales para cada paso de tiempo. El objetivo de los proyectos de
simulación de campo es usualmente comparar alternativas con el propósito de ayudar a tomar
decisiones. Un caso base es corrido usualmente, el cual representa operaciones actuales
continuas. Entonces otros casos son corridos, los cuales representan operaciones alternativas,
tales como perforación de nuevos pozos, aumentar los compresores del campo, estimulación de
pozos, inyección de fluidos (no es común esta discusión para yacimientos de gas seco, excepto
para yacimientos almacenadores de gas), etc. Decisiones de operación se efectúan con base en
los pronósticos de comportamiento y económicos.

Un modelo quizás sea considerado razonable si es que no viola ningún límite

físico conocido. En muchos casos el modelo es aceptado si es que es razonable,
en otras situaciones no solo los límites físicos deben ser satisfactorios, sino
que también los pasos del proceso deben ser aprobados para que el modelo
sea aceptado, de otra manera si la solución está basada en una técnica
innovadora que es razonable pero que no ha sido aprobada el modelo no será

Simulación de Reservorio
92
aceptado. Los modeladores pueden usar un método que está en la literatura
pero si el comisionado del estudio tiene una objeción empírica o filosófica no
lo aprueba. El modelado de un área de ventana es un buen ejemplo de un
método que puede ser razonable pero no es aceptado por qué falla al describir
adecuadamente el flujo a través los bordes del área de ventana y su
rendimiento son resultados pobres. Si alguien en el cargo de autoridad o de
influencia tuvo una mala experiencia con métodos de modelado, este quizás se
rehúse aceptar los resultados de ese modelo. Simultáneamente los
modeladores necesitan estar al tanto de que algunos métodos de modelado no
son universalmente aceptados.

El objetivo principal del manejo de reservorio es determinar las condiciones

óptimas de operación necesitadas para maximizar la recuperación económica
de hidrocarburos. Esto es logrado en principio ordenando los recursos
accesibles para:

 Mejorar la recuperación de un reservorio.

 Minimizar la inversión de capital y gastos operacionales.

El impacto comercial del estudio de simulación es la preparación de una

predicción de fluido neto desde el desarrollo del campo proyectado.
De esta manera el estudio es a menudo completado haciendo predicciones del
desarrollo del campo.

5.5.1.- Desarrollo de las Predicción 20

El desarrollo de las predicciones es valioso para una variedad de propuestas.

Las predicciones pueden ser usadas para interpretar y comprender mejor el
comportamiento del reservorio y estas proveen un significado que determina
modelos sensitivos a cambios en los datos de entrada. Este análisis sensitivo
puede guiar la adquisición de datos adicionales para mejorar el manejo
(administración) del reservorio.

Las predicciones permiten a las personas estimar la vida de un proyecto por

la predicción de recuperación vs. Tiempo. La vida del proyecto no solo depende
del comportamiento del fluido del reservorio, sino también de los problemas
comerciales. Los modelos dejan que el usuario imponga una variedad de
constreñimientos económicos, en el futuro desarrollo del reservorio durante el
proceso de estimación de la vida del proyecto. Estos contrastes son claramente
importantes para el futuro de un proyecto y con problemas técnicos es a
menudo necesario comparar procesos de recuperación diferentes como parte
del estudio, puesto que hay solo un campo es difícil de creer que muchos

20SCHLUMBERGER Publications, “PRODUCTION FORECASTING MANUAL”, 2005, J.

Robert, P. Mc Donald, Pág.4-5.

Simulación de Reservorio 93
procesos diferentes de recuperación puedan ser evaluados en el campo incluso
como pequeña escala de proyecto piloto.

Los proyectos pilotos tienden a ser sustancialmente más costosos de correr

que los estudios de simulación. En algunos casos, sin embargo vale la pena
confirmar el estudio de simulación con un proyecto piloto. Esto es
especialmente verdadero con procesos costosos como el fluido térmico y
químico. Aun otro uso para el modelo de predicción es la preparación de un
plan de manejo de reservorio. Los planes de manejo del reservorio han sido
discutidos en secciones previas. Su preparación es a menudo la motivación
más importante para el desarrollo del estudio de simulación.

Los procesos de predicción comienzan con la calibración del modelo. Es

usualmente necesario asegurar continuidad en el caudal del pozo cuando los
modeladores de control del caudal cambian durante el ajuste del historial para
presionar el control durante la fase de predicción del estudio. Una
discontinuidad en el caudal es observada entre el fin de la historia y el
principio de la predicción. La diferencia de caudal usualmente es producida
por que el índice de productividad actual del pozo, especialmente el efecto skin;
no es adecuadamente modelado por el modelo. Se necesita hacer un ajuste al
modelo para el caudal final del histórico ajuste con el caudal citado
inicialmente.

El siguiente paso es preparar un caso base de predicción. El caso base de

predicción es una condición de proyección de operación existente aplicada. Por
ejemplo: el caso base para un campo nuevamente mejorado que esta
encaminado a una depleción primaria debería ser un caso de depleción
primaria que extienda el uso de un limite de uso especifico. Por contraste si el
campo ha sido circulado con agua el flujo de agua debe ser el caso base y una
alternativa puede incluir inyección de gas y WAG (agua alternada con gas).

Simulación de Reservorio
94
BIBLIOGRAFIA

1.- "Fundamentals of Numerical Reservoir Simulation" - Donald Peaceman - Elsevier Scientific

Publishing Company - 1977; 173 pag.
2.- "Reservoir Simulation" - Calvin C. Mattax and Robert L. Dalton - SPE Monograph Volume
13 - 1990; 161 pag.
3.- "Modern Reservoir Engineering - A Simulation Approach" - Henry B. Crichlow - Prentice
Hall Inc. - 1977 - 354 pag.
4.- “Petroleum Reservoir Simulation”. K. Azis and A. Setari. Elsevier
5.- “Principles of Aplied Reservoir Simularion”. J. Fanchi. Gulf Publishing Company, 1997.
6.- “Principles of Hydrocarbon Reservoir Simulation”. G. W. Thomas. International Human
Resources Development Corporation, 1982.
7.- “Notas sobre simulación numérica de yacimientos”. Gildardo Osorio Gallego, Ph D.
Universidad Nacional, Sede Medellín. Febrero de 2002.

Simulación de Reservorio
95
Capítulo VI Fundamentos básicos de la simulación de reservorio

6.1.-Introduccion

La simulación matemática ayuda a describir, con cierta precisión, el comportamiento de

procesos físicos que ocurren en los yacimientos. Para ello el ingeniero debe identificar dichos
procesos y formular las ecuaciones matemáticas que las gobiernan. Sin embargo, esta tarea no
es nada fácil ya que el flujo de fluidos en medios porosos es un fenómeno muy complejo y para
representarlo se deben considerar ecuaciones que describan el flujo de los fluidos en 1, 2, ó 3
fases, a través de “canales de flujo” que presentan variaciones de uno o varios órdenes de
magnitud en donde los fluidos pueden ser tratados como incompresibles, ligeramente
compresibles o compresibles. Además para representar el sistema de flujo pueden considerarse
1, 2 e incluso 3 dimensiones, incluyendo si se desea, heterogeneidad en las propiedades
petrofísicas, efectos gravitacionales, efectos capilares y transferencia de masa entre las fases.
Con lo anterior es fácil suponer que la habilidad para predecir el comportamiento de un
yacimiento estará en función, primero, de la habilidad que tenga el ingeniero para identificar el
proceso físico que se presenta en el yacimiento y después, para representar dicho proceso en
forma matemática. Las ecuaciones que se emplean en la simulación de yacimientos se obtienen
de la combinación de varios principios físicos como son:

a) Ley de la conservación de masa

b) Ley de la conservación de momento
c) Ley de la conservación de energía (1ª Ley de la termodinámica)
d) Ecuación de flujo (Ley de Darcy) e) Ecuaciones de Estado

6.2.- Energías de los Fluidos en Medios Porosos

Las energías de los fluidos en medios porosos están relacionado con tres tipos de potencial de
energía que son21:

a) Energía gravitacional,
b) Energía de presión,
c) Energía cinética

6.2.1.- Energía Gravitacional

La fuerza gravitacional está definida Fg  mg donde el cambio de energía

gravitacional está dada por:

dEg  Fg dZ  mgdZ Integrando Eg  Eo  mgZ  Z o  E g  mgZ

La fuerza necesaria para levantar el cuerpo a una distancia Z es el peso mg

21
Simulación de Yacimientos principios y conceptos Sepulveda y Escobar

Simulación de Reservorio 96
6.2.2.- Energía de Presión

Para un proceso de desplazamiento de un fluido viene dado como; P= F*A

La diferencial de trabajo está dada: dW  fuerza * dis tan cia  area * presion * dis tan cia
Dónde: Distancia* Área =Volumen dW  VdP

Como la energía disponible de un fluido es equivalente al trabajo

dW  VdP m p dP dP
E p  Eo  m E p  m
p

V p0  p0 

6.2.3.- Energía Cinemática

Describe el movimiento de fluido en el reservorio, el cual toma como base las ecuaciones de
flujo propuesta por Darcy. La cual puede aplicare para expresar la velocidad volumétrica de un
fluido en sus tres dimensiones o ejes

K xp K y p K z p
Vx   Vy   Vz  
ux uy uz

Si ocurre un flujo multifásico, la ley de Darcy puede ser aplicada para expresar la velocidad
volumétrica de flujo en cada fase
K p K p K p
Vxo   x o Vyo   y o Voz   z o  o g
ux uy uz
En forma análoga se aplica para flujo Radial

kkropo
Vxo  
uor

El potencial está definido como:   P   / 144Z esto es estrictamente valido para

fluidos incompresibles. Por lo tanto la gradiente de potencial debería definiré
para fluidos compresibles   P   / 144Z

6.2.3.1.- Ley General de Darcy (Anisotropía)

Si se tiene permeabilidad direccional (Anisotropía) la expresión general es :

k

u  0.00633 
u

Simulación de Reservorio
97
La permeabilidad se convierte en tensor, k, el cual posee forma matricial. Este tensor tiene nueve
valores pero la matriz es simétrica (i.e., k31 = k13), de modo que la matriz tiene solo seis valores
distintos. Para flujo en la dirección x,

k11  k12  k13 

ux    
u x u y u z

Esta expresión puede simplificarse orientando las coordenadas del sistema de flujo (coordenadas
del simulador) a lo largo de los ejes de permeabilidad. Estos ejes son ortogonales entre sí y se
alinean con las permeabilidades máximas y mínimas. Cuando las coordenadas se orientan de
esta manera, entonces cada dirección tiene su propia permeabilidad y el tensor solo tiene tres
elementos diferentes de cero.

El flujo en la dirección x se simplifica a:

K x K x
Vx   Vx  0.00633
ux ux
6.2.3.2.- Flujo Multifásico

La ley de Darcy se extiende a flujo multifásico. La presión, gravedad y efectos capilares se

incluyen en estas ecuaciones. El primer paso es definir los potenciales de las tres fases.
Arbitrariamente, se usa la presión del petróleo como la presión de referencia. Esto conduce a la
inclusión de los términos capilares en las ecuaciones de potencial de gas y agua.

o  g 
o  Po  Z  P o Z g  Pg  Z  P  g Z  Pcog
144 144 144 144

w 
w  Pw  Z  P  wg Z  Pcow
144 144
La presión de petróleo en las ecuaciones anteriores es remplazada por los
términos de presión capilar las cuales son mostradas en la siguiente forma.

 ko   
uo  0.00633  P  o z 
uo  144 
kg    kw   

ug  0.00633  P  g z  Pcog  
uw  0.00633  P  wg z  Pcow 
ug  144  uw  144 

Simulación de Reservorio
98
6.2.4.- Efectos de la Presión Capilar

Los efectos de la presión capilar tiene similar idad con los efectos gravitacionales sobre un
fluido, pero el efecto capilar es menos directo. Es importante saber cuándo existen efectos
capilares en un simulador. El uso más común de los datos de presión capilar está en la
determinación de las saturaciones originales a condiciones iniciales. El equilibro
capilar/gravedad se representan por la traducción de la curva de presión capilar bajo condiciones
de drenaje contra la saturación de agua. La relación entre elevación y presión capilar está dada
por:

 w  o
Pcow  H o  H 
144
Donde H es la elevación del contacto de agua libre

A menudo se nota que es importante simular los efectos capilares en el desplazamiento en un

modelo a escala de campo. Si se ignoran los efectos capilares, se encuentra que el simulador no
proporciona el frente brusco que da la solución de Welge. En vez de ésto, un error de
Modelamiento llamado dispersión numérica resulta en un frente disperso similar al causado por
la inclusión del efecto de la presión capilar. Note que la dispersión física en el yacimiento es
más importante por rocas de baja permeabilidad pero estarán causando efectos de dispersión de
las heterogeneidades.

6.2.5.- Fuerza de Fluido en el Yacimiento

A menudo, analizamos el comportamiento del flujo en términos de fuerzas implicadas por la ley
de Darcy. Nos referimos entonces a las fuerzas viscosas, gravitacionales y capilares. Estas se
expresan usualmente en términos del gradiente de presión:

u 
Fuerza  Vis cos idad  u, psi / pie Fuerza  Gravitacional  , psi / pie
0.0063k 144
dP
Fuerza  Capilares  PC  C S , psi / pie
dS

Simulación de Reservorio
99
Es común analizar los términos de flujo en el yacimiento en términos de relaciones. Recordemos
el número de Reynolds, para flujo en tuberías como una relación las fuerzas inerciales con las
fuerzas viscosas. Un ejemplo de una relación de fuerzas para flujo en yacimientos es el número
gravitacional: la relación entre fuerzas gravitacionales con las fuerzas de viscosidad

 w  o 0.00633k 1
N gravity 
144 u u
6.3.- Ecuación de Continuidad

La ecuación de continuidad en un medio poroso es la relación de la derivada de la porosidad por

la densidad en función del tiempo:

 
   u   
t
El caudal de flujo másico es   n.u ds   u.n ds para obtener la masa total
 

de fluido
 
 
flujo  total.masa.salida   . u.n ds flujo  total.masa. perdidad .dV   dV
s s t
 
.masa.total. perdida   dV Igualando las ecuaciones
v
t
Masa total que sale = masa total perdida

 
. u  n  ds   dV Aplicando el teorema de divergencia de
s
v
t
Gauss
 
 
. u  n  ds     u  dV    u dV

  dV
s
v v v
dt
A medida que dV tiende a cero, los integrando se igualan dando la forma final
de la ecuación de continuidad.
 
 
 u   
t
6.4.- Aplicación Regla de la cadena

Si w=f(u1, u2, …um) y uj=f(x1,x2,…xn)

w w u1 w u2 w um

 *  *  ...  *
xi u1 xi u2 x1 um xi

Simulación de Reservorio 10
0
Simulación de Reservorio 10
1