Simulacion Triple Porosidad

UNIVERSIDAD NACIONAL AUTONOMA DE MEXICO
PROGRAMA DE MAESTRIA Y DOCTORADO EN

INGENIERIA
FACULTAD DE INGENIERIA
SIMULACIN DE YACIMIENTOS NATURALMENTE

FRACTURADOS CON TRIPLE POROSIDAD
T E S I S
QUE PARA OPTAR POR EL GRADO DE:
MAESTRO EN INGENIERIA
CAMPO DE CONOCIMIENTO - CAMPO DISCIPLINARIO
P R E S E N T A :
ERIKA FABIOLA CORTES RUBIO

TUTOR:
DR. VCTOR HUGO ARANA

2008
NDICE
RESUMEN............................................................................................................... i
LISTA DE FIGURAS...............................................................................................
ii
LISTA DE TABLAS................................................................................................. iii
CAPTULO 1. INTRODUCCIN.. 1
CAPTULO 2. FORMULACIN MATEMTICA DEL PROBLEMA.....
2.1 Consideraciones Geolgicas.........................................................................
2.1.1 Porosidad.
2.1.2 Compresibilidad.
2.1.3 Presin Capilar...
2.1.4 Permeabilidad.....
6
7
9
11
2.2 Modelo Conceptual y Matemtico....
13
2.2.1 Modelo Conceptual.... 16

2.2.2 Modelo Matemtico... 19
2.2.2.1 Ecuaciones Diferenciales para el Medio Poroso 1...........
19
20
21
2.3 Condiciones Iniciales y de Frontera....
23
2.3.1 Condiciones Iniciales.......

2.3.2 Condiciones de Frontera.....
23
25
CAPTULO 3. SOLUCIN NUMRICA DEL PROBLEMA....................................
26
3.1.
Ecuaciones de Flujo en Diferencias Finitas..................................
26
3.1.1. Ecuaciones para el Medio 1................................................

26
30
31
3.2 Solucin del Sistema de Ecuaciones.................................................
32
3.2.1 Mtodo Iterativo de Newton-Raphson................................... 32

3.2.2 Reduccin Matricial del Sistema de Ecuaciones................. 40
3.2.3 Solucin del Sistema Lineal de Ecuaciones........................ 43
3.3 Implementacin del Modelo Conceptual en el Simulador.....
47
CAPTULO 4. VALIDACIN DEL MODELO.......
49
4.1 Solucin para el Flujo Radial Monofsico en un Medio....................

4.2 Solucin para Flujo Radial Monofsico en Dos Medios...................
4.3 Solucin Flujo Radial Monofsico en Tres Medios...........................
4.4 Resultados de la Simulacin de Flujo Monofsico en Uno,
Dos y Tres Medios..............................................................................
49
52
54
57
CAPTULO 5. RESULTADOS........
60
5.1 Caso 1. Simulacin de un Sistema de Triple Porosidad por

Medio del Modelo de Doble Porosidad ........................ 60
5.2 Caso 2. Simulacin de Flujo Multifsico en un Sistema de
Triple Porosidad ...................................................... 63
5.4 Caso 3. Simulacin de Flujo Multifsico en un Sistema de
Triple Porosidad ............................................ 74
CAPTULO 6. CONCLUSIONES...........
81
NOMENCLATURA.................................................................................................. 83
REFERENCIAS....................................................................................................... 86
APNDICE A. OBTENCIN DE LAS ECUACIONES DE FLUJO

MULTIFSICO...........
88
APNDICE B. APROXIMACIN EN DIFERENCIAS FINITAS....
92
APNDICE C. ECUACIONES Y DERIVADAS PARA UN MODELO DE

TRIPLE POROSIDAD....
97
RESUMEN
Aunque el modelado de los yacimientos naturalmente fracturados ha sido desarrollado

a lo largo de ms de cinco dcadas, sigue y seguir habiendo una intensa investigacin
para obtener un modelamiento prctico y robusto de este tipo de yacimientos. El
desarrollo ptimo de un yacimiento naturalmente fracturado requiere de un
entendimiento adecuado de la fsica del flujo en el yacimiento y de su implementacin
correcta en un simulador numrico de yacimientos. Algunos ajustes de historia
realizados por medio de la simulacin numrica han mostrado que la clsica
formulacin de doble porosidad utilizada desde finales de los sesentas, no siempre
muestra buenos resultados. Cuando un yacimiento naturalmente fracturado presenta
ms de dos medios porosos afectando el flujo de fluidos, la formulacin de doble
porosidad no puede ser usada. Los medios debern ser agrupados de alguna manera
para reproducir dos medios y de esta forma predecir el comportamiento.
El objetivo del presente trabajo es presentar un modelo simple y prctico para simular
numricamente el flujo multifsico en yacimientos petroleros de triple porosidad, e
incorporar el mismo en un simulador multifsico-multidimensional (11), para llevarlo de
doble a triple porosidad, incrementando as la versatilidad del mismo.
El algoritmo desarrollado en el presente trabajo es una extencin a flujo multifsico de
los modelos propuestos por Camacho y cols (8) y por Rodrguez y cols (27). Debido a
que el modelo matemtico es aplicable a cualquier medio poroso, los tres medios
porosos que representan la triple porosidad pueden ser la combinacin de varias
escalas de fracturas, y/o vgulos, y /o matrices. Se asume que el modelo de flujo
pseudo-estacionario (modelo de Warren y Root 35) que se utiliza, describe el flujo
multifsico de fluidos en estos medios porosos. En este trabajo, el sistema ms
externo es considerado continuo y, por lo tanto, es el nico que se encuentra
conectado a los pozos productores.
Una vez establecidas las ecuaciones que describen el flujo multifsico en un sistema de
triple porosidad y despus de aplicar el mtodo de Newton al sistema discretizado de
ecuaciones no lineales, considerando un ordenamiento normal, se obtiene una
estructura matricial tri-diagonal, la cual se reduce aplicando cclicamente el
complemento de Schur.(40), economizando tiempo y memoria de cmputo. El sistema
matricial obtenido es similar al que se obtiene en el flujo multifsico de aceite, gas y
agua en un yacimiento de una porosidad-una permeabilidad. El trabajo est
estructurado de la manera siguiente: el captulo uno presenta una introduccin sobre
el desarrollo y evolucin de los modelos de simulacin, el captulo dos presenta la
formulacin matemtica del modelo propuesto as como sus condiciones iniciales y de
frontera, el captulo tres presenta la solucin numrica del problema y la
implementacin del modelo en un simulador multifsico y multidimensional, el captulo
cuatro presenta la validacin del modelo usando la solucin analtica publicada por
Camacho y cols (8) y por Rodrguez y cols (27), el captulo cinco presenta la
aplicacin del modelo en tres casos sintticos, y finalmente, el captulo seis presenta
las conclusiones del presente trabajo.
LISTA DE TABLAS
TABLA
PGINA
Tabla 2.1 Posibles Configuraciones de un sistema de Triple Porosidad.

Tabla 4.1 Datos Generales Utilizados para la Validacin del Modelo
Propuesto.
Tabla 4.2 Datos PVT Utilizados para la Validacin del Modelo Propuesto.
Tabla 5.1 Datos Generales Utilizados para la Simulacin del Caso 1.
Tabla 5.2 Datos Utilizados para la Simulacin del Modelo de Triple
Porosidad.
Tabla 5.3 Datos Utilizados para la Simulacin del Modelo de Doble
Porosidad; Promedio Aritmtico.
Porosidad; Promedio Armnico.
Porosidad; Promedio Geomtrico.
Tabla 5.7 Datos PVT para el Aceite Utilizados para la Simulacin del Caso 2.
Tabla 5.8 Datos PVT para el Gas Utilizados para la Simulacin del Caso 2.
Tabla 5.9 Datos PVT para el Agua Utilizados para la Simulacin del Caso 2.
17
58
iii
58
60
60
61
61
61
63
63
64
64
74
LISTA DE FIGURAS
FIGURA
PGINA
Figura 2.1 a) Espacio Poroso Consolidado (matriz), b) Representacin

Simplificada del Espacio Poroso de Vgulos y Fracturas.
Figura 2.2 Curvas de Presin Capilar para los Procesos de Drene e
Imbibicin.
Figura 2.3 Curvas de Presin Capilar si el Promedio de los Dimetros de
Poros es: a) Uniforme, b) Disperso.
Figura 2.4 Permeabilidad Relativa de un Sistema Agua-Aceite (curvas
de imbibicin).
Figura 2.5 Permeabilidad Relativa de un Sistema Gas-Aceite (curvas de
drene).
Figura 2.6 Idealizacin de un Yacimiento Naturalmente Fracturado
(modelo de Warren y Root).
Figura 2.7 Idealizacin de un Sistema de Triple Porosidad.
Figura 2.8 Representacin Grfica del Modelo de Una Porosidad Una
Permeabilidad.
Figura 2.9 Representacin Grfica del Modelo de Doble Porosidad
Una Permeabilidad.
Figura 2.10 Representacin Grfica del Modelo de Doble Porosidad
Doble Permeabilidad.
Figura 2.11 Representacin Grfica del Modelo de Triple Porosidad
Doble Permeabilidad.
Triple Permeabilidad.
Una Permeabilidad.
Figura 2.14 Distribucin de Presiones en el Yacimiento.
Figura 3.1 Diagrama de flujo del Simulador Camalen (11).
Figura 4.1 Resultados Obtenidos para las Soluciones Numrica y
Analtica.
Figura 5.1 Resultados de la Simulacin para los Modelos de Triple y
Doble Porosidad.
Figura 5.2 Grficas de Presin Capilar y Permeabilidad Utilizados en la
Simulacin del Caso 2.
Figura 5.3 Grfica de pwf contra Tiempo para la Simulacin de Flujo
Multifsico en un Sistema de Triple Porosidad con y sin Presiones
Capilares.
Figura 5.4 Grfica de Qo contra Tiempo para la Simulacin de Flujo
Capilares.
Figura 5.5 Grfica de Qg contra Tiempo para la Simulacin de Flujo
Capilares.
iv
10
11
13
13
15
15
17
17
17
18
18
18
24
48
59
62
65
67
67
68
Figura 5.6 Grfica de RGA contra Tiempo para la Simulacin de Flujo

Capilares.
Figura 5.7 Grfica de Qw contra Tiempo para la Simulacin de Flujo
Capilares.
Figura 5.8 Resultados de Presin de Aceite para el Medio 1 con y sin Pc.
Figura 5.9 Resultados de Saturacin de Agua para el Medio 1 con y sin
Pc.
Figura 5.10 Resultados de Saturacin de Gas para el Medio 1 con y sin
Pc.
Figura 5.11 Resultados de Saturacin de Aceite para el Medio 1 con y
sin Pc.
Figura 5.12 Malla Utilizada en la Simulacin del Caso 3.
Figura 5.13 Grficas de Permeabilidad Utilizadas en la Simulacin del
Caso 3.
Figura 5.14 Grficas de Presin Capilar Utilizadas en la Simulacin del
Caso 3.
Figura 5.15 Grfica de pwf contra Tiempo para la Simulacin de Flujo
Multifsico en un Sistema de Triple Porosidad (Caso 3).
Figura 5.16 Grfica de Qo contra Tiempo para la Simulacin de Flujo
Figura 5.17 Grfica de Qg contra Tiempo para la Simulacin de Flujo
Figura 5.18 Grfica de Qw contra Tiempo para la Simulacin de Flujo
Figura 5.19 Grfica de RGA contra Tiempo para la Simulacin de Flujo
Figura 5.20 Distribucin de Presiones de Aceite, po (Caso 3).
Figura 5.21 Distribucin de Saturaciones de Aceite, So (Caso 3).
Figura 5.22 Distribucin de Saturaciones de Gas, Sg (Caso 3).
Figura 5.23 Distribucin de Saturaciones de Agua, Sw (Caso 3).
68
69
70
71
72
73
74
75
75
76
76
77
78
78
79
79
80
80
Captulo 1
Introduccin
CAPTULO 1
INTRODUCCIN
La modelacin de Yacimientos Naturalmente Fracturados (YNF) ha sido, desde la

formulacin propuesta por Barenblat.4, un reto permanente para los investigadores
involucrados en la explotacin de yacimientos de este tipo. La descripcin correcta de
los YNF, junto con una descripcin sensata de los procesos fsicos del flujo de fluidos
que se llevan a cabo, representan la base para obtener pronsticos de produccin lo
ms cercanos a la realidad.
Desde los aos 60s hasta la actualidad, han aparecido diversas formulaciones para
conceptualizar un YNF. Las aproximaciones matemticas desarrolladas para modelar el
flujo de fluidos en yacimientos fracturados, se basan en aproximaciones continuas que
consisten en: el establecimiento de un modelo conceptual, la incorporacin de la
informacin geomtrica de un determinado sistema matriz-fractura, del
establecimiento de ecuaciones de conservacin de masa y energa para los dominios
de la matriz y la fractura y, finalmente, de la solucin del sistema de ecuaciones
discretas no lineales.
Los formulaciones comnmente utilizadas para modelar el flujo de fluidos en rocas
fracturadas, incluyen:1) Red de fracturas discretas, 2) Modelo de doble medio
continuo, incluyendo doble y multiporosidad, doble permeabilidad , el mtodo ms
general de Mltiple Interaccin Continua (MINC) 25, y 3) Mtodo de Continuidad
Efectiva (ECM), pasando por sistemas equivalentes y continuos, hasta la
representacin de fracturas mediante distribuciones fractales [Diodato 10].
De los mtodos anteriormente mencionados el ms utilizado ha sido el de doble medio
continuo, debido a que es menos demandante computacionalmente que otras
aproximaciones. Bsicamente los diferentes modelos pueden distinguirse por la forma
como es considerado el almacenamiento y la continuidad de los fluidos.
Todas las formulaciones poseen ventajas y desventajas entre ellas, las cuales
dependen de la aplicacin en cuestin. En lo que respecta a la aplicacin en ingeniera
petrolera, la conceptualizacin que ha mostrado una mejor aplicacin a escala de
campo, ha sido la formulacin de doble medio continuo [Warren y Root 35]. Lo
anterior no presupone que esta formulacin represente correctamente todos los
procesos dinmicos que se llevan acabo en los YNF.
En general, un YNF ha sido conceptualizado como la interaccin de dos medios o
sistemas porosos. Un medio est representado por la matriz; es decir, la porcin de la
roca que no presenta fracturas y el otro medio est representado por las fracturas
que intercambian fluidos de alguna manera con la matriz y que abstractamente
coinciden en espacio y tiempo. Los bloques de matriz estn rodeados de fracturas, no
existe comunicacin entre los bloques de matriz y la red de fracturas es continua.
Captulo 1
Introduccin
En la simulacin de los yacimientos fracturados, nos encontramos con el hecho de que

los bloques de matriz pueden contener alrededor del 90% de la reserva total de aceite.
El principal problema para la recuperacin de aceite en un yacimiento fracturado es
precisamente la extraccin de aceite de estos bloques de matriz. Por lo anterior, es
importante comprender los mecanismos que se llevan a cabo en los bloques de matriz
y simularlos lo ms precisamente posible [Tadjar 21].
Algunos YNF muy probablemente presentan ms de dos sistemas porosos. Sin
embargo, debido a que no existen actualmente algoritmos para simular YNF con tres
sistemas porosos de diferentes propiedades, de alguna manera se agrupan los
sistemas identificados para seguir aplicando el modelo de doble porosidad y de esta
manera realizar predicciones.
De acuerdo con la teora de multi-porosidad [Aifantis 2], cualquier medio que exhibe
discontinuidades finitas en la porosidad, se considera que posee propiedades multiporosidad, esto es, si un sistema no es homogneo, entonces puede presentar dos o
ms porosidades.
Para poder establecer cual es el modelo ms adecuado para representar a un
yacimiento de hidrocarburos y obtener una exitosa prediccin del flujo de fluidos en el
yacimiento, es necesario definir las caractersticas de permeabilidad y porosidad de los
tipos de roca que existen, ya que, como se mencion anteriormente, los diferentes
modelos se distinguen por la forma como es considerado el almacenamiento y la
continuidad de los fluidos. Por ejemplo, el modelo de una porosidad-una
permeabilidad, se aplica en general a yacimientos no fracturados con porosidad y
permeabilidad uniformes. Un yacimiento fracturado con permeabilidad de matriz
relativamente baja pero alto almacenamiento puede ser caracterizado por medio del
modelo de doble porosidad-una permeabilidad. Un modelo comnmente aceptado es el
de doble porosidad-doble permeabilidad, en el cual, las porosidades y las
permeabilidades de la matriz y de la fractura son significativamente diferentes y
continuas.
Adems del concepto tradicional de doble porosidad, diversos modelos de triple
porosidad o triple continuidad han sido propuestos para describir el flujo a travs de
rocas fracturadas. En general, estos mtodos se han enfocado en manejar las
heterogeneidades de la matriz subdividiendo esta en dos o ms subdominios, con
diferentes propiedades para cada uno de los medios.
Actualmente tambin se han propuesto modelos multiporosidad-multipermeabilidad
[Aifantis 2] para representar diversos yacimientos fracturados con: permeabilidad
moderada (modelo de triple porosidad-doble permeabilidad) o permeabilidad alta
(modelo de triple porosidad-triple permeabilidad).
Entre los modelos multiporosidad-multipermeabilidad que se han publicado en la
literatura se pueden mencionar los siguientes:
Captulo 1
Introduccin
En 1985, Pruess y Narasimhan 25 desarrollaron un mtodo de Multiple Interaccin

Continua, MINC, el cual manejaba el flujo transiente multifsico y multidimensional,
tanto en las fracturas como en los bloques de matriz, por medio de una aproximacin
numrica. La suposicin principal del mtodo MINC es que las condiciones
termodinmicas en la matriz dependen principalmente de la distancia entre esta y la
fractura ms cercana. En dicho mtodo, el dominio de flujo poda ser particionado en
elementos de volumen de tal forma que todas las interfases entre los elementos de
volumen dentro de la matriz, fueran paralelos a la fractura ms cercana. Refinando los
bloques de matriz con base en la distancia a las fracturas, se tiene como resultado un
patrn de elementos de volumen anidados. Cada elemento de volumen tiene un estado
termodinmico asignado al mismo. El concepto bsico del MINC de particionar los
bloques de matriz de acuerdo a la distancia de las caras de la fractura puede
extenderse a tamaos y formas irregulares.
Un modelo de triple porosidad-una permeabilidad fue propuesto por Abdassah y
Ershaghi 1. Estos autores consideraron un modelo de flujo no estacionario entre el
sistema de fracturas, con dos tipos de bloques de matriz, y flujo primario solo a travs
del sistema de fracturas. Ellos consideraron flujo paralelo entre un sistema de
fracturas, que tienen propiedades homogneas, e interactan con dos grupos de
bloques de matriz separados, que tienen distintas permeabilidades y porosidades.
En 2002, Camacho y cols 8 presentaron un estudio para modelar porosidades
secundarias, principalmente porosidad vugular, en yacimientos naturalmente
fracturados. Presentaron soluciones para dos diferentes casos: cuando no existe flujo
primario a travs de los vgulos, lo cual es una extensin del modelo de Warren y
Root 33, y para el caso en que el proceso de disolucin ha creado un sistema
interconectado de vgulos. En ambos casos, existe una interaccin entre la matriz, los
vgulos y el sistema de fracturas.
En 2004, Wu, Y. S. y cols 37, presentaron un modelo conceptual de triple porosidad y
triple permeabilidad, para simular los procesos de flujo en rocas fracturadas. Ellos
conceptualizan el sistema matriz-fractura formado por una matriz y dos tipos de
fracturas: fracturas grandes y fracturas pequeas; y extienden el concepto de doble
permeabilidad por medio de la adicin de una conexin ms (va fracturas pequeas),
entre las fracturas grandes y los bloques de matriz.
En 2004, Rodrguez, F. y cols 27 presentaron una solucin analtica para caracterizar
la porosidad secundaria en yacimientos naturalmente fracturados. Ellos desarrollaron
un modelo de triple porosidad anidada-una permeabilidad para el transiente de presin
de un pozo produciendo en un yacimiento naturalmente fracturado. Se consider un
sistema de triple porosidad actuando a diferentes escalas: matriz, porosidad
secundaria de pequea y gran escala, donde el flujo del fluido a travs de estos medios
es en serie de la siguiente manera: la matriz intercambia fluidos con la porosidad
secundaria de pequea escala, la cual a su vez alimenta la porosidad secundaria de
gran escala.
Captulo 1
Introduccin
En 2006 Kang, Z. y cols 16, presentaron e implementaron un modelo conceptual

multi-continuo. El modelo conceptual, propuesto para investigar el flujo y
desplazamiento multifsico a travs de yacimientos vugulares naturalmente
fracturados, se basa en la observacin y anlisis de datos geolgicos. En este modelo
conceptual, la roca vugular naturalmente fracturada es considerada un medio de triple
continuidad compuesta por: (1) fracturas de gran escala, de alta permeabilidad y bien
conectadas; (2) roca matriz de baja o nula permeabilidad; y (3) vgulos o cavidades
de varios tamaos. Los vgulos pueden estar: (1) directamente conectados a las
fracturas grandes, (2) indirectamente conectados a las fracturas grandes a travs de
las fracturas pequeas o microfracturas, o (3) aislados de las fracturas grandes por la
roca matriz.
En 2006, Tajdar y cols 21 utilizaron el modelo generalizado de doble porosidad
denominado de Mltiple Interaccin continua, (MINC) y lo implementaron en un
simulador en paralelo para modelar el flujo de fluidos en matriz y en fractura. De
acuerdo al modelo MINC, los bloques de matriz fueron discretizados en subdominios,
tanto en direccin vertical como horizontal para ofrecer una precisin mayor del flujo
transiente en los bloques de matriz. En dicho modelo se asumi la existencia de
equilibrio termodinmico local dentro de cada medio a cualquier nivel de tiempo y,
basndose en esa suposicin de equilibrio local, definieron variables termodinmicas
tales como: presin, saturaciones, concentraciones y temperatura.
En 2006 Yu-Shu y cols 38 desarrollaron un modelo conceptual de mltiple medio
continuo y lo implementaron en un simulador de yacimientos tridimensional y trifsico
utilizando el modelo generalizado de triple continuidad. El modelo conceptual
consideraba las rocas vugulares fracturadas como un medio de triple o mltiple
continuidad, constituido por: 1) fracturas de alta permeabilidad, 2) matriz de baja
permeabilidad y 3) vgulos de diversos tamaos. En dicho trabajo tambin proponen
una metodologa para el clculo de la permeabilidad de los bloques que contienen
vgulos.
En 2007 Rubin 30 desarroll un modelo de Doble subdominio-Doble permeabilidad,
para simular el drene gravitacional debido al flujo directo entre los bloques de matriz
(continuidad capilar) y a la reinfiltracin en la matriz de la fractura. Rubin aplic su
modelo en yacimientos fracturados con capa secundaria de gas.
De los trabajos anteriormente mencionados, es notable el esfuerzo realizado por tratar
de comprender los fenmenos que gobiernan el flujo de fluidos en los yacimientos
naturalmente fracturados. Con el fin de contribuir con dichos esfuerzos, el presente
trabajo presenta un modelo conceptual para yacimientos que presentan triple
porosidad. Los objetivos del presente trabajo son:
1) Desarrollar un modelo de triple porosidad-una permeabilidad, para simular
numricamente el flujo multifsico y multidimensional en yacimientos
petroleros que presentan triple porosidad.
2) Validar el modelo desarrollado contra soluciones analticas.
Captulo 1
Introduccin
3) Implementar el modelo propuesto en un simulador trifsico, tridimensional, de

doble porosidad.
4) Utilizar el simulador en casos de aplicacin (sintticos).
Una vez incorporado el
caractersticas siguientes:
modelo
conceptual
al
simulador,
ste
presenta
las
Tres Dimensiones
Totalmente implcito
Tres Fases
Triple Porosidad
En el modelo propuesto para ser implementado en el simulador, se considera un

sistema compuesto por tres medios, en donde el medio uno posee propiedades bien
definidas e interacta con otros dos medios que poseen diferentes porosidades y
permeabilidades. El medio uno se considera continuo, mientras que los medios dos y
tres proporcionan almacenamiento pero no contribuyen con el transporte y flujo de los
fluidos. Para el modelo propuesto, un sistema con triple porosidad est compuesto por
tres medios: uno continuo y dos discontinuos anidados, similar a los modelos
propuestos por Camacho y cols (8) y por Rodrguez y cols (27).
El trabajo presente est estructurado de la manera siguiente:
En una primera parte se presenta una introduccin sobre el desarrollo y

evolucin de los modelos de simulacin.
En una segunda parte se establece la formulacin matemtica del problema y

se asignan condiciones iniciales y de frontera.
En una tercera parte se realiza la solucin numrica del problema que consiste
en convertir las ecuaciones diferenciales en ecuaciones en diferencias finitas, se
aplica el mtodo de Newton para generar un sistema algebraico de ecuaciones
lineales, se realiza una reduccin matricial del sistema de ecuaciones aplicando
el mtodo del Complemento de Schur y se resuelve el sistema de ecuaciones
lineales mediante la aplicacin de un algoritmo para la solucin de matrices
dispersas. En una tercera parte se incorpora el modelo conceptual a un
simulador trifsico-tridimensional.
En una cuarta parte se valida el modelo con soluciones analticas.
En una quinta parte se utiliza el simulador en casos prcticos.
En una sexta parte se realizan las conclusiones del presente trabajo.
Finalmente cabe mencionar que aunque implcitamente el algoritmo presente algunas

desventajas en el sentido de adquisicin de informacin, tambin presenta valiosas
ventajas: hay un mayor grado de libertad para poder realizar ajustes de historia en
YNF con triple porosidad.
Captulo 2
Formulacin Matemtica del Problema

CAPTULO 2
FORMULACIN MATEMTICA DEL PROBLEMA
2.1 Consideraciones Geolgicas

2.1.1 Porosidad
Las caractersticas petrofsicas de una roca carbonatada estn controladas bsicamente
por dos redes de poros: una red de poros interparticular y una red de poros vugular.
La Porosidad Interparticular puede definirse como el espacio poroso localizado entre las
partculas de la roca, el cual no es significativamente mayor al tamao propio de las
partculas. La porosidad interparticular entre granos multicristalinos se denomina
porosidad intergranular y entre partculas de cristal simple se denomina porosidad
intercristalina. La Porosidad Vugular puede definirse como ese espacio poroso de
mayor o igual tamao que las partculas de la roca y comnmente se presenta como
fracturas o grandes cavidades irregulares. Los vgulos se encuentran interconectados
generalmente de dos formas: (1) a travs de redes de poros interparticulares (vgulos
separados), y (2) por contacto directo con otro vgulo (vgulos conectados), [Lucia
20
].
A la Porosidad Interparticular se le denomina porosidad primaria y es tpico de las
arenas y las calizas. A la Porosidad Vugular se le denomina porosidad secundaria o,
cuando se refiere solo a vgulos o fracturas, porosidad vugular o porosidad de fractura
(figura 2.1).
Vgulos y Fracturas
Figura 2.1 a) Espacio Poroso Consolidado (matriz), b) Representacin Simplificada del Espacio Poroso de Vgulos y Fracturas
La porosidad secundaria generalmente se encuentra en rocas compactas, quebradizas

y con porosidad intergranular relativamente baja, tales como calizas compactas. La
porosidad secundaria se origina normalmente por el fracturamiento de la roca, el
agrietamiento y la disolucin debido al agua circulante [Van Golf-Racht 34].
Para estimar los parmetros petrofsicos por medio de examen visual en rocas con
ausencia de vgulos conectados, debe determinarse la cantidad de porosidad
interparticular, el tamao de partcula y la cantidad de porosidad vugular separada.
Tambin es importante determinar la presencia o ausencia de vgulos conectados. Si
Captulo 2
estn presentes, la permeabilidad puede ser mucho ms alta y las propiedades

capilares muy diferentes a las que pudieran predecirse si no hubiera vgulos
conectados [Lucia 20]. Adems de lo anterior, si existe presencia de vgulos
conectados, el modelo a utilizar sera el de doble permeabilidad.
Sistema de Mltiple Porosidad
En un yacimiento heterogneo, la porosidad total (t), es el resultado de la simple
adicin de las porosidades primarias y secundarias, esto es:
t = 1 + 2 + ... + n .
(2.1)
Esta porosidad total es equivalente a la definicin esttica de almacenamiento de roca

o volumen poroso total.
Evaluacin Cuantitativa de la Porosidad Secundaria
En general, la evaluacin de la porosidad total se realiza por medio de un
procedimiento indirecto (registros geofsicos), o por medio de un procedimiento directo
(anlisis de ncleos).
Los valores de porosidad secundaria pueden variar si lo que se encuentra presente
son vgulos o fracturas. Resultados diferentes se esperan tambin cuando las fracturas
son de un tipo o de otro; la ms importante subdivisin de fracturas relaciona dos
categoras: macrofracturas y microfracturas. Las macrofracturas son fracturas
extendidas con grandes aperturas, las cuales se desarrollan a travs de varias capas;
mientras que las microfracturas (o fisuras), son fracturas con aperturas estrechas y de
extensin limitada, a menudo limitadas a una sola capa [Van Golf 34].
2.1.2 Compresibilidad
En los yacimientos heterogneos, la compresibilidad tiene un papel muy importante,
especialmente si existe gran contraste entre la porosidad de cada uno de los medios.
La compresibilidad, se define en forma general, como el cambio de volumen, V , por
unidad de volumen, V , para una presin aplicada, p ; esto es:
c=
1 V
.
V p
(2.2)
Captulo 2
Dependiendo del volumen que se utilice como referencia, la compresibilidad puede ser
la propiedad de un cierto volumen de roca sometido a compresin, tal como el
volumen total ( VB ), el volumen poroso ( V p ) o volumen de fluido ( V f ).
Compresibilidad de la Roca
La compresibilidad de la roca en un yacimiento convencional se refiere al volumen
total, de roca y poroso; esto es:
Vb = Vr + V p
;
(2.3)
V p = Vr .
y se expresan por medio de su definicin bsica, como:
Compresibilidad Total, cb =
1 dVb
Vb dp
Compresibilidad de Roca, c r =
1 dVr
;
Vr dp
Compresibilidad de Poro, c p =
1 dV p
V p dp
(2.4)
En el caso de yacimientos heterogneos, adems de la porosidad primaria de la matriz,

es necesario tomar en consideracin la porosidad secundaria de las fracturas (f),
cavernas (c) y vgulos (v). La compresibilidad de la roca se expresa entonces por
medio de:
c p = cm + c cc + f c f + v cv .
(2.5)
Compresibilidad Total (roca y fluido)

La compresibilidad del sistema, incluyendo la roca y los fluidos que saturan los poros,
puede expresarse como:
cT = So co + S wcw + S g cg + c p .
(2.6)
Captulo 2
2.1.3 Presin Capilar

La presin capilar existe siempre que los poros estn saturados con dos o ms fases.
En un sistema multifsico, la presin capilar es, por definicin, la presin de la fase no
mojante menos la presin de la fase mojante.
Ya que las fuerzas interfaciales, manifestadas como Presiones Capilares, son las
mayores fuerzas dentro del volumen representativo del medio poroso, la presin
capilar es una caracterstica bsica del sistema roca-fluido en el flujo multifsico.
Entre los factores que afectan esta variable, se encuentran: la tensin interfacial entre
los fluidos, la mojabilidad de la roca, la geometra del espacio poroso, as como la
historia de saturaciones (drene o imbibicin).
Para sistemas gas-aceite, la expresin para determinar la presin capilar est dada
por:
Pcgo ( S g ) = pg po .
(2.7)
Para sistemas agua-aceite, la expresin para determinar la presin capilar est dada
por:
Pcwo ( S w ) = po pw .
(2.8)
Curvas de Presin Capilar

Si una roca del yacimiento est completamente saturada por la fase mojante, el
desplazamiento de esta fase por medio de una fase no mojante, estar relacionado con
la magnitud de la presin capilar y el correspondiente decremento de la saturacin de
la fase mojante. Este desplazamiento se conoce como desplazamiento de drene, y la
relacin entre la presin capilar y la saturacin es conocida como curva de presin
capilar de drene.
La presin capilar de imbibicin se desarrolla cuando la roca del yacimiento est
saturada con una fase no mojante, la cual debe ser desplazada por medio de una fase
mojante. Para una distribucin de poro convencional, la presin capilar de imbibicin
es aproximadamente la mitad en el orden de magnitud de la presin de drene y, por lo
tanto, la curva mostrar una histresis considerable entre los desplazamientos de
drene e imbibicin para una misma roca (figura 2.2).
Presin Capilar, Pc
Captulo 2
Drene
Imbibicin
Saturacin de la Fase Mojante
100%
Figura 2.2 Curvas de Presin Capilar para los Procesos de Drene e Imbibicin
La forma de la curva capilar del proceso de drene refleja la homogeneidad del tamao
de poro. Si en un gran intervalo de saturacin de agua la Pc permanece constante
(figura 2.3a), esto prueba la homogeneidad de la dimensin del poro y viceversa, la
variacin de Pc contra saturacin (figura 2.3b), prueba la heterogeneidad de la
dimensin del poro.
Por otro lado, solo las curvas de imbibicin muestran la existencia de una saturacin
particular en la saturacin de la fase no mojante (en el caso de agua-aceite, S ro ). Esta
saturacin residual ocurre durante el desplazamiento de imbibicin, cuando la fase no
mojante queda atrapada llenando poros individuales o pequeos grupos de poros. Este
entrampamiento es el resultado de la interrupcin parcial de la comunicacin entre
varios poros que son evitados durante el desplazamiento de imbibicin. Es fsicamente
evidente que hay mayores probabilidades de evitar ms poros si existe una
distribucin mayor en el tamao de poro, donde las variaciones de las dimensiones de
los poros, reducen la inter-comunicacin. Por lo tanto, ocurrir una mayor saturacin
residual de la fase no mojante durante el desplazamiento de imbibicin si la
distribucin de poro (no normal), es similar al caso de la figura 2.3b.
10
Captulo 2
Swi
Swi
Sor
Sor
b)
a)
Figura 2.3. Curvas de Presin Capilar si el Promedio de los Dimetros de Poros es: a) Uniforme, b) Disperso
2.1.4 Permeabilidad
La permeabilidad es una propiedad del medio poroso y es una medida de la capacidad
del medio para transmitir fluidos. En presencia de dos o ms medios porosos, la
permeabilidad puede definirse como: permeabilidad del medio 1, permeabilidad del
medio 2,, permeabilidad del medio n.
Permeabilidad Relativa
Originalmente, la ley de Darcy se deriv para el flujo de una fase donde la
permeabilidad del medio poroso al fluido es la permeabilidad absoluta del medio.
Cuando dos o ms fluidos fluyen simultneamente a travs del medio poroso, la ley de
Darcy puede modificarse para calcular el gasto de cada fase. Estas modificaciones
incluyen: el uso de la permeabilidad efectiva, el potencial y la viscosidad de cada fase.
Con estas modificaciones, el gasto de aceite y agua en un sistema bifsico y en la
direccin x , puede escribirse como:
qox =
kox Ax d o
o dx
(2.9)
y,
11
Captulo 2
qwx =
kwx Ax d w
,
w dx
(2.10)
donde las permeabilidades efectivas,
kox y kwx dependen de la saturacin. La
permeabilidad efectiva del aceite, puede expresarse como:
k
kox = k x ox = k x krow ,
kx
donde
(2.11)
k x =permeabilidad absoluta del medio poroso en la direccin
krow =permeabilidad relativa al aceite. Similarmente para el agua, se tiene:

kwx = k x krw .
(2.12)
Sustituyendo las ecuaciones (2.11) y (2.12) en (2.9) y (2.10), se tiene:
qox = k x Ax
krow d o
,
o dx
(2.13)
krw d w
.
w dx
(2.14)
y,
qwx = k x Ax
Similarmente para el flujo de aceite y gas, se tiene:
qox = k x Ax
krog d o
o dx
(2.15)
y,
qgx = k x Ax
krg d g
g dx
(2.16)
krow y krw en funcin de la

y krg como funciones de la
Las figuras 2.4 y 2.5 muestran las permeabilidades relativas

saturacin,
S w , para un sistema agua-aceite y, krog
saturacin,
S g , para un sistema gas-aceite, respectivamente. Ntese que la suma de
12
Captulo 2
krow y krw para cualquier S w (y krog y krg para cualquier S g ), es siempre menor o
igual a uno. Este ltimo efecto es causado por la tensin interfacial entre los dos
fluidos residentes.
Saturacin de Gas
Saturacin de Agua
Figura 2.5 Permeabilidad Relativa de un Sistema

Gas-Aceite (curvas de drene).
Figura 2.4 Permeabilidad Relativa de un Sistema

Agua-Aceite (curvas de imbibicin).
En resumen, para flujo multifsico en medios porosos, los gastos de flujo de las
diferentes fases, estn dados por:
qox = k x Ax
qgx = k x Ax
qwx = k x Ax
donde
kro d o
;
o dx
krg d g
(2.17)
(2.18)
krw d w
.
w dx
(2.19)
g dx
krw , krg y kro
son las permeabilidades relativas del agua, gas y aceite
respectivamente, en un sistema trifsico.
13
Captulo 2
2.2 Modelo Conceptual y Matemtico

Para un mejor entendimiento del modelo de Triple Porosidad, es necesario definir el
modelo de doble porosidad.
Modelo de Doble Porosidad
El modelo de doble porosidad fue introducido por Barenblatt y cols 4, y ms tarde
extendido por Warren y Root 35. Los modelos de doble porosidad estn basados en
un medio de flujo idealizado (figura 2.6), que consta de una porosidad primaria creada
por depositacin y litificacin, y una porosidad secundaria creada por fracturamiento
[Warren y Root 35]. La base de estos modelos es la observacin de que la masa de
roca no fracturada (matriz), presenta gran parte de la porosidad del medio
(almacenamiento) pero poca permeabilidad (flujo). Por otro lado, la fractura puede
presentar poco almacenamiento pero alta permeabilidad. La matriz y la fractura estn
idealizados como dos medios separados pero que se encuentran interactuando en
espacio y tiempo, donde la transferencia de fluidos ocurre de acuerdo al potencial del
fluido entre los dos medios.
En el modelo de doble porosidad, un yacimiento fracturado se concibe formado por dos
sistemas: uno de fracturas, que constituye un medio continuo; y un sistema de
bloques de matriz, que constituye un medio discontinuo. Adicionalmente se realizan las
suposiciones siguientes con respecto a la estructura del yacimiento y a la formulacin
matemtica.
1. El sistema continuo de fracturas proporciona la trayectoria principal de flujo del
fluido en el yacimiento. El fluido que es desplazado de los bloques de matriz
fluye cerca de la interfase matriz-fractura, y entonces fluye a travs del sistema
de fracturas hacia el pozo. Lo anterior significa que los pozos se encuentran
ubicados nicamente en el sistema de fracturas.
2. La ley de Darcy es aplicable para describir el flujo multifsico en la fractura.
3. Los sistemas de fracturas y de bloques de matriz ocupan una misma malla
computacional, y por lo tanto, para cada conjunto de soluciones para el sistema
de fracturas, se obtiene un conjunto de soluciones para el sistema de bloques
de matriz.
La esencia de esta aproximacin es que en el medio poroso fracturado, las fracturas se
caracterizan por tener mayor difusividad (de aqu los pequeos tiempos de respuesta),
que la roca matriz. Por lo tanto, la respuesta temprana del sistema es controlada por
las fracturas, mientras que la respuesta tarda, es influenciada por la matriz.
14
Captulo 2
Vgulos
Matriz
Fracturas
Matriz
Fracturas
Figura 2.6 Idealizacin de un Yacimiento Naturalmente Fracturado (Modelo de Warren y Root)
Modelo de Triple Porosidad

Una vez definido el modelo de doble porosidad, considrese ahora un sistema donde el
medio uno tiene propiedades homogneas e interacta con otros dos medios
separados, que tienen diferentes porosidades y permeabilidades (figura 2.7).
Figura 2.7 Idealizacin de un Sistema de Triple Porosidad (11)
Este sistema se refiere como un sistema de triple porosidad y es una representacin

ms realista de los yacimientos fracturados que el modelo tradicional de doble
porosidad.
15
Captulo 2
2.2.1 Modelo Conceptual

Para poder establecer cual es el modelo ms adecuado para representar a un
yacimiento en particular y de esta forma obtener una prediccin exitosa del flujo de
fluidos en el yacimiento, es necesario definir las caractersticas de permeabilidad y
porosidad del tipo de roca en particular, ya que los diferentes modelos se distinguen
por la forma como consideran el almacenamiento y la continuidad de los fluidos. Por
ejemplo, el modelo de una porosidad-una permeabilidad (figura 2.8), se aplica en
general a yacimientos no fracturados. Un yacimiento fracturado que presenta dos
diferentes porosidades, de las cuales solo una de ellas es continua; esto es, contribuye
con el transporte y flujo de los fluidos, puede caracterizarse por medio del modelo de
doble porosidad-una permeabilidad (figura 2.9). Un modelo comnmente aceptado es
el de doble porosidad-doble permeabilidad (figura 2.10), en el cual, existen dos
porosidades y ambas son continuas. Finalmente, se ha propuesto un modelo
multiporosidad-multipermeabilidad en la literatura [Aifantis.2], para representar
diversos yacimientos fracturados con: permeabilidad moderada (modelo de triple
porosidad-doble permeabilidad, figura 2.11) o permeabilidad alta (modelo de triple
porosidad-triple permeabilidad, figura 2.12).
En el modelo propuesto para la construccin del simulador se considera un sistema
compuesto por tres medios. El flujo de fluidos entre estos tres medios es en serie; esto
es, el medio tres intercambia fluidos con el medio dos, el cual a su vez intercambia
fluidos con el medio uno. El medio uno se considera continuo, mientras que los medios
dos y tres proporcionan almacenamiento pero no contribuyen con el transporte y flujo
de los fluidos (figura 2.13). Bsicamente, el modelo propuesto toma las ideas
planteadas por Camacho y cols (8) y por Rodrguez y cols (27) y las extiende a flujo
multifsico.
Para el modelo propuesto, un sistema con triple porosidad puede estar compuesto por
tres diferentes sistemas porosos con diferentes propiedades petrofsicas, que afectarn
de alguna manera el flujo de fluidos dentro del yacimiento. En la tabla 2.1 se
presentan las posibles configuraciones de un sistema de triple porosidad.
A diferencia del modelo MINC (25) que considera que las diferentes porosidades estn
ordenas en serie y en sentido creciente, donde la porosidad menor ( 1 ), se encuentra
en contacto con una porosidad mayor ( 2 ) y as sucesivamente, de forma tal que
<
(i +1) ;
(i )
en el modelo propuesto las porosidades no tienen un ordenamiento; esto es,
puede estar en contacto con
o con cualquier otra porosidad que integre el
sistema, con la nica restriccin de que la transferencia de fluidos entre los medios sea
en serie. Lo anterior se hizo considerando que un YNF con triple porosidad puede estar
compuesto por fracturas, vgulos y matriz, o tambin por fracturas grandes, fracturas
pequeas y matriz.
Tabla 2.1 Posibles Configuraciones de un sistema de Triple Porosidad
16
Captulo 2
Sistema 1
F
F
F
Sistema 2
f
m
m
Sistema 3
m
v
m
donde:
F
f
v
m
= fracturas grandes
= micro fracturas
= vgulos
= matriz
Esto es, un YNF con triple porosidad puede estar compuesto por fracturas grandes
(sistema 1), fracturas pequeas (sistema 2) y matriz (sistema 3); o bien por fracturas
(sistema 1), matriz (sistema 2) y vgulos (sistema 3), etc.
En el trabajo presente se supone que un sistema multifsico isotrmico en un
yacimiento naturalmente fracturado est compuesto de tres fases: aceite, gas y agua.
Se considera que cada fase fluye en respuesta a las fuerzas de presin, gravitacionales
y capilares, de acuerdo a la ley de Darcy.
i-1
i+1
Figura 2.8 Representacin Grfica del Modelo de Una Porosidad Una Permeabilidad.
i-1
i+1
Figura 2.9 Representacin Grfica del Modelo de Doble Porosidad Una Permeabilidad.
i-1
i+1
Figura 2.10 Representacin Grfica del Modelo de Doble Porosidad Doble Permeabilidad.
17
Captulo 2
i-1
i+1
Figura 2.11 Representacin Grfica del Modelo de Triple Porosidad Doble Permeabilidad.
i-1
i+1
Figura 2.12 Representacin Grfica del Modelo de Triple Porosidad Triple Permeabilidad.
1
i-1
1i
1
i+1
Figura 2.13 Representacin Grfica del Modelo Utilizado en el Presente Trabajo (Triple Porosidad Una Permeabilidad).
18
Captulo 2
2.2.2 Modelo Matemtico
Las ecuaciones que describen el flujo trifsico en un yacimiento naturalmente

fracturado se desarrollan en el Apndice A, y comprenden el conjunto siguiente de
ecuaciones en cada uno de los medios: medio 1, medio 2 y medio 3.
2.2.2.1 Ecuaciones Diferenciales para el Medio Poroso 1
Ecuaciones de Flujo
La descripcin del flujo de fluidos en un medio poroso matemticamente se modela a
travs de ecuaciones que dominan los principios de conservacin de masa, de cantidad
de movimiento (aproximada por la ley de Darcy) y de una ecuacin de estado,
obtenindose as 3 ecuaciones de flujo (una para cada fase). De acuerdo al modelo
propuesto, el flujo de fluidos entre los tres medios es en serie; esto es, el medio tres
intercambia fluidos con el medio dos, el cual a su vez intercambia fluidos con el medio
uno; por lo tanto, no existe transferencia directa entre los medios tres y uno. De
acuerdo a lo anterior, las ecuaciones que describen el flujo de fluidos en el medio uno,
el cual es considerado continuo, estn definidas por las siguientes expresiones:
Aceite:
kk
ro ( po oD ) + bo qo + o 21 = ( bo So ) 1 .
1
t
o Bo
1
(2.20)
Gas:
kk
kkrg
ro Rs ( po oD ) +
p g g D )
(
g Bg
o Bo
1
(2.21)
+ Rsbo qo + bg q g + [ Rso ]21 + g = ( Rsbo So + bg S g )

21
1
1
1
t
Agua:
kk
rw ( pw wD ) + bw q w + [w ]21 = [ bw S w ]1 .
1
t
Bw w
1
19
(2.22)
Captulo 2
En las ecuaciones (2.20) a (2.22), :o21, g21 y w21, representan el ritmo de intercambio
de masa entre el sistema 2 y 1 de las fase aceite, gas y agua respectivamente, a
condiciones de yacimiento por unidad de volumen de roca.
Ecuaciones Auxiliares
Se tienen dos ecuaciones de presin capilar:
Para sistemas gas-aceite, la expresin para determinar la presin capilar est dada
por:
Pc go ,1 ( S g1 ) = pg ,1 po ,1 .
(2.23)
Para sistemas agua-aceite, la expresin para determinar la presin capilar est dada
por:
Pc wo ,1 ( S w1 ) = po ,1 pw,1 .
(2.24)
Ecuaciones de Restriccin
Se tiene una ecuacin de restriccin para las saturaciones de las fases:
(2.25)
So1 + S g1 + S w1 = 1 .
2.2.2.2
Ecuaciones Diferenciales para el Medio Poroso 2
En este caso, las ecuaciones de transferencia no presentan trminos de flujo debido a

que es un medio discontinuo; por lo tanto, nicamente aparecen los trminos de
transferencia entre los medios 1 y 2, y entre los medios 2 y 3.
Ecuaciones de Transferencia
Una para cada fase (agua, aceite y gas):
Aceite:
[o ]21 + [o ]32 =
( bo So ) 2 .
t
(2.26)
20
Captulo 2
Gas:
[ Rso ]21 g + [ Rso ]32 + g =

21
32

( Rsbo So + bg S g ) .
2
t
(2.27)
Agua:
[w ]21 + [w ]32 =
[ bw Sw ]2 .
t
(2.28)
En las ecuaciones (2.26) a (2.28), : o32, g32 y w32, representan el ritmo de

intercambio de masa entre el sistema 3 y 2 de las fase aceite, gas y agua
respectivamente, a condiciones de yacimiento por unidad de volumen de roca.
Ecuaciones Auxiliares.
Expresiones de presin capilar para describir la relacin entre las presiones de las
fases;
Pc go,2 ( S g 2 ) = pg ,2 po,2 .
(2.29)
Pc wo ,2 ( S w 2 ) = po ,2 pw,2 .
(2.30)
Ecuacin de Restriccin.
Una ecuacin de restriccin para las saturaciones de las fases:
So 2 + S g 2 + S w 2 = 1 .
2.2.2.3
(2.31)
Ecuaciones Diferenciales para el Medio Poroso 3
Al igual que en el medio poroso 2, las ecuaciones de transferencia no presentan

trminos de flujo debido a que es un medio discontinuo; por lo tanto, nicamente
aparecen el trmino de transferencia entre los medios 2 y 3.
21
Captulo 2
Ecuaciones de Transferencia.
Una para cada fase (agua, aceite y gas):
Aceite:
o 32 =
( bo So ) 3 .
t
(2.32)
Gas:
[ Rso ]32 g
3,2

( Rsbo So + bg S g ) .
3
t
(2.33)
Agua:
[w ]32 =
[bw Sw ]3 .
t
(2.34)
Ecuaciones de Presiones Capilares.

Expresiones de presin capilar para describir la relacin entre las presiones de las
fases:
Pc go ,3 ( S g 3 ) = pg ,3 po ,3 ;
(2.35)
Pc wo ,3 ( S w3 ) = po ,3 pw,3 .
(2.36)
Ecuacin de Restriccin.
Una ecuacin de restriccin para las saturaciones de las fases:
So 3 + S g 3 + S w3 = 1 .
(2.37)
22
Captulo 2
2.3
Condiciones Iniciales y de Frontera
Adems de las ecuaciones diferenciales en la formulacin de un problema de flujo de

fluidos en un medio poroso, se requiere definir las condiciones iniciales y de frontera.
La condicin inicial es aquella que representa la condicin del sistema a un tiempo
especfico, y las condiciones de frontera indican las condiciones que se tienen en los
lmites del sistema.
2.3.1 Condiciones Iniciales
En la mayora de los estudios de simulacin, las condiciones iniciales se obtienen
suponiendo condiciones de equilibrio gravitacional/capilar inicial, lo que implica cero
flujo de las fases a t=0. En esta forma de inicializacin, la distribucin de presiones se
obtiene especificando una presin de referencia a una profundidad determinada y,
utilizando el potencial del fluido, se determinan las presiones a las dems
profundidades. La distribucin inicial de saturaciones se obtiene a partir de los
contactos de los fluidos y las relaciones de presin capilar. El procedimiento descrito
anteriormente se puede expresar matemticamente mediante el algoritmo siguiente:
1. Se calcula la presin de la fase por medio de la siguiente ecuacin:
po = po, ref + o ( po ) z zref .

Donde
(2.38)
po es la presin de la fase aceite a cualquier nivel del yacimiento; po , ref es la
presin de la fase aceite al nivel de referencia fijado, y

fase aceite a una presin
es la densidad relativa de la
po .
Ya que la ecuacin (2.38) contiene implcitamente a po, (a travs de o ), la solucin de

po puede obtenerse en forma iterativa por medio del mtodo de Newton-Raphson.
2. Si pi < pb , existe gas libre en el yacimiento y, por lo tanto,
(2.39)
p p , ref = pb ,
Donde
pb , es la presin de saturacin.
Entonces, la presin del gas en el contacto gas-aceite, se puede calcular por medio de
la expresin siguiente:
23
Captulo 2
pg , zcgo = Pc go , e + po , zcgo .
(2.40)
De forma similar para el agua:

(2.41)
pw, zcwo = po , zcwo Pc wo , e .
3. La distribucin de presiones del gas y del agua se calculan, de forma similar a la del
aceite, por medio de las expresiones siguientes:
p g = p g , zcgo + g p g z zcgo
pw = pw, zcwo + w pw [ z zcwo ]
(2.42)
(2.43)
4. Conocida la distribucin de presiones en el yacimiento (figura 3.9), se calcula la

distribucin de presiones capilares, con las expresiones siguientes:
Pcgo ( S g ) = pg ( z ) po ( z ) ;
(2.44)
Pcwo ( S w ) = po ( z ) pw ( z ) .
(2.45)
5. Por medio de interpolacin inversa y haciendo uso de las presiones capilares, se

obtiene la distribucin inicial de saturaciones.
Pcgo ( z ) = pg ( z ) po ( z )
Pcwo ( z ) = po ( z ) pw ( z )
Figura 2.14 Distribucin de Presiones en el Yacimiento
24
Captulo 2
2.3.2 Condiciones de Frontera

Las ecuaciones diferenciales a resolver por el simulador requieren que se especifiquen
todas las fronteras, esto incluye tanto las fronteras exteriores (lmites del yacimiento,
x=L), como fronteras interiores (pozos productores e inyectores, x=0). Para las
fronteras externas, x=L, generalmente se supone que no existe flujo en la frontera del
yacimiento. Matemticamente esto puede expresarse de la siguiente forma:
( p
p D )
x=L
(2.46)
=0,
donde:
L = Frontera del yacimiento.
p = fases agua, gas y aceite
Para fronteras internas, x=0, pueden especificarse, ya sea gastos del pozo o presiones
de fondo fluyendo. Si se especifica el gasto del pozo, la ley de Darcy puede utilizarse
para generar una condicin de frontera tipo Neumann:
( p
p D )
x=0
qw B
.
kA
(2.47)
Si se especifica la presin para el radio de pozo, entonces se puede obtener una

condicin de frontera tipo Dirichlet,
x =0
(2.48)
= pwf .
25
Captulo 3
Solucin Numrica del Problema

CAPTULO 3
SOLUCIN NUMRICA DEL PROBLEMA
La solucin numrica consiste en obtener una representacin aproximada de las

ecuaciones en derivadas parciales en puntos predeterminados del dominio, en espacio
y tiempo, mediante el empleo de mtodos de discretizacin en diferencias finitas.
El conjunto de ecuaciones que describen el comportamiento de flujo de fluidos en el
yacimiento, ecuaciones (2.20) a (2.22), (2.26) a (2.28) y (2.32) a (2.34), para los
medios 1, 2 y 3, son no lineales, y por consiguiente no se pueden resolver por
mtodos analticos. Se recurre a mtodos numricos para resolverlas mediante su
aproximacin en diferencias finitas, lo que genera
un sistema algebraico de
ecuaciones no lineales, que puede ser resuelto mediante el mtodo iterativo de
Newton-Raphson. En cada iteracin, el mtodo de Newton-Rapshon resulta en un
sistema algebraico de ecuaciones lineales que puede resolverse mediante algn
algoritmo aplicable a la solucin de sistemas de ecuaciones con matrices dispersas (15).
En este captulo se presentan las ecuaciones aproximadas de flujo para el modelo de
triple porosidad y se revisa la aplicacin del mtodo de Newton-Raphson a la solucin
de los sistemas de ecuaciones algebraicas no lineales correspondientes. Tambin se
aprovecha el esquema matricial resultante para reducirlo mediante operaciones
matriciales (27).
3.1.
Ecuaciones de Flujo en Diferencias Finitas.
Los sistemas de ecuaciones no lineales (2.20) a (2.22), (2.26) a (2.28) y (2.32) a

(2.34), para los medios 1, 2 y 3, se resuelven numricamente. El carcter continuo de
estas ecuaciones en espacio y tiempo, se cambia por un carcter discreto mediante
su aproximacin en diferencias finitas. Los trminos de flujo de las ecuaciones del
medio 1, se aproximan mediante diferencias centrales y los trminos de acumulacin
en los medios 1, 2 y 3,
mediante
diferencias regresivas. Este proceso de
discretizacin da como resultado un sistema de ecuaciones algebraicas no lineales en
cada etapa de tiempo.(15). Como se muestra ms adelante, las ecuaciones, que son las
que constituyen el modelo, se escriben en trminos de operadores en diferencias. La
aproximacin en diferencias finitas de los sistemas de ecuaciones (2.20) a (2.22),
(2.26) a (2.28) y (2.32) a (2.34), se desarrollan en detalle en el Apndice B.
3.1.1. Ecuaciones para el Medio 1.
Las ecuaciones de flujo en el medio 1 (establecidas anteriormente), aproximadas
mediante diferencias finitas para el caso de flujo tridimensional son las siguientes :
Para el aceite:
V
n +1
n +1
n +1
To ( po o D ) 1, ijk + bo qo + [o 21 ]ijk = r , ijk t ( bo So ) 1,ijk .
1, ijk
t
26
(3.1)
Captulo 3
Para el gas:
To Rs ( po o D ) + Tg ( pg g D )
+ Rsbo qo
n +1
1, ijk
+ bg q g
n +1
1, ijk
n +1
1, ijk
Vr , ijk
n +1
n +1
+ [ Rso ]21 + g =
t ( Rsbo So + bg S g )
21
1, ijk
t
. (3.2)
Para el agua :
n +1
n +1
Vr , ijk
n +1
Tw ( pw w D ) 1, ijk + bw q w + w 21 =
t [ bw S w ]1, ijk .
1, ijk
1, ijk
t
(3.3)
Las ecuaciones de presin capilar se discretizan de la siguiente manera:
Pc go1,n +ijk1 ( S g ) = pg po ;
1, ijk
n +1
(3.4)
n +1
Pcwo1,n +ijk1 ( S w ) = po pw .
1, ijk
(3.5)
La ecuacin de restriccin restante, es la siguiente:
(S
+ Sg + Sw )
n +1
1, ijk
= 1.
(3.6)
En las ecuaciones (3.1) a (3.3), se tiene que la transmisibilidad de las fases en el

medio 1, se define como:
b k
To = o ro ,
o
(3.7)
bg krg
Tg =

g
(3.8)
k
Tw = rw ,
Bw w
(3.9)
27
Captulo 3
Donde , es un factor geomtrico que depende del sistema de coordenadas

geomtricas elegidas.
Para evaluar las transmisibilidades de los fluidos en el medio 1 (fronteras de las
celdas), se utiliza el concepto corriente arriba; es decir, las propiedades dentro de los
parntesis se evalan a las condiciones de presin, y saturacin de la celda de mayor
presin (28).
Las aproximaciones para la transferencia de fluidos entre el medio 2 y el medio 1,
mostradas en las ecuaciones (3.1), (3.2) y (3.3), son las siguientes (32):
n +1
( o 21 )ijk
= To 21 ( po 2 po1 ) ijk ,
(3.10)
( )
= Tg 21 ( pg 2 pg1 )
(3.11)
n +1
n +1
g 21 ijk
( w21 )ijk
n +1
n +1
ijk
n +1
= Tw 21 ( pw 2 pw1 ) ijk .
Ntese que
p 21 = Vr ,ijkp 21
(3.12)
, p = o, g, w. Las transmisibilidades de las fases entre los
medios 1 y 2 se evalan corriente arriba; esto es, se calcula el potencial entre los
medios y las propiedades se evalan a las condiciones del medio de mayor presin.
Esto es, si el medio 2 es el de mayor presin, entonces las propiedades estn
definidas como:
b k
To 21 = Vr k2 o ro ,
o 21
(3.13)
b k
Tg 21 = Vr k2 g rg

g
,
21
(3.14)
k
Tw 21 = Vr k2 rw .
Bw w 21
(3.15)
De lo contrario, si el medio 1 es el de mayor presin (ej. inyeccin de fluidos en el

yacimiento), entonces las propiedades se evalan a las condiciones de este medio.
28
Captulo 3
Ntese que las ecuaciones (3.4) a (3.9) pueden acoplarse en las ecuaciones (3.1) a
(3.3) para reducir el nmero de ecuaciones y de incgnitas del sistema.
Las presiones de las fases gas y agua, p g y p w respectivamente, pueden eliminarse
usando las relaciones de presin capilar, dadas en las ecuaciones (3.4) y (3.5):
(p
= po + Pcgo )
n +1
1, ijk
(3.16)
( pw = po Pc wo )1,ijk .
n +1
(3.17)
La saturacin de aceite, S o , puede reemplazarse usando la relacin expresada por la

ecuacin (3.6):
(S
= 1 S g Sw )
n +1
1, ijk
(3.18)
La porosidad de cada medio en el nivel de tiempo n+1, se puede escribir en funcin de

la compresibilidad de la roca, cr , considerada como constante :
n +1 = n 1 + cr ( pon +1 pon ) .
(3.19)
Acoplando las expresiones (3.16) a (3.19) en (3.1) a (3.3), el sistema reducido de

ecuaciones es el siguiente:
Para el aceite:
n +1
n +1
n +1
To ( po o D ) 1, ijk + bo qo + To 21 ( po 2 po1 ) ijk

1, ijk
=
Vr , ijk
t
.
t bo (1 S g S w )
1, ijk
(3.20)
Para el gas:
Tg ( po + Pcgo g D )
n +1
1, ijk
n +1
+ Tg 21 ( po 2 po1 + Pcgo 2 Pcgo1 ) + RsTo 21 ( po 2 po1 )

=
Vr , ijk
t
n +1
+ RsTo ( po o D ) 1, ijk + bg qg + Rs bo qo
1, ijk
t bg S g + bo Rs (1 S g S w )
1, ijk
29
n +1
ijk
(3.21)
Captulo 3
Para el agua:
Tw ( po Pcwo w D ) 1, ijk + [bw qw ]1, ijk

n +1
n +1
n +1
+ Tw 21 ( po 2 po1 Pcwo 2 + Pcwo1 ) ik =
3.1.2.
Vr , ik
t
t [ bw S w ]1,ijk .
(3.22)
Ecuaciones para el Medio 2.
Siguiendo un procedimiento similar para los medios 2 y 3, el sistema de ecuaciones

diferenciales que modelan el flujo mutifsico en el medio 2, expresadas por las
ecuaciones (2.26) a (2.28), aproximadas mediante diferencias finitas, son las
siguientes:
Para el aceite:
[ o 21 ]ijk + [ o 32 ]ijk =
n +1
n +1
Vr , ijk
t
t bo (1 S g S w )
(3.23)
2, ijk
Para el gas:
g 21 + Rs o 21
n +1
ijk
+ g 32 + Rs o 32
n +1
ijk
Vr , ijk
t
.
2, ijk
(3.24)
Para el agua:
[ w 21 ]ijk + [ w 21 ]ijk =
n +1
3.1.3.
n +1
Vr , ijk
t
(3.25)
Ecuaciones para el Medio 3.
Para el aceite:
[ o 32 ]ijk =
n +1
Vr , ijk
t
t bo (1 S g S w ) .
3, ijk
30
(3.26)
Captulo 3
Para el gas:
g 32 + Rs o 32
n +1
ijk
Vr , ijk
t
(3.27)
3, ijk
Para el agua:
[ w32 ]ijk =
n +1
Vr , ijk
t
Los trminos de transferencia,
(3.28)
o 21 , o 32 , g 21 , g 32 , w 21
w32 ,
en los medios 1, 2 y 3,
se definen de la forma siguiente:

Aceite:
[ o 21 ]ijk
= To 21 ( po 2 po1 ) ijk ,
[ o32 ]ijk
= To 32 ( po 3 po 2 ) ijk .
n +1
n +1
n +1
n +1
Gas:
g 21 + Rs o 21 = Tg 21 ( po 2 po1 + Pcgo 2 Pcgo1 ) + ( RsTo )21 ( po 2 po1 ) ,

ijk
ijk
n +1
n +1
g 32 + Rs o 32 = Tg 32 ( po 3 po 2 + Pcgo 3 Pcgo 2 ) + ( RsTo )32 ( po 3 po 2 ) .

ijk
ijk
n +1
n +1
Agua:
[ w21 ]ijk
= Tw 21 ( po 2 po1 Pcwo 2 + Pcwo1 ) ijk ,
[ w32 ]ijk
= Tw32 ( po 3 po 2 Pcwo 3 + Pcwo 2 ) ijk .
n +1
n +1
n +1
n +1
31
Captulo 3
Concluyendo, el sistema algebraico de ecuaciones en diferencias finitas para los

medios 1, 2 y 3, consiste en cada nivel de tiempo n+1, y en cada celda ijk de la malla
de clculo, de un conjunto de 9 ecuaciones con el mismo nmero de incgnitas.
3.2 Solucin del Sistema de Ecuaciones
El conjunto de ecuaciones en diferencias finitas, que describe el comportamiento de
flujo de fluidos en el yacimiento, ecuaciones (3.20) a (3.22), (3.23) a (3.25) y (3.26) a
(3.28) constituyen un sistema algebraico de ecuaciones no lineales. Debido a esto, su
solucin se obtiene mediante el mtodo iterativo de Newton-Raphson, lo que genera
en cada iteracin un sistema lineal de ecuaciones.
3.2.1 Mtodo Iterativo de Newton Raphson
La aplicacin del Mtodo de Newton-Raphson comienza con definir las funciones de
residuos siguientes:
Medio 1.
Aceite:
n +1
n +1
n +1
1
Ron1,+ijk
= To ( po o D ) 1, ijk + bo qo + (To )21 ( po 2 po1 )
1, ijk
1, ijk
Vr1,ijk
t
t bo (1 S g S w ) = 0 .
1ijk
(3.29)
Gas
Rgn1,+1ijk = To Rs ( po o D ) 1, ijk + Tg ( po + Pcgo g D ) + bg qg + bo qo Rs

1, ijk
1, ijk
n +1
n +1
n +1
+ (To Rs )21 ( po 2 po1 ) + (Tg ) ( po 2 po1 + Pcgo 2 Pcgo1 )
1, ijk
21
n +1
Vr1, ijk
t
t bo (1 S g S w ) Rs + bg S g = 0 .
1, ijk
(3.30)
Agua:
Rwn1,+1ijk = Tw ( po Pcwo wD ) 1, ijk + [bw qw ]1, ijk

n +1
n +1
n +1
+ Tw 21 ( po 2 po1 Pcwo 2 + Pcwo1 ) 1,ijk
32
Vr1, ijk
t
t [ bw S w ]1, ijk = 0 .
(3.31)
Captulo 3
Medio 2.
Aceite:
V
n +1
n +1
=0
Ron2,+1ijk = (To )21 ( po 2 po1 ) + (To )32 ( po 3 po 2 ) r 2, ijk t bo (1 S g S w )
ijk
ijk
2, ijk
t
. (3.32)
Gas:
Rgn2,+1ijk = (To Rs )21 ( po 2 po1 ) + (Tg ) ( po 2 po1 + Pcgo 2 Pcgo1 )

21
ijk
n +1
+ (To Rs )32 ( po 3 po 2 ) + (Tg ) ( po 3 po 2 + Pcgo 3 Pcgo 2 )

32
ijk
n +1
Vr 2, ijk
t
t bo (1 S g S w ) Rs + bg S g
=0.
2, ijk
(3.33)
Agua:
n +1
n +1
Rwn 2,+1ijk = Tw 21 ( po 2 po1 Pcwo 2 + Pcwo1 ) ijk + Tw32 ( po 3 po 2 Pcwo 3 + Pcwo 2 ) ijk
Vr 2, ijk
t
t [ bw S w ]2, ijk = 0 .
(3.34)
Medio 3.
Aceite:
V
n +1
Ron3,+1ijk = (To )32 ( po 3 po 2 ) r 3, ijk t bo (1 S g S w ) = 0 .
ijk
3, ijk
t
(3.35)
Gas:
Rgn3,+1ijk = (To Rs )32 ( po 3 po 2 ) + (Tg ) ( po 3 po 2 + Pcgo 3 Pcgo 2 )

32
ijk
n +1
Vr 3, ijk
t
t bo (1 S g S w ) Rs + bg S g = 0 .
3, ijk
33
(3.36)
Captulo 3
Agua:
V
n +1
Rwn 3,+1ijk = Tw32 ( po 3 po 2 Pcwo 3 + Pcwo 2 ) ijk r 3, ijk t [ bw S w ]3, ijk = 0 .
t
(3.37)
Una vez expresadas las funciones de residuos, se establece la dependencia de las

mismas, para esto, se definen los siguientes vectores de incgnitas, en la celda ijk y
para cada uno de los medios:
U1, ijk = ( po , S g , S w )
1, ijk
U 2, ijk = ( po , S g , S w )
U 3, ijk = ( po , S g , S w )
2, ijk
T
3, ijk
,
,
.
En forma general, la dependencia de las funciones de residuos de las incgnitas de

cada medio, sern:
Medio 1:
Ro1,ijk = Ro1, ijk ( U1, ijk 1 , U1, ij 1k , U1, i 1 jk , U1, ijk , U 2, ijk , U1, i +1 jk , U1, ij +1k , U1, ijk +1 ) = 0 ,
Rg1, ijk = Rg1, ijk ( U1, ijk 1 , U1, ij 1k , U1, i 1 jk , U1, ijk , U 2, ijk , U1, i +1 jk , U1, ij +1k , U1, ijk +1 ) = 0 ,
Rw1, ijk = Rw1, ijk ( U1, ijk 1 , U1, ij 1k , U1, i 1 jk , U1,ijk , U 2,ijk , U1, i +1 jk , U1, ij +1k , U1, ijk +1 ) = 0 .
Medio 2:
Ro 2, ijk = Ro 2, ijk ( U1, ijk , U 2, ijk , U 3, ijk ) = 0 ,

Rg 2, ijk = Rg 2, ijk ( U1, ijk , U 2, ijk , U 3, ijk ) = 0 ,
Rw 2,ijk = Rw 2, ijk ( U1, ijk , U 2, ijk , U 3, ijk ) = 0 .
Medio 3:
Ro3,ijk = Ro3, ijk ( U 2, ijk , U 3, ijk ) = 0 ,
34
Captulo 3
Rg 3,ijk = Rg 3, ijk ( U 2, ijk , U 3, ijk ) = 0 ,

Rw3, ijk = Rw3, ijk ( U 2, ijk , U 3, ijk ) = 0 .
El proceso iterativo se establece expandiendo las funciones de residuo mediante una
serie de Taylor truncada, alrededor del nivel de iteracin ( ), de la que solo se
conservan los trminos de menor orden (29), esto es:
Medio 1.
Fase aceite o gas:
Rp,ijk +1 Rp,ijk +1 Rp,ijk +1 Rp,ijk +1 Rp,ijk +1

po,klm +
Sg ,klm +
Sw,klm +
po2,ijk +
Sg 2,ijk
S
klm
o
klm
g
klm
w
klm
o
ijk
g
ijk
,
,
,
2,
2,
Rp,ijk +1
Sw2,ijk = Rp1,ijk .
+
Sw2,ijk
(3.38)
Para agua:
klm
Rp,ijk +1 Rp,ijk +1 Rp,ijk +1 Rp,ijk +1

po,klm +
Sw,klm +
po2,ijk +
Sw2,ijk = Rp1,ijk ,
Sw,klm
Sw2,ijk
po,klm
po2,ijk
(3.39)
donde:
p=o, g, w.
klm = ijk 1, ij 1k , i 1 jk , ijk , i + 1 jk , ij + 1k , ijk + 1

Medio 2.
Fase aceite o gas:
Rp,ijk +1 Rp,ijk +1 Rp,ijk +1 Rp,ijk +1 Rp,ijk +1 Rp,ijk +1

po1,ijk +
Sg1,ijk +
Sw1,ijk +
po2,ijk +
Sg 2,ijk +
Sw2,ijk
po1,ijk
Sg1,ijk
Sw1,ijk
po2,ijk
Sg 2,ijk
Sw2,ijk
Rp,ijk +1 Rp,ijk +1 Rp,ijk +1
+
po3,ijk +
Sg3,ijk +
Sw3,ijk = Rp2,ijk
po3,ijk
Sg3,ijk
Sw3,ijk
. (3.40)
35
Captulo 3
Para agua:

po1,ijk +
Sw1,ijk +
po2,ijk +
Sw2,ijk +
po3,ijk +
Sw3,ijk = Rw2,ijk
po1,ijk
Sw1,ijk
po2,ijk
Sw2,ijk
po3,ijk
Sw3,ijk
. (3.41)
Donde:
p=o, g, w.
Medio 3.
Fase aceite o gas:

po2,ijk +
Sg 2,ijk +
Sw2,ijk +
po3,ijk +
Sg3,ijk +
Sw3,ijk = Rp3,ijk
po2,ijk
Sg 2,ijk
Sw2,ijk
po3,ijk
Sg3,ijk
Sw3,ijk
.(3.42)
Para agua:
Rp,ijk +1 Rp,ijk +1 Rp,ijk +1 Rp,ijk +1

po2,ijk +
Sw2,ijk +
po3,ijk +
Sw3,ijk = Rp3,ijk .
po2,ijk
Sw2,ijk
po3,ijk
Sw3,ijk
(3.43)
En las ecuaciones (3.38) a (3.43), se tiene un sistema de ecuaciones donde las

incgnitas son los cambios iterativos de las presiones y saturaciones en cada una de
las celdas y en cada uno de los medios; esto es:
y
( p , S
o
, Sw )
+1
definiciones de
3, ijk
( p , S , S )
+1
w 1, ijk
, donde i = 1, 2, ..., I , j = 1, 2, ..., J ,
( p , S
o
, Sw )
+1
2, ijk
k = 1, 2, ... , K. Las
po , S g , S w , para el medio poroso 1 son las siguientes:
1
1
po1,+ijk
= po1,+ijk
po1, ijk ,
(3.44)
Sg1,+1ijk = Sg1,+1ijk Sg1,ijk ,
(3.45)
36
Captulo 3

(3.46)
Sw1,+1ijk = Sw1,+1ijk Sw1,ijk .

Las definiciones para los medios 2 y 3, son similares.
Ntese que con el fin de simplificar la escritura, el superndice correspondiente al nivel

de tiempo n+1 fue eliminado de las expresiones (3.38) a (3.43).
El proceso iterativo del nivel de tiempo n+1 se inicia comnmente con la estimacin de
la solucin siguiente:
Para el medio 1:
( p ,S
o
, Sw )
1, ijk
( p ,S
o
, Sw )
1, ijk
(3.47)
Para el medio 2:
( p ,S
o
, Sw )
0
2, ijk
= ( po , S g , S w )
n
2, ijk
(3.48)
Para el medio 3:
( p ,S
o
, Sw )
0
3, ijk
= ( po , S g , S w )
n
3, ijk
(3.49)
Y termina cuando los cambios iterativos de las incgnitas en cada uno de los medios
son menores a una tolerancia estipulada,
+1
, por ejemplo: p o ,ijk
< . Las incgnitas
restantes se definen similarmente.

El sistema lineal de ecuaciones obtenido mediante la aplicacin del mtodo iterativo de
Newton-Raphson en la iteracin
( + 1) , expresiones (3.38) a (3.43), puede escribirse
en forma compacta como sigue :
[J ]
(3.50)
U +1 = R ,
37
Captulo 3
Donde
[J ]

se conoce como la matriz Jacobiana,
es el vector de incgnitas y
R es el vector de funciones de residuos.

La estructura matricial del subsistema de ecuaciones para la celda ijk, es el siguiente:
38
Captulo 3
U1
U
+
2
U3 i 1, jk
i 1, jk
+1
G11
[ 0]
[ 0]
+1
A11 A12
A21 A22
[ 0] A32
E11
[ 0]
[ 0]
[0]
[0]
[0]
[0]
[0]
[0]
C11
[ 0]
[ 0]
[0]
[0] U1
[0] U2 +
[0]ij,k 1 U3 ij,k 1
B11
[ 0]
[ 0]
D11
[ 0]
[ 0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
+1
U1
U
+
2
U3 i, j 1, k
i , j 1, k
[0]
A23
A33 ijk
+1
U1
U
+
2
U3 i, j , k
F11 [ 0] [ 0]
U1
[0] [0] U1
[0] [0] U2 + [0] [0] [0] U2 +
[0] [0]ij,k +1 U3 ij,k +1 [0] [0] [0]i, j +1,k U3 i, j +1,k
+1
R1
[0] [0] U1
[0] [0] U2 = R2
[0] [0]i +1, jk U3 i +1, jk R3 ijk
+1
2539
+1
(3.51)
Captulo 3
Donde E11 , G11 , C11 , A11 , A12 , B11 , F11 , D11 , A22 , A23 , A33 , A32 son submatrices de orden (3)(3).
Las submatrices
A11 y A12 , contienen las derivadas de las funciones de residuos del
medio 1 con respecto a las incgnitas del medio 1 y del medio 2 en la celda ijk,
respectivamente. Las submatrices
A21 , A22 y A23 , contienen las derivadas de las
funciones de residuos del medio 2 con respecto a las incgnitas del medio 1, del medio
2 y del medio 3 en la celda ijk, respectivamente. Las
submatrices A32 y A33 ,
contienen las derivadas de las funciones de residuos del medio 3 con respecto a las
incgnitas del medio 2 y del medio 3 en la celda ijk, respectivamente. Las submatrices
E11 , G11 , C11 , B11 , F11 , y D11 contienen las derivadas de las funciones de residuos
del medio 1 con respecto a las incgnitas del medio 1 en las celdas i-1,j,k; i,j-1,k;
i,j,k-1 ; i,j,k+1 ; i,j+1,k e i+1,j,k respectivamente.
3.2.2 Reduccin Matricial del Sistema de Ecuaciones
Primera Reduccin para la formulacin TI.
La estructura matricial del sistema lineal de ecuaciones (3.51), permite reducirlo a un
sistema de 3 ecuaciones con 3 incgnitas por celda, similar al que se obtiene en el
caso de yacimientos homogneos de aceite negro. La reduccin del sistema de
ecuaciones se realiza mediante la aplicacin del complemento de Schur.(40)(21), de la
forma siguiente: en una primera etapa se acoplan las ecuaciones del medio 3 en las
ecuaciones del medio 2 y posteriormente se acoplan las del medio 2 en las del medio
1. Lo anterior se muestra paso a paso a continuacin, para esto, el subsistema de
ecuaciones 3.51, puede escribirse de la siguiente manera:
Ecuacin del medio 1:
C11
U1 +1 + G11 U1 +1 + E11 U1 +1 + A11 U1 +1 + A12
U 2 +1 + D11 U1 +1 + F11 U1 +1 + B11 U1 +1 =
R1
. (3.52)
(3.53)
A21 U1 +1 + A22 U 2 +1 + A23 U 3 +1 = R2 .

(3.54)
A32
U 2 +1 + A33
U 3 +1 = R3 .
40
Captulo 3
Ahora bien, la expresin (3.54) se resuelve para las incgnitas del medio 3,
trminos de las incgnitas del medio 2,
U 2 +1
U 3 +1 , en
, como sigue :
(3.55)
U 3 +1 = R3* A32* U 2 +1 .
Donde:
1
A32
A32 = A33
.
(3.56)
(3.57)
R = A33 R3 .
La expresin (3.55), representa la solucin de las incgnitas en el medio 3, en funcin
de las incgnitas del medio 2.
Sustituyendo la ecuacin (3.55), en la ecuacin del medio 2, expresin (3.53), se tiene
lo siguiente:
*
A21 U1 +1 + A22 U 2 +1 + A23 R3* A32
U 2 +1 = R2 .
(3.58)
Realizando las operaciones correspondientes:

*
A21 U1 +1 + A22 U 2 +1 A23 R3* A23 A32
U 2 +1 = R2 .
Agrupando:
*
= R2 + A23 R3* .
A21 U1 +1 + U 2 +1 A22 A23 A32
Finalmente:
(3.59)
*
A21 U1 +1 + A22
U 2 +1 = R2* ,
Siendo:
41
Captulo 3

(3.60)
22
A = A22 A23 A
*
32 ,
(3.61)
R2 = R2 A23 R3* .
Ahora bien, de la expresin (3.59), se resuelve para las incgnitas del medio 2,
U 2 +1 , en trminos de las incgnitas del medio, U1 +1
, como sigue :
(3.62)
*
U 2 +1 = R2** A21
U1 +1 ,
Donde:
(3.63)
A21 ,
= A22
A21
(3.64)
*
R2* .
R2** = A22
La expresin (3.62), representa la solucin de las incgnitas del medio 2, en funcin

de las incgnitas del medio 1.
Sustituyendo la ecuacin (3.62) en la ecuacin del primer medio, expresin (3.52), se
tiene lo siguiente:
C11
U1 +1 + G11 U1 +1 + E11 U1 +1 + A11 U1 +1 + A12 R2** A21* U1 +1 +
D11
U1 +1 + F11 U1 +1 + B11 U1 +1 = R1 .
(3.65)
Simplificando:
C11
U1 +1 + G11 U1 +1 + E11 U1 +1 + A11 U1 +1 A12 R2** A12 A21* U1 +1 +
D11
U1 +1 + F11 U1 +1 + B11 U1 +1 = R1 .
Agrupando:
C11
U1 +1 + G11 U1 +1 + E11 U1 +1 + U1 +1 A11 A12 A21* +
D11
U1 +1 + F11 U1 +1 + B11 U1 +1 = R1 + A12 R2** .
42
(3.66)
Captulo 3
Finalmente:
C11
U1 +1 + G11 U1 +1 + E11 U1 +1 + A11* U1 +1 + D11 U1 +1 + F11 U1 +1 + B11 U1 +1 = R1* .
(3.67)
Donde:
(3.68)
*
A11* = A11
A12
A21 ,
(3.69)
R1* = R1 A12
R2** .
Ntese que el subsistema de ecuaciones en (3.67), consiste de 3 ecuaciones con 3

incgnitas por cada celda de la malla de clculo. El subsistema (3.67) de ecuaciones
contiene slo incgnitas del medio 1, o sea,
U1 +1 .
La estructura matricial del subsistema de ecuaciones reducido (3.67), posee una

estructura matricial similar a la que se obtiene para un yacimiento de aceite negro de
una porosidad.
Una vez que se resuelven las incgnitas en el medio 1, en el nivel de tiempo n+1, es
posible resolver las incgnitas del medio 2 (ecuacin 3.62) y del medio 3 (ecuacin
3.55), de manera directa en cada nodo de la malla de clculo (5).
3.2.3 Solucin del Sistema Lineal de Ecuaciones
La solucin del sistema lineal se puede obtener mediante mtodos directos o mtodos
iterativos.
Mtodos Directos
Los mtodos directos realizan un nmero finito de operaciones para resolver un
sistema de ecuaciones.
En los mtodos directos, el trabajo computacional y el requerimiento de memoria estn
directamente relacionados con el nmero de incgnitas a resolver. Entre mayor es el
nmero de ecuaciones, el ancho de banda es generalmente mayor y la cantidad de
memoria de cmputo y de tiempo de procesamiento son tambin mayores.
El nmero de ecuaciones y de incgnitas, para un problema dado, depende del nmero
total de bloques o de celdas, en que se discretiza el dominio de inters y del nmero
de componentes considerados.
Los mtodos directos se basan en el Mtodo de Eliminacin Gaussiana, y constan de
dos pasos:
43
Captulo 3
1. En un barrido hacia delante se eliminan todos los elementos diferentes de cero

ubicados debajo de la diagonal principal, esto es, se transforma la matriz de
coeficientes en una matriz triangular.
2. En un barrido hacia atrs, se obtiene la solucin de las incgnitas.
Dentro de los mtodos directos, los Algoritmos de Banda y las Tcnicas para resolver
sistemas de ecuaciones con Matrices Dispersas, son de inters en la Simulacin
Numrica de Yacimientos.
Las tcnicas de matrices dispersas reconocen la existencia de un gran nmero de ceros
en la matriz de coeficientes y evitan su almacenamiento y operaciones innecesarias.
Dentro de este tipo de mtodos se encuentran: El Mtodo de Eliminacin Gaussiana y
el Algoritmo de Thomas.
Mtodos Iterativos
En los mtodos iterativos se aplican algoritmos cclicos; se busca que a medida que las
iteraciones progresen, el mtodo converja a la solucin.
En el presente trabajo, el sistema lineal de ecuaciones se resuelve mediante el mtodo
iterativo GMRES, el cual es un mtodo tipo Krylov y se describe a continuacin:
Mtodos de Krylov
Estos mtodos son utilizados comnmente en combinacin con el mtodo de
Descomposicin de Dominio para resolver problemas cuya dimensin requiere del uso
de cmputo paralelo, as como en la solucin de problemas no lineales.
Un algoritmo iterativo resulta una opcin adecuada para la solucin de sistemas de
gran escala, sin embargo las tcnicas usuales como los mtodos de tipo SOR, resultan
poco efectivos en estos casos, entre los algoritmos iterativos ms eficientes se
encuentran los mtodos de Krylov, los cuales han sido objetos de estudio en diversas
reas de la matemtica aplicada.
Estos algoritmos se basan en la idea de buscar una aproximacin a la solucin del
sistema dentro de un subespacio de dimensin pequea, esto es descrito a
continuacin.
Sea
r0 = b Ax0
donde
A R nn
con
y
x0 una aproximacin inicial a x* , solucin del sistema Ax = b
b R n , se desea determinar
xk x0 + K k
(4.57)
b Axk Lk
(4.58)
44
Captulo 3
tal que
xk sea una aproximacin razonable de x * y k << n , con K k := { A j r0 }
k 1
j =0
Lk un subespacio dependiente del sistema.

La eleccin del subespacio Lk conduce a diferentes mtodos de aproximacin.
Gradiente Conjugado
Uno de los mtodos mas utilizados debido a su simplicidad, eficiencia y bajo costo
computacional, es el mtodo de Gradiente Conjugado (CG), en este caso se considera
Lk = K k , con esta eleccin el proceso lleva a buscar, en cada iteracin, reducir el
tamao del error medido en una norma dependiente de la matriz
iteracin de CG se resuelve el problema:
min x x*
xx0 + K k
donde
A , as, en cada
v A :=
con
v, Av
el producto escalar usual.
Desafortunadamente este algoritmo requiere que la matriz de coeficientes sea
A = At ) y la convergencia solo es garantizada si A es definida positiva, es

t
n
decir, xAx > 0 x R , esto no se presenta en el problema que deseamos resolver,
simtrica (
sin embargo existen diversas tcnicas y que permiten el uso de ste mtodo en
problemas mas generales, adems de diversas variantes del mtodo para matrices no
simtricas.
El siguiente algoritmo representa la forma ms simple del Gradiente Conjugado.
Algoritmo GC
Dados
A, b, x0
r0 = b Ax0 , v0 = r0
Para
k = 1, 2,...
Determinar
Si
( b Ax
Construir
xk
k 2
< tol )
vk K k +1
Salir
tal que
vi , v j = 0
45
para
j = 0,..., k 1
Captulo 3
GMRES
Dadas las caractersticas del sistema, es necesario buscar un algoritmo con requisitos
mnimos sobre la matriz de coeficientes, el mtodo de Residuo Mnimo Generalizado
(GMRES), es uno de los algoritmos ms representativos de los mtodos de Krylov y
uno de los ms eficientes para matrices generales, en este proceso se elige
Lk = AK k .
En cada aproximacin de GMRES, al igual que en gradiente conjugado, se resuelve un

problema de minimizacin, sin embargo, en este caso se busca reducir la norma del
residual asociado a cada aproximacin, as, en cada paso, GMRES resuelve el
problema:
min b Ax
xx0 + K k
donde
v 2 :=
v, v
Es importante notar que si n = k entonces
xk = x* es decir, la solucin exacta ser
alcanzada en a lo mas n iteraciones, en aritmtica exacta.

Para llevar al cabo la solucin del problema, se requiere de una base ortogonal de
es decir, se requiere
{v }
j
k 1
j =0
{ }
tal que K k := span v j
k 1
j =0
Kk ,
vi , v j = 0 para i j por lo
que en cada paso debe obtenerse un elemento de dicha base y almacenar los
anteriores.
Algoritmo GMRES
Dados
A, b, x0
r0 = b Ax0 , v0 = r0
Para
k = 1,...
Determinar
Si
( b Ax
Construir
xk
k 2
< tol )
Salir
vk K k +1 tal que vi , v j = 0
para
j = 0,..., k 1
Debido a la necesidad de almacenar los vectores generados en cada iteracin, el

consumo de memoria al aplicar GMRES puede ser muy elevado.
Una estrategia para controlar el consumo de memoria, es el reinicio del algoritmo si
ste no ha alcanzado una aproximacin aceptable despus de un nmero de
46
Captulo 3
iteraciones previamente establecido ( m ), el algoritmo as obtenido es conocido como

GMRES con reinicio (GMRES( m )) .
Algoritmo GMRES( m )
Dados
A, b, x0
r0 = b Ax0 , v0 = r0
Para
k = 1,..., m 1
Determinar xk
Si
( b Ax
Construir
Si
k 2
< tol )
Salir
vk K k +1 tal que vi , v j = 0
( b Ax
k 2
> tol )
para
j = 0,..., k 1
entonces
x0 = x k
r0 = b Ax 0
v 0 = r0
reiniciar el ciclo
Esta versin del algoritmo requiere nicamente del almacenamiento de m vectores
ortogonales, sin embargo la eleccin del parmetro de reinicio no es sencilla, en
general, no existen criterios que permitan elegir el valor de m a partir de las
caractersticas del sistema.
3.3 Implementacin del Modelo Conceptual en el Simulador
El simulador en el cual fue implementado el simulador es un simulador multifsico y
tridimensional capaz de modelar el flujo de fluidos en yacimientos naturalmente
fracturados. Este simulador fue desarrollado previamente y la contribucin del
presente trabajo fue implementar inicialmente el modelo de triple porosidad, para
posteriormente expandirlo a mltiple porosidad.
A continuacin se presenta el algoritmo de programacin del simulador, en el cual se
incluyen las principales tareas realizadas en el cdigo. Para este simulador existen
varios mdulos los cuales realizan tareas especficas entre las que se encuentran:
lectura de datos, creacin de la malla, clculos PVT, generacin del jacobiano, solucin
del sistema de ecuaciones e impresin de resultados [Pea, O y Resendiz, J.T.,
2006 24].El diagrama de flujo del simulador a grandes rasgos es el siguiente:
47
Captulo 3
INICIO
Inicializacin
de variables
Lectura
de datos
Verificacin de
Convergencia
Tsim<=Tempo
no
si
Petrofsica
TolVar = NO_CONV .or.

MaxVar = NO_CONV
Alojamiento
de memoria
Cambio de
unidades
AllTolVars =
NO_CONV
Ciclo while
de tiempo
Inicializacin
con ceros
no
Aux1 = Tempo
+ DT
MaxVar =
NO_CONV
no
si
AllTolVars =
SI_CONV
no
Aux1>=Tsim
si
AllTolVars =
NO_CONV
DT =tsim
- tempo
Niter = niter+1
var_prim(+1)=var_prim()
si
Coord.
Cilndricas
NTS=NTS+1
si
Niter > ItMax
Niter = 0
no
Tipo == K
DT=FMDT*DT
si
NTS=0
Tipo == C
no
NITER = 0
no
Niter == 1
no
Niter = 0
Petrofsica
si
si
Coord.
Cartesianas
var(+1)=var()
dt = dt / FMDT
.
Formacin del
Jacobiano
Inicializacin
Po, Sg, Sw
Llamado del
SOLVER
Clculo vol.
Inicial poroso
no
si
AllTolvars =
NO_CONV
var_prim(+1)=var_prim()
var_prim(n+1)=var_prim(n)
Impresin
en pantalla
TolVar =
SI_CONV
AllTolVars=SI_CONV
and.NITER>=itmin
Figura 3.1 Diagrama de flujo del Simulador Camalen
48
40
(11)
FIN
EXIT
Captulo 4
Validacin del Modelo

CAPTULO 4
VALIDACIN DEL MDELO
La validacin del modelo propuesto se realiz mediante la simulacin de flujo radial

monofsico en uno, dos y tres medios porosos. Posteriormente se graficaron los
resultados de Pwf vs log (t), para observar su comportamiento.
Los resultados de la simulacin del flujo radial en uno, dos y tres medios, se validaron
mediante soluciones analticas de la siguiente forma: la validacin para el flujo en un
solo medio se realiz mediante el uso de la solucin de lnea fuente, para dos medios
mediante el uso de la solucin de Warren y Root (35), y para tres medios mediante la
solucin propuesta por Rodriguez y cols (27), para un modelo de triple porosidad-una
permeabilidad. A continuacin se describen las soluciones analticas utilizadas y
posteriormente se muestran los resultados de la simulacin:
4.1 Solucin para el Flujo Radial Monofsico en un Medio
La ecuacin que gobierna la transmisin de presin en un medio poroso conteniendo
un fluido ligeramente compresible, est dado por:
2 p 1 p 2 p 2 p c p
.
+
+
+
=
k r t
r 2 r r 2 z 2
(4.1)
La ecuacin (4.1), es obtenida bajo las siguientes suposiciones:

a. Flujo laminar o Darcy.
b. Difusividad hidrulica constante.
c. Medio poroso homogneo e isotrpico.
d. Compresibilidad del fluido pequea.
e. Gradientes de presin ( p r ) pequeos, de tal forma que
f.
( p
Flujo monofsico.
g. Efectos de gravedad y trmicos no considerados.

Si se considera nicamente flujo radial, la ecuacin (4.1) se reduce a:
49
r ) 0 .
2
Captulo 4
2 p 1 p c p
.
+
=
k r t
r 2 r r
(4.2)
Esta ecuacin es conocida como la ecuacin de Difusividad, la cual es la base para la

teora de pruebas de presin.
Solucin para Yacimiento Infinito.
Despus del periodo de almacenamiento, el transiente de presin que se podra medir
en la cara del pozo refleja la transmisin de la presin a travs del yacimiento. A
tiempos largos, la respuesta est bajo la influencia de las fronteras del yacimiento,
pero antes de estos efectos, la respuesta de presin no se entera de las fronteras y el
yacimiento acta como si fuera infinito en extensin. Aunque existen varios tipos de
flujos identificados durante este periodo de yacimiento infinito, uno de los ms
fcilmente identificados es el flujo radial. El flujo radial para yacimientos
comportndose infinitamente es la base de la interpretacin de las tcnicas de pruebas
de presin.
En ausencia de almacenamiento y dao, el transiente de presin (respuesta de presin
que se podra medir en la cara del pozo) bajo un comportamiento infinito de un flujo
radial, considerando al pozo como una lnea fuente, source line y gasto constante,
est dada por:
1 r2
p D = Ei D ,
2 4t D
(4.3)
donde:
pD =
2kh
( pi p(r ,t ) )
qB
(4.4)
tD =
kt
crw2
(4.5)
rD =
r
rw
(4.6)
50
Captulo 4
Sustituyendo ecuaciones (4.4) a (4.6) en la ecuacin (4.3), se tiene lo siguiente:
p( r , t ) = pi
q B r 2
Ei
,
4 kh 4 t
(4.7)
donde:
es la Difusividad hidrulica
p(r ,t ) es la presin en un punto r al tiempo t.

La integral exponencial, Ei, en la ecuacin (4.7), puede ser aproximada como:
Ei ( x ) = ln x + 0.5772 cuando x < 0.01
donde x =
(4.8)
rD2
r 2c
r2
r2
=
=
=
4t D 4t 4kt
4kt
c
Cuando la ecuacin 4.3 es evaluada en la cara del pozo (disparos), o sea, cuando
r=rw, se tiene que:
rD =
rw
1
=1 y x =
4t D
rw
Haciendo
ln x = ln
1
= ln (1) ln ( 4t D ) = 0 ln ( 4 ) + ln ( t D ) ,
4t D
entonces la ecuacin 4.3 se reduce a:
pwD =
1
[ ln tD + 0.80907] + s ,
2
donde s es el dao de la formacin.
51
Captulo 4
4.2 Solucin para Flujo Radial Monofsico en Dos Medios

Los efectos de doble porosidad estn descritos en trminos de dos parmetros
adimensionales, y , que relacionan las propiedades de la matriz con las
propiedades de la fractura. depende de la relacin de transmisibilidades matrizfractura, mientras que relaciona la capacidad de almacenamiento de la fractura con
la capacidad de almacenamiento del sistema completo. y estn definidas por las
ecuaciones siguientes:
km 2
rw ,
kf
f ctf
f ctf + m ctm
(4.9)
(4.10)
es un factor de forma que depende de la geometra de flujo entre la matriz y las

fracturas.
Las expresiones en forma adimensional que describen el flujo en medio de doble
porosidad obtenidas por Warren y Root (35), fueron:
Para el medio continuo:
2 PfD 1 PfD
PfD
PmD
+
=
+
1
(
)
rD2 rD rD
tD
tD
(4.11)
Para el medio discontinuo:
(1 )
PmD
= ( PfD PmD ) .
tD
(4.12)
La solucin de estas ecuaciones para un yacimiento infinito y uno finito, con un pozo
productor a gasto constante [Prats 41], son:
52
Captulo 4
Yacimiento Infinito:
PwD =
t D
t D
1
Ln t D + 0.80908 + Ei
Ei
+ S .
2
1
1
(
)
(
)
(4.13)
Yacimiento finito:
PwD
tD
2
3r 4 4r 4 Ln r 2r 2 1
1)
2 1
(
(1 )
D
D
D
D
= 2
+
t
+
1
.
D
2
2
4
1
rD 1 4
(D )
(4.14)
Donde rD y tD son el radio y tiempo adimensionales, definidos por:
rD =
tD =
r
,
rw
(4.15)
0.0002637 k f t
( ct ) + ( ct ) rw2
f
m
(4.16)
Donde S de representa el factor de dao y el termino Ei(x) representa la funcin

Integral Exponencial. La presin real puede ser obtenida de la definicin de presin
adimensional
PD =
k f h Pi P ( r , t )
141.2qB
(4.17)
Una grfica semi-logartmica de la presin vs. tiempo revela la presencia de dos lneas
rectas paralelas en lugar de una, como lo obtenido generalmente para el caso de un
yacimiento convencional, Figura 5.1. Dado que el sistema de fracturas es de alta
transmisibilidad y esta directamente conectado al pozo, este responde primero; la
matriz no esta conectada de forma directa al pozo y es de baja transmisibilidad, por lo
tanto responde despus. La separacin entre las dos rectas depende de la capacidad
de almacenamiento relativa, , mientras que el periodo de transicin de la primera a
la segunda lnea ser una funcin de la relacin de transmisibilidades, .
53
Captulo 4
4.3 Solucin Flujo Radial Monofsico en Tres Medios

Las ecuaciones que describen el flujo de un fluido ligeramente compresible, con todas
las suposiciones usuales, en un medio de triple porosidad y una permeabilidad, de
acuerdo a Rodrguez y cols (27), se describen a continuacin:
El flujo radial en el medio con porosidad secundaria de gran escala, hacia un pozo
produciendo a gasto constante en un yacimiento infinito, se describe por:
1 p DF
rD
rD rD
rD
p
+ fF ( p Df p DF ) = F DF .
t D
(4.18)
Las condiciones iniciales y de frontera estn dadas por:

(4.19)
p DF (rD ,0) = 0 ,
rD
p DF
rD
= 1 ,
(4.20)
rD =1
y,
(4.21)
Lim p (r , t ) = 0 .
rD DF D D
La ecuacin de flujo para el medio con porosidad secundaria de pequea escala es:
fF ( p Df p DF ) + mf ( p Dm p Df ) = f
p Df
t D
(4.22)
con la siguiente condicin inicial:
p Df (rD ,0) = 0 .
(4.23)
54
Captulo 4
La ecuacin de flujo para la matriz es:
mf ( p Dm p Df ) = m
p Dm
.
t D
(4.24)
con la siguiente condicin inicial:
p Dm (rD ,0) = 0 .
(4.25)
Las variables adimensionales estn definidas por:
p DF = 2k F h( pi p F ) / qB ,
(4.26)
p Df = 2k F h( pi p f ) / qB ,
(4.27)
p Dm = 2k F h( pi p m ) / qB ,
(4.28)
t D = k F t / (ct )rw2 ,
(4.29)
rD = r rw ,
(4.30)
mf = k m mf rw2 k F
(4.31)
fF = k f fF rw2 k F
(4.32)
F = (ct ) F
(ct ) ,
(4.33)
F = (ct ) F
(ct ) ,
(4.34)
55
Captulo 4
m = 1 (ct ) F (ct ) f
(4.35)
(ct ) = (ct ) F + (ct ) f + (ct ) m .
(4.36)
La solucin matemtica al problema descrito por las ecuaciones (4.18) a (4.25), es:
pwDF (t D ) =
1
1 t 1 t 1 t 1 t
Ln ( t D ) + 0.80908 E1 D + E1 D E1 D + E1 D ,
2
2 a1 2 a2 2 a3 2 a4
(4.37)
a j , j=1,2,,4, depende del almacenamiento y de los parmetros de interflujo
donde
del medio de triple porosidad.

A partir de la ecuacin (4.37), se puede obtener una expresin para describir la
derivada de la presin (27), esto es:
tD (
dpwD
) = 1 (1 +
2
dt D
tD
a1
tD
a2
tD
a3
Donde:
a1 =
1
c1 c12 4c3 ,
2
a2 =
1
c2 c2 2 4c4
2
a3 =
1
c1 + c12 4c3 ,
2
a4 =
1
c2 + c2 2 4c4
2
y,
56
tD
a4
).
(4.38)
Captulo 4
c1 =
+ F ) m
mf
c2 =
m ( m + f )
+
,
mf
fF
c3 =
m f F
mf fF
c4 =
m f
mf fF
+ f ) F
fF
4.4 Resultados de la Simulacin de Flujo Monofsico en Uno, Dos y Tres

Medios
La tabla 4.1 muestra los datos utilizados en la validacin del modelo. Cabe mencionar
que para poder realizar la validacin, el valor de la porosidad se mantuvo constante en
ambos casos (analtica y numrica), ya que las soluciones analticas utilizadas hacen
esta suposicin.
La figura 4.1 muestra los resultados obtenidos para la simulacin de flujo monofsico
en uno, dos y tres medios, graficados con los resultados de las soluciones analticas
descritas anteriormente. En esta figura puede observarse que ambos resultados
(analticos vs numricos), se ajustan perfectamente.
Como puede observarse en esa misma figura, cuando se tiene flujo en un solo medio,
se presenta un comportamiento lineal; cuando se trata de dos medios, se presenta un
cambio de pendiente y; cuando se trata de tres medios, dos cambios de pendiente. Lo
anterior se debe a que el medio uno es de alta transmisibilidad y esta directamente
conectado al pozo, por lo tanto, este responde primero; los medios 2 y 3 no estn
conectados de forma directa al pozo y son de baja transmisibilidad, por lo tanto
responden despus.
En la figura 4.1 tambin puede observarse que al inicio de la simulacin el
comportamiento de los resultados de los tres casos (uno, dos y tres medios), es el
mismo; esto es porque al inicio de la simulacin no se sienten los efectos del segundo
y tercer medio, estos se van sintiendo conforme avanza el transiente. Finalmente, en
los resultados obtenidos mediante la solucin numrica, al final de la simulacin
pueden apreciarse los efectos de frontera, situacin que no se presenta en los
resultados de las soluciones analticas, ya que suponen yacimiento infinito.
57
Captulo 4
Pi
Qo
rw
h
k
Ct
Bo
Unidades Prcticas
4351.1[psi]
3773.8[bpd]
0.3281[ft]
164.0[ft]
1000/100/1[mD]
0.01/0.1/0.3[Frac.]
4.0E-05/4.0E-05/4.0E-05[psi-1]
0.16/0.16/0.16[Frac.]
9.29E-02/9.29E-02[ft-2]
0.695[cp]
2.0[bl/bl]
2.44E-02/2.44E-01/7.32E-01
1.00E-03/1.00E-05
Unidades del SI
30.0E+06[Pa]
600[m3/Seg]
0.1[m]
50.0[m]
9.87E-13/9.87E-14/9.87E-16[m2]
0.01/0.1/0.3[Frac.]
5.8E-09/5.8E-09/5.8E-09[Pa-1]
0.16/0.16/0.16[Frac.]
1.0/1.0[m-2]
[Pa6.95E-04Seg]
2.0[m3/m3]
Tabla 4.1 Datos Generales Utilizados para la Validacin del Modelo Propuesto
Presin (Mpa) Bo (m3/m3)

0.101
1.062
1.825
1.150
3.548
1.207
6.996
1.295
13.891
2.000
17.338
2.000
20.786
2.000
27.680
2.000
62.154
2.000
o (cp)
1.040
0.975
0.910
0.830
0.695
0.695
0.695
0.695
0.695
o (kg/m3)
740.754
697.511
677.378
656.820
624.634
613.571
605.129
593.420
593.420
Rs (m3/m3)
0.178
16.118
32.057
66.073
113.267
138.022
165.627
226.180
226.180
Tabla 4.2 Datos PVT Utilizados para la Validacin del Modelo Propuesto
58
Captulo 4
Pwf vs Log (t)

Sol. Analtica
Triple Porosidad
4355.0
Simulador Triple
Porosidad
1=0.01
k1=1000
Pwf (psi)
4345.0
Sol. Analtica
Doble Porosidad
Simulador Doble
Porosidad
2=0.06
k2=100 mD
4335.0
Sol. Analtica
Simple Porosidad
3=0.1
k3=1 mD
4325.0
Simulador Simple
Porosidad
4315.0
Efectos de
Frontera
4305.0
1.0E-07
1.0E-06
1.0E-05
1.0E-04
1.0E-03
1.0E-02
1.0E-01
1.0E+00
1.0E+01
t (horas)
Figura 4.1 Resultados Obtenidos para las Soluciones Numrica y Analtica
40
59
1.0E+02
1.0E+03
Captulo 5
Resultados
CAPTULO 5
RESULTADOS
El objetivo del presente captulo es mostrar el funcionamiento del simulador una vez
que se ha incorporado el modelo de triple porosidad, por medio de la simulacin de
tres casos sintticos. Los casos simulados son:
5.1 Caso 1. Simulacin de un Sistema de Triple Porosidad por Medio del
Modelo de Doble Porosidad.
El objetivo de este caso es mostrar la importancia de la caracterizacin y modelado
correctos de los yacimientos de mltiple porosidad. Para lo anterior se utiliz un caso
sinttico de un yacimiento compuesto de tres medios porosos, el cual es explotado por
medio de un pozo productor ubicado al centro de la malla. En cada uno de los
extremos de la malla se cuenta adems con un pozo inyector, cuyo objetivo es
mantener la presin en el yacimiento. La malla utilizada para este caso fue de 20x20x3
celdas. Los datos utilizados para la simulacin del modelo de triple porosidad se
muestran en las tablas 5.1 y 5.2.
Propiedad
Unidades Prcticas
Unidades del SI
4206.094[psi]
4014.644[psi]
3144.9, 3144.9, 3144.9[bpd]
6289.8, 6289.8[bpd]
16568.24[ft]
16568.24[ft]
164.0[ft]
Pi
Pb
Qw(iny)
Qo(prod)
Hx
Hy
Hz
29.00E+06[Pa]
27.68E+06[Pa]
5.78E-3, 5.78E-3, 5.78E-3[m3/seg]
11.575E-3, 11.575E-3[m3/Seg]
5050[m]
5050[m]
50.0[m]
Tabla 5.1. Datos Generales Utilizados para la Simulacin del Caso 1.
Propiedades para el modelo de triple porosidad:
Propiedad
Unidades Prcticas
300/50/10
0.03/0.04/0.04
Ct
Swi
Unidades del SI
[mD]
[Frac.]
-1
4.35E-4/4.35E-4/4.35E-4 [psi ]
0.0/0.16/0.16
[Frac.]
2.96E-13/4.93E-14/9.87E-16 [m2]
0.03/0.04/0.04
[Frac.]
6.31D-08/6.31E-08/6.31E-08 [Pa-1]
0.0/0.16/0.16
[Frac.]
Tabla 5.2. Datos Utilizados para la Simulacin del Modelo de Triple Porosidad.
Una vez realizada la simulacin del yacimiento tomando en cuenta las propiedades de
los tres medios, se busc la forma de simular este yacimiento por medio del modelo de
60
Captulo 5
Resultados
doble porosidad, para lo cual se fusionaron los medios dos y tres, tomando en cuenta
las propiedades de ambos. Para el caso de la porosidad, nicamente se sumaron los
valores de ambos medios, pero para el caso de la permeabilidad, fue necesario probar
con tres tipos de promedios (aritmtico, armnico y geomtrico) de ambos valores y
simular los tres casos para observar cual de los tres valores arrojaba mejores
resultados. Los datos utilizados para la simulacin de los modelos de doble porosidad
se muestran en las tablas 5.3, 5.4 y 5.5.
Propiedades para los modelos de doble porosidad:
a) Promedio Aritmtico
Propiedad
Unidades Prcticas
300/30
Unidades del SI
[mD]
2.96E-13/2.96E-14
0.03/0.08
[Frac.]
0.03/0.08
Ct
4.35E-04/4.35E-04
[psi-1]
6.31D-08/6.31E-08
Swi
0.0/0.16
[Frac.]
0.0/0.16
[m2]
[Frac.]
[Pa-1]
[Frac.]
Tabla 5.3. Datos Utilizados para la Simulacin del Modelo de Doble Porosidad. Promedio Aritmtico.
b) Promedio Armnico
Propiedad
Unidades Prcticas
300/16.6
Unidades del SI
[mD]
2.96E-13/1.64E-14
0.03/0.08
[Frac.]
0.03/0.08
Ct
4.35E-04/4.35E-04
[psi-1]
6.31D-08/6.31E-08
Swi
0.0/0.16
[Frac.]
0.0/0.16
[m2]
[Frac.]
[Pa-1]
[Frac.]
Tabla 5.4. Datos Utilizados para la Simulacin del Modelo de Doble Porosidad. Promedio Armnico.
c) Promedio Geomtrico
Propiedad
Unidades Prcticas
Unidades del SI
300/22.36
[mD]
0.03/0.08
[Frac.]
Ct
4.35E-04/4.35E-04
Swi
0.0/0.16
[psi
2.96E-13/2.2E-14
-1
[Frac.]
0.03/0.08
6.31D-08/6.31E-08
0.0/0.16
[m2]
[Frac.]
[Pa-1]
[Frac.]
Tabla 5.5. Datos Utilizados para la Simulacin del Modelo de Doble Porosidad. Promedio Geomtrico.
La figura 5.1 muestra los resultados de las cuatro simulaciones realizadas. Como se
puede observar en dicha figura, aunque en los primeros aos de la simulacin los
resultados de ambos modelos (triple y doble porosidad) son muy similares, al final de
la simulacin la respuesta de presin obtenida con el modelo de triple porosidad (lnea
61
Captulo 5
Resultados
rosa), es muy diferente a la respuesta obtenida con los modelos de doble porosidad, lo
que conducira a toma de decisiones equivocadas.
p wf vs t
4500
Productor 1
(10,10,1)
4000
3500
Doble
Porosidad
(Aritmtico)
pwf (psi)
3000
2500
2000
Doble
Porosidad
(Armnico)
1500
1000
Doble
Porosidad
(Geomtrico)
500
0
0
400
800
1200
1600
2000
2400
2800
3200
3600
Tiempo (das)
Figura 5.1. Resultados de la Simulacin para los Modelos de Triple y Doble Porosidad.
62
Captulo 5
5.2
Resultados
Caso 2. Simulacin de Flujo Multifsico en un Sistema de Triple

Porosidad.
Para este caso se utiliz una malla de 101x101x1 con 2 pozos productores y 3 pozos
que inyectan agua despus de 500 y 750 das de produccin. Todos los pozos fueron
posicionados de tal forma que la distribucin de presiones y saturaciones fuera
simtrica. Los datos generales utilizados para la simulacin de este caso se muestran
en la tabla 5.6.
Propiedad
Pi
Pb
Qw(iny)
Qo(prod)
Hx
Hy
Hz
k
Ct
Swi
Unidades Prcticas
Unidades del SI
4206.094[psi]
29.00E+06[Pa]
4014.644[psi]
27.68E+06[Pa]
12579.6, 12579.6, 12579.6[bpd]
23.15E-3, 23.15E-3, 23.15E-3[m3/seg]
12579.6, 12579.6[bpd]
23.15E-3, 23.15E-3[m3/Seg]
16568.24[ft]
5050[m]
16568.24[ft]
5050[m]
164.0[ft]
50.0[m]
1.48E-13/4.93E-14/9.87E-15[m2]
150/50/10[mD]
0.03/0.04/0.08[Frac.]
0.03/0.04/0.08[Frac.]
4.35D-09/4.35E-09/4.35E-09[Pa-1]
3.0E-05/3.0E-05/3.0E-05[psi-1]
0.0/0.18/0.18[Frac.]
0.0/0.18/0.18[Frac.]
Los datos PVT para este caso se muestran en las tablas 5.7 a 5.9.
Tabla 5.7 Datos PVT para el Aceite Utilizados para la Simulacin del Caso 2.
Po (MPa)
0.10135
1.82504
3.54873
6.99611
13.89086
17.33824
20.78562
27.68037
62.15415
Bo (m3/m3)
1.062
1.150
1.207
1.295
1.435
1.500
1.565
1.695
1.579
o (cp)
1.04000
0.97500
0.91000
0.83000
0.69500
0.64100
0.59400
0.51000
0.74000
63
o (kg/m3)
740.753
697.511
677.378
656.820
624.634
613.571
605.192
593.420
637.015
Rso(m3/m3)
0.178093
16.117453
32.056812
66.072651
113.267402
138.022385
165.626862
226.178618
226.178618
Captulo 5
Resultados
Tabla 5.8 Datos PVT para el Gas Utilizados para la Simulacin del Caso 2.
Pg (MPa)
0.10135
1.82504
3.54873
6.99611
13.89086
17.33824
20.78562
27.68037
62.15415
Bg (m3/m3)
0.936
0.068
0.035
0.018
0.009
0.007
0.006
0.005
0.002
g (cp)
0.00800
0.01000
0.01100
0.01400
0.01900
0.02100
0.02300
0.02700
0.04700
g (kg/m3)
1.036
14.282
27.528
54.022
107.007
133.469
159.915
212.958
447.432
rsg(m3/m3)
0
0
0
0
0
0
0
0
0
Tabla 5.9 Datos PVT para el Agua Utilizados para la Simulacin del Caso 2.
Pw (MPa)
0.10135
1.82504
3.54873
6.99611
13.89086
17.33824
20.78562
27.68037
62.15415
w (cp)
0.310
0.310
0.310
0.310
0.310
0.310
0.310
0.310
0.310
Bw (m3/m3)
1.041
1.041
1.040
1.038
1.035
1.034
1.032
1.029
1.029
w (kg/m3)
996.4848
997.2328
997.9819
999.4834
1002.5001
1004.0131
1005.5350
1008.5884
1008.5884
La simulacin de este caso se llev a cabo en dos etapas: en una primera etapa se
realiz la simulacin con valores de presin capilar nicamente para el medio 3 y en
una segunda etapa se realiz la simulacin sin valores de presin capilar; esto con el
fin de observar el efecto entre simular con y sin valores de presin capilar. Los datos
de presin capilar y permeabilidad utilizados se muestran en la figura 5.2. Cabe
mencionar que por ser un caso sinttico se decidi usar curvas de permeabilidad
relativa en forma de x para la fractura, aunque en la realidad estas curvas pueden
tener un comportamiento diferente (39).
64
Captulo 5
Resultados
Pcwo vs Sw
Pcgo vs Sg
5.0
20
4.5
18
16
3.5
14
3.0
12
2.5
Pcw o
Pcgo
Medio 1
Medio 3
4.0
Medio 2
1.5
1.0
Medio 3
10
2.0
Medio 1
Medio 2
0.5
0.0
0
0.0
0.2
0.4
0.6
0.8
1.0
0.0
0.2
0.4
Sg
0.8
1.0
Sw
Permeabilidades Agua-Aceite
Permeabilidad Gas-Aceite
1.0
1.0
Krg (medio 1)
0.9
krw (medio 1)
0.9
0.8
0.8
krow (medio 1)
Krog (medio 1)
0.7
0.7
0.6
0.6
krg (medio 3)
K ro w
K rog
0.6
0.5
0.4
Krog (medio 3)
Krw (medio 2)
0.5
Krow (medio 2)
0.4
0.3
0.3
Krw (medio 3)
Krg (medio 2)
0.2
0.2
0.1
0.1
krow (medio 3)
Krog (medio 2)
0.0
0.0
0.0
0.2
0.4
0.6
0.8
0.0
1.0
Sg
0.2
0.4
0.6
0.8
Sw
Figura 5.2 Grficas de Presin Capilar y Permeabilidad Utilizados en la Simulacin del Caso 2.
65
1.0
Captulo 5
Resultados
En la figura 5.3 se observa el comportamiento de la presin de fondo fluyendo tanto en

los pozos productores como en los inyectores. Al inicio de la simulacin el
comportamiento en ambos casos (con y sin presin capilar) es similar, pero a partir de
los 1000 das de simulacin se observa que la cada de presin para el caso en que se
utilizaron datos de presin capilar (Pc0), fue mayor que en el caso donde no se
utilizaron. La similitud al inicio de la simulacin se debe a que los medios uno y dos
tienen pc=0, y es hasta que se sienten los efectos del medio 3 (con Pc0), que se
presenta una diferencia en los resultados de pwf, ya que para que el aceite fluya se
tienen que vencer los efectos de las fuerzas capilares, las cuales se encuentran
relacionadas con la tensin interfacial, la mojabilidad y el tamao de poro, lo que
provoca una cada adicional de la presin.
En la figura tambin se observa que a pesar de que a los 500 das de simulacin se
inicia la inyeccin de agua en uno de los pozos, esto no es suficiente para mantener la
presin en el yacimiento y esta sigue disminuyendo, por lo que, a los 750 das se inicia
la inyeccin en dos pozos ms, y no es sino hasta aproximadamente l000 das de
simulacin cuando se observa un cambi en el comportamiento de la presin, la cual
tiene un comportamiento relativamente lineal hasta aproximadamente 2000 das de
simulacin, cuando nuevamente empieza a disminuir.
En la figura 5.4 se encuentra graficado el gasto de aceite vs tiempo. En esta figura
puede observarse que el gasto de aceite se mantiene constante durante todo el tiempo
de simulacin. Lo anterior se debe a que en el simulador se encuentra especificado
como constante el gasto, en dicho caso y de acuerdo a la Ley de Darcy, si el gasto se
mantiene constante, la que vara es la presin. Como es lgico, en los pozos de
inyeccin no se tiene produccin de aceite, por lo tanto, estos presentan un gasto de
aceite igual a cero.
En la figura 5.5 se encuentra graficado el gasto de gas vs tiempo. En esta figura se
observa que al inicio y hasta aproximadamente 1000 das de simulacin se tiene un
comportamiento ascendente en el gasto de gas, esto se debe a que al ir disminuyendo
la presin la presencia de gas aumenta dentro del yacimiento, ya que el gas que
inicialmente se encontraba disuelto en el aceite, se libera. De 1000 a 2000 das de
simulacin se presenta un comportamiento descendente, esto debido al mantenimiento
de presin que se presenta con la inyeccin de agua, finalmente, de los 2000 das y
hasta el final de la simulacin se presenta nuevamente un comportamiento
ascendente, nuevamente provocado por el depresionamiento del yacimiento. En esta
figura tambin se observa que el gasto de gas para el caso en el que se tienen datos
de presin capilar, es mayor debido a que la cada de presin tambin es mayor (figura
5.3), por lo tanto, se libera mayor cantidad de gas.
En la figura 5.6 se encuentra graficada la RGA vs tiempo. Los cambios en el
comportamiento de la RGA estn directamente relacionados con la presencia de gas
dentro del yacimiento. Al inicio de la simulacin, la RGA aumenta, ya que la presencia
de gas aumenta en el yacimiento (figura 5.5), para posteriormente disminuir y volver
aumentar, de acuerdo a la presencia de gas en el yacimiento.
En la figura 5.7 se encuentra graficado el gasto de agua vs tiempo, donde se puede
observar que el agua en los pozos productores se presenta hasta aproximadamente
2000 das de produccin, esto es, 1500 das despus de que se inicia la inyeccin de
agua.
66
Captulo 5
Resultados
Pwf vs Tiempo
Pwf (psi)
4500
Productor 1 (Con Presin Capilar)
4000
Productor 1 (Sin Presin Capilar)
3500
3000
2500
Inyector 1 (Con Presin Capilar)
Inyector 1 (Sin presin capilar)
2000
1500
1000
500
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
Tiempo (das)
Figura 5.3 Grfica de Pwf vs Tiempo para la Simulacin de Flujo Multifsico en un Sistema de Triple
Porosidad con y sin Presiones Capilares.
Qo vs Tiempo
2500
2000
Qo (m3/dia)
1500
1000
500
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
Tiempo (das)
Figura 5.4 Grfica de Qo vs Tiempo para la Simulacin de Flujo Multifsico en un Sistema de Triple
67
Captulo 5
Resultados
Qg vs Tiempo
45
Qg (m3/dia)
40
35
30
25
20
15
10
5
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
Tiempo (das)
Figura 5.5 Grfica de Qg vs Tiempo para la Simulacin de Flujo Multifsico en un Sistema de Triple
RGA vs Tiempo
2500
2000
1500
RGA
1000
500
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
Tiempo (das)
Figura 5.6 Grfica de RGA vs Tiempo para la Simulacin de Flujo Multifsico en un Sistema de Triple
68
Captulo 5
Resultados
Qw vs Tiempo
2500
2000
Qw (m 3/dia)
1500
1000
500
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
Tiempo (das)
Figura 5.7 Grfica de Qw vs Tiempo para la Simulacin de Flujo Multifsico en un Sistema de Triple
Finalmente, en las figuras 5.8 a 5.11, se muestran los grficas de contorno para las
propiedades de Po, Sw, Sg y So, obtenidas en el ltimo paso de tiempo, donde pueden
apreciarse claramente las diferencias entre los resultados obtenidos para la simulacin
con tablas de presin capilar y sin ellas.
69
Captulo 5
Resultados
Po (Sin Pc)
100
3200
3400
3600
3800
4000
90
80
70
Y Data
60
50
40
30
20
10
10
20
30
40
50
60
70
80
90
100
X Data
Po (Con Pc)
100
3200
3400
3600
3800
4000
90
80
70
Y Data
60
50
40
30
20
10
10
20
30
40
50
60
70
80
90
100
X Data
Figura 5.8. Resultados de Presin de Aceite para el Medio 1 con y sin Pc.
70
Captulo 5
Resultados
Sw (Sin Pc)
100
0.0
0.2
0.4
0.6
0.8
1.0
90
80
70
Y Data
60
50
40
30
20
10
10
20
30
40
50
60
70
80
90
100
X Data
Sw (Con Pc)
100
0.0
0.2
0.4
0.6
0.8
1.0
90
80
70
Y Data
60
50
40
30
20
10
10
20
30
40
50
60
70
80
90
100
X Data
Figura 5.9. Resultados de Saturacin de Agua para el Medio 1 con y sin Pc.
71
Captulo 5
Resultados
Sg (Sin Pc)
100
0.00
0.02
0.04
0.06
0.08
0.10
0.12
0.14
0.16
90
80
70
Y Data
60
50
40
30
20
10
10
20
30
40
50
60
70
80
90
100
X Data
Sg (Con Pc)
100
0.00
0.02
0.04
0.06
0.08
0.10
0.12
0.14
0.16
90
80
70
Y Data
60
50
40
30
20
10
10
20
30
40
50
60
70
80
90
100
X Data
Figura 5.10. Resultados de Saturacin de Gas para el Medio 1 con y sin Pc.
72
Captulo 5
Resultados
So (Sin Pc)
100
0.0
0.2
0.4
0.6
0.8
1.0
90
80
70
Y Data
60
50
40
30
20
10
10
20
30
40
50
60
70
80
90
100
X Data
So (Con Pc)
100
90
0.0
0.2
0.4
0.6
0.8
1.0
80
70
Y Data
60
50
40
30
20
10
10
20
30
40
50
60
70
80
90
100
X Data
Figura 5.11. Resultados de Saturacin de Aceite para el Medio 1 con y sin Pc.
73
Captulo 5
Resultados
5.3 Caso 3. Simulacin de Flujo Multifsico en un Sistema de Triple Porosidad.

Para este caso se utiliz una malla de 20x20x3 en la cual se tienen 5 pozos: un
productor ubicado en el centro y cuatro inyectores ubicados en cada uno de los
extremos de la malla, los cuales empiezan a inyectar agua a partir de los 750 das de
simulacin. La malla utilizada y la ubicacin de los pozos se muestran en la figura
5.12.
Figura 5.12. Malla Utilizada en la Simulacin del Caso 3.
Propiedad
Pi
Pb
Qw(iny)
Qo(prod)
Hx
Hy
Hz
k
Ct
Swi
Unidades Prcticas
Unidades del SI
4206.094[psi]
29.09E+06[Pa]
4014.644[psi]
27.68E+06[Pa]
17.36E-3, 17.36E-3, 17.36E-3
9434.7, 9434.7, 9434.7,
9434.7[bpd]
17.36E-3[m3/seg]
15724.15[bpd]
28.93E-3[m3/Seg]
6561.67[ft]
2000.0[m]
6561.67[ft]
2000.0[m]
492.12[ft]
150.0[m]
300/50/10[mD]
2.96E-13/4.93E-14/9.87E-15[m2]
0.03/0.04/0.04[Frac.]
0.03/0.04/0.04[Frac.]
-1
3.0E-06/3.0E-06/3.0E-06[psi ]
4.35E-10/4.35E-10/4.35E-10[Pa-1]
0.0/0.16/0.16[Frac.]
0.0/0.16/0.16[Frac.]
74
Captulo 5
Resultados
El propsito de este caso es observar el efecto de las fuerzas gravitacionales, las

cuales se reflejan en la distribucin de los fluidos dentro del yacimiento. Los datos
generales, de presin capilar y permeabilidad utilizados en este caso, se presentan en
la tabla 5.10 y las figuras 5.13 y 5.14. Los datos PVT para el aceite, gas y agua, son
los mismos que en el caso anterior.
Permeabilidad Gas-Aceite
Permeabilidad Agua-Aceite
1.0
1.0
Krg
(medio 1)
0.9
0.8
0.8
Krow
(medio 1)
Krog
(medio 1)
0.7
0.6
0.7
Krg
(medio 2)
K ro w
Krog
(medio 2)
0.4
Krw
(medio 2)
0.6
0.5
0.5
Krow
(medio 2)
0.4
0.3
0.3
Krg
(medio 3)
Krw
(medio 3)
0.2
0.2
Krog
(medio 3)
0.1
Krow
(medio 3)
0.1
0.0
0.0
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
0.0
0.1
0.2
0.3
0.4
Sg
0.5
0.6
0.7
0.8
0.9
1.0
Sw
Figura 5.13 Grficas de Permeabilidad Utilizadas en la Simulacin del Caso 3.

Presin Capilar Gas-Aceite
Presin Capilar Agua-Aceite
2.5
70.0
60.0
Pcgo
(medio 1)
2.0
Pcwo
(medio 1)
50.0
1.5
Pcw o
40.0
Pcgo
(medio 2)
Pcgo
K ro g
Krw
(medio 1)
0.9
Pcwo
(medio 2)
30.0
1.0
20.0
0.5
Pcgo
(medio 3)
Pcwo
(medio 3)
10.0
0.0
0.0
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
0.0
1.0
0.1
0.2
0.3
0.4
0.5
0.6
Sw
Sg
Figura 5.14 Grficas de Presin Capilar Utilizadas en la Simulacin del Caso 3.
75
0.7
0.8
0.9
1.0
Captulo 5
Resultados
En las figuras 5.15 y 5.16, se encuentran graficados los resultados de pwf, vs t y de Qo

vs t, respectivamente. En la figura 5.16 puede observarse que el gasto de aceite al
inicio y hasta aproximadamente 1250 das de simulacin es de 2500 m3/da, para
posteriormente disminuir hasta menos de 500 m3/da. Lo anterior se debe a la cada
de presin que se observa en la figura 5.15, esto es, el pozo produce al gasto
especificado mientras la presin en el yacimiento no llega a un mnimo (1000 psi) y,
una vez que cae hasta dicho mnimo, produce nicamente la cantidad que alcanza con
dicha presin. Posteriormente el gasto vuelve a incrementarse una vez que se siente el
efecto de los pozos inyectores de agua.
Pwf vs t
4500
4000
Productor 1
(10,10,1)
3500
Inyector 1
(2,2,1)
Pwf (psi)
3000
2500
Inyector 2
(19,19,1)
2000
Inyector 3
(2,19,1)
1500
1000
Inyector 4
(19,2,1)
500
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
Tiempo (das)
Figura 5.15 Grfica de Pwf vs Tiempo para la Simulacin de Flujo Multifsico

en un Sistema de Triple Porosidad (Caso 3).
Qo vs t
Qo (m3/da)
3000
2500
Productor 1
(10,10,1)
2000
Inyector 1
(2,2,1)
1500
Inyector 2
(19,19,1)
1000
Inyector 3
(2,19,1)
500
Inyector 4
(19,2,1)
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
Tiempo (das)
Figura 5.16 Grfica de Qo vs Tiempo para la Simulacin de Flujo Multifsico

76
Captulo 5
Resultados
En la figura 5.17 se muestran los resultados de gasto de gas vs tiempo. En esta figura
se observa un aumento en el gasto de gas ocasionado por el depresionamiento del
yacimiento, para posteriormente volver a disminuir al sentirse el efecto de los pozos
inyectores.
Qg vs t
350
300
Productor 1
(10,10,1)
250
Qg (m3/da)
Inyector 1
(2,2,1)
200
Inyector 2
(19,19,1)
150
Inyector 3
(2,19,1)
100
Inyector 4
(19,2,1)
50
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
Tiempo (das)
Figura 5.17 Grfica de Qg vs Tiempo para la Simulacin de Flujo Multifsico

En la figura 5.18 se observa que el agua nicamente se encuentra presente en los

pozos inyectores y nunca se presenta en el pozo productor.
Finalmente la figura 5.19 contiene graficados los resultados de RGA vs t. De forma
similar al gasto de aceite y gas, es evidente que el comportamiento de la RGA se
encuentra influenciado por los cambios de presin dentro del yacimiento.
Las figuras 5.20 a 5.23 muestran la distribucin de presiones de aceite y de
saturaciones de aceite, gas y agua al final de la simulacin. En dichas figuras puede
apreciarse como se distribuyen los fluidos en el yacimiento debido a las fuerzas
gravitacionales. Esto es, el agua que es el fluido de mayor densidad se distribuye en el
fondo del yacimiento, posteriormente el aceite y finalmente el gas que es el fluido con
menor densidad.
77
Captulo 5
Resultados
Qw vs t
1800
1600
Productor 1
(10,10,1)
1400
Inyector 1
(2,2,1)
Qw (m3/da)
1200
1000
Inyector 2
(19,19,1)
800
600
Inyector 3
(2,19,1)
400
Inyector 4
(19,2,1)
200
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
Tiempo (das)
Figura 5.18 Grfica de Qw vs Tiempo para la Simulacin de Flujo Multifsico

RGA vs t
9000
8000
Productor 1
(10,10,1)
7000
Inyector 1
(2,2,1)
6000
5000
RGA
Inyector 2
(19,19,1)
4000
Inyector 3
(2,19,1)
3000
2000
Inyector 4
(19,2,1)
1000
0
0
500
1000
1500
2000
2500
3000
3500
4000
Tiempo (das)
Figura 5.19 Grfica de RGA vs Tiempo para la Simulacin de Flujo Multifsico

78
4500
Captulo 5
Resultados
Figura 5.20 Distribucin de Presiones de Aceite, Po (Caso 3).
Figura 5.21 Distribucin de Saturaciones de Aceite, So (Caso 3).
79
Captulo 5
Resultados
Figura 5.22 Distribucin de Saturaciones de Gas, Sg (Caso 3).
Figura 5.23 Distribucin de Saturaciones de Agua, Sw (Caso 3).
80
Captulo 6
Conclusiones
CAPTULO 6
CONCLUSIONES
Se desarroll un modelo de triple porosidad-una permeabilidad para describir el

flujo de fluidos en Yacimientos Naturalmente Fracturados que presentan triple
porosidad.
Se incorpor el modelo propuesto a un simulador trifsico-multidimensional de

doble porosidad, extendindolo de doble a triple porosidad, considerando
medios inelsticos.
Se valid el modelo con la solucin analtica obtenida por Camacho y cols (8) y
por Rodrguez y cols (27), para un sistema de triple porosidad-una
permeabilidad, bajo las condiciones de: 1) flujo radial monofsico, 2) fluido
ligeramente compresible, 3) gasto constante y, 4) yacimiento infinito;
obtenindose un buen ajuste entre la solucin numrica y la analtica.
Se aplic el modelo en tres casos sintticos:

o
El primer caso consisti en simular el comportamiento de un yacimiento

compuesto por tres medios con los modelos de doble y triple porosidad,
respectivamente. Para poder realizar la simulacin con el modelo de
doble porosidad se agruparon las porosidades de los medios dos y tres,
y se realizaron diferentes promedios (aritmtico, armnico y
geomtrico), entre las permeabilidades de dichos medios. Los resultados
mostraron que todas las respuestas obtenidas con el modelo de doble
porosidad presentan diferencias con respecto a la respuesta del modelo
de triple porosidad. Dichas diferencias van desde un 9% (promedio
geomtrico), hasta un 79% (promedio armnico).
El segundo caso consisti en simular el comportamiento de un

yacimiento que presenta triple porosidad. Dicho yacimiento cuenta con
dos pozos productores de aceite y tres pozos inyectores de agua. El
objetivo de este caso bsicamente fue probar el funcionamiento del
simulador. Como complemento de este caso se realizaron corridas con y
sin curvas de presin capilar para observar el efecto de las mismas en la
respuesta de presin. Los resultados obtenidos mostraron que se
presenta una mayor cada de presin en el caso en que se utilizaron las
curvas de presin capilar. Lo anterior se debe a que para que el aceite
fluya necesita vencer los efectos de las fuerzas capilares, lo que provoca
dicha diferencia en el comportamiento de presin.
Finalmente, el caso tres consisti en simular el comportamiento de un

yacimiento que presenta triple porosidad y que cuenta con cinco pozos:
un productor y cuatro inyectores. El objetivo de este caso fue probar el
funcionamiento del simulador bajo condiciones de flujo multifsico
(aceite, gas y agua). Los resultados obtenidos fueron satisfactorios ya
que el simulador fue capaz de representar el comportamiento de las tres
fases, en cada uno de los medios presentes en el yacimiento.
81
Captulo 6
Conclusiones
Se mostr la importancia de la correcta caracterizacin y modelado de los

yacimientos de mltiple porosidad.
Para una correcta caracterizacin de este tipo de yacimientos se recomienda el

uso de pruebas de presin, las cuales pueden apoyarse en las grficas
desarrolladas por Rodrguez y cols (27). para la caracterizacin de yacimientos
con triple porosidad anidada. La informacin obtenida por medio de las pruebas
de presin puede complementarse con la petrofsica, ya que muchas veces es
posible observar la presencia de mltiples porosidades en el anlisis de ncleos.
82
REFERENCIAS
1. Abdassah, D. y Ershaghi, I.: Triple Porosity Systems for Representing Naturally
Fractured Reservoirs, SPEFE, 1, 113-127, Abril, 1986.
2. Aifantis, E. C.: On the Problem of Diffusion in Solids, Acta Mechanics, 1980.
3. Ausbrooks, R. y Hurley, N.: Pore-Size Distributions in Vuggy Carbonates From Core
Images, NMR, and Capillary Pressure, SPE 56506, Octubre 1999.
4. Barenblatt, G.I., Zheltov, I.P., y Kochina, I.N.: Basic Concepts in the Theory of
Seepage of Homogenous Liquid in Fissured Rocks, PMM, Sov. Appl. Math. Mech. 24-5,
1960.
5. Barenblatt, G.I.: The Matematical Model of Nonequilibrium Effects in Water-Oil
Displacement, SPE 75169, 2002.
6. Blaskovich, F. T. y cols.: Numerical Simulation of Naturally Fractured Reservoirs, SPE
15627, 1983.
7. Bourdet, D.: Pressure Behavior of Layered Reservoirs with Crossflow, paper SPE
13628 presentado en la SPE California Regional Meeting, Bakersfield, California, 27-29
Marzo, 1985.
8. Camacho, V. R., Vsquez, C. M., Castrejn A. R. y Arana, O. V. : Pressure
Transient and Decline Curve Behaviors in Naturally Fractured Vuggy Carbonate
Reservoirs, SPE 77689, 2002.
9. De la Porte, J.J., Kossack, C.A. y Zimmerman, R.W.: The Effect of Fracture Relative
Permeabilities and Capillary Pressures on the Numerical Simulation of Naturally
Fractured Reservoirs, SPE 95241, Octubre, 2005.
10. Diodato, D.M.: A Compendium of Fracture Flow Models, Training Center,
Hohenfels, Germany, 1994.
11. Grupo de Simulacin Numrica de Yacimientos: Desarrollo de un Simulador
Multifsico Multipropsito", PEMEX-UNAM, 2006.
12. Ertekin, T. : Basic Applied Reservoir Simulation, SPE Textbook Series, USA, 2001.
13. Hill, A.C. y Thomas, G.W.: A New Approach for Simulating Complex Fractured
14. Iffly, R., Rousselet, D.C. y Vermeulen, J.L.: Fundamental Study of Imbibition in
Fissured Oil Field, SPE 4102, 1962.
15. Infante del Rio, J.A.: Metodos Numericos. Teoria, Problemas y Practicas con
Matlab, Piramide, Mexico, 1999.
16. Kang, Z. y Wu, Y.S.: Modeling Multiphase flow in Naturally Fractured Vuggy
Petroleum Reservoirs, SPE 102356, 2006.
17. Kazemi, H. : Pressure-Transient Analysis of Naturally Fractured Reservoirs with
Uniform Fracture Distribution, SPEJ 451-462, Diciembre, 1969.
18. Kazemi, H. y Gilman, J.: Improvements in Simulation of Naturally Fractured
19. Kossack, C. A., y Gurpinar, O.: A Methodology for Simulation of Vuggy and
Fractured Reservoirs, SPE 66366, 2001.
20. Lucia, F.J.: Petrophysical Parameters Estimated From Visual Descriptions of
Carbonate Rocks: A Field Classification of Carbonate Pore Space, SPE 10073, 1983.
86
21. Naimi-Tajdar, R., C. Han, K. Sepehrnoori, T.J. Arbogast and M.A. Miller: A Fully
Implicit, Compositional, Parallel Simulator for IOR Processes in Fractured Reservoirs,
SPE 100079, 2006.
22. Nelson, R.: Geological Analysis of Naturally Fractured Reservoirs, Gulf Professional
Publishing, USA, 2001.
23. Prats, M.: Thermal Recovery Monograph, Society of Petroleum Engineers, USA,
1986.
24. Pea, O. y Resendiz, J.T.: Simulacin Multifsica Tridimensional en Yacimientos
Naturalmente Fracturados, F.I. UNAM, 2006.
25. Pruess, K. y Narasimhan, T. N.: A Practical Method for Modeling Heat and Flow in
Fractured Porous Media, SPEJ, 1985.
26. Reiss, L.H.: Reservoir Engineering in Fractured Reservoirs, 1976.
27. Rodriguez, F., Arana, O.V. y Cinco, L. H.: Well Test Characterization of Small and
Large Scale Secondary Porosity in Naturally Fractured Reservoirs, SPE 90287, 2004.
28. Rodrguez de la Garza, F.: Fundamentos de Simulacin Numrica de Yacimientos,
notas de la material, DEPFI, UNAM, 2005.
29. Rodrguez de la Garza, F.: Apuntes del Curso Simulacin Avanzada de yacimientos,
DEPFI, UNAM, (2005).
30. Rubin, B.: Simulating Gravity Drainage and Reinfiltration with a Subdomain-DualPermeability Hybrid Fracture Model, SPE 106191-MS, Febrero, 2007.
31. Siavoshi, J.: Triple Porosity Systems Bring New Insight to Naturally Fractured
Reservoirs, work presented in the Eleventh Oil, Gas and Petrochemical congress, Oct.
2001, Teheran, Iran.
32. Thomas, L.K., Dixon, T.N. y Pierson, R.G. : Fractured Reservoir Simulation , SPEJ
(1983).
33. Thomas, G.W., A New Approach for Simulating Complex Fractured Reservoirs, SPE
13537, 1985.
34. Van Golf-Racht, T: Fundamentals of Fractured Reservoir Engineering, Elsevier, New
York, 1982.
35. Warren, J.E. y Root, P.J.: The Behavior of Naturally Fractured Reservoirs, SPEJ,
245-255, Septiembre, 1963.
36. Wu, Y. S. y cols.: A Physically Based Approach for Modeling Multiphase FractureMatriz Interaction in Fractured Porous Media, Advances in Water Resources, Elsevier,
2004.
37. Wu, Y. S. y cols.: A Triple-Continuum Approach for Modeling Flow and Transport
Processes in Fractured Rock, Journal of Contaminant Hydrology, Elsevier, Enero,
2004.
38. Wu, Y. S. y cols.: Analysis of Flow Behavior in Fractured Lithophysal Reservoirs,
Journal of Contaminant Hydrology, Lawrence Berkeley National Laboratory, 2006.
39. Yamamoto, R.H., y cols.: Compositional Reservoir Simulator for Fissured Systems
The Single Block Model, SPEJ (1971).
40. Zhang, F.: The Schur Complement and its Aplications, Numerical Methods and
Algorithms, Springer Science, USA, 2005.
87
Apndice A
APNDICE A
OBTENCIN DE LAS ECUACIONES DE FLUJO MULTIFSICO
En este apndice se tratan los principios fundamentales de la fsica relacionados con
la transferencia de masa, y se determina como se relacionan las representaciones
matemticas de stos para desarrollar las ecuaciones diferenciales que representan
los balances de masa dentro de un yacimiento.
En el trabajo presente se utiliza el modelo de flujo de fluidos tipo Beta, donde el gas
y aceite obtenidos despus de separadores son los pseudo-componentes que forman
las fases aceite y gas en el yacimiento. Lo anterior se ejemplifica mediante el
diagrama siguiente:
COMPONENTE
FASE
Aceite del separador

Gas del separador
Agua
Aceite
Gas
Agua
Se considera que:
Existen tres fases mviles posibles en un yacimiento
Las fases aceite y gas se encuentran en equilibrio en espacio y tiempo
La fase agua no intercambia masa con las otras fases
Para derivar las ecuaciones del flujo multifsico Tipo Beta Modificado, se parte de un
volumen elemental representativo del medio poroso (figura A1), en donde existe
flujo unidimensional.
( o vox ) x + x
( o vox ) x
Figura A.1. Volumen elemental representativo del medio poroso

Con base en el principio de conservacin de materia, se puede establecer el balance
de materia siguiente en el volumen elemental en la direccin x.
88
Apndice A
Ritmo de entrada de masa al volumen de control:
( Ayz o vox ) .
(A.1)
Ritmo de salida de masa del volumen de control:
( Ayz o vox )
x + x
(A.2)
Donde:
Ayz = rea transversal del elemento en el plano yz expuesto al flujo.
o = Densidad de la fase aceite
vox = Velocidad macroscpica de la fase aceite en la direccin x.

El ritmo de acumulacin de masa en el volumen de control, est dado por el ritmo de
cambio de masa, contenidos en las fases aceite y gas en el espacio poroso, esto es:
Ritmo de acumulacin de masa del elemento en el volumen de control:
Ax
( S o o ) .
t
(A.3)
El ritmo de produccin/inyeccin de masa en el volumen de control, definiendo q p =

gasto volumtrico de la fase p=o,g,w a condiciones de yacimiento por unidad de
volumen de roca, es:
Ax ( o qo ) .
(A.4)
Substituyendo las ecuaciones (A.1) a (A.4) en el postulado de conservacin de

materia, dividiendo la expresin resultante entre el volumen de control, Ax, rearreglando y tomando lmites cuando x 0, se llega a:
( o vox ) + o qo = [ ( o S o )].
x
t
(A.5)
La ecuacion (A.5) puede extenderse al caso de flujo multidimensional. En notacin

del operador diferencial nabla , estas ecuaciones pueden escribirse en forma ms
general como:
89
Apndice A
( o vo ) + o qo =
( o So ) .
t
(A.6)
Para flujo laminar en el yacimiento, podemos expresar la velocidad de cada fase

mediante la ecuacin de Darcy:
vp =
kk rp
p
( p
p D .
(A.7)
Donde: p = o, g, w.
Sustituyendo la ecuacin (A.7) en la ecuacin (A.6), se tiene lo siguiente:
kk
o ro (po o D ) + o qo = [ ( o S o )] .
o
t
(A.8)
Para el caso del aceite se considera que:
@ c. y. =
@ c.s.
Bo
(A.9)
Substituyendo la ecuacin (A.9) en la ecuacin (A.8), se tiene que:
@ c.s. kkro
@ c.s.
@ c.s.

( po oD ) +
qo =
( So ) .
t Bo
Bo o
Bo
Se tiene que
@ c.s. =cte y que bo =
1
; entonces:
Bo
b kk
o ro (p o o D ) + bo q o + omf = [ (bo S o )] .
t
o
(A.10)
(A.11)
Ntese que el operador diferencial nabla , se define en el sistema de coordenadas

cartesianas como:
90
Apndice A
i+
j+ k.
x
y
z
(A.12)
y el vector de velocidades de la fase p = o,g,w,
v p , es :
v p = v px i + v py j + v pz k .
(A.13)
Un tratamiento similar para la fase gas considerando el gas disuelto que se tiene en
el aceite se conduce a lo siguiente:
kk rg
kk
R s bo ro ( p o o D ) + bg
p g g D )
(
o
g

+ R s bo q o + bg q g + R s omf + gmf =
R s bo S o + bg S g .
(A.14)
Un tratamiento similar para la fase agua conduce a lo siguiente:
kk
rw (p w w D ) + bw qw + wmf = [bw S w ].
t
Bw w
(A.15)
Las ecuaciones (A.11), (A.14) y (A.15) describen el flujo multifsico en un
yacimiento bajo condiciones isotrmicas y laminares.
91
Apndice A
92
Apndice B
APNDICE B
APROXIMACIN EN DIFERENCIAS FINITAS
A continuacin se presenta la aproximacin de las ecuaciones diferenciales en
ecuaciones de diferencias, obtenidas a travs de un esquema en diferencias finitas.
Realizando el anlisis para el medio uno y partiendo de la ecuacin (2.20), solo para
el aceite, considerando una sla direccin en el nivel de tiempo (n+1) y empleando
un esquema implcito, se tiene:
n +1
n +1
n +1
kkro po
D
n +1
o
+ bo qo 1, i + [o 21 ]1, i = ( bo So ) 1, i
x o Bo x
x 1,i
t
(B.1)
Considerando,
D
p
u = ox o o
x
x
(B.2)
donde :
ox =
kkro
o Bo
(B.3)
Ahora, sustituyendo la ecuacin (B.2) en el trmino de flujo de la ecuacin (B.1), se

tiene:
n +1
po
D
ox
o
x x
x i
uin +1
xi
(B.4)
Empleando diferencias centrales en la ecuacin (B.4), se tiene:
uin+1
x i
u n+11 u n+11
=
i+
(B.5)
xi
92
Apndice B
donde :
Di +1 Di
n +1 po , i +1 po , i
n +1
u 1 = ( ox )i + 1
1
i+
o,i +
2
x 1
x 1
2
2
+
i
i+
2
2
n +1
(B.6)
Di Di 1
n +1 po , i po , i 1
n +1
u 1 = ( ox )i 1
1
i
o,i
2
x 1
x 1
2
2
i
i
2
2
n +1
(B.7)
Sustituyendo las ecuaciones (B.6) y (B.7) en (B.5), se tiene:
n +1
D D
u
1
n +1 p
=
( ox )i + 1 o,i +1 o,i o,i + 1 i +1 i
2
xi
x 1
xi
2
xi + 12
i+

2
n +1
i
(B.8)
n +1
D D
1
n +1 p
( ox )i 1 o,i o,i 1 o,i 1 i i 1

2
xi
x 1
2
xi 12
i

2
Discretizando el trmino de acumulacin empleando diferencias regresivas en

tiempo, se tiene:
n +1
1
n +1
n
( bo So ) i = n +1 ( bo So ) ( bo So )
i
t
t
(B.9)
Sustituyendo las ecuaciones (B.8) y (B.9) en (B.1), se tiene:
n +1
n +1

p
p
1
D D
D D
1
n +1 p
n +1 p
( ox )i + 1 o,i +1 o,i o,i + 1 i +1 i ( ox )i 1 o,i o,i 1 o,i 1 i i 1
2
2
xi
x 1 xi
x 1
2
2
xi + 12
xi 12
i+
i

2
2
+ bo qo
n +1
i
+ [o 21 ]i
n +1
1
n +1
n
bo So ) ( bo So )
n +1 (
i
t
(B.10)
93
Apndice B
Multiplicando la ecuacin (B.10) por el volumen de roca de la celda i, Vri = Axi , se

tiene:
n +1
ox A
po , i +1 po , i o , i + 1 ( Di +1 Di )
2
xi i + 1
n +1
n +1
ox po , i po , i 1 1 ( Di Di 1 )
o,i
2
xi i 1
n +1
+ bo qo
n +1
i
+ [ o 21 ]i
n +1
Vri
n +1
n
bo So ) ( bo So )
n +1 (
i
t
(B.11)
Introduciendo el concepto de transmisibilidad para la fase aceite:
Tox =
ox A
x i
(B.12)
Las transmisibilidades de la fase gas y aceite se definen de forma similar.

Sustituyendo la ecuacin (B.12) en (B.11), se tiene:
n +1
(To )i + 1 ( po o D )i + 1 (To )i 1 ( po o D )i 1 + bo qo i
n +1
n +1
n +1
n +1
+ [ o 21 ]i
n +1
Vri
t ( bo So ) i
t
(B.13)
La ecuacin (B.13) se puede expresar de una manera ms compacta, empleando

operadores en diferencias, como se muestra a continuacin:
n +1
To ( po o D ) i
+ bo qo
n +1
i
+ [ o 21 ]i
n +1
Vri ,1
t
t ( bo So ) i
(B.14)
Ntese que en (B.14) se ha incluido los operadores diferencias centrales, ui
n+1
diferencias regresivas,
t uin+1,
definidos como:
uin+1 = u n+11 u n+11

i+
, y
(B.15)
t uin+1 = uin+1 uin
(B.16)
94
Apndice B
Expandiendo a tres dimensiones y procediendo en forma similar para las ecuaciones

(2.21) y (2.22), se tienen las siguientes expresiones:
Para el aceite:
Vr
n +1
n +1
n +1
To ( po o D ) 1, ijk + bo qo + [ o 21 ]ijk = 1, ijk t ( bo So ) 1, ijk
1, ijk
t
(B.17)
Para el gas:
RsTo ( po o D ) 1,ijk + Tg ( pg g D ) + Rsbo qo + bg qg + Rs o 21 + g 21

ijk
1, ijk
1, ijk
n +1
n +1
n +1
Vr1, ijk
t
n +1
t ( Rsbo So + bg S g )
1, ijk
(B.18)
Para el agua:
Vr
n +1
n +1
n +1
Tw ( pw w D ) 1, ijk + bw qw + w 21 = 1, ijk t [ bw S w ]1, ijk
1, ijk
ijk
t
(B.19)
Para el medio 2, ecuaciones (2.26) a (2.28), debido a su conceptualizacin,

nicamente se discretiza el trmino de acumulacin, presentando las siguientes
ecuaciones:
Para el aceite:
o 21
n +1
ijk
+ o 32
n +1
ijk
Vr2, ijk
t
t ( bo So ) 2, ijk
(B.20)
Para el gas:
Rs o 21 + g 21
n +1
ijk
+ Rs o 32 + g 32
n +1
ijk
Vr2, ijk
t
95
t ( Rs bo So + bg S g )
2, ijk
(B.21)
Apndice B
96
Apndice B
Para el agua:
[ w 21 ]ijk + [ w32 ]ijk =

n +1
n +1
Vr2,ijk
t
t [ bw S w ]2, ijk
(B.22)
Para el medio 3, ecuaciones (2.32) a (2.34), debido a su conceptualizacin,

nicamente se discretiza el trmino de acumulacin, presentando las siguientes
ecuaciones:
Para el aceite:
o 32
n +1
ijk
Vr3, ijk
t
t ( bo So ) 3, ijk
(B.23)
Para el gas:
Rs o 32 + g 32
n +1
ijk
Vr3, ijk
t
t ( Rs bo So + bg S g )
3, ijk
(B.24)
Para el agua:
[ w32 ]ijk =
n +1
Vr3, ijk
t
t [ bw S w ]3, ijk
Ntese el cambio de
(B.25)
cp ,mf y q p , por p ,mf y q p , respectivamente, donde p = o,
g, w.
96
Apndice C
APNDICE C
ECUACIONES Y DERIVADAS PARA UN MODELO DE TRIPLE POROSIDAD
ACEITE:
ECUACIN PARA EL MEDIO 1:
Trmino
Fuente
Trmino de flujo
Trmino
Transferencia
2-1
Trmino de acumulacin
Vr
n +1
n +1
n +1
To ( po o D ) 1, ijk + bo qo + [ o 21 ]ijk = 1, ijk t bo (1 S g S w )
1, ijk
t
Trmino
Trmino
Transferencia Transferencia
2-1
3-2
[ o 21 ]ijk
n +1
+ [ o 32 ]ijk
n +1
Vr2, ijk
t
t bo (1 S g S w )
Trmino
Transferencia
3-2
[ o 32 ]ijk
n +1
Vr3, ijk
t
t bo (1 S g S w )
3, ijk
97
2, ijk
1, ijk
Apndice C
Donde:
bo kkro
( po 2 po1 ) = To21 ( po 2 po1 )
o 21
o 21 =
bo kkro
( po 3 po 2 ) = To32 ( po 3 po 2 )
o 32
o 32 =
FUNCIN DE RESIDUO DEL ACEITE PARA EL MEDIO 1:
Vr
n +1
n +1
n +1
Fo1,ijk = To ( po o D ) 1, ijk + bo qo + [ o 21 ]ijk 1, ijk t bo (1 S g S w )
1, ijk
t
=0
1, ijk
Expansin de la funcin de residuo del medio 1:
Fo1,ijk
n +1
= To 1
i + 2 , j,k
po
po
( o D )
1
i, j,k
i +1, j , k
i + , j,k
2
To
i 12 , j , k
1
To
i , j + 12 , k
po
po
( o D )
1
i, j,k
i , j +1, k
i, j + ,k
2
To
i , j 12 , k
1
T
o i, j,k + 1
2
p
po
( o D )
i, j,k
1
o i , j , k +1
i, j,k +
T
o
ii , j , k 1
2
n +1
n +1
n +1
n +1
p
po
( o D )
1
i , j 1, k
oi , j , k
i, j ,k
p
po
( o D )
1
i , j , k 1
oi , j , k
i, j ,k
bo qo + To 21 ( po 2 po1 )
ijk
1, ijk
Vp1, ijk
1 + Cr ( pon +1 pon ) bo (1 S w S g )
n +1
b 1 S S
o(
w
g)
)}
n
=0
1, ijk
Derivada de la funcin de residuo del aceite con respecto a la presin de

aceite en el medio 1, en ijk:
Fo1,ijk
po1, ijk
{T.Flujo}
po1, ijk
{T. Fuente}
po1, ijk
{FT21}
po1, ijk
98
n +1
po
po
( o D )
1
i 1, j , k
i, j,k
i , j,k
2
{T.Acum.}
po1, ijk
n +1
n +1
Apndice C
Trminos de flujo:
{T.Flujo}
po1, ijk
= Tox 2
+ Toy 2
+ Toz 2
Potox 2
T
P
T
+ Potox 2 ox 2 Tox1 otox1 Potox1 ox1
po1, ijk
po1, ijk
po1, ijk
po1, ijk
Potoy 2
po1, ijk
+ Potoy 2
Toy 2
po1, ijk
Toy1
Potoy1
po1, ijk
Potoy1
Toy1
po1, ijk
Potoz 2
T
P
T
+ Potoz 2 oz 2 Toz1 otoz1 Potoz1 oz1
po1, ijk
po1, ijk
po1, ijk
po1, ijk
Trminos fuente:
{T. Fuente}
po1, ijk
= boijk
qoijk
po1, ijk
+ qoijk
boijk
po1, ijk
Funcin de transferencia 2-1:
{FT21}
po1, ijk
P
T
= To 21 oto 21 + Poto 21 o 21
po1,ijk
po1,ijk
Trmino de acumulacin:
{T.Acum.}
po1, ijk
Vr

bo (1 S w S g ) + bo (1 S w S g )
= 1, ijk

po1, ijk
t po1,ijk
n +1
ijk
Vp1, ijk
n +1 b
n +1
n
o
+ Cr bo (1 S w S g )
1 + Cr ( po po ) (1 S w S g )
t
po1, ijk
Derivadas de potenciales:
oijk
Potox 2

1
=
po i +1, jk po ijk o i + 1 , jk ( Di +1, jk Di , jk ) = 1 ( Di +1, jk Dijk )
po1, ijk po1, ijk
po1, ijk
2
99
n +1
ijk
Apndice C
oijk
Potox1

1
=
po ijk po i 1, jk o i 1 , jk ( Dijk Di 1, jk ) = 1 ( Dijk Di 1 jk )
po1, ijk po1, ijk
po1, ijk
2
2
Derivadas de las transmisibilidades:
To = F
bo kro
1 bo
1 o
To
= To
b p
po1, ijk
o o1, ijk o po1, ijk
Derivada de la funcin de residuo del aceite con respecto a la saturacin de

gas en el medio 1, en ijk:
Fo1, ijk
S g1, ijk
{T.Flujo}
S g1, ijk
{T. Fuente}
S g1,ijk
{FT21}
S g1, ijk
{T.Acum.}
S g1, ijk
Trminos de flujo:
{T.Flujo}
S g1,ijk
= Potox 2
Toy 2
Toy1
Tox 2
Tox1
Toz 2
Toz1
Potox1
+ Potoy 2
Potoy1
+ Potoz 2
Potoz1
S g1,ijk
S g1,ijk
S g1, ijk
S g1,ijk
S g1, ijk
S g1,ijk
{FT21}
S g1,ijk
To 21ijk
= Poto 21
S g1,ijk
Trminos Fuente:
( bo qo )ijk
n +1
S g1,ijk
= boi , j , k
qoi , j ,k
S g1,ijk
100
Apndice C
Trmino de Acumulacin:
{T.Acum.}
S g1, ijk
=
=
Vr1, ijk
bo (1 S g S w )
t S g1, ijk
ijk
1 + Cr ( pon +1 pon ) ( bon +1 )
ijk
t
Vp1, ijk
Potox 2
=0
S g1, ijk
Potox1
=0
S g1, ijk
Derivadas de transmisibilidades:
To
b kro
=F o
o S g1,ijk
S g1, ijk
agua en el medio 1, en ijk:
Fo1,ijk
S w1, ijk
{T.Flujo}
S w1, ijk
{T. Fuente}
S w1, ijk
{FT21}
S w1, ijk
{T.Acum.}
S w1, ijk
Trminos de flujo:
{T.Flujo}
S w1, ijk
= Potox 2
Toy 2
Toy1
Tox 2
Tox1
Toz 2
Toz1
Potox1
+ Potoy 2
Potoy1
+ Potoz 2
Potoz1
S w1, ijk
S w1, ijk
S w1, ijk
S w1, ijk
S w1,ijk
S w1, ijk
Trminos fuente:
101
Apndice C
{T. Fuente}
S w1, ijk
=0
Trmino de transferencia:
{ FT 21}
S w1, ijk
= Poto 21
To 21
S w1, ijk
{T.Acum.}
S w1, ijk
=
=+
Vr1, ijk
bo (1 S g S w )
t S w1, ijk
ijk
1 + Cr ( pon +1 pon ) bon +1
1,ijk
t
Vp1, ijk
Derivadas de transmisibilidades:
To
b kro
=F o
o S w1,ijk
S w1, ijk
Derivada de la funcin de residuo del aceite en el medio 1, en i-1,j,k:
Fo1
Potox1
Tox1
= Tox1
+ Potox1
po1,i 1 jk
po1,i 1 jk
po1,i 1 jk
Fo1
Tox1
= Potox1
S g1,i 1 jk
S g1,i 1 jk
Fo1
Tox1
= Potox1
S w1,i 1 jk
S w1,i 1 jk
Derivada de la funcin de residuo del aceite en el medio 1, en i,j-1,k:
Potoy1
Toy1
Fo1
= Toy1
+ Potoy1
po1,ij 1, k
po1,ij 1, k
po1,ij 1, k
102
Apndice C
Toy1
Fo1
= Potoy1
S g1,ij 1, k
S g1,ij 1, k
Toy1
Fo1
= Potoy1
S w1ij 1, k
S w1ij 1, k
Derivada de la funcin de residuo aceite en el medio 1, en el nodo i,j,k-1:
Fo1
Potoz1
Toz1
= Toz1
+ Potoz1
po1,ijk 1
po1,ijk 1
po1,ijk 1
Fo1
Toz1
= Potoz1
S g1,ijk 1
S g1,ijk 1
Fo1
Toz1
= Potoz1
S w1,ijk 1
S w1,ijk 1
Derivada de la funcin de residuo del aceite en el medio 1, en el nodo
i+1,j,k:
Fo1
Potox 2
Tox 2
= Tox 2
+ Potox 2
po1,i +1 jk
po1,i +1 jk
po1,i +1 jk
Fo1
Tox 2
= Potox 2
S g1,i +1 jk
S g1,i +1 jk
Fo1
Tox 2
= Potox 2
S w1,ijk
S w1,ijk
Derivada de la funcin de residuo del aceite en el medio 1, en el nodo
i,j+1,k:
Potoy 2
Toy 2
Fo1
= Toy 2
+ Potoy 2
po1,ij +1k
po1,ij +1k
po1,ij +1k
Toy 2
Fo1
= Potoy 2
S g1,ij +1k
S g1,ij +1k
103
Apndice C
Toy 2
Fo1
= Potoy 2
S w1,ij +1k
S w1,ij +1k
Derivada de la funcin de residuo de la fase aceite el medio 1, en el nodo
i,j,k+1:
Fo1
Potoz 2
Toz 2
= Toz 2
+ Potoz 2
po1,ijk +1
po1,ijk +1
po1,ijk +1
Fo1
Toz 2
= Potoz 2
S g1,ijk +1
S g1,ijk +1
Fo1
Toz 2
= Potoz 2
S w1,ijk +1
S w1,ijk +1
n +1
n +1
Fo 2,ijk = To 21 ( po 2 po1 ) ijk + To 32 ( po 3 po 2 ) ijk
Vr2, ijk
t
t bo (1 S g S w )
=0
2, ijk

Fo 2, ijk
po 2, ijk
{FT21}
po 2, ijk
{FT32}
po 2,ijk
{T.Acum.}
po 2, ijk
{FT21}
po 2,ijk
P
T
= To 21 oto 21 + Poto 21 o 21
po 2,ijk
po 2,ijk
{FT32}
po 2,ijk
P
To 32
= To 32 oto 32 + Poto 32
po 2,ijk
po 2,ijk
104
Apndice C
Derivada del trmino de acumulacin con respecto a la presin en el medio 2:
{T.Acum.}
po 2, ijk
Vr

bo (1 S w S g ) + bo (1 S w S g )
= 2, ijk

po 2, ijk
t po 2, ijk
n +1
2, ijk
Vp2, ijk
n +1 b
n +1
n
o
1 + Cr ( po po ) (1 S w S g )
t
po 2, ijk

Fo 2, ijk
S g 2, ijk
{FT21}
S g 2, ijk
{FT32}
S g 2, ijk
{T.Acum.}
S g 2,ijk
{FT21}
S g 2,ijk
To 21ijk
= Poto 21
S g 2, ijk
{FT32}
S g 2,ijk
To 32ijk
= Poto 32
S g 2,ijk
Derivadas del trmino de acumulacin con respecto a la saturacin de gas en el

medio 2:
105
n +1
2,ijk
Apndice C
{T.Acum.}
S g 2, ijk
=
=
Vr2,ijk
bo (1 S g S w )
t S g 2, ijk
2, ijk
1 + Cr ( pon +1 pon ) ( bon +1 )
2, ijk
t
Vp2, ijk

Fo 2, ijk
S w 2, ijk
{FT21}
S w 2, ijk
{FT32}
S w 2, ijk
{T.Acum.}
S w 2, ijk
{ FT 21}
S w 2, ijk
= Poto 21
To 21
S w 2, ijk
{ FT 32}
S w 2, ijk
= Poto 32
To 32
S w 2, ijk
{T.Acum.}
S w 2, ijk
=
=+
Vr2, ijk
bo (1 S g S w )
t S w 2, ijk
Vp2,ijk
t
{1 + C ( p
r
n +1
o
pon ) bonijk+1
2, ijk
2, ijk
n +1
Fo3,ijk = To 32 ( po 3 po 2 ) ijk
Vr3, ijk
t
t bo (1 S g S w )
106
3, ijk
=0
Apndice C

Fo 3, ijk
po 3,ijk
{FT32}
po 3, ijk
{T.Acum.}
po 3, ijk
{FT32}
po 3,ijk
P
To 32
= To 32 oto 32 + Poto 32
po 3,ijk
po 3,ijk
Derivada del trmino de acumulacin con respecto a la presin en el medio 3:
{T.Acum.}
po 3,ijk
Vr3, ijk

=
bo (1 S w S g ) + bo (1 S w S g )

p
t po 3,ijk
o 3, ijk
n +1
3, ijk
Vp3, ijk
n +1 b
n +1
n
o
1 + Cr ( po po ) (1 S w S g )
t
po 3, ijk

Fo 3, ijk
S g 3, ijk
{FT32}
S g 3,ijk
{T.Acum.}
S g 3, ijk
{FT32}
S g 3,ijk
To 32ijk
= Poto 32
S g 3, ijk
Derivadas del trmino de acumulacin con respecto a la saturacin de gas en el

medio 3:
107
n +1
3, ijk
Apndice C
{T.Acum.}
S g 3, ijk
=
=
Vr3, ijk
bo (1 S g S w )
t S g 3, ijk
3, ijk
1 + Cr ( pon +1 pon ) ( bon +1 )
3, ijk
t
Vp3, ijk

Fo 3, ijk
S w3, ijk
{FT32}
S w3, ijk
{T.Acum.}
S w3, ijk
{ FT 32}
S w3, ijk
= Poto 32
To 32
S w3, ijk
{T.Acum.}
S w3, ijk
=
=+
Vr3, ijk
bo (1 S g S w )
t S w3, ijk
3, ijk
Vp3, ijk
t
{1 + C ( p
r
n +1
o
pon ) bonijk+1
3, ijk
Donde:
Derivadas de los potenciales con respecto a la presin en el medio 1:
Poto 21
= 1
p o 1,ijk
Poto 21
=1
po 2,ijk
108
Apndice C
Poto 32
= 1
po 2,ijk
Poto 32
=1
po 3,ijk
Derivadas de las transmisibilidades con respecto a la presin en el medio 1:
Si poto21 0
To 21
=0
po1,ijk
Si poto21 < 0
1 bo 1 o
To 21
= To 21
po1,ijk
bo po o po 1
Si poto21 0
1 bo 1 o
To 21
= To 21
po 2,ijk
bo po o po 2
Si poto21 < 0
To 21
=0
po 2,ijk
Si poto32 0
To 32
=0
po 2,ijk
Si poto32 < 0
1 bo 1 o
To 32
= To 32
po 2,ijk
bo po o po 2
109
Apndice C
Si poto32 0
1 bo 1 o
To 32
= To 32
po 3,ijk
bo po o po 3
Si poto32 < 0
To 32
=0
po 3,ijk
Derivadas de los potenciales con respecto a la saturacin de gas en el medio 1:
Poto 21
=0
S g1,ijk
Poto 21
=0
S g 2,ijk
Poto 32
=0
S g 3,ijk
Derivadas de las transmisibilidades con respecto a la saturacin de gas en el medio
1:
Si poto21 0
To 21
=0
S g1,ijk
Si poto21 < 0
To 21
b kro 2
= F21 o 2
S g 2,ijk
o 2 S g 2,ijk
2:
Si poto21 0
110
Apndice C
To 21
b kro 2
= F21 o 2
S g 2,ijk
o 2 S g 2,ijk
Si poto21 < 0
To 21
=0
S g 2,ijk
Si poto32 0
To 32
=0
S g 2,ijk
Si poto32 < 0
To 32
b kro 2
= F32 o 2
S g 2,ijk
o 2 S g 2,ijk
3:
Si poto32 0
To 32
b kro 3
= F32 o 3
S g 3,ijk
o3 S g 3,ijk
Si poto32 < 0
To 32
=0
S g 3,ijk
Derivadas de los potenciales con respecto a la saturacin de agua en el medio 1:
Poto 21
=0
S w1,ijk
Poto 21
=0
S w 2,ijk
111
Apndice C
Poto 32
=0
S w3,ijk
Derivadas de las transmisibilidades con respecto a la saturacin de agua en el medio
1:
Si poto21 0
To 21
=0
S w1,ijk
Si poto21 < 0
To 21
b kro 2
= F21 o 2
S w 2,ijk
o 2 S w 2,ijk
2:
Si poto21 0
To 21
b kro 2
= F21 o 2
S w 2,ijk
o 2 S w 2,ijk
Si poto21 < 0
To 21
=0
S w 2,ijk
Si poto32 0
To 32
=0
S w 2,ijk
Si poto32 < 0
To 32
b kro 2
= F32 o 2
S w 2,ijk
o 2 S w 2,ijk
3:
Si poto32 0
112
Apndice C
To 32
b kro 3
= F32 o 3
S w3,ijk
o 3 S w3,ijk
Si poto32 < 0
To 32
=0
S w3,ijk
GAS:
Trminos de flujo
Tg ( pg + g D )
n +1
1, ijk
Trmino
Transferencia
2-1
Trminos Fuente
n +1
n +1
n +1
+ RsTo ( po + oD ) 1,ijk + bg qg + Rsbo qo + Rs o 21 + g 21

ijk
1, ijk
=
Vr1,ijk
t
t bo Rs (1 Sg Sw ) + bg Sg
1, ijk
Trmino Transferencia Trmino Transferencia

2-1
3-2
RS o 21 + g 21
n +1
ijk
+ RS o 32 + g 32
n +1
ijk
Vr2, ijk
t
t bo RS (1 S g S w ) + bg S g
113
2ijk
Apndice C
Trmino Transferencia
3-2
RS o 32 + g 32
n +1
ijk
Vr3, ijk
t
t bo RS (1 S g S w ) + bg S g
3ijk
Donde:
bo kkro
( po 2 po1 ) = To21 ( po 2 po1 )
o 21
o 21 =
bg kkrg
( pg 2 pg1 ) = Tg21 ( po 2 po1 + Pcgo 2 Pcgo1 )
g 21
g 21 =
bo kkro
( po 3 po 2 ) = To32 ( po 3 po 2 )
o 32
o 32 =
bg kkrg
( pg 3 pg 2 ) = Tg32 ( po 3 po 2 + Pcgo 3 Pcgo 2 )
g 32
g 32 =
FUNCIN DE RESIDUO DEL GAS PARA EL MEDIO 1:
Fg1,ijk = Tg ( pg + g D )
Vr1, ijk
t
n +1
1, ijk
n +1
n +1
n +1
+ RsTo ( po + o D ) 1, ijk + bg qg + Rs bo qo + g 21 + Rs o 21
1, ijk
ijk
t bg S g + Rsbo (1 S g S w )
1, ijk
=0
114
Apndice C
Fg1,ijk
= Tg 1
i + 2 , jk
po
po + Pcgo
Pcgoijk ( g D )
ijk
i +1, jk
i+1, jk
1
i + , jk
Tg 1 po po
+ Pcgo Pcgo
( g D )
i 1, jk
ijk
i 1, jk
1
i 2 , jk ijk
i , jk
2
n +1
( RsTo )i + 1 , jk poi+1, jk poijk ( o D ) 1

i + , jk
2
n +1
n +1
Tg 1 po po
+ Pcgo Pcgo
( g D )
i , j 1, k
ijk
i , j 1, k
1
i , j 2 ,k ijk
i , j ,k
2
n +1
n +1
( RsTo )i , j 1 , k po po
( o D )
ijk
i , j 1, k
1
i , j ,k
2
n +1
n +1
Tg 1 po po + Pcgo Pcgo ( g D )
ijk
ijk +1
ijk
1
ijk + 2 ijk +1
ijk +
2
Tg 1 po po + Pcgo Pcgoijk 1 ( g D )
ijk 1
ijk
1
ijk 2 ijk
ijk
2
( RsTo ) po po ( o D )
ijk +1
ijk
1
1
ijk +
ijk +
2
2
n +1
Tg 1 po
po + Pcgo
Pcgoijk ( g D )
ijk
i , j +1, k
1
i , j + 2 ,k i , j +1,k
i , j + ,k
2
( RsTo )i , j + 1 , k poi , j +1,k poijk ( o D ) 1

i , j + ,k
2
n +1
( RsTo )i 1 , j , k po po ( o D )
ijk
i 1 jk
1
i , jk
2
n +1
n +1
( RsTo ) po po ( o D )
ijk
ijk 1
1
1
ijk
ijk
2
2
1
n +1
bg qg + Rs bo qo + Tg 21 ( po 2 po1 + Pcgo 2 Pcgo1 ) + RsTo 21 ( po 2 po1 )
ijk
1, ijk
n +1
n +1
Vp1, ijk
1 + Cr ( pon +1 pon ) bg S g + Rs bo (1 S w S g )
t
115
n +1
b S + R b 1 S S
g g
s o(
w
g)
)}
n
1, ijk
=0
Apndice C
Derivada de la funcin de residuo del gas con respecto a la presin de aceite

en el medio 1, en ijk:
Fg1, ijk
po1, ijk
{T.Flujo}
po1, ijk
{T. Fuente}
po1, ijk
{FT21}
po1, ijk
{T.Acum.}
po1,ijk
Trminos de flujo:
{T.Flujo}
po1, ijk
= Tgx 2
Potgx 2
po1, ijk
( RsTo ) x 2
Tgy 2
Potgy 2
po1, ijk
+ Potgy 2
Potoy 2
po1, ijk
Potgz 2
po1, ijk
( RsTo ) z 2
po1, ijk
Tgx1
Potgx1
po1, ijk
Potgx1
Tgx1
po1, ijk
( RsTo ) x 2
( RsTo ) x1
Potox 2
P
+ Potox 2
( RsTo ) x1 otox1 Potox1
+
po1, ijk
po1, ijk
po1, ijk
po1, ijk
( RsTo ) y 2
Tgz 2
Tgx 2
+ Potgx 2
Tgy 2
po1, ijk
+ Potoy 2
+ Potgz 2
Potgy1
Tgy1
po1, ijk
( RsTo ) y 2
po1, ijk
Tgz 2
po1, ijk
Tgz1
Potgy1
( RsTo ) y1
Potgz1
po1, ijk
Potgz1
Tgy1
po1, ijk
Potoy1
po1, ijk
Tgz1
po1, ijk
Potoy1
( RsTo ) y1
po1, ijk
( RsTo ) z 2
( RsTo ) z1
Potoz 2
P
+ Potoz 2
( RsTo ) z1 otoz1 Potoz1
po1, ijk
po1, ijk
po1, ijk
po1, ijk
Trminos fuente:
{T. Fuente}
po1, ijk
= bgijk
qgijk
po1, ijk
+ qgijk
bgijk
po1,ijk
qo
1 boijk
1 Rsijk
ijk
+ ( Rsbo )ijk
+ qoijk
+
bo po1, ijk Rs po1, ijk
po1, ijk
ijk
ijk
{FT21}
po1, ijk
P
T
= Tg 21 otg 21 + Potg 21 g 21
po1, ijk
po1, ijk

( RsTo )21
P
+ ( RsTo )21 oto 21 + Poto 21
po1, ijk
po1, ijk

116
Apndice C
{T.Acum.}
po1, ijk

bg S g + Rsbo (1 S g S w ) + bg S g + Rs bo (1 S g S w )
t

p
t
ijk
o1, ijk
po1, ijk
Vr , ijk
bg
Vp1, ijk
n +1
bo
Rs
n +1
n
=
+ (1 S w S g ) Rs
+ bo
1 + Cr ( po po ) S g
t
po1, ijk
po1,ijk
po1, ijk
Cr bg S g + Rs bo (1 S w S g )
n +1
ijk
gijk

1
po i +1 jk po ijk + Pcgoi +1 jk Pcgoijk + g i + 1 jk ( Di +1 jk Dijk ) = 1 + ( Di +1 jk Dijk )
po1,ijk
po1,ijk
2
gijk

1
po ijk po i 1 jk + Pcgoijk Pcgoi 1 jk + g i 1 jk ( Dijk Di 1 jk ) = 1 + ( Dijk Di 1 jk )
po1,ijk
po1,ijk
2
Potgx 2
po1,ijk
Potgx1
po1,ijk
oijk
Potox 2

1
=
po i +1 jk po ijk + o i + 1 jk ( Di +1 jk Dijk ) = 1 + ( Di +1 jk Dijk )
po1,ijk po1,ijk
po1,ijk
2
2
oijk
Potox1

1
=
po ijk po i 1 jk + o i 1 jk ( Dijk Di 1 jk ) = 1 + ( Dijk Di 1 jk )
po1,ijk po1,ijk
po1,ijk
2
2
Tg = F
bg kkrg
To = F
bo kkro
g
o
1 b
To
1 o
o
= To
bo po1,
po1,ijk
o po1,ijk
ijk
117
n +1
Apndice C
Tg
po1,ijk
1 bg
1 g
= Tg
bg po1,
g po1,ijk
ijk
RsTo
To
Rs
= Rs
+ To
po1,ijk
po1,ijk
po1,ijk
Derivada de la funcin de residuo del gas con respecto a la saturacin de
Fg1, ijk
S g1, ijk
{T.Flujo}
S g1, ijk
{T. Fuente}
S g1,ijk
{FT21}
S g1, ijk
{T.Acum.}
S g1, ijk
Trminos de flujo:
{T.Flujo}
S g1,ijk
Potgx 2
= Tgx 2
Tgy 2
Tgz 2
S g1,ijk
Potgy 2
S g1, ijk
Potgz 2
S g1,ijk
+ Potgx 2
+ Potgy 2
+ Potgz 2
Tgx 2
S g1,ijk
Tgy 2
S g1, ijk
Tgz 2
S g1, ijk
Tgx1
Tgy1
Tgz1
Potgx1
S g1,ijk
Potgy1
S g1,ijk
Potgz1
S g1, ijk
Potgx1
Potgy1
Potgz1
Tgx1
S g1, ijk
Tgy1
S g1, ijk
Tgz1
+ Potox 2
+ Potoy 2
+ Potoz 2
S g1,ijk
( RsTo ) x 2
S g1,ijk
( RsTo ) y 2
S g1,ijk
( RsTo ) z 2
S g1, ijk
{FT21}
S g1, ijk
( RsTo )21,
ijk
= Poto 21
S g1, ijk

Tg 21ijk
P
+ Tg 21ijk otg 21 + Potg 21

S g1,ijk
S g1,ijk

Trminos Fuente:
( bg qg + bo Rs qo )
S g1, ijk
n +1
ijk
= bg ijk
qgijk
S g1, ijk
bg ijk Pcgoijk
+ qgijk
po1, S g1,
ijk
ijk
118
qoijk
+ ( bo Rs )ijk
S g1, ijk
Potox1
Potoy1
Potoz1
( RsTo ) x1
S g1,ijk
( RsTo ) y1
S g1,ijk
( RsTo ) z1
S g1,ijk
Apndice C
{T.Acum.}
S g1, ijk
Vr1, ijk
bg S g + Rsbo (1 S g S w )
t S g1, ijk
ijk

bg Pcgo
Vp1, ijk
n +1
n
Rsbon +1
1 + Cr ( po po ) bg + S g
t
po1,ijk S g1, ijk
ijk
Derivadas de los potenciales con respecto a la saturacin de gas:
Potox 2
=0
S g1, ijk
Potox1
=0
S g1, ijk
Potgx 2
S g1 ijk
Po tg x 1
S g 1 ijk
Pcgo
S g1 ijk
Pc g o
gijk Pcgo
1
Di +1 jk Dijk )
(
2
po1ijk S g1 ijk
1
+ ( D i + 1 jk D ijk
2
S g 1 ijk
1
jk
2
Pc g o
p o 1ijk S g 1 ijk
Derivadas de transmisibilidades con respecto a la saturacin de gas:
Tg
S g1 ijk
= Fk
krg bg
S g1 ijk
1 bg Pcgo
1 krg
1 g Pcgo
+
= Tg
bg po1,ijk S g1 ijk krg S g1 ijk g po1,ijk S g1 ijk
Reescribiendo la ec. anterior, tenemos:

Tg
S g1 ijk
Fbg krg
g S g1
ijk
Tg
Pcgo
po1,ijk S g1 ijk
119
Apndice C
( RsTo )
S g1 ijk
To
S g1 ijk
= Rs
To
Fkbo kro
=
S g1 ijk
o S g1 ijk
Fg1, ijk
S w1, ijk
{T.Flujo}
S w1, ijk
{T. Fuente}
S w1, ijk
{FT21}
S w1, ijk
{T.Acum.}
S w1, ijk
Trminos de flujo:
{T.Flujo}
S w1, ijk
= Potox 2
Potoz 2
( RsTo ) x 2
S w1,ijk
( RsTo ) z 2
S w1, ijk
Potox1
Potoz1
( RsTo ) x1
S w1, ijk
( RsTo ) z1
S w1, ijk
Trminos fuente:
{T. Fuente}
S w1, ijk
=0
Trmino de transferencia:
{ FT 21}
S w1, ijk
= Poto 21
( RsTo )21
S w1, ijk
120
+ Potoy 2
( RsTo ) y 2
S w1, ijk
Potoy1
( RsTo ) y1
S w1, ijk
Apndice C
{T.Acum.}
S w1, ijk
=
=+
Vr1, ijk
t S w1, ijk
ijk
Vp1, ijk
t
{1 + C ( p
r
n +1
o
pon ) ( Rs bon +1 )
ijk
Derivadas de transmisibilidades con respecto a la saturacin de agua:
( RsTo )
S w1, ijk
= Rs
To
S w1, ijk
To
F k bo k ro
=
o S w 1, ijk
S w 1, ijk
Derivada de la funcin de residuo del pseudo componente gas en el nodo i1,j,k:
Fg1
po1,i 1 jk
Fg1
S g1,i 1 jk
Fg1
S w1,i1 jk
( RsTo ) x1
Potgx1
Tgx1
P
= Tgx1
+ Potgx1
+ ( RsTo ) x1 otox1 + Potox1
po1,i1 jk
po1,i1 jk
po1,i1 jk
po1,i 1 jk
= Tgx1
Potgx1
S g1,i 1 jk
= Potox1
Potgx1
Tgx1
S g1,i1 jk
Potox1
( RsTo ) x1
S g1,i 1 jk
( RsTo ) x1
S w1,i1 jk
Derivada de la funcin de residuo del pseudo componente gas en el nodo i,j1,k:
Fg1
po1,ij 1,k
Fg1
S g1ij 1,k
( RsTo ) y1
Potgy1
Tgy1
Potoy1
= Tgy1
+ Potgy1
+ ( RsTo ) y1
+ Potoy1
po1,ij 1,k
po1,ij 1,k
po1,ij 1,k
po1,ij 1,k
= Tgy1
Potgy1
S g1ij 1,k
Potgy1
Tgy1
S g1ij 1,k
Potoy1
( RsTo ) y1
S g1ij 1,k
121
Apndice C
Fg1
S w1,ij 1k
= Potoy1
( RsTo ) y1
S w1,ij 1k
Derivada de la funcin de residuo del pseudo componente gas en el nodo

i,j,k-1:
Fg1
po1,ijk 1
Fg1
S g1,ijk 1
Fg1
S w1ij 1k
( RsTo ) z1
Potgz1
Tgz1
P
= Tgz1
+ Potgz1
+ ( RsTo ) z1 otoz1 + Potoz1
po1,ijk 1
po1,ijk 1
po1,ijk 1
po1,ijk 1
= Tgz1
Potgz1
S g1,ijk 1
= Potoy1
Potgz1
Tgz1
S g1,ijk 1
Potoz1
( RsTo ) z1
S g1,ijk 1
( RsTo ) y1
S w1ij 1k
Derivada de la funcin de residuo del pseudo componente gas con respecto

a la presin del aceite en el nodo i+1,j,k
Fg1
po1,i+1, jk
Fg1
S g1,i+1 jk
Fg1
S w1i+1, jk
= Tgx 2
= Tgx 2
Potgx 2
po1,i+1, jk
Potgx 2
S g1,i+1 jk
= Potox 2
+ Potgx 2
+ Potgx 2
Tgx 2
po1,i+1, jk
Tgx 2
S g1,i+1 jk
+ ( RsTo ) x 2
+ Potox 2
( RsTo ) x 2
Potox 2
+ Potox 2
po1,i+1, jk
po1,i+1, jk
( RsTo ) x 2
S g1,i+1 jk
( RsTo ) x 2
S w1i +1, jk

a la presin del aceite en el nodo i,j+1,k
Fg1
po1,ij +1,k
= Tgy 2
Potgy 2
po1,ij +1,k
+ Potgy 2
Tgy 2
po1,ij +1,k
+ ( RsTo ) y 2
122
Potoy 2
po1,ij +1,k
+ Potoy 2
( RsTo ) y 2
po1,ij +1,k
Apndice C
Fg1
S g1,ij +1,k
Fg1
S w1,ij +1,k
Potgy 2
= Tgy 2
S g1,ij +1,k
= Potoy 2
+ Potgy 2
Tgy 2
S g1,ij +1,k
+ Potoy 2
( RsTo ) y 2
S g1,ij +1,k
( RsTo ) y 2
S w1,ij +1,k

a la presin del aceite en el nodo i,j,k+1
Fg1
po1,ijk +1
Fg1
S g1,ijk +1
Fg1
S w1,ijk +1
= Tgz 2
= Tgz 2
Potgz 2
po1,ijk +1
Potgz 2
S g1,ijk +1
= Potoz 2
+ Potgz 2
Tgz 2
po1,ijk +1
+ Potgz 2
Tgz 2
S g1,ijk +1
+ ( RsTo ) z 2
+ Potoz 2
( RsTo ) z 2
Potoz 2
+ Potoz 2
po1,ijk +1
po1,ijk +1
( RsTo ) z 2
S g1,ijk +1
( RsTo ) z 2
S w1,ijk +1
Fg 2,ijk = Tg21 ( po 2 po1 + Pcgo 2 Pcgo1 ) + ( RsTo )21 ( po 2 po1 )

ijk
n +1
+ Tg32 ( po 3 po 2 + Pcgo 3 Pcgo 2 ) + ( RsTo )32 ( po 3 po 2 )
Vr2, ijk
t
2, ijk
n +1
ijk
=0

Fg 2, ijk
po 2, ijk
{FT21}
po 2,ijk
{FT32}
po 2, ijk
{T.Acum.}
po 2, ijk
123
Apndice C
{FT21}
po 2,ijk
P
T
( RsTo )21
P
= Tg 21 otg 21 + Potg 21 g 21 + ( RsTo )21 oto 21 + Poto 21
po 2, ijk
po 2,ijk
po 2, ijk
po 2, ijk
{FT32}
po 2, ijk
= Tg 32
Potg 32
po 2, ijk
+ Potg 32
Tg 32
po 2, ijk
+ ( RsTo )32
( RsTo )32
Poto 32
+ Poto 32
po 2,ijk
po 2, ijk
{T.Acum.}
po 2, ijk
bg S g + Rs bo (1 S g S w ) + bg S g + Rs bo (1 S g S w )
t
t
po 2, ijk
po 2, ijk
ijk
Vr2, ijk
bg
Vp
n +1
bo
Rs
= 2, ijk 1 + Cr ( pon +1 pon ) S g
+ (1 S w S g ) Rs
+ bo
t
po 2, ijk
po 2, ijk
po 2,ijk
n +1
ijk

Fg 2, ijk
S g 2, ijk
{FT21}
S g 2, ijk
{FT32}
S g 2, ijk
{T.Acum.}
S g 2,ijk
{FT21}
S g 2,ijk
( RsTo )21,
ijk
= Poto 21
S g 2, ijk

Tg 21ijk
P

S g 2, ijk
S g 2, ijk

124
n +1
Apndice C
{FT32}
S g 2, ijk
( RsTo )32,
ijk
= Poto 32
S g 2, ijk

Tg 32ijk
P

S g 2, ijk
S g 2, ijk

{T.Acum.}
S g 2, ijk
Vr2, ijk
t S g 2, ijk
ijk
Vp2, ijk
bg Pcgo
n +1
n
n +1
C
p
p
b
S
R
b
1
r( o
o )
g
s o
g
p
S
t
o
g
2,
2,
ijk
ijk
ijk

Fg 2, ijk
S w 2, ijk
{FT21}
S w 2, ijk
{FT32}
S w 2, ijk
{T.Acum.}
S w 2, ijk
Trmino de transferencia 2-1:
{ FT 21}
S w 2, ijk
( RsTo )21
= Poto 21
S w 2,ijk
{ FT 32}
S w 2, ijk
= Poto 32
( RsTo )32
S w 2, ijk
125
Apndice C
{T.Acum.}
S w 2, ijk
Vr2, ijk
t S w 2, ijk
ijk
=+
Vp2, ijk
t
{1 + C ( p
n +1
o
pon ) ( Rsbon +1 )
ijk
Fg 3,ijk = Tg32 ( po 3 po 2 + Pcgo 3 Pcgo 2 ) + ( RsTo )32 ( po 3 po 2 )

ijk
n +1
Vr3, ijk
t
3, ijk
=0

Fg 3, ijk
po 3, ijk
{FT32}
po 3, ijk
{T.Acum.}
po 3, ijk
{FT32}
po 3, ijk
= Tg 32
Potg 32
po 3,ijk
+ Potg 32
Tg 32
po 3, ijk
+ ( RsTo )32
126
( RsTo )32
Poto 32
+ Poto 32
po 3, ijk
po 3, ijk
Apndice C
{T.Acum.}
po 3,ijk

bg S g + Rs bo (1 S g S w ) + bg S g + Rs bo (1 S g S w )
t

p
t
po 3,ijk
o 3, ijk
ijk
Vr3, ijk
bg
Vp
n +1
bo
Rs
= 3, ijk 1 + Cr ( pon +1 pon ) S g
+ (1 S w S g ) Rs
+ bo
t
po 3,ijk
po 3, ijk
po 3, ijk
n +1
ijk

Fg 3, ijk
S g 3, ijk
{FT32}
S g 3,ijk
{T.Acum.}
S g 3, ijk
{FT32}
S g 3, ijk
( RsTo )32,
ijk
= Poto 32
S g 3, ijk

Tg 32ijk
P

S g 3, ijk
S g 3, ijk

{T.Acum.}
S g 3, ijk
Vr3,ijk
1
b
S
R
b
S
S
+
(
)
g g
s o
g
w
t S g 3, ijk
ijk
Vp3, ijk
bg Pcgo
n +1
n
n +1
C
p
p
b
S
R
b
1
+
(
)
r
o
o
g
s o
g
po 3,ijk S g 3, ijk
t
ijk

Fg 3, ijk
S w3, ijk
{FT32}
S w3, ijk
{T.Acum.}
S w3, ijk
127
n +1
Apndice C
{ FT 32}
S w3, ijk
= Poto 32
( RsTo )32
S w3, ijk
{T.Acum.}
S w3, ijk
=
=+
Vr3, ijk
t S w3, ijk
ijk
Vp3, ijk
t
{1 + C ( p
r
n +1
o
pon ) ( Rsbon +1 )
ijk
Donde:
Potg 21
po1,ijk
= 1
Poto 21
= 1
po1ijk
Potg 21
p o 2 ,ijk
=1
Poto 21
=1
po 2ijk
Potg 32
p o 2,ijk
= 1
Poto 32
= 1
po 2ijk
128
Apndice C
Potg 32
=1
p o 3,ijk
Poto 32
=1
po 3ijk
Si poto21 0 y potg21 0
Tg 21
po1ijk
=0
( RsTo )21
po1ijk
=0
Si poto21 < 0 y potg21 < 0
Tg 21
po1,ijk
1 bg1
1 g1
= Tg 21
bg1 po1,
g1 po1,ijk
ijk
( RsTo )21
po1,ijk
= Rs 21
To 21
R
+ To 21 s 21
po1,ijk
po1,ijk

Tg 21
po 2,ijk
1 bg 2
1 g 2
= Tg 21
bg 2 po 2,
g 2 po 2,ijk
ijk
( RsTo )21
po 2,ijk
= Rs 21
To 21
Rs 21
+ To 21
po 2,ijk
po 2,ijk
1 b
To 21
1 o 2
o2
= To 21
bo 2 po 2,
po 2,ijk
o 2 po 2,ijk
ijk
129
Apndice C
Tg 21
po 2,ijk
=0
( RsTo )21
po 2,ijk
=0
Tg 32
po 2,ijk
=0
( RsTo )32
po 2,ijk
=0
Tg 32
po 2,ijk
1 bg 2
1 g 2
= Tg 32
bg 2 po 2,
g 2 po 2,ijk
ijk
( RsTo )32
po 2,ijk
= Rs 32
To 32
Rs 32
+ To 32
po 2,ijk
po 2,ijk
1 b
To 32
1 o 2
o2
= To 32
bo 2 po 2,
po 2,ijk
o 2 po 2,ijk
ijk

Tg 32
po 3,ijk
1 bg 3
1 g 3
= Tg 32
bg 3 po 3,
g 3 po3,ijk
ijk
( RsTo )32
po 3,ijk
= Rs 32
To 32
Rs 32
+ To 32
po 3,ijk
po 3,ijk
1 b
To 32
1 o 3
o3
= To 32
bo 3 po 3,
po 3,ijk
o 3 po3,ijk
ijk
130
Apndice C
Tg 32
po 3,ijk
=0
( RsTo )32
=0
po 3,ijk
Potg 21
S g 1,ijk
Pcgo1
S g 1,ijk
Poto 21
=0
S g1,ijk
Poto 21
=0
S g 2,ijk
Potg 21
S g 2,ijk
Pcgo 2
S g 2,ijk
Poto 32
=0
S g 2,ijk
Potg 32
S g 2,ijk
Pcgo 2
S g 2,ijk
Poto 32
=0
S g 3,ijk
Potg 32
S g 3,ijk
Pcgo 3
S g 3,ijk
Derivadas de las transmisibilidades con respecto a la saturacin de gas en el medio 1:

131
Apndice C
Tg 21
S g1,ijk
=0
( RsTo )21
S g1,ijk
=0
Tg 21
S g1,ijk
= F21
kbg1 krg1 Tg 21 pcgo1
g1 S g1, po1, S g1,

ijk
( RsTo )21
S g1,ijk
= Rs 21
ijk
ijk
To 21
S g1,ijk
To 21
kb kro
= F21 o1
S g1,ijk
o1 S g1,ijk
Tg 21
S g 2,ijk
= F21
kbg 2 krg 2 Tg 21 pcgo 2
g 2 S g 2, po 2, S g 2,
ijk
( RsTo )21
S g 2,ijk
= Rs 21
ijk
ijk
To 21
S g 2,ijk
To 21
kb kro 2
= F21 o 2
S g 2,ijk
o 2 S g 2,ijk
Tg 21
S g 2,ijk
=0
( RsTo )21
S g 2,ijk
=0
132
Apndice C
Tg 32
S g 2,ijk
=0
( RsTo )32
S g 2,ijk
=0
Tg 32
S g 2,ijk
= F32
g 2 S g 2, po 2, S g 2,
ijk
( RsTo )32
S g 2,ijk
= Rs 32
ijk
ijk
To 32
S g 2,ijk
To 32
kb kro 2
= F32 o 2
S g 2,ijk
o 2 S g 2,ijk
Tg 32
S g 3,ijk
= F32
g 3 S g 3, po 3, S g 3,
ijk
( RsTo )32
S g 3,ijk
= Rs 32
ijk
ijk
To 32
S g 3,ijk
To 32
kb kro 3
= F32 o 3
S g 3,ijk
o 3 S g 3,ijk
Tg 32
S g 3,ijk
=0
133
Apndice C
( RsTo )32
S g 3,ijk
=0
AGUA:
Trmino
Fuente
Trmino de flujo
Trmino
Transferencia
2-1
n +1
Tw ( pw + w D ) 1, ijk + bw qw + [ w 21 ]ijk
1, ijk
n +1
n +1
Vr1, ijk
t
t [ bw S w ]1, ijk
Trmino de
acumulacin
Trmino
Trmino
Transferenci Transferencia
a
3-2
[ w 21 ]ijk + [ w32 ]ijk =

n +1
n +1
Trmino de
acumulacin
Vr 2, ijk
t
t [ bw S w ]2, ijk
134
Apndice C
Trmino
Transferenci
a
Trmino de
acumulacin
[ w32 ]ijk =
n +1
Vr 3, ijk
t
t [ bw S w ]3, ijk
Donde:
bw kkrw
( pw 2 pw1 ) = Tw21 ( po 2 po1 Pcwo 2 + Pcwo1 )
w 21
w 21 =
bw kkrw
( pw3 pw 2 ) = Tw32 ( po 3 po 2 Pcwo 3 + Pcwo 2 )
w 32
w32 =
FUNCIN DE RESIDUO DEL AGUA PARA EL MEDIO 1:
V
n +1
n +1
n +1
Fw1ijk = Tw ( pw wD ) 1, ijk + bw qw + [ w 21 ]ijk r1, ijk t [ bw S w ]1, ijk
1, ijk
t
135
Apndice C
Fw1,ijk
= Tw 1 po
po Pcwo
+ Pcwo ( wD )
i +1, jk
ijk
i +1, jk
ijk
1
i + , jk
i+ 2 , jk
2
Tw 1 poijk poi1, jk Pcwoijk + Pcwoi1, jk ( w D ) 1

i , jk
i , jk
2
2
n +1
+
1
n +1
Tw 1 poi , j +1,k poijk Pcwoi , j +1,k + Pcwoijk ( w D ) 1

i , j + ,k
i , j + 2 ,k
2
Tw 1 poijk poi , j 1,k Pcwoijk + Pcwoi , j 1,k ( w D ) 1

i , j ,k
i , j ,k
2
+
1
n +1
n +1
Tw 1 poijk +1 poijk Pcwoijk +1 + Pcwoijk ( w D ) 1

ijk +
ijk + 2
2
Tw 1 poijk poijk 1 Pcwoijk + Pcwoijk 1 ( w D ) 1

ijk
ijk
2
2
+ bwqw
n +1
1, ijk
+ [ w 21 ]ijk
n +1
Vp1,ijk
t
{1 + C ( p
r
n +1
o
n +1
n +1
n +1
n
pon ) ( bw S w ) ( bw S w )
=0
1, ijk
Derivada de la funcin de residuo del agua con respecto a la presin de

Fw1, ijk
po1, ijk
{T.Flujo}
po1, ijk
{T. Fuente}
po1, ijk
{FT21}
po1, ijk
{T.Acum.}
po1, ijk
Trminos de flujo:
{T.Flujo}
P
T
P
T
= Twx 2 otwx 2 + Potwx 2 wx 2 Twx1 otwx1 + Potwx1 wx1
po1, ijk
po1, ijk
po1, ijk
po1, ijk
po1, ijk
Potwy 2
T
P
T
+ Potwy 2 wy 2 Twy1 otwy1 + Potwy1 wy1 +
Twy 2
po1, ijk
po1, ijk
po1, ijk
po1, ijk
P
T
P
T
Twz 2 otwz 2 + Potwz 2 wz 2 Twz1 otwz1 + Potwz1 wz1
po1, ijk
po1, ijk
po1, ijk
po1, ijk
136
Apndice C
Trmino fuente:
{T. Fuente}
po1, ijk
= bw
qw ijk
po1, ijk
+ qw ijk
bw ijk
po1, ijk
{FT21}
po1, ijk
= Tw 21
Potw 21
T
+ Potw 21 w 21
po1, ijk
po1, ijk
{T.Acum.}
po1, ijk
=
=
[ bw Sw ]ijk
t po1, ijk
Vr , ijk

Vpijk
bw
n +1
n
1
S
C
p
p
b
C
+
+
w
r( o
o )
w r
t
po1, ijk

ijk
po i+1, jk po ijk Pcwoi+1, jk + Pcwoijk wi + 1 , jk ( Di +1, jk Dijk )

2
wijk
1
= 1 ( Di +1, jk Dijk )
2
po1,ijk
Potwx 2
=
po1,ijk po1,ijk
po ijk po i1, jk Pcwoijk + Pcwoi1, jk wi 1 , jk ( Dijk Di 1 jk )
2
wi , j ,k
1
= 1 ( Dijk Di 1, jk )
po1,ijk
2
Potwx1
=
po1,ijk po1,ijk
137
Apndice C
Tw = F
bw krw
1 bw
1 w
Tw
= Tw
b p
po1, ijk
w o1, ijk w po1, ijk
Derivada de la funcin de residuo del agua con respecto a la saturacin de

Fw1, ijk
S g1, ijk
{T.Flujo}
S g1, ijk
{T. Fuente}
S g1, ijk
{FT21}
S g1, ijk
{T.Acum.}
S g1, ijk
Trminos de flujo:
{T.Flujo}
S g1, ijk
=0
Funcin de transferencia matriz fractura (FTMF):
{FT21}
S g1, ijk
=0
Trmino Fuente:
{T. Fuente}
S g1, ijk
=0
{T.Acum.}
S g1, ijk
=0

138
Apndice C
Fw1, ijk
S w1, ijk
{T.Flujo}
S w1, ijk
{T. Fuente}
S w1, ijk
{FT21}
S w1, ijk
{T.Acum.}
S w1, ijk
Trminos de flujo:
{T.Flujo}
S w1, ijk
P
T
P
Twx1
= Twx 2 otwx 2 + Potwx 2 wx 2 Twx1 otwx1 + Potwx1
+
S
w
1,
ijk
w
1,
ijk
w
1,
ijk
w
1,
ijk
Twy1
Potwy 2
T
P
+ Potwy 2 wy 2 Twy1 otwy1 + Potwy1
Twy 2
+
S w1, ijk
S w1, ijk
S w1, ijk
S w1, ijk
P
Twz 2
P
Twz1
Twz 2 otwz 2 + Potwz 2
Twz1 otwz1 + Potwz1
S
S
S
S
w
1,
ijk
w
1,
ijk
w
1,
ijk
w
1,
ijk
{FT21}
S w1,ijk
P
T
= Tw 21 otw 21 + Potw 21 w 21
S w1,ijk
S w1,ijk
Trmino fuente:
{T. Fuente}
S w1, ijk
= bw
qw
bw Pc wo
qw
S w1, ijk
po1,ijk S w1,ijk
{T.Acum.}
S w1, ijk
Vr1, ijk
( bw Sw )
t S w1, ijk
ijk
Vp1, ijk
bw Pcow
n +1
n
1 + Cr ( po po ) bw S wijk
t
po S w1, ijk

ijk
Derivadas de potenciales con respecto a la saturacin de agua:
139
Apndice C
Potwx 2
=
S w1,ijk S w1,ijk
=
Pcwoijk
S w1,ijk
po i+1, jk po ijk Pcwoi+1, jk + Pcwoijk wi + 1 , jk ( Di +1, jk Dijk )

2
wijk Pcwoijk
1
+ ( Di +1, jk Dijk )
2
po1,ijk S w1,ijk
Potwx1

=
po ijk po i1, jk Pcwoijk + Pcwoi1, jk wi 1 , jk ( Dijk Di 1, jk )
S w1,ijk S w1,ijk
2
Pcwoijk 1
wijk Pcwoijk
=
+ ( Dijk Di 1, jk )
S w1,ijk 2
po1,ijk S w1,ijk
Derivadas de transmisibilidades con respecto a la saturacin de agua:
1 k
Tw
1 w Pcwo
1 bw Pcwo
rw
= Tw
+
S w1,ijk
krg S w1,ijk w po1,ijk S w1,ijk bw po1,ijk S w1,ijk
Derivada de la funcin de residuo del agua en el nodo i-1,j,k:
Fw1
P
Twx1
= Twx1 otwx1 + Potwx1
po1,i 1 jk
po1,i1 jk
po1,i1 jk
Fw1
P
Twx1
= Twx1 otwx1 + Potwx1
S w1,i 1 jk
S w1,i 1 jk
S w1,i1 jk
Derivada de la funcin de residuo del agua en el nodo i,j-1,k:
P
Twy1
Fw1
= Twy1 otwy1 + Potwy1
po1,ij 1k
po1,ij 1k
po1,ij 1k
P
Twy1
Fw1
= Twy1 otwy1 + Potwy1
S w1,ij 1k
S w1,ij 1k
S w1,ij 1k
Derivada de la funcin de residuo del agua en el nodo i,j,k-1:
140
Apndice C
Fw1
P
T
= Twz1 otwz1 + Potwz1 wz1
poijk 1
poijk 1
poijk 1
Fw1
P
T
= Twz1 otwz1 + Potwz1 wz1
S wijk 1
S wijk 1
S wijk 1
Derivada de la funcin de residuo del pseudo componente aceite en el nodo

i+1,j,k:
Fw1
P
Twx 2
= Twx 2 otwx 2 + Potwx 2
po1,i+1 jk
po1,i+1 jk
po1,i+1 jk
Fw1
P
Twx 2
= Twx 2 otwx 2 + Potwx 2
S w1,i +1 jk
S w1,i+1 jk
S w1,i+1 jk

i,j+1,k:
P
T
Fw1
= Twy 2 otwy 2 + Potwy 2 wy 2
poij +1k
poij +1k
poij +1k
P
Twy 2
Fw1
= Twy 2 otwy 2 + Potwy 2
S w1,ij +1k
S w1,ij +1k
S w1,ij +1k

i,j,k+1:
Fw1
P
Twz 2
= Twz 2 otwz 2 + Potwz 2

po1,ijk +1
po1,ijk +1
po1,ijk +1
Fw1
P
Twz 2
= Twz 2 otwz 2 + Potwz 2

S w1,ijk +1
S w1,ijk +1
S w1,ijk +1

n +1
n +1
Fw 2 = Tw21 ( po 2 po1 Pcwo 2 + Pcwo1 ) + Tw32 ( po 3 po 2 Pcwo 3 + Pcwo 2 )

ijk
ijk
Vr 2, ijk
t
t [ bw S w ]2, ijk = 0

141
Apndice C
Fw 2, ijk
po 2, ijk
{FT21}
po 2, ijk
{FT32}
po 2, ijk
{T.ACUM.}
po 2, ijk
Funcin de Transferencia 2-1:
{FT21}
po 2, ijk
P
T
= Tw 21 otw 21 + Potw 21 w 21
po 2, ijk
po 2,ijk
{FT32}
= Tw32
po 2, ijk
Potw32
T
+ Potw32 w32
po 2, ijk
po 2, ijk
{T.ACUM.}
po 2, ijk
=
=
Vr 2, ijk
t po 2, ijk
[bw Sw ]2,ijk

Vp2, ijk
bw
n +1
n
+ bwCr
S w 1 + Cr ( po po )
t
po 2, ijk
2,ijk

Fw 2, ijk
S g 2, ijk
{FT21}
S g 2, ijk
{FT32}
S g 2, ijk
{T.ACUM.}
S g 2, ijk
{FT21}
S g 2, ijk
=0
{FT32}
S g 2,ijk
=0
142
Apndice C
{T.ACUM.}
S g 2, ijk
=0

Fw 2, ijk
S w 2, ijk
{ FT 21}
S w 2, ijk
{ FT 32}
S w 2, ijk
{T . ACUM .}
S w 2, ijk
{FT21}
S w 2, ijk
P
Tw 21
= Tw 21 otw 21 + Potw 21
S w 2, ijk
S w 2,ijk
{FT32}
S w 2,ijk
= Tw32
Potw32
T
+ Potw32 w32
S w 2,ijk
S w 2,ijk
V
{T . ACUM .}
= r ,ijk
S w,ijk
t
=
m
( bw S w )
S wijk
ijk
bwm Pcwom
n +1
n
1 + Cr ( pom
p
)
om bwmijk S wmijk
t
pomijk S wmijk
V p ijk
V
n +1
Fw3 = Tw32 ( po 3 po 2 Pcwo 3 + Pcwo 2 ) r 3,ijk t [ bw S w ]3, ijk = 0
ijk
t

Fw3, ijk
po 3, ijk
{FT32}
po 3, ijk
{T.ACUM.}
po 3, ijk
143
Apndice C
{FT32}
P
T
= Tw32 otw32 + Potw32 w32

po 3, ijk
po 3, ijk
po 3, ijk
{T.ACUM.}
po 3, ijk
=
=
[bw Sw ]3,ijk
t po 3, ijk
Vr 3, ijk

Vp3, ijk
bw
n +1
n
+ bwCr
S w 1 + Cr ( po po )
t
po 3, ijk
3,
ijk

Fw3, ijk
S g 3, ijk
{FT32}
S g 3, ijk
{T.ACUM.}
S g 3,ijk
{FT32}
=0
S g 3, ijk
{T.ACUM.}
S g 3, ijk
=0

Fw3, ijk
S w3, ijk
{ FT 32}
S w3, ijk
{T . ACUM .}
S w3, ijk
144
Apndice C
{FT32}
S w3,ijk
P
T
= Tw32 otw32 + Potw32 w32

S w3,ijk
S w3,ijk
{T . ACUM .}
S w3, ijk
Vr 3, ijk
( bw S w )
t S w3, ijk
3, ijk
V p 3, ijk
bw Pcwo
n +1
n
C
p
p
1
(
)
r
o
o bw S w
t
po 3,ijk S w3, ijk
Donde:
Derivadas del potencial con respecto a la presin en el medio 1:
Potw 21
= 1
p o 1,ijk
Potw 21
=1
po 2,ijk
Potw32
= 1
po 2,ijk
Potw32
=1
po 3,ijk
Si potw21 0
Tw 21
=0
po1,ijk
Si poto21 < 0
145
3, ijk
Apndice C
1 b
1 w
Tw 21
w
= Tw 21
bw po1,
w po1,ijk
po1,ijk
ijk

Si potw21 0
1 b
1 w
Tw 21
w
= Tw 21
bw po 2,
w po 2,ijk
po 2,ijk
ijk
Si potw21 < 0
Tw 21
=0
po 2,ijk
Si potw32 0
Tw32
=0
po 2,ijk
Si potw32 < 0
1 b
1 w
Tw32
w
= Tw32
bw po 3,
w po 3,ijk
po 2,ijk
ijk

Si potw32 0
1 b
1 w
Tw32
w
= Tw32
w po3,ijk
po 3,ijk
b p
w o 3,ijk
Si potw32 < 0
Tw32
=0
po 3,ijk
146
Apndice C
Potw 21 Pcwo1
=
S w1,ijk S w1,ijk
Potw 21
P
= cwo 2
S w 2,ijk
S w 2,ijk
Potw32 Pcwo 2
=
S w 2,ijk S w 2,ijk
Potw32
P
= cwo 3
S w3,ijk
S w3,ijk
Derivadas de las transmisibilidades con respecto a la saturacin de agua en el medio 1:
Si potw21 0
Tw 21
=0
S w1,ijk
Si potw21 < 0
Tw 21
=F
S w1, ijk
S w1, ijk
krwbw
P
Tw 21
bw krw
cwo
=F
w S w1,ijk S w1,ijk po1,ijk
w

Si potw21 0
Tw 21
krwbw
Pcwo Tw 21
bw krw
=F
=F
S w 2, ijk
S w 2, ijk w
w S w 2,ijk S w 2,ijk po 2,ijk
Si potw21 < 0
Tw 21
=0
S w 2, ijk
Si potw32 0
147
Apndice C
Tw32
=0
S w 2, ijk
Si potw32 < 0
Tw32
krwbw
Pcwo Tw32
bw krw
=F
=F
S w 2, ijk
S w 2, ijk w
w S w 2,ijk S w 2,ijk po 2,ijk
Si potw32 0
Tw32
krwbw
Pcwo Tw32
bw krw
=F
=F
S w3, ijk
S w3, ijk w
w S w3,ijk S w3,ijk po3,ijk
Si potw32 < 0
Tw32
=0
S w3, ijk
148

Simulacion Triple Porosidad

Cargado por

Copyright:

Formatos disponibles

Simulacion Triple Porosidad

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Simulacion Triple Porosidad

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL AUTONOMA DE MEXICO

PROGRAMA DE MAESTRIA Y DOCTORADO EN

SIMULACIN DE YACIMIENTOS NATURALMENTE

ERIKA FABIOLA CORTES RUBIO

DR. VCTOR HUGO ARANA

LISTA DE TABLAS................................................................................................. iii

CAPTULO 2. FORMULACIN MATEMTICA DEL PROBLEMA.....

2.1 Consideraciones Geolgicas.........................................................................

2.2 Modelo Conceptual y Matemtico....

2.2.1 Modelo Conceptual.... 16

2.3 Condiciones Iniciales y de Frontera....

2.3.1 Condiciones Iniciales.......

CAPTULO 3. SOLUCIN NUMRICA DEL PROBLEMA....................................

Ecuaciones de Flujo en Diferencias Finitas..................................

3.1.1. Ecuaciones para el Medio 1................................................

3.2 Solucin del Sistema de Ecuaciones.................................................

3.2.1 Mtodo Iterativo de Newton-Raphson................................... 32

CAPTULO 4. VALIDACIN DEL MODELO.......

4.1 Solucin para el Flujo Radial Monofsico en un Medio....................

5.1 Caso 1. Simulacin de un Sistema de Triple Porosidad por

APNDICE A. OBTENCIN DE LAS ECUACIONES DE FLUJO

APNDICE B. APROXIMACIN EN DIFERENCIAS FINITAS....

APNDICE C. ECUACIONES Y DERIVADAS PARA UN MODELO DE

Aunque el modelado de los yacimientos naturalmente fracturados ha sido desarrollado

Tabla 2.1 Posibles Configuraciones de un sistema de Triple Porosidad.

Figura 2.1 a) Espacio Poroso Consolidado (matriz), b) Representacin

Figura 5.6 Grfica de RGA contra Tiempo para la Simulacin de Flujo

La modelacin de Yacimientos Naturalmente Fracturados (YNF) ha sido, desde la

En la simulacin de los yacimientos fracturados, nos encontramos con el hecho de que

En 1985, Pruess y Narasimhan 25 desarrollaron un mtodo de Multiple Interaccin

En 2006 Kang, Z. y cols 16, presentaron e implementaron un modelo conceptual

3) Implementar el modelo propuesto en un simulador trifsico, tridimensional, de

En el modelo propuesto para ser implementado en el simulador, se considera un

En una primera parte se presenta una introduccin sobre el desarrollo y

En una segunda parte se establece la formulacin matemtica del problema y

En una cuarta parte se valida el modelo con soluciones analticas.

En una quinta parte se utiliza el simulador en casos prcticos.

En una sexta parte se realizan las conclusiones del presente trabajo.

Finalmente cabe mencionar que aunque implcitamente el algoritmo presente algunas

Formulacin Matemtica del Problema

2.1 Consideraciones Geolgicas

La porosidad secundaria generalmente se encuentra en rocas compactas, quebradizas

Formulacin Matemtica del Problema

estn presentes, la permeabilidad puede ser mucho ms alta y las propiedades

Esta porosidad total es equivalente a la definicin esttica de almacenamiento de roca

Formulacin Matemtica del Problema

En el caso de yacimientos heterogneos, adems de la porosidad primaria de la matriz,

Compresibilidad Total (roca y fluido)

Formulacin Matemtica del Problema

2.1.3 Presin Capilar

Curvas de Presin Capilar

Formulacin Matemtica del Problema

Saturacin de la Fase Mojante

Formulacin Matemtica del Problema

Formulacin Matemtica del Problema

donde las permeabilidades efectivas,

kox y kwx dependen de la saturacin. La

permeabilidad efectiva del aceite, puede expresarse como:

k x =permeabilidad absoluta del medio poroso en la direccin