DESCRIPCIÓN BREVE

Libro de texto de la materia de Modelado de

Yacimientos para Licenciatura y Maestría en
Ingeniería Petrolera.

Sharon Harumi Hernández Vite

MODELADO DE Isabel Melissa Suárez Cruz

José Isidro Cruz Ramírez
Edna Idalia Herrera Martínez
Jesús Eduardo Negrín Núñez
YACIMIENTOS Miguel Alberto Bonfil Oramas

Universidad Olmeca

Revisó: Mtro. Tomás Izaguirre Longoria

 Prefacio

Este libro de texto es la base de la materia de modelado de yacimientos que en su integración

permite de una manera rápida la consulta de la información básica requerida para la
construcción de un modelado de yacimientos petroleros, esto, tomando en cuenta a uno de los
simuladores con mayor impacto y mayormente utilizado en México, eclipse 100 y 300.

1
 Agradecimientos
A La Universidad Olmeca

Son muchos los docentes que han sido parte de nuestro camino universitario, y a
todos ellos les queremos agradecer por transmitirme los conocimientos
necesarios para poder estar aquí. Sin ustedes los conceptos serían solo palabras,
y las palabras ya sabemos quién se las lleva, el viento y por último agradecer a la
universidad que nos ha exigido tanto, pero al mismo tiempo nos permitió obtener
el tan ansiado certificado. Agradecemos a cada directivo por su trabajo y por su
gestión, sin lo cual no estarían las bases ni las condiciones para adquirir
conocimientos. Un trabajo de investigación es siempre fruto de ideas, proyectos
y esfuerzos previos que corresponden a otras personas. En este caso nuestro más
sincero agradecimiento al Ingeniero Tomas Izaguirre Longoria, de la Universidad
Olmeca, con cuyo trabajo estaremos siempre en deuda. Gracias por su
amabilidad, su valioso tiempo y sus ideas para la recopilación de este trabajo.

Agradecimientos a Viviana

Le agradecemos muy profundamente a Viviana por su dedicación y paciencia, sin

sus palabras y correcciones precisas no habríamos podido lograr llegar a esta
instancia tan anhelada. Gracias por su guía y todos sus consejos, los llevaremos
grabados para siempre en la memoria de nuestro futuro profesional.

Sharon Hernández Vite

Isabel Melissa Suárez Cruz
José Isidro Cruz Ramírez
Edna Idalia Herrera Martínez
Jesús Eduardo Negrín Núñez
Miguel Alberto Bonfil Oramas

2
 Antecedentes

Derivado de la problemática para la rápida consulta de bibliografía en la materia de Modelado

de Yacimientos, surge la necesidad de generar un libro de texto integrando los temas de mayor
relevancia e impacto para la materia, este libro de texto, se realizó con la integración de la
investigación de cada uno de los alumnos del tercer semestre grupo “A” de la Maestría en
Ingeniería Petrolera a diciembre de 2022.

Objetivo
Eficientar la explotación de un yacimiento a través del modelado integral del yacimiento para
hacerlo más rentable y evitar errores durante el desarrollo de la perforación, delimitación,
reparaciones Ma y Me, barrido y explotación del yacimiento.

Sírvase como una herramienta de rápida consulta para la integración del modelo de simulación
para efectuar el mismo y analizar los resultados con la finalidad de ajustar los planes de
explotación de los yacimientos.

3
 Contenido
Aspectos Generales ................................................................................................................................... 8
Definición y Objetivos ........................................................................................................................... 8
Utilidad de la Simulación....................................................................................................................... 8
Ajuste de Simulador con la Historia del Yacimiento. ........................................................................... 9
Resultados de una Simulación. ........................................................................................................... 10
Etapas Para Desarrollar Un Modelo........................................................................................................ 10
Como trabaja un modelo......................................................................................................................... 11
Información Requerida Para Utilizar Un Simulador. .............................................................................. 13
Introducción......................................................................................................................................... 13
Descripción física del yacimiento. ....................................................................................................... 14
Mecanismos De Desplazamiento. ....................................................................................................... 14
Propiedades Petrofísicas. .................................................................................................................... 15
Información Adicional. ........................................................................................................................ 16
Datos de Producción y de Relación de Flujo....................................................................................... 16
Estado Mecánico de los Pozos. ........................................................................................................... 17
Aspecto Económico ............................................................................................................................. 18
Mapas................................................................................................................................................... 18
Obtención de Permeabilidades Relativas. .......................................................................................... 20
Geometría de los campos. ....................................................................................................................... 20
Caracterización geológica. ................................................................................................................... 22
Caracterización Petrofísica. ................................................................................................................. 24
Los modelos básicos. ........................................................................................................................... 29
Características del flujo Estacionario. ............................................................................................. 30
Ecuación de continuidad ................................................................................................................. 31
Teorema de existencia y unicidad....................................................................................................... 33
Construcción de las ecuaciones de estado. ........................................................................................ 34
Ecuación de presión. ........................................................................................................................ 35
Ecuación de saturación. ................................................................................................................... 36
Simplificación de los modelos. ........................................................................................................ 37
Fluidos incompresibles dentro de una roca.................................................................................... 38

4
 Formulación matemática de un fluido incompresible. .................................................................. 39
Flujo incompresible en dos fases. ................................................................................................... 40
Conservación de la cantidad de movimiento. ................................................................................ 40
Ecuación de conservación de la masa. ............................................................................................ 41
Ecuación de la presión capilar. ........................................................................................................ 41
Modelo físico final. .......................................................................................................................... 42
Modelo black oil. ............................................................................................................................. 42
El Modelo De Dos Fases Compresible ..................................................................................................... 44
El método de elementos finitos para flujos incompresibles.................................................................. 44
Descripción del método MEF. ......................................................................................................... 44
MEF para flujos incompresibles ...................................................................................................... 45
Resumen del método MEF. ............................................................................................................. 46
Planteamiento matemático de MEF. .............................................................................................. 47
El método de elementos finitos para flujos en dos fases. ................................................................. 49
Flujo en dos fases. ........................................................................................................................... 49
Flujo en dos fases, patrones de flujo. ............................................................................................. 50
Ecuación de Bernoulli. ..................................................................................................................... 51
Bernoulli, velocidad y presión. ........................................................................................................ 52
Construcción de un modelo simple de simulación. ................................................................................ 53
¿Qué es simular? ................................................................................................................................. 53
¿Cuándo simular? ................................................................................................................................ 53
Planeación de un estudio de simulación de yacimientos .................................................................. 55
Revisión de la Información Disponible de Campo. ............................................................................. 56
Modelo PVT de los Fluidos. ................................................................................................................. 56
Revisión de la Información Disponible de Campo. ......................................................................... 56
Modelo PVT de los Fluidos .............................................................................................................. 57
Modelo Petrofísico. ............................................................................................................................. 59
Construcción de la malla de simulación. ............................................................................................ 60
Integración del modelo numérico (data) ............................................................................................ 64
Manual simulador Eclipse. ...................................................................................................................... 68
Objetivo. .............................................................................................................................................. 68
Modelos de Simulación. ...................................................................................................................... 69
Definición de complejos mecanismos de los fluidos del yacimiento. ............................................... 70

5
Tratamiento con petróleo pesado. ..................................................................................................... 71
Estructura de la malla de simulación. ................................................................................................. 71
Diseño de la malla de simulación........................................................................................................ 71
Eclipse 100 ........................................................................................................................................... 72
Secciones del archivo principal o “data file”. ..................................................................................... 75

6
 INTRODUCCIÓN

El presente material, dedicado al modelado de Yacimientos, se efectuó por los estudiantes de 3er
semestre de la Maestría en ingeniería Petrolera, con la intención de documentar lo concerniente al
Modelado Numérico de yacimientos, es un texto dirigido a estudiantes de pregrado y postgrado en
Ingeniería de Petróleos que reúne las experiencias y los conocimientos que los autores han adquirido a
lo largo de sus carreras docentes y de sus estudios de postgrado. El texto, dirigido a estudiantes de
ingeniería y posgrado en Ingeniería de Petróleos, que tratan de resaltar los aspectos más importantes
de la Simulación de Yacimientos, dada la necesidad de describir lo más exactamente posible el
comportamiento de un yacimiento de hidrocarburos. Se presentan los conceptos fundamentales de
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, en el texto, podrán encontrar una unidad dedicada a el software
eclipse, donde en forma de guía, menciona paso a paso como debe operarse el simulador y la
información que se requiere.

7
 Aspectos Generales

Definición y Objetivos

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,
acopladas con ecuaciones representativas de flujo de fluidos, temperatura y/o la concentración
de estos fluidos a través de medios porosos. Dichas ecuaciones 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 de las ecuaciones obedece a lo siguiente:

• 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 un

determinado yacimiento y con base a los resultados obtenidos, optimizar ciertas condiciones
para aumentar la recuperación. Para ello se requiere de la experiencia y buen juicio del
ingeniero para decidir cuando es preciso utilizar un modelo y que tipo de modelo es el más
conveniente en cada caso; así como evaluar de una manera apropiada tanto los datos que se
van a utilizar en la simulación como los resultados que se obtengan de ella. La selección del
modelo a utilizar, además del aspecto económico, está en función de lo que se desea 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.

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
caracterizados y el modelo matemático de simulación ha sido probado y calibrado

8
adecuadamente. Mientras que físicamente el yacimiento puede producirse una sola vez y lo
más probable 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 aceite.

• Tener una buena idea del movimiento de los fluidos dentro del yacimiento.
• Determinar el comportamiento de un campo de aceite bajo diversos mecanismos de
desplazamiento, como puede ser: la inyección de agua, la inyección de gas, el
depresionamiento natural o el uso de algún método de recuperación mejorada.
• Determinar la conveniencia de inyectar agua en un yacimiento de aceite 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 la rata 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.

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 de calidad, el
objetivo de la simulación tenderá a satisfacerse y la predicción del comportamiento será mejor.
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 que se vaya obteniendo mayor información del
yacimiento.

9
Lo primero que se hace para ajustar el simulador con la historia del yacimiento, es calcular el
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, el 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. De esta forma, a lo que en realidad llega el ingeniero es a una combinación
de estas variables (que da como resultado un ajuste), la cual no es única, por lo que dicha
combinación no puede representar de una manera precisa las condiciones del 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 que se tenga en dicho yacimiento. Sin embargo, a medida que el
periodo ajustado sea mayor, la predicción que se haga será más 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.

Resultados de una Simulación.

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

presiones 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-aceite y gas-aceite 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.

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

10
• Predecir el comportamiento futuro
El proceso iterativo mencionado se puede observar en la Fig. 1.1, ya que al ir avanzando en
las 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:
• 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.

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 va a 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 halla
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 el Diagrama 1 se puede observar 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 familiarizando 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.

11
Diagrama 1. Etapas para desarrollar un modelo.

12
 Diagrama 2, Flujo de trabajo de un modelo.

Información Requerida Para Utilizar Un Simulador.

Introducción.

Para que todo trabajo de ingeniería de yacimientos tenga éxito se debe contar con 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 del mismo y de los
tipos de mecanismo por medio de los cuales va a producir. Los resultados que se obtengan de

13
la simulación serán tan buenos como los datos que se hayan empleado para realizarla y, el
tiempo que se pueda perder en preparar esta información es un tiempo bien 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 cambio
significativo 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 se varía
dicha propiedad dentro de este rango y se obtienen resultados parecidos, se puede tomar como
buena una de las predicciones, o bien, relegar a segundo término esfuerzos adicionales para
medir con precisión dicha propiedad. Si, por el contrario, variando esa propiedad se alteran los
resultados considerablemente, es necesario redoblar esfuerzos para obtener 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.

Descripción física del yacimiento.

Para obtener una descripción física del yacimiento es necesario llevar a cabo un estudio
geológico de detalle que proporcione un conocimiento estratigráfico, estructural y
petrográfico, que permita de esta manera caracterizar al yacimiento perfectamente. Dicho
estudio geológico se completa con métodos geofísicos. La información de este tipo que
interesa 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

Mecanismos De Desplazamiento.

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

14
• 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, dando

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

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

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

de echado muy pronunciado) que tengan valores de permeabilidad altos en el sentido vertical
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 y puede

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

Propiedades Petrofísicas.

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

obtenidos del yacimiento y que se procuran sean representativos. Para asegurar una mayor
precisión en estos datos se puede obtener información complementaria de estas
propiedades. Dicha información la proporcionan los registros eléctricos y los análisis de
prueba de presión. Además, existen correlaciones publicadas para la obtención de estas
propiedades 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, aceite y gas, Sw, So, Sg
• Presión capilar entre diferentes interfases, Pcw-o, Pcg-o, Pcg-w
• Permeabilidad relativa al agua, aceite y al gas, krw, kro, krg
• Compresibilidad de la formación, cr

15
propiedades PVT de los fluidos

Las propiedades de los fluidos son también obtenidas en el laboratorio por medio de muestras
sacadas de los pozos. Para que los valores que se obtengan sean aceptables (lo mismo ocurre
con las propiedades petrofísicas), se requiere que las mediciones se hagan lo más
cuidadosamente posible y tratando de acercar al máximo las condiciones del laboratorio a las
condiciones existentes en el yacimiento. Estas propiedades de los fluidos que se requieren en
un trabajo de simulación son:

• Factores de volumen del agua, del aceite y del gas, Bw, Bo, Bg
• Relación de solubilidad en el aceite y en el agua, Rs, Rsw
• Viscosidades del agua, del aceite y del gas, μw, μo, μg
• Compresibilidad del agua, del aceite y del gas, cw, co, cg
• Comportamiento de fases
• Presión de Saturación

Información Adicional.

Además de la información fundamental que se menciona con anterioridad, existen otros

datos irrelevantes cuando se realiza una simulación. Dicha información corresponde casi
en su totalidad como se verá a continuación, a características de los pozos.

Datos de Producción y de Relación de Flujo.

Cuando se trata de hacer un ajuste del modelo con la historia del yacimiento, se requieren
conocer los ritmos de producción y la declinación de la presión. Estos datos de producción
que se necesitan para cada pozo, se pueden desglosar en los siguientes puntos:

• Flujo de aceite vs. tiempo

• Flujo de gas vs. tiempo
• Flujo de agua vs. tiempo
• Cualquier presión media vs. Tiempo

Además, es preciso contar con los índices de productividad y si es el caso, con los índices
de inyectividad de los pozos que integran el yacimiento.

En la práctica generalmente se cuenta con un registro completo de la rata de producción de

aceite de cada pozo, pero no pasa lo mismo con las ratas de producción de gas y de agua,

16
cuya información la mayoría de las veces es limitada. Por ello se necesita que con los datos
disponibles se elabore una gráfica como la que se presenta en la figura 1 que permita
interpolando, obtener una información más completa.

Figura 1. Obtención de los gastos de producción de gas y agua con información incompleta.

Estado Mecánico de los Pozos.

Según lo visto hasta el momento, al parecer para llevar a cabo una simulación, cualquier
información sobre el estado mecánico de los pozos que integran el yacimiento carecería de
interés, pues, aunque los pozos forman parte integral del sistema, la influencia que puedan
tener en él parece haber sido considerada ya en los datos de producción. Además, si la
simulación es un estudio a nivel del yacimiento, ¿para qué sirve entonces la información sobre
el estado mecánico de los pozos? Un avance muy significativo en simulación es acoplar el
comportamiento que tiene los fluidos dentro del yacimiento al que presenta a lo largo de las
tuberías de producción en su camino hacia la superficie. Para ello se requiere contar con el
método de flujo multifásico que entre como subrutina en el simulador. Es de suponer,
lógicamente, que un trabajo de esta naturaleza requiere de las características mecánicas de los
pozos, Fig. 2.2. Existe una gran cantidad de correlaciones que tratan sobre el comportamiento
de los fluidos en las tuberías de producción. El uso de dichas correlaciones, al igual que los
estudios de simulación, está sujeto a ciertas consideraciones importantes.

El estado mecánico de los pozos lo comprende la información siguiente:

• Profundidad máxima del pozo, indicando si es vertical, direccional o desviado, y
horizontal

17
• Diámetro interior del pozo
• Características de las tuberías de revestimiento: diámetro, profundidad, peso y grado
• Características del aparejo de producción. Tipo de terminación: diámetro, grado y peso
de las tuberías de producción. Equipo para métodos artificiales de producción.

Figura 2. Acoplamiento del flujo de fluidos

Aspecto Económico

En todo trabajo de ingeniería debe ser considerado como un punto primordial el aspecto
económico. En simulación de yacimientos la información de este tipo que se debe tomar en
cuenta es la siguiente:

• Precio del barril de aceite

• Costo del pozo
• Límite económico
• Máxima relaciones agua-aceite y gas-aceite con que se piensa trabajar
• Mínima presión de fondo fluyendo
• Precio de gas
• Gastos de operación

Mapas

Al preparar la información que se necesita para realizar una simulación, se elaboran los
siguientes mapas:
• Mapa estructural
• Mapa isópaco
• Mapa de isoporosidades

18
• Mapa de isopermeabilidades

Los mapas estructurales sirven para determinar a través de las curvas de nivel, las
profundidades de los pozos, efectos geológicos del subsuelo como fallas, así como la vista de
planta del yacimiento, límites del mismo, contactos agua-aceite, gas-aceite y/o gas-agua.

Al mapa isópaco lo componen líneas que unen puntos en el yacimiento de igual espesor. Entre
otras cosas sirve para cuantificar volumétricamente el volumen original de aceite y/o el
volumen original de gas. Se comentó al tratar sobre la información petrofísica requerida la
importancia que tenía ésta y la forma de obtenerse. Así pues, las porosidades y las
permeabilidades se conocen en localizaciones discretas del yacimiento y el simulador requiere
un conocimiento de estas propiedades en todos y cada uno de los puntos del mismo. Con este
fin se construyen los mapas de isoporosidades e isopermeabilidades.

En ocasiones se elaboran mapas en los cuales se encuentra la distribución de

combinaciones o productos de propiedades como, por ejemplo:

• Porosidad-espesor, φh
• Porosidad-saturación-espesor, φSoh

Permeabilidades Relativas.

Se ha comentado a lo largo de este capítulo la gran importancia que tiene el contar con una
buena información, haciendo ver que la clave para obtener buenos resultados se basa en buena
información de entrada disponible. Pues bien, sin lugar a duda que la información crítica que se
emplea dentro de toda la información que se requiere al efectuar una simulación, son las
permeabilidades relativas, dado que una relación determinada de ellas define los resultados
que entrega el modelo. Si en el transcurso de un estudio de simulación el yacimiento cambia de
un mecanismo de desplazamiento a otro, el ingeniero debe determinar si este cambio afecta el
desplazamiento de los fluidos o la recuperación dentro de las celdas en que ha sido dividido el
yacimiento y, de ser así, lo que es casi seguro los datos de permeabilidades relativas deben ser
modificados para reflejar dicho cambio. Sin embargo, cambiar el juego de permeabilidades
relativas dependiendo del mecanismo de desplazamiento que opera en el yacimiento por un
tiempo determinado no es nada sencillo, si se considera que la información más difícil de
obtener son precisamente estas curvas de permeabilidades relativas, por la dificultad que existe
de llevar a cabo pruebas de desplazamiento en el laboratorio bajo las mismas condiciones en
que se encuentra el sistema roca-fluidos en el yacimiento. Esto es cierto incluso para
permeabilidades relativas en dos fases. Con lo anterior no se trata de decir que éste sea un
problema particular de la simulación, sino se trata de un problema que se presenta en todo
trabajo que se realiza sobre ingeniería de yacimientos y por el cual es común tener errores
considerables en los resultados obtenidos.

19
Obtención de Permeabilidades Relativas.

Las permeabilidades relativas pueden obtenerse por medio de las siguientes maneras:
• Pruebas de desplazamiento en el laboratorio
• Cálculos por medio de datos de presión capilar
• Ajuste con los datos del campo
• Correlaciones

Geometría de los campos.

La caracterización de yacimientos consiste en detectar y evaluar las características y los

parámetros de la formación que afectan el comportamiento de flujo, entre los cuales se tiene:

Diagrama 4. Parámetros de la formación

El objetivo principal es conocer las características de las rocas productoras de hidrocarburos y

de los fluidos presentes en el yacimiento, tanto cualitativa como cuantitativamente, por lo que,
desde este punto de vista, se consideran dos tipos de caracterización: cualitativa y cuantitativa.

20
Un estudio completo de caracterización de yacimientos involucra diversas etapas y
actividades, entre las que sobresalen las siguientes:
a) Recopilación y validación de la información necesaria.
b) Creación de una base de información o datos (organización).
c) Establecer el modelo geológico-petrofísico estático del yacimiento: Definir el Marco
Estratigráfico-Estructural geológico. Definir el Marco Estratigráfico-Estructural
geofísico (sismológico). Realizar la Evaluación Petrofísica. o Establecer el Modelo
Estático Integral 3D del Yacimiento.
d) Establecer el Modelo Dinámico del Yacimiento: Obtener datos de producción. Realizar
el diagnóstico de geometrías de flujo. Detectar interferencia entre pozos. Estimar
parámetros del yacimiento.

La caracterización estática de yacimientos se puede definir como la disciplina que se encarga

de analizar e integrar la geometría y los límites físicos, convencionales, el contacto agua–
aceite original, así como las variaciones internas verticales y laterales de sus características
petrofísicas, con el fin de precisar la estimación del volumen original de hidrocarburos,
reservas probadas, probables y posibles de los yacimientos.

La información que se procesa en la caracterización estática es:

a) Datos geofísicos: Levantamientos geofísicos, particularmente sísmicos.
b) Datos geológicos: Aspectos geológicos de la roca almacenadora; geometría del cuerpo;
la heterogeneidad del medio poroso; la forma, el tamaño y la orientación del cuerpo
sedimentario.
c) Datos petrofísicos: Los registros geofísicos proporcionan litología, porosidad,
saturación de fluidos y permeabilidad, zonas productoras, determinar profundidad y

21
 espesor de zonas, distinguir entre aceite, gas o agua en el yacimiento y para estimar
las reservas de hidrocarburos.

Figura 2. Información de caracterización.

Caracterización geológica.

Los insumos se pueden obtener principalmente de registros geofísicos de pozos, de muestras

de roca, de pruebas de variación de presión y de datos sísmicos. Con esta información se
obtienen características geológicas del yacimiento, como son: litología, tipo de trampa, tipos de
límites, geometría, profundidad, relieve y cierre estructural.
Comprende tres aspectos principales:
1. Sedimentológicos
2. Estratigráficos
3. Estructurales.

22
 Diagrama 5. Aspectos de la caracterización geológica.

La caracterización sísmica en las áreas de modelación y caracterización de yacimientos se

requiere de una interpretación geológica del yacimiento a mayor detalle, contribuye a la
explotación optima del yacimiento.
La resolución sísmica depende de varios factores, tales como; la longitud de onda, la frecuencia,
la fase y la amplitud, además de que con los datos sísmicos se pueden obtener algunos atributos
sísmicos que se relacionan con propiedades físicas de interés como las fracturas.
Durante el desarrollo de la interpretación sísmica, y debido a que es un método indirecto, se
tiene incertidumbre en el margen de error; desde los parámetros utilizados en el levantamiento
sísmico, el procesado de datos sísmicos, recolección de la información, la identificación de
horizontes (estrato geológico a estudiar), el modelo de velocidad, la conversión tiempo a
profundidad, etc., que se refleja en el resultado final de la estructura geológica y su volumen de
hidrocarburos.

23
 Figura 3. Imagen Sísmica

Caracterización Petrofísica.

Se realiza mediante la integración de datos que se obtienen de núcleos, registros

geofísicos de pozos y pruebas de presión, principalmente. Una correcta
correlación de esta información aporta datos para determinar propiedades físicas
del medio poroso, como son: la porosidad, permeabilidad absoluta, la
compresibilidad, así como también aspectos relacionados con la textura de la
roca, como: la forma, el tamaño, la clasificación y el empacamiento de los granos.

24
 Diagrama 6. Caracterización de la roca.

Una de las características que debe conocerse en un yacimiento, es su habilidad para permitir
el flujo de fluidos a través de él. La cuantificación de esta característica es de suma importancia
para estudios relativos a la explotación de un yacimiento. Esta propiedad del medio poroso
recibe el nombre de "permeabilidad".
LEY DE DARCY Enuncia que "la velocidad de un fluido homogéneo en un medio poroso es
proporcional al gradiente de presión e inversamente proporcional a la viscosidad del fluido “.

La ley de Darcy es una ley estadística que promedia el comportamiento de muchos canales
porosos.

25
Donde:

V= Velocidad aparente (cm/seg)

μ= Viscosidad del fluido (cP)
dp/ds= Gradiente de presión (Amt/cm)
K= permeabilidad

La unidad de permeabilidad es el "Darcy". Se dice que una roca tiene la permeabilidad de un

Darcy cuando un "fluido con una viscosidad de un centipoise avanza a una velocidad de un
centímetro por segundo, bajo un gradiente de presión de una atmósfera por centímetro“.
es:

A continuación, se presenta la obtención de una ecuación que describe el flujo de un fluido

bifásico en un medio poroso.
Ley de Darcy: es una relación empírica entre el flujo de fluidos a través de un medio poroso y
el gradiente de potencial. A condiciones de flujo laminar se define:

Donde:

Vα= Velocidad de la fase α

μα= Viscosidad dinámica del fluido
Pα= Presión del fluido
Krα= Densidad relativa de la fase α
K = Tensor de la permeabilidad.
ρf= Densidad del fluido
g= Gravedad
h= Altura

26
Este desarrollo se hizo tomando en cuenta las siguientes consideraciones:
1. El fluido es homogéneo, en una fase y Newtoniano.
2. No hay reacciones químicas entre el fluido y el medio poroso.
3. Prevalecen condicione de flujo laminar.
4. La permeabilidad es una propiedad del medio poroso y es independiente de la presión,
la temperatura y del flujo de fluidos.
5. No hay efecto de desplazamiento (efecto Klinkenberg).
6. No hay efecto electrocinético.

Ley de Darcy sistemas lineales: la ecuación

se puede integrar cuando se conoce la geometría del sistema a través del cual se mueve el
fluido. La integración se resuelve de la siguiente manera:

El gasto volumétrico de flujo es constante para líquidos, pues la densidad no cambia

significativamente con la presión.

Figura 4. Dirección del flujo.

Esta ecuación es la convencional de flujo lineal utilizada en los cálculos del flujo de fluidos

27
Ley de Darcy sistemas lineales: la ley de Darcy generalizada
pueda emplearse para describir el flujo en cualquier medio poroso donde la geometría del
sistema sea fácil de integrar. La siguiente figura se ilustra el flujo radial que ocurre en la vecindad
de un pozo productor.

Para este modelo de flujo radial la expresión de Darcy es

Figura 5. Expresión de Darcy

Donde q es la tasa de flujo en el yacimiento en cm3/seg; k es la permeabilidad absoluta en

Darcy, h es el espesor en cm, re es radio de drenaje en cm, r es el radio del pozo en cm, Pc es la
presión en el radio de drenaje en atm, Pwf la presión de flujo en el fondo del pozo en atm y la
m es la viscosidad en cp.
Esta última ecuación supone que el yacimiento es homogéneo y está completamente saturado
con una sola fase liquida.

28
Los modelos básicos.

Ecuación de continuidad: es una expresión del principio de conservación de masa. En la

derivación se considera que fluye un fluido en una fase a través del medio poroso.
Es una ecuación de balance de materia en un volumen de control de un sistema
(masa que entra) – (masa que sale) = (acumulación).
La ecuación de continuidad es una expresión del principio de conservación de masa.
El proceso general que se sigue en la derivación de las ecuaciones de continuidad se muestra a
continuación:

a) Se selecciona un volumen elemental, de acuerdo con la geometría de flujo de interés,

puede ser: coordenadas cartesianas, cilíndricas y esféricas.
b) Se expresan todos los gastos que entran y salen del volumen elemental en un periodo
fijo.

c) Para expresar la ecuación de continuidad debe suponerse un medio poroso como se

muestra en las figuras, se muestra un flujo en dirección x, y, z.

Figura 6. Direcciones de flujo.

29
Ecuación de continuidad en coordenadas cartesiana (x, y, z)

Ecuación de difusión, es de movimiento que permite analizar el comportamiento de la presión

y tiempo en un medio poroso. Una ecuación diferencial parcial que describe la variación, en el
tiempo y el espacio, de una cantidad física que es regida por la difusión. La ecuación de difusión
proporciona un buen modelo matemático para la variación de la temperatura a través de la
conducción del calor y de la propagación de las ondas electromagnéticas en un medio
altamente conductivo. La ecuación de difusión es una ecuación diferencial parcial parabólica
cuya forma característica relaciona la primera derivada parcial de un campo con respecto al
tiempo con sus segundas derivadas parciales con respecto a las coordenadas espaciales. Está
íntimamente relacionada con la ecuación de ondas.

Comportamiento cualitativo de la solución en estado de la no difusión y sin presión capilar.

Características del flujo Estacionario.

a) Flujo másico que entra = Flujo másico que sale

b) El sistema responde instantáneamente a los cambios en presión o gasto en cualquier
parte del sistema
c) La presión y el gasto permanecen constantes en algún punto del yacimiento como una
función del tiempo

Difusión: es el movimiento de átomos en un material con el fin de eliminar las diferencias de

concentración y llegar a una composición homogénea y uniforme.

30
La primera Ley de Fick se utiliza en estado estacionario de difusión, es decir, cuando la
concentración dentro del volumen de la difusión no cambia con respecto al tiempo. Esta Ley
afirma que la densidad de corriente de partículas es proporcional al gradiente de concentración.

La constante de proporcionalidad se denomina coeficiente de difusión D, J es el flujo de difusión

en dimensiones de distancia-tiempo. Mide la cantidad de sustancia que pasara por un área por
un intervalo de tiempo.

Para obtener las ecuaciones correspondientes para el flujo en de aceite y agua, donde el agua
es la fase mojante, de adecuar la ecuación siguiente para la fase w y o.

Ecuación de continuidad

Se tiene las siguientes ecuaciones obtenidas

Sustituyendo la ley de Darcy para cada fase y expandiendo las derivadas

31
Donde:
Ø=Porosidad
Vα= Velocidad de la fase α
ρα= Densidad de la fase
Sα=Saturación de la fase
λα= Movilidad de la fase
K = Tensor de la permeabilidad.
Pα= Presión de la fase α
g= Gravedad
qα= Gasto de la fase α

A continuación, se dividen las ecuaciones por sus respectivas densidades

Donde el Qt es el gasto total

Y la compresibilidad de las fases

32
Se define la presión promedio

Definir las movilidades de fase, ⁄o y ⁄w, que serán funciones de la saturación

Sustituyendo las ecuaciones y reacomodando, así como simplificando la ecuación en términos

de la velocidad total definida:

La ecuación de presión en términos de la velocidad total de los fluidos es

Teorema de existencia y unicidad.

Para una clase de problemas de valor inicial, la existencia y unicidad de una solución puede ser
demostrado.
El teorema de picard-Lindelöf garantiza una solución única en e intervalo que contiene algunos
to si f y sus derivadas parciales df/dy son continuas en una región que contiene to e Yo.
El teorema construye una secuencia de funciones que convergen a la solución integral de la
ecuación, por lo tanto, la solución del problema de valor inicial. También se denomina el
método de a proximidades sucesivas.

33
El teorema de Picard-Lindelöf (muchas veces llamado simplemente teorema de Picard, otro
teorema de Cauchy-Lipschitz o teorema de existencia y unicidad) es un resultado matemático
de gran importancia dentro del estudio de las ecuaciones diferenciales ordinarias (EDO).
Establece en qué condiciones puede asegurarse la existencia y unicidad de solución de una EDO
dado un problema de Cauchy (problema de valor inicial).

Construcción de las ecuaciones de estado.

Las ecuaciones de estado son muy importantes en los estudios termodinámicos de diferentes
sistemas. Estas ecuaciones permiten estimar propiedades termodinámicas de compuestos y
sustancias, así como, evaluar los estados de sustancias puras y de mezclas, ya que son modelos
matemáticos que no requieren el uso de diagramas termodinámicos, eso sí, requieren un
cálculo numérico preciso.

Fig. 7. Ecuación de estado de los gases ideales.

Desde un punto de vista matemático, se puede decir que las ecuaciones de estado son
expresiones matemáticas que permiten modelar el comportamiento de una sustancia.
Las ecuaciones de estado tienen la siguiente forma general:

Dicha expresión es función de la presión (P), la temperatura (T) y el volumen (V).

34
Es posible expresar la forma general en las siguientes formas:

Las principales ecuaciones de estado utilizadas en la industria petrolera son:

Fluido incompresible. 𝝆 = 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕𝒆

Fluido de compresibilidad constante. 𝝆 = 𝝆𝟎 ∗ 𝒆𝒄(𝑷−𝑷𝟎 )

Fluido ligeramente compresible. 𝝆 = 𝝆𝟎 [𝟏 + 𝒄(𝑷 − 𝑷𝟎 )]
Gas ideal. 𝑷𝑽 = 𝒏𝑹𝑻
Gas real. 𝑷𝑽 = 𝒁𝒏𝑹𝑻

Ecuación de presión.

Para obtener las ecuaciones correspondientes para el flujo en de aceite y agua, donde el agua
es la fase mojante. Se tiene las siguientes ecuaciones obtenidas de adecuarla ecuación de
CONTINUIDAD para la fase w y o.

Sustituyendo y simplificando la ecuación obtenemos:

Donde:
Qt = Gasto total
Ut = Velocidad total
Esta es la ecuación de presión en términos de la velocidad total ut de los fluidos.

35
Ecuación de saturación.

Al desarrollar una ecuación para la saturación podemos enfocarla ya sea en la fase mojante o
no mojante.
Primero se obtiene la velocidad de fase mojante en términos de la velocidad total.

Combinando las ecuaciones se obtiene:

Se establecen las funciones de saturación:

Sustituyendo obtenemos:

Asumiendo incompresibilidad, de manera que ρo y ρw son constantes, así como la porosidad,

la ecuación se convierte en:

36
Simplificación de los modelos.

Considerando:
1. La porosidad es constante a lo largo de los medios de comunicación.
2. Los fluidos son incompresibles e inmiscibles.
3. No hay fuentes ni sumideros.
4. No se consideran los efectos de la presión capilar, ni de la fuerza de gravedad.
las ecuaciones que gobiernan se convierten en:

Donde:

37
Usando las correlaciones de Stone para el flujo de dos fases las movilidades de fase están
dadas por:

Aquí µ = µw/µo es la relación de movilidad, que es un número entre 0 y 1 que indica la

facilidad con la que la fase no mojante puede fluir. Con las expresiones anteriores y la
movilidad total de fases se tiene:

Estas son las ecuaciones empleadas para el flujo de fases.

Fluidos incompresibles dentro de una roca.

Para identificar los fluidos del yacimiento se toma en cuenta el coeficiente de compresibilidad
isotérmica. En general, los fluidos se clasifican en tres grupos:

a) Fluidos compresibles. 𝜌 = 𝜌0 ∗ 𝑒 𝑐(𝑃−𝑃0 )

b) Fluidos ligeramente compresibles. 𝜌 = 𝜌0 [1 + 𝑐(𝑃 − 𝑃0 )]

c) Fluidos incompresibles. 𝜌 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒

Un fluido incompresible es aquel cuyo volumen (o densidad) no cambia con la presión.

38
Para un fluido incompresible dentro del yacimiento, el flujo es igual en cualquier parte del
yacimiento, es decir, el flujo es constante.

Fig. 8. Comportamiento de P vs V de los diferentes fluidos

Formulación matemática de un fluido incompresible.

La expresión matemática de un fluido incompresible está dada por la ecuación:

ρ=constante

La cual nos indica que la densidad del fluido (peso/volumen) no se ve afectada por
incrementos de presión.
La variación de volumen del fluido respecto a la variación de presión es 0.
𝜕𝑉
=0
𝜕𝑃

La variación de densidad del fluido respecto a la variación de presión es 0.

𝜕𝜌
=0
𝜕𝑃

39
Flujo incompresible en dos fases.

El modelamiento y simulación del flujo bifásico agua-aceite en un medio poroso homogéneo.

Modelamientos de este tipo son necesarios para ayudar a comprender ciertas técnicas como la
recuperación secundaria de hidrocarburos. Los dos fluidos considerados son inmiscibles e
incompresibles. Se utiliza la ley de Darcy generalizada para el modelamiento del flujo de los
fluidos.
En el contexto de la mecánica de fluidos, un medio poroso puede ser definido como un material
en el que existen cavidades interconectadas por canales, en los cuales fluyen los fluidos;
estas cavidades son llamadas poros.
Consideremos un medio poroso homogéneo e incompresible, donde se encuentran presentes
dos fases: una fase agua denotada w y otra fase aceite.

Conservación de la cantidad de movimiento.

La ecuación de conservación de la cantidad de movimiento se traduce por la ley de

Darcy generalizada:

40
Ecuación de conservación de la masa.

Teniendo en cuenta que el medio poroso es incompresible y el flujo de los fluidos es también
incompresible, la ecuación de conservación de la masa para cada fase se escribe:

Ecuación de la presión capilar.

41
Modelo físico final.

Resumiendo, se tiene que el sistema de ecuaciones que gobierna el flujo bifásico

incompresible en un medio poroso e incompresible es:

Modelo black oil.

Este modelo maneja las tres fases petróleo, gas y el agua como componentes de masas donde
sólo permiten al gas para disolverse en el petróleo y el agua. Esta solubilidad de gas es descrita
en el petróleo y el agua por los factores de solubilidad de gas (o proporciones de gas petróleo
en solución), Rso y Rsw, respectivamente. Estas cantidades son dependientes de la presión y
esto es incorporado en el modelo de petróleo negro.
Los procesos de yacimiento que pueden ser modelados usando el modelo black oíl incluyen:

 Recuperación por expansión del fluido, mecanismo de gas en solución (agotamiento

primario).

42
  Inyección de agua incluyendo fuerzas viscosas, capilares y de gravedad (recuperación
secundaria).
 Inyección Inmiscible de gas.
 Tres procesos de recuperación de fase como agua inmiscible alternando con gas.
 Procesos de imbibición capilares.
El modelo composicional puede ser considerado como una generalización del modelo black oíl.
Se requieren un modelo de simulación de yacimiento composicional cuando efectos de
transferencia de masa en la interfaz son significativos que ocurren en el proceso de
desplazamiento del fluido. Este modelo por lo general define tres fases (gas, petróleo y agua)
pero las composiciones reales de las fases del petróleo y gas explícitamente son reconocidas
debido a su comportamiento PVT más complicado. Es decir, los componentes separados (C1,
C2, C3, etc.) en las fases del petróleo y de gas explícitamente son rastreados.
Toda conservación de masas es aplicada a cada componente (el petróleo, el gas y agua) como
en el modelo black oíl. Los ejemplos de los procesos de yacimiento que pueden ser modelados
usando un modelo composicional incluyen:

 Inyección de gas con el petróleo movible por primer contacto o desarrollado (multi-
contacto) miscibilidad (por ejemplo, inyección de CO2).
 La simulación de inyección de gas cerca en yacimientos que se encuentran cerca de
condiciones críticas.
 Procesos de recirculación de gas en yacimientos condensados.
Este modelo se basa en la composición del fluido proveniente del yacimiento tomando en
cuenta todas las fracciones de petróleo tanto livianas como pesadas, y la aplicación de
correlaciones empíricas para similar las fases presentes a determinadas condiciones de presión
y temperatura basadas en las ecuaciones termodinámicas de equilibrio para las fases líquidas y
la fase de gas.

Este modelo Composicional es un método más exacto, pero a su vez es más complicado, este
modelo es justificado para cálculo de problemas de fluidos volátiles necesitando un cálculo de
transferencia de calor riguroso.

43
 Fig. 9. Fases del modelo composicional.

El Modelo De Dos Fases Compresible

El método de elementos finitos para flujos incompresibles

Descripción del método MEF.

El método de los
elementos finitos es
un método de
aproximación de Nodo
problemas
s
continuos, de tal
forma que:
Element
os
El continuo se divide
en un número finito
de partes,
“elementos”, cuyo
comportamiento se

44
especifica mediante un número finito de parámetros asociados a ciertos puntos característicos
denominados “nodos”. Estos nodos son los puntos de unión de cada elemento con sus
adyacentes.

• Las incógnitas del problema dejan de ser funciones matemáticas y pasan a ser el valor
de estas funciones en los nodos.
• El comportamiento en el interior de cada elemento queda definido a partir del
comportamiento de los nodos mediante las adecuadas funciones de interpolación o
funciones de forma.

El MEF está pensado para ser usado en computadoras y permite resolver ecuaciones
diferenciales asociadas a un problema físico o ingenieril sobre geometrías complicadas, así
como en la simulación de sistemas físicos complejos.

MEF para flujos incompresibles

El desarrollo de técnicas eficientes, son varios aspectos a considerar cuando planteamos un

modelo de simulación numérica de flujos de fluidos incompresibles mediante el Método de los
Elementos Finitos. El estudio teórico de las ecuaciones relevantes que gobiernan los flujos de
fluidos incompresibles, esencialmente las de Stokes y Navier-Stokes, fue realizado y completado
en los años 60 y primeros 70.
La resolución numérica de flujos de fluidos incompresibles presenta dos restricciones de
estabilidad de relevancia, bien conocidas por la comunidad de investigadores y usuarios en
Mecánica de Fluidos computacional. Por una parte, la condición de incompresibilidad impone
condiciones de compatibilidad bastante estrictas entre los espacios de discretización de
velocidad y presión empleados.
La simulación de tales flujos mediante el uso directo del método de discretización Elementos
Finitos produce en general soluciones numéricas y los espacios de discretización de velocidad y
presión son compatibles.

45
Resumen del método MEF.

El MEF (método de elementos finitos), por tanto, se basa en transformar un cuerpo de

naturaleza continua en un modelo discreto aproximado, esta transformación se denomina
discretización del modelo.
El conocimiento de lo que sucede en el interior de este modelo del cuerpo aproximado, se
obtiene mediante la interpolación de los valores conocidos en los nodos.
Es por tanto una aproximación de los valores de una función a partir del conocimiento de un
número determinado y finito de puntos.

Se trata de introducir espacios de discretización diferentes para velocidad y presión, de forma

que satisfagan una condición de compatibilidad específica, conocida como de Brezzi-Babuska,
con objeto de eliminar la presión (ef. Griffits (1979), Hecht (1981)).
“Ver soluciones de las ecuaciones de Stokes, así como sobre su aproximación mediante el MEF (método de
elementos finitos)”.

46
La clave para conseguir estabilizar Elementos Finitos "incompresibles" es debilitar la condición
de divergencia nula. Ello es realizado en los métodos clásicos enriqueciendo el espacio de
velocidades en grados de libertad frente al espacio de presiones. Otra forma de lograrlo es
relajar explícitamente la condición:

𝛻∙𝑢 =0 Cambiando a

Siendo gh una función "pequeña" conveniente

𝛻 ∙ 𝑢ℎ = 𝑔ℎ
Con ello se consigue cierto control sobre el gradiente de presión (y no sólo
sobre la presión), impidiéndose la aparición de modos de presión. Apoyarse del MEF para
encontrare soluciones a un modelo de flujos incompresibles.

Planteamiento matemático de MEF.

Un elemento finito e viene definido por sus nodos (i,j,m) y por su contorno formado por líneas
que los unen. Los desplazamientos u de cualquier punto del elemento se aproximan por un
vector columna u.

Figura 10. Coordenadas nodales (i, j, k) y desplazamientos de los nodos

47
N son funciones de posición dadas (funciones de forma) y ae es un vector formado por los
desplazamientos nodales de los elementos considerados.

Las funciones Ni, Nj, Nm, han de escogerse de tal forma que al sustituir en (5.1) las coordenadas
nodales, se obtengan los desplazamientos nodales.
Conocidos los desplazamientos de todos los puntos del elemento, se pueden determinar las
deformaciones (ε) en cualquier punto.
El MEF se puede entender, desde un punto de vista estructural, como una generalización del
cálculo matricial de estructuras al análisis de sistemas continuos. De hecho, el método nació
por evolución de aplicaciones a sistemas estructurales a la a actualidad a aplicaciones de ciencia
e ingeniería.

48
El método de elementos finitos para flujos en dos fases.

Flujo en dos fases.

Por definición, el flujo multifásico es el flujo interactivo de dos o más fases distintas con
interfaces comunes en un conducto. Cada fase, que representa una fracción de volumen (o
fracción de masa) que tiene sus propias propiedades, velocidad y temperatura.

Figura 11. Flujo en tuberías.

La influencia de la fuerza del cuerpo gravitacional en los flujos multifásicos es de mucha mayor
importancia que en el caso de los flujos monofásicos.

 La velocidad del sonido cambia drásticamente para los materiales que experimentan un
cambio de fase, y puede ser de diferente magnitud. Esto influye significativamente en
un flujo del fluido.
 El cambio de fase significa que las caídas de presión inducidas por el flujo pueden causar
un cambio de fase adicional.
 Hay muchos tipos de inestabilidades en el flujo multifásico.

Uno de los aspectos más desafiantes de lidiar con el flujo de dos fases es el hecho de que puede
tomar muchas formas diferentes. Las distribuciones espaciales y las velocidades de las fases
líquida y de vapor en el canal de flujo es un aspecto muy importante.
Las estructuras de flujo observadas se definen como patrones de flujo de dos fases y tienen
características de identificación particulares. Estos diferentes patrones de flujo han sido
categorizados de acuerdo con la dirección del flujo en relación con la aceleración gravitacional.

49
Flujo en dos fases, patrones de flujo.

Los principales regímenes de flujo en tubos verticales se muestran en la tabla. Debe tenerse en
cuenta que los valores de calidad y caudal dependen del fluido y la presión.

Figura 12. patrones de flujo.

El método de elementos finitos (FEM) ha tenido importantes avances para poder resolver con
éxito los desafíos de la mecánica de fluidos, pero todavía para problemas complejos

50
relacionados a flujos multifásicos o fuertemente turbulentos, el método de volúmenes finitos
es en la práctica la única solución.

Los métodos los cuales son la base de la modelación matemática para la mecánica del continuo,
por definición, son respetadas por las ecuaciones discretas deducidas por el método de
volúmenes finitos. El método no se limita sólo a problemas de mecánica de fluidos, y de forma
general envuelve los siguientes pasos:
1. Descomponer el dominio en volúmenes de control
2. Formular las ecuaciones integrales de conservación para cada volumen de control
3. Aproximar numéricamente las integrales
4. Aproximar los valores de las variables en las caras y las derivadas con la información de
las variables nodales
5. Ensamblar y resolver el sistema algebraico obtenido

Ecuación de Bernoulli.

El principio de Bernoulli es un enunciado que parece ir en contra de la intuición, acerca de cómo

la velocidad de un fluido se relaciona con la presión del fluido.
Dentro de un flujo horizontal de fluido, los puntos de mayor velocidad del fluido tendrán menor
presión que los de menor velocidad.
Así que dentro de una tubería horizontal de agua que cambia de diámetro, las regiones donde
el agua se mueve más rápido se encontrarán a menor presión que las regiones donde se mueve
más lento.

51
 Figura 12. Flujo a través de un reductor.

Básicamente es la suposición de que toda la energía cinética del fluido está contribuyendo
directamente al proceso de avance del flujo del fluido -. Ello debería hacer evidente que la
existencia de turbulencias o cualquier movimiento caótico del fluido implicaría que algo de la
energía cinética no está contribuyendo al avance del fluido a través del tubo.

Bernoulli, velocidad y presión.

Los fluidos incompresibles tienen que aumentar su velocidad cuando alcanzan una sección más
estrecha para mantener el volumen de flujo constante. Por esta razón, una boquilla estrecha
en una manguera causa que el agua salga más rápido. Puede ser que algo te esté molestando
sobre este fenómeno: si el agua se acelera en la constricción, también gana energía cinética.

Esta relación inversa entre la presión y la velocidad en un punto en un fluido se llama el principio
de Bernoulli. en puntos a lo largo de una línea horizontal de flujo, las regiones de mayor presión

52
tienen una menor velocidad del fluido, y las regiones de menor presión tienen una mayor
velocidad del fluido.

Construcción de un modelo simple de simulación.

¿Qué es simular?

Representar un fenómeno físico, mediante un conjunto de programas de computación que

resuelve los modelos reales, mediante algoritmos matemáticos y obtiene soluciones
aproximadas de las condiciones de entorno real. Se ingresa las ecuaciones y relaciones
fundamentales que modela el hecho físico. En el caso de reservorios, es fundamental incluir
ecuaciones de masa, energía y continuidad.

¿Cuándo simular?

No existe un tiempo generalizado para empezar una simulación, pero es recomendable hacerlo
luego de un lapso de tiempo prudencial de producción, que nos permita obtener una data
confiable del yacimiento para de esta forma garantizar los resultados.

Figura 13, Cuándo simular

53
Tabla 1. Simuladores más usados.

54
Planeación de un estudio de simulación de yacimientos

Las actividades más importantes que se llevan a cabo durante un estudio de simulación de
yacimientos son:

Diagrama 7. Actividades de simulación.

55
Revisión de la Información Disponible de Campo.

La información recopilada debe ser revisada para garantizar su confiabilidad y detectar

inconsistencias, además de organizarse de manera que permita determinar las características
del modelo que se va a construir.

Para la construcción del modelo de simulación, se contó con la siguiente información:

• Modelo ajustado de un PVT composicional.
• Mapa de cimas y bases.
• Mapa de pozos.
• Definición de límites (físicos y convencionales).
• Valores promedio de porosidad, permeabilidad y relación de espesores (NTG).
• Modelo petrofísico.
• Historia de presión-producción.

Modelo PVT de los Fluidos.

Para modelar el comportamiento termodinámico de los fluidos del yacimiento, se utilizó un

modelo PVT composicional, el cual se obtuvo de ajustar los principales parámetros volumétricos
reportados en un análisis PVT experimental realizado a una muestra tomada en uno de los
pozos, la cual se consideró representativa de los fluidos ya que cuando se tomó dicha prueba,
la presión del yacimiento se encontraba por arriba de la presión de saturación.

Revisión de la Información Disponible de Campo.

La información recopilada debe ser revisada para garantizar su confiabilidad y detectar

inconsistencias, además de organizarse de manera que permita determinar las características
del modelo que se va a construir.
Para la construcción del modelo de simulación, se contó con la siguiente información:
 Modelo ajustado de un PVT composicional.
 Mapa de cimas y bases.

56
  Mapa de pozos.
 Definición de límites (físicos y convencionales).
 Valores promedio de porosidad, permeabilidad y relación de espesores (NTG).
 Modelo petrofísico.
 Historia de presión-producción.
 Actualmente, se cuenta con diferentes programas capaces de aplicar técnicas de
mallado y escalamiento para modelos de simulación de yacimientos; entre los de
mayor uso se puede mencionar algunos como: GRID©, FloGrid© y los módulos de
PETREL©

Modelo PVT de los Fluidos

 Para modelar el comportamiento termodinámico de los fluidos del yacimiento, se
utilizó un modelo pVT composicional, el cual se obtuvo de ajustar los principales
parámetros volumétricos reportados en un análisis pVT experimental realizado a una
muestra tomada en uno de los pozos, la cual se consideró representativa de los fluidos
ya que cuando se tomó dicha prueba, la presión del yacimiento se encontraba por
arriba de la presión de saturación.

Figura 14. Diagrama de fases de una mezcla de 6 pseudocomponentes.

57
58
Modelo Petrofísico.

El modelo petrofísico empleado, surgió del análisis de la distribución de saturación del agua con respecto a la
profundidad, y con la definición de presión capilar, se obtuvo la distribución de ésta con respecto a la profundidad.

59
Construcción de la malla de simulación.

Normalmente, para construir una malla de simulación, se toma como base un modelo estático que posee las
heterogeneidades geológicas, la descripción estructural, estratigráfica y distribución de propiedades petrofísicas
del yacimiento. Para la construcción de la malla de simulación, se delimitó el área de interés de la configuración
estructural del horizonte de la Brecha Paleoceno del Cretácico Superior (BTPKS) definida por los límites físicos y
convencionales asociados con el contacto agua-aceite, fallas presentes, cambio de facies y criterios inferidos.

60
Presenta el mapa de cimas de la formación BTPKS, de fallas, y el mapa de pozos, el cual contiene únicamente a
dos, H-1 y H-2

61
62
Fig. 3.10 Construcción de la malla

63
Integración del modelo numérico (data)

El archivo DATA, también llamado “archivo de datos”, es un archivo de texto ASCII que contiene
toda la información para el modelado del yacimiento. Está subdividido en secciones
encabezadas por “palabras clave” que identifican los datos de entrada y datos de salida de cada
sección. A continuación, se describe las secciones del archivo DATA y la información que
contiene cada una de ellas:

64
 • RUNSPEC

• GRID
Contiene los archivos de salida provenientes del preprocesador GRID, como son los sistemas
coordenados (COORD y ZCOORD), la distribución de porosidad, celdas activas y celdas inactivas,
relación de espesores netos y brutos (NTG), permeabilidades, factores que involucran el
tamaño del bloque de matriz y el control de flujo entre la matriz y la fractura, lo que se
denomina como SIGMA.

65
66
 • EDIT
En la sección EDIT, se pueden modificar algunos parámetros que permitan ajustar el modelo,
como el volumen poroso y las transmisibilidades en (i, j, k). En este modelo, no se utilizó la
sección EDIT debido a que el volumen poroso se modificó en la sección GRID con la instrucción
EQUALS, que asigna o reemplaza el valor de una propiedad a un grupo de celdas dentro de la
malla. Fig. 3.14-2

• PROPS
Contiene las tablas de propiedades con los datos de presión y saturación de la roca y fluidos del
yacimiento en función de la densidad, viscosidad, permeabilidad relativa, presión capilar, datos
pVT, compresibilidades

• REGIONS
Define las diferentes regiones correspondientes a las funciones de saturación y volumen en
base a los contactos agua-aceite (regiones de equilibrio)

67
 Manual simulador Eclipse.

Objetivo.

Manual de uso simulador eclipse

La Simulación diferencial finita de yacimientos es una de las herramientas más usadas para
apoyar las decisiones de la Gerencia de yacimientos, desde la determinación de las reservas, a
través de una planificación inicial de la etapa de producción y el diseño de las instalaciones de
superficie, así como el diagnóstico de problemas.

El simulador de yacimientos ECLIPSE es un software que permite a los ingenieros predecir y

manejar eficientemente el flujo de fluidos creando modelos de una manera más certera.

68
 Simulación de yacimientos es posible:
Determinar el comportamiento de un yacimiento sometido Determinar los efectos de localización y espaciamiento de
a inyección de agua/o gas, bajo agotamiento natural los pozos

Estimar el efecto de la tasa de producción sobre la Estimar líneas de drenaje en yacimientos heterogéneos de
recuperación petróleo y gas

Calcular la producción total de un yacimiento para un

numero dado de pozos en localizaciones especificas, entre
otros.

ECLIPSE permite plantearse cuestionamientos sobre la viabilidad económica de la explotación

de un yacimiento, como por ejemplo:

- Cuan cierto es la estimación de reservas? Puedo minimizar la incertidumbre?

- Como es la distribución de los hidrocarburos en la formación?
- Que tan difícil puede ser la recuperación de los hidrocarburos y que técnicas de producción y
tasas son las mejores?
- Donde debo localizar las próximos pozos para maximizar el recobro y recuperar la inversión?,
etc.

Modelos de Simulación.

ECLIPSE ofrece tres modelos de simulación, blackoil, composicional y termal.

El modelo Blackoil se utiliza para situaciones donde necesito conocer el efecto de una de las
fases en las reservas y el recobro. Este modelo asume que los fluidos del yacimiento consisten

69
de tres “componentes” (agua, petróleo y gas) donde las composiciones de petróleo y gas son
constantes en el tiempo y a cada paso de presión.

La solución blackoil permite modelar extensivamente el control de pozos, soportar

eficientemente la planificación de las operaciones del yacimiento, incluyendo la inyección de
gas, debido a que proporciona una simulación de las tres fases en 3D.

Al simular yacimientos altamente fracturados se requiere un apoyo para la inhibición

gravitacional y el drenaje, difusión molecular y el desplazamiento de los fluidos viscosos a través
de la matriz de la roca. La solución blackoil ofrece amplias opciones duales para la porosidad y
permeabilidad.

Definición de complejos mecanismos de los fluidos del yacimiento.

Cuando el petróleo consiste de dos o más hidrocarburos y estos exhiben diferentes fases y
cambios en la composición, relativos a la temperatura y presión, y el simulador requiere
complejas definiciones del comportamiento de los fluidos para la predicción del recobro, se
utiliza un simulador composicional.

En ECLIPSE Composicional los fluidos pueden ser representados por ecuaciones de estado.
Este modelo es una buena opción para estudiar los crudos volátiles o condensados, los
programas de inyección de gas y el estudio del recobro secundario. El Conocimiento del
comportamiento de la composición es importante para planificar y diseñar las instalaciones de
superficie.

La simulación composicional es ideal para:

- Incrementar la inyección de gas o mantener la presión del yacimiento.

- Inundaciones miscibles debido a la inyección de gas
- Flujos de dióxido de carbono provenientes del gas en solución.
- Yacimientos con fluidos composicionales cercanos al punto de burbujeo.
- Altas presiones y temperaturas del yacimiento.
- Modelaje de yacimientos fracturados.

70
Tratamiento con petróleo pesado.

Debido a que los fluidos con altas densidades y viscosidades presentan cambios masivos "en
sitio" con respecto a la temperatura del yacimiento y la energía calorífica, avanzados métodos
de recobro termal son típicamente requeridos para petróleos pesados, extrapesados y
yacimientos bituminosos.

La simulación termal en los sistemas de producción usa una implícita jerarquía de multinivel
para generar modelos de pozos de flujo cruzado que incrementa el realismo. Esto permite
analizar el petróleo, gas, agua y el comportamiento de los sólidos en los pozos con modelos de
petróleo-vivo o petróleo-muerto.

Estas soluciones incorporan las buenas-prácticas en las tecnologías para que clientes,
investigadores y universidades del mundo simulen una amplia variedad de procesos de
recobro termal.

Estructura de la malla de simulación.

Los diferentes parámetros para la construcción de la malla incluyen elevaciones o topes de estructuras,
permeabilidad en tres direcciones, porosidad, fallas, espesores, fracturas y acuíferos. La resolución del
modelo depende de la resolución de la malla. Una malla fina divide el yacimiento en pequeños bloques
y numéricamente más exactos. La malla gruesa tiene menor cantidad de bloques, pero la malla gruesa
es más larga que los bloques de la malla fina en volumen. Esto quiere decir, que la malla gruesa tiene
menos problemas en realizar corridas que la malla fina, pero menos exactos.

Diseño de la malla de simulación.

Se puede diseñar la malla del yacimiento de diferentes formas. Puede ser construida en dos o
tres dimensiones y en coordenadas cartesianas o cilíndricas. Existen técnicas disponibles para
aproximar la distribución de los fluidos en modelos en 2D y 3D como modificando curvas de
permeabilidades relativas y presiones capilares.

71
ECLIPSE es un software, una herramienta para la simulación de yacimientos.
Presenta varias opciones: modelamiento de petróleo negro (ECLIPSE 100) y
modelamiento composicional (ECLIPSE 300).

Eclipse 100

para propósitos generales con la opción para simular condensado de gas, define a las fases
petróleo y gas como unidades cuya composición es constante a cualquier presión y tiempo.
ECLIPSE es construido en lenguaje de programación FORTRAN77 y puede ser ejecutado en
cualquier computador con un compilador ANSI-Standard y suficiente memoria.

72
¿Qué es FORTRAN 77?

El FORTRAN es uno de los lenguajes de programación más ampliamente utilizados,

especialmente entre la comunidad cientıfica. Fué el primer lenguaje de alto nivel que fué
aceptado. Su nombre deriva de FORmula TRANslation (traducción de fórmulas) y ha tenido
varias revisiones.

El American National Standards Institute (Instituto Nacional Estadounidense de Estándares),

cuyo acrónimo es ANSI, tiene por objetivo crear directrices que sean ampliamente aceptadas
en una amplia gama de sectores. Su objetivo es crear uniformidad para las mercancías y los
procesos de la industria.

Nuevas Facilidades:

ECLIPSE 100
• RUNSPEC section
• GRID section
• PROPS section
• REGIONS section
• SUMMARY section

Es versión de Eclipse es la opción para simular los fluidos del reservorio con sus variables;
composición fraccional de componentes y la simulación de estos mediante ecuaciones cúbicas
de estado, valores de permeabilidad dependientes de la presión y tratamiento de petróleo
negro.

La elección del tipo de simulador está determinada por las características del reservorio
respecto de su punto crítico o presión de burbuja, cuando se debe modelar un sistema que se
encuentre cerca del punto crítico donde los más pequeños cambios de presión y temperatura
afectan considerablemente al comportamiento composicional de los diferentes fluidos es
necesario el uso de un simulador de modelamiento composicional.

El programa accede a la información desde un archivo plano único de datos que contiene la
descripción completa del modelo. Este archivo es una recopilación de comandos conocidos con
el nombre de “keywords” que indican que tipo de información ingresada y/o generada se
empleará, el programa la procesa y genera resultados que son leídos y visualizados con otros
keywords destinados para el efecto. El archivo único “data file” es diseñado para ser leído en el
programa por secciones, cada sección está diseñada para un propósito específico.

Las secciones del archivo son definidas y determinadas para ejecutar y almacenar información
de acuerdo a los pasos y descripción del modelo, las secciones que representa los pasos que
sigue el simulador para la generación de un modelo dinámico son las siguientes.

73
Un archivo de entrada de datos ECLIPSE se divide en secciones, cada una de las cuales se
introduce mediante una palabra clave. una lista de todos Las palabras clave del encabezado de
sección se proporcionan a continuación, junto con una breve descripción del contenido de cada
sección. A Se puede encontrar un desglose más detallado de los contenidos de la sección en los
resúmenes de sección que siguen. inmediatamente después de esta descripción general.

74
Secciones del archivo principal o “data file”.

RUNSPEC (especificaciones de la corrida)

Esta sección está diseñada para determinar o especificar varios componentes como la cantidad
requerida de memoria RAM, también la distribución de la memoria entre los componentes de
la simulación que son: pozos, datos tubulares, la malla de simulación, vectores solución, en
otras palabras las dimensiones del problema. También es necesaria para la especificación de las
características básicas del modelo, como es el tipo de simulador que se empleará, el título de
la simulación, la fecha de inicio, el número de fases empleadas, el tamaño de la malla de
simulación, el tipo de unidades que se desea emplear en la simulación, entre otras.
ECLIPSE genera la simulación ejecutando las soluciones para cada unidad de la malla y
consecutivamente a la siguiente unidad de modo ortogonal, por lo que requiere de rápido
acceso a muchas soluciones previas, este proceso es conocido como apilamiento de las
soluciones o “Solver Stack”.
La palabra clave CPR se ha agregado para activar el solucionador lineal CPR, que emplea el
Restringido Método Residual de Presión para resolver las ecuaciones lineales. A menudo, esto
puede ser útil para tareas más difíciles. problemas, cuando el solucionador lineal ORTHOMIN
estándar toma una gran cantidad de iteraciones lineales para converger.
La palabra clave MENSAJES ahora tiene un elemento número 13 que controla el número
máximo de veces que algunos los mensajes se envían al archivo PRT.

GRID (malla de simulación)

La sección GRID está diseñada para proveer al simulador de la información necesaria para el
cálculo del volumen poroso y de la transmisibilidad en todas las direcciones con el propósito de
determinar los cálculos del flujo de fluidos.

75
En esta sección se determina la geometría básica de la malla de simulación y propiedades de la
roca, porosidad, permeabilidad absoluta, espesores netos, en cada celda. Esta información es
importada de pre-procesadores como Gocad.

Upscaling (Escalamiento)
En la definición del modelo se determina las medidas de la malla de simulación y de cada bloque
o celda. Una celda representa una unidad volumétrica del reservorio, esta debe tener relación
con los objetivos que se desean lograr por lo que se debe tomar en cuenta la relación de
escalamiento, medidas y forma de las celdas.
Dentro de cada celda se define a sus propiedades como las mismas en cada punto dentro de los
límites (I1,I2) x (J1,J2) x (K1,K2). Para efectos de cálculos y visualización las propiedades de la
celda están determinadas para su centro:
((I1+I2)/2; (J1+J2)/2; (K1+K2)/2).
La palabra clave MINNNCT se ha ampliado para permitir al usuario especificar una difusividad
mínima y transmisibilidad térmica para una conexión no vecina, así como (o en lugar de) un
mínimo transmisibilidad.
La palabra clave ROCKFRAC, anteriormente admitida solo en ECLIPSE 300, se ha agregado a
ECLIPSE 100 para simplificar la entrada de fracciones de volumen de celda que representan roca
en metano de lecho de carbón o gas de esquisto carreras modelo.
El elemento 215 de la palabra clave OPTIONS ahora se puede usar para reemplazar los valores
enumerados no establecidos con el valor de una celda adyacente. Esto puede ser útil en
modelos gruesos donde una celda representativa es inactivo en el modelo original.
El tipo de geometría es muy importante, ECLIPSE presenta la posibilidad de escoger de entre
los siguientes tipos:
Corner Point: utilizado en estructuras geológicas difíciles.
Block Centered: celdas rectangulares con horizontes paralelos y perpendiculares
Radial: Utilizado para estudios en las cercanías de los pozos.
PEBI (unstructured): se utiliza una malla no-estructurada en secciones falladas y donde se
requiere de mayor detalle.
No siendo posible mezclar varios tipos en el mismo estudio, pero con la ayuda de la acción
refinamiento es posible redefinir ciertos sectores en los cuales se puede subdividir las celdas y
cambiar el tipo de geometría.

76
EDIT (edición)
Aunque esta sección es opcional es utilizada cuando se requiere modificar el volumen poroso,
profundidad del centro de cada bloque, tansmisibilidades, difusividad y otros datos generados
en la sección GRID.

PROPS (propiedades)
Esta sección esta diseñada para la importación, edición y manejo de datos de los fluidos,
requeridos para determinar la distribución de los mismos en el reservorio en la etapa inicial y
durante el tiempo de flujo, para realizar cálculos de balance de materiales en cada celda en
cada intervalo de tiempo; estos datos son: PVT, funciones de saturación (permeabilidad
relativa y presión capilar) y datos de compresibilidad de la roca.
• Los datos PVT requeridos son: Factor Volumétrico, Viscosidad, GOR.
• La información de fluidos son: PVT en función de la presión, Densidad o gravedad.
• Información requerida de la roca es: permeabilidad relativa en función de la
saturación, presión capilar en función de la saturación, compresibilidad de la roca en
función de la presión.

REGIONS (regiones)
Esta sección es opcional pero está diseñada para dividir al reservorio en regiones de equilibrio
de similares características PVT. Esta acción es realizada para la mejor caracterización del
reservorio de acuerdo a variaciones zonales en las propiedades y/o para propósitos de reportes.
Generalmente se asigna una región de equilibrio para cada contacto agua petróleo existente.
Se añaden las palabras clave SURFSTES y ESSNODE para activar el tensioactivo dependiente de
la salinidad. modelo que permite que la tensión superficial agua-petróleo varíe con la salinidad
así como con el surfactante concentración. Estas tablas se pueden enviar al archivo PRINT con
el argumento RPTPROPS SURF.
La palabra clave LANGMPL se puede utilizar para escalar la presión utilizada para calcular las
capacidades de adsorción para un modelo de metano de lecho de carbón.
Las palabras clave de escalado de punto final LSWL, LSWLPC, LSWCR, LSWU, LSOWCR, LSOGCR,
LKRW, LKRWR, Se agregaron LKRO, LKRORW, LKRORG y LPCW para permitir escalar la baja
salinidad (mojado por aceite) puntos finales de la tabla de saturación de agua y petróleo cuando
las opciones LOWSALT y ENDSCALE han sido seleccionadas. especificado.

77
Las palabras clave de escalado de punto final HWSWL, HWSWLPC, HWSWCR, HWSWU,
HWSOWCR, HWSOGCR, Se agregaron HWKRW, HWKRWR, HWKRO, HWKRORW, HWKRORG y
HWPCW para permitir escalar los puntos finales de la tabla de saturación de aceite y agua
mojada por agua del surfactante (alta salinidad) cuando el SURFACTW y se han especificado las
opciones ENDSCALE.

SCHEDULE
Las palabras clave PLYVISC y PLYVISCS pueden usarse para actualizar los multiplicadores de
viscosidad del polímero. previamente especificado en la sección PROPS.
La palabra clave CSKIN se ha agregado para permitir que los factores de máscara de conexión
se actualicen utilizando un usuario- argumento definido.
Los argumentos POLY y POLYVM se agregaron a las palabras clave RPTSCHED y RPTRST para
generar concentración de polímero y multiplicador de viscosidad de polímero en estudios de
inundación de polímero, respectivamente.
Los argumentos ESALTP y ESALTS se agregaron a las palabras clave RPTSCHED y RPTRST para
valores de salinidad efectivos de salida calculados usando entradas a la palabra clave SALTEFF
en polímero y estudios de inundación de surfactante, respectivamente.

Palabras claves
Las palabras clave en el archivo de datos de entrada (incluidas las palabras clave de encabezado
de sección) tienen hasta 8 caracteres cada una. longitud y debe comenzar en la columna 1.
Todos los caracteres hasta la columna 8 son significativos. Cualquier carácter en el misma línea
que una palabra clave desde la columna 9 en adelante se tratará como un comentario.

Data de palabras claves

Los datos de una palabra clave deben seguir a esa palabra clave en una nueva línea. Los datos
se pueden dividir en espacios en blanco o nuevas líneas, y la posición en la línea no es
significativa.

Repetir recuentros
En los datos que siguen a una palabra clave, se pueden usar asteriscos para indicar recuentos
repetidos. Una cantidad de datos puede ser repetido un número requerido de veces

78
precediéndolo con el número requerido y un asterisco. Por ejemplo, los dos casos siguientes
son equivalentes:

Valores predeterminados
La descripción de la palabra clave indica cuándo se pueden aplicar los valores predeterminados.
Hay dos formas de establecer las cantidades en sus valores predeterminados. En primer lugar,
al terminar un registro de datos prematuramente con una barra inclinada (/) las cantidades
restantes sin especificar se establecen en sus valores predeterminados.
En segundo lugar, las cantidades seleccionadas colocadas antes de la barra inclinada se pueden
establecer por defecto ingresando n* donde n es la número de cantidades consecutivas a ser
predeterminadas.

Comentarios
Todas las líneas que comienzan con los dos caracteres ' -- ' se tratan como comentarios y
ECLIPSE las ignora. Las líneas de comentarios (y también las líneas en blanco) se pueden insertar
en cualquier parte del archivo de datos. Los comentarios también pueden ser agregado al final
de las líneas de datos comenzando el comentario con los dos caracteres '--', pero en este caso
los comentarios no deben contener comillas.

Palabras claves Global

Algunas palabras clave pueden aparecer en cualquier sección del archivo de entrada. Estas
palabras clave y breves descripciones de sus funciones se enumeran a continuación.

Leer y repetir el archivo de entrada

Las palabras clave ECHO y NOECHO activan y desactivan el eco del archivo de entrada en el
archivo de impresión. La inicial por defecto está haciendo eco. La palabra clave INCLUDE
permite leer datos de otro archivo con nombre. Le sigue el nombre de un archivo del que se

79
tomará la entrada. Una vez leído, el archivo se cierra y se reanuda la entrada desde el archivo
principal, a partir de la palabra clave después de INCLUDE.

Ejemplo:

La palabra clave SKIP hace que se omitan (ignoren) todas las palabras clave subsiguientes hasta
que aparezca una palabra clave ENDSKIP. encontrado. La palabra clave SKIP100 omite las
palabras clave intermedias cuando se ejecuta ECLIPSE 100. De manera similar, la palabra clave
SKIP300 omite las palabras clave intermedias cuando se ejecuta ECLIPSE 300.
La palabra clave COLUMNS se puede utilizar para restablecer los márgenes izquierdo y derecho
para leer el archivo de datos de entrada. Solamente los datos dentro de estos márgenes son
leídos y procesados. La configuración predeterminada para los márgenes izquierdo y derecho
es 1 y 132 respectivamente. La palabra clave FORMFEED establece el carácter de avance de
formulario en el archivo de impresión. El valor predeterminado es usar estándar. Control de
carro FORTRAN.

Configuración de límites de impresión y detención de mensajes

ECLIPSE tiene un sistema de manejo de mensajes que tiene seis niveles de severidad:
1.Mensajes
2. Comentarios
3. Advertencias
4. Problemas
5. Errores
6. Errores

Control de depuración
Las palabras clave DEBUG y DEBUG3 controlan la salida al archivo de depuración en ECLIPSE
100 y ECLIPSE 300 respectivamente. Cada uno de los enteros que siguen a la palabra clave

80
controla la salida de una opción de programa o activa una prueba. Esto está destinado
principalmente a fines de desarrollo de programas.

Finalización de archivos de entrada

La palabra clave END finaliza la lectura de datos antes del final real de un archivo de entrada. el
programa hace no repetir ni procesar datos después de esto. END puede usarse en un archivo
INCLUDE. Para finalizar la lectura de un INCLUYE archivo antes del final real del archivo, y
devolver el control al archivo de entrada principal, se puede utilizar ENDINC. Fin y ENDINC se
generan automáticamente al final real de los archivos relevantes, por lo que normalmente no
es necesario usó.

Lista de palabras claves global

SOLUTION (solución)
Esta sección esta diseñada para ejecutar la simulación y definir una solución inicial. Se
determina las condiciones iniciales del reservorio como son: la distribución de presiones,
caudales iniciales gas-petróleo o petróleo-gas para cada celda, dependencia de las propiedades

81
de los fluidos del reservorio con respecto a la profundidad, condiciones iniciales de acuífero. La
inicialización puede ser realizada de tres maneras distintas:
Equilibración: el software calcula la presión y saturación inicial utilizando los datos ingresados.
Como resultado se tiene un “RESTART FILE” con lascondiciones iniciales determinadas para el
campo y desde el cual parte el ajuste del modelo y las predicciones para el campo.
Reinicio: la solución inicial es obtenida de un archivo de reinicio creado en una corrida anterior.

Enumeración: el usuario es quién define de acuerdo a un amplio conocimiento del reservorio

las condiciones iniciales, especificando explícitamente la solución inicial para cada celda de la
malla.

Antes de la producción de un campo, los fluidos se encuentran en equilibrio hidrostático, para

la recreación de esta etapa ECLIPSE requiere de datos de profundidad, contactos agua–petróleo
y/o gas–petróleo y las presiones iniciales a una profundidad de referencia. Las funciones de
saturación toman los datos ingresados en la sección PROPS para definir la distribución de los
fluidos en el reservorio.
El modelamiento de acuíferos es determinante al generar un modelo, debe ser modelado
numérica y analíticamente, en la sección GRID se realiza la especificación numérica y en
SOLUTION la determinación analítica. El número de acuíferos y el máximo número de celdas
que están conectados son especificados en la sección RUNSPEC.

SUMMARY (vectores resultado)

Esta sección es opcional pero es de especial interés cuando se requiere la especificación de
variables que tienen que ser escritas en el o los archivos de salida luego de cada simulación. Si
no hay una sección SUMMARY, ECLIPSE no crea un archivo de resumen. Esta sección es
básicamente usada para la solicitud de información definida sea por bloque, campo o pozo.
SCHEDULE
Esta es la última sección y es usada para especificar los “timesteps” o periodos de tiempo en
que se dividirá la simulación. También es la sección en donde se ingresa los datos reales de
presión y producción para ser comparados con los calculados por el programa.
Esta sección es diseñada para ser utilizada en dos modos:
• Ajuste Histórico
• Predicciones

82
El ajuste histórico se realiza al especificar las mediciones de: presión, producción e inyección
históricas registradas en los pozos y las mediciones obtenidas como resultado del proceso de
simulación. ECLIPSE simula el comportamiento de producción, inyección y presión tomando los
datos ingresados al programa (modelo estático e inicialización) y herramientas matemáticas.
El comportamiento del reservorio es reproducido por el simulador y comparado con los datos
reales históricos del campo para determinar su ajuste. El proceso de ajuste continúa con una
serie de cambios bien sea globales (campo o reservorio) o locales (pozos) para determinar las
condiciones de flujo y reproducir el historial de caudales de producción y presiones. En el
proceso de ajuste se debe realizar todas las corridas necesarias con cada uno de los cambios
planeados, dependiendo de la complejidad del modelo y número de pozos éste proceso puede
tomar semanas o incluso meses de dedicación. Dependiendo de este ajuste se determinará la
exactitud de las predicciones.
El ajuste se realiza al variar en un rango tolerable la información ingresada del campo, ésta no
refleja necesariamente la realidad del campo por ser datos generalmente puntuales obtenidos
de pozos. Las variables del reservorio que generalmente son aceptadas como verdaderas por
ser consideradas de baja incertidumbre, se muestran a continuación.
• Porosidad.
• Espesor total de la formación.
• Espesor neto.
• Propiedades de los fluidos, solo si son datos obtenidos en laboratorio.
• Compresibilidad de la roca, solo cuando se obtienen de datos experimentales.
• Presión capilar.
• Presión ajustada al punto de referencia.
• Contacto original de los fluidos.
• Caudales de producción.

Las variables del reservorio que llevan un grado mayor de incertidumbre son:

• Permeabilidad direccional.
• Tranmisibilidad en las fallas.
• Relación Kv/Kh.
• Curvas de permeabilidad relativa.
• Propiedades del acuífero.
• Propiedades de la roca y fluidos cuando son obtenidas de correlaciones.
• Parámetros de flujo del pozo.

Para el ajuste histórico se recomienda seguir lineamientos generales pero no

determinativos, las etapas para el ajuste del modelo se presentan a continuación:

(1) Ajustar la presión (campo o presión del reservorio)

(2) Ajustar la relación gas petróleo, corte de agua, etc.

83
(3) Ajustar presiones de flujo e índice de productividad.

El proceso de ajuste es individual y requiere de soluciones particulares para cada caso, pero se puede
determinar varias recomendaciones generales:

• Mantener el modelo lo más simple.

• Mantener el modelo físicamente posible.
• Proceder con el ajuste desde el comportamiento del reservorio en general para continuar con el ajuste
individual de pozos.

El ajuste de la presión del reservorio se realiza con las siguiente acciones:

 Especificar la producción total del reservorio para cada pozo. Para asegurar el volumen de
caudal total del reservorio es necesario ajustar al modo de control a “Liquid Rate” (LRAT).

 El ajuste de la presión se realiza al variar parámetros como: volumen poroso, compresibilidad

total, tamaño del acuífero y transmisibilidad.

 Otras propiedades que son consideradas comúnmente para el ajuste son: tamaño del acuífero,
transmisibilidad e IP del acuífero, compresibilidad total, porosidad y tamaño del reservorio.

El ajuste de los caudales de producción se obtiene al hacer cambios en la curva de permeabilidad relativa
del petróleo y agua.

Cuando se tiene un corte de agua muy alto pero un tiempo de ruptura correcto, disminuir la relación
Krw/Kro y el tamaño del acuífero.

Por el contrario cuando se tiene el tiempo de ruptura adelantado al histórico, considerar el ajuste de
lutitas y barreras alrededor del pozo, la permeabilidad vertical en el pozo.

Predicciones

El procedimiento de predicción es la etapa diseñada para la generación de escenarios con el fin de

determinar el comportamiento futuro del reservorio.

Con las corridas de simulación se predice el comportamiento del reservorio (campo o pozos) a partir del
ajuste de los datos históricos. De esta manera, la simulación se convierte en una herramienta de
predicción para el desarrollo del campo (perforación, work overs, mantenimiento de presión, etc).

ECLIPSE se presenta por un icono en el escritorio del computador, que al ser activado es posible visualizar
el recopilador o Launcher en donde se despliegan sus aplicaciones.

84
 Bibliografía

1. “Caracterización Integrada de Yacimientos Petroleros”, México, D.F. Octubre 2014

2. Modelos de simulación de yacimientos basados en leyes de conservación
hiperbólica, aplicando el método de seguimiento de frentes sobre líneas de
corriente. México, D.F., Ciudad Universitaria 2016.
3. Caracterización Estática de Yacimientos Petroleros a Partir de Análisis de Muestras
de Roca. México, D.F. Cd. Universitaria Noviembre 2004.
4. Flujo de fluidos en medios porosos, primera edición. Facultad de ciencias de la
ingeniería.
5. Crichlow, H.: “Modern Reservoir Engineering: A Simulation Approach”, Prentice
Hall Inc., Englewood Cliffs, N.J., 1977.
6. Rivera, V. G.: “Simulación de Yacimientos”, Seminario de Ingeniería de Petróleos,
Universidad de América, Santafé de Bogotá, 1981.
7. Thomas, G.W.: “Principles of Hydrocarbon Reservoir Simulation”, International
Human Resources Development Corporation, Boston, 1982.

85