Simulacion Streamline

SIMULACIN STREAMLINE: ESTADO DEL ARTE Y APLICACIONES
FUNDAMENTALES EN EL MODELAMIENTO Y ESTUDIO DE YACIMIENTOS

ALTAMENTE HETEROGNEOS
ADRIANA MARCELA MNDEZ BOHRQUEZ

HERNANDO ABRIL PREZ
UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE INGENIERAS FISICO-QUMICAS
ESCUELA DE INGENIERA DE PETRLEOS
BUCARAMANGA
2005
SIMULACIN STREAMLINE: ESTADO DEL ARTE Y APLICACIONES

ALTAMENTE HETEROGNEOS
ADRIANA MARCELA MNDEZ BOHRQUEZ

HERNANDO ABRIL PREZ
Tesis para optar al ttulo de

Ingeniero de Petrleos
Director
ING. OLGA PATRICIA ORTZ CANCINO
UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE INGENIERAS FISICO-QUMICAS
ESCUELA DE INGENIERA DE PETRLEOS
BUCARAMANGA
2005
DEDICATORIA
A Dios por concederme la vida, por sus gracias, bendiciones y ser la fuerza que
me impulsa a alcanzar mis ideales.
A mi madre, por el don de la vida, por su incansable entrega, amor y esfuerzo,

ya que con su ayuda y consejo hoy puedo disfrutar de este logro. Gracias
Mami.
A mi hermano Vctor por su compaa y cario y por ensearme a su corta

edad que lo ms importante es tener sueos.
A Evangelista, por su apoyo y ayuda en los duros y gratos momentos y por la

plena confianza que ha depositado en m.
A toda mi familia, por los ms cercanos y los ms distantes, porque siempre

encontrar una mano amiga en cada uno de ellos.
A mis luces celestiales, quienes desde el ms all me aportaron su ayuda

especial en todo instante.
A mis amigos, compaeros de estudios y aventuras, a ellos gracias por su

amistad.
ADRIANA
DEDICATORIA
A mi madre porque me dio el regalo ms preciado que es la vida.
A mis hermanos Jos Angel, Arnulfo y Henry por su apoyo y ayuda en los
momentos ms difciles.
A mis tos Hector Melndez y Samuel Gmez por su ayuda y colaboracin ya

que sin ellos este sueo de estar en una universidad y culminar nunca se
hubiese llevado a cabo.
A Dios, que en vez de darme un padre me dio dos. A ellos muchas gracias por
sus bendiciones.
A mis tas Hilda y Alcira por su ayuda y enseanza que cada vez son ms tiles
para mi vida.
A Marlene por su amor y apoyo incondicional sin esperar nada a cambio.
A mis primos en especial a Joan Sebastin y Eduardo, por aceptarme como un

hermano ms.
A mi compaera, por su fuerza de voluntad, ya que sin ella este trabajo no se

hubiese llevado a cabo. Gracias Marce.
HERNANDO
AGRADECIMIENTOS
Los autores expresan sus agradecimientos a:
La ingeniera Olga Patricia Ortiz Cancino, por su invaluable aporte, direccin y
colaboracin en el desarrollo de este proyecto.
Los ingenieros Nicols Santos y Fernando Calvete, por la atencin y

orientacin brindada.
Los ingenieros Edwin Javier Carrero, William Navarro y Roy Elicer Sandoval
por su aporte intelectual y ayuda incondicional en el manejo de los simuladores
empleados.
A la Escuela de Ingeniera de Petrleos, en cabeza de su directora, Dra. Zuly

Himelda Caldern, por su amplia colaboracin y disposicin de los medios
necesarios requeridos para el desarrollo del proyecto.
A los ingenieros Jorge Mantilla V. y Jos Arnobio Vargas Medina, quienes

desde la distancia nos aportaron un conocimiento ms detallado del progreso
de esta nueva tecnologa en nuestro pas.
CONTENIDO
pg.
INTRODUCCION
22
1. GENERALIDADES
23
1.1 PLANTEAMIENTO DEL PROBLEMA
23
1.2 DESPLAZAMIENTOS DOMINADOS POR CONVECCIN
23
1.3 SIMULACIN STREAMLINE BASADA EN EL FLUJO
24
1.3.1 Contexto histrico.
24
1.3.2 Formulacin general de la simulacin streamline.
26
1.3.3 Simulacin streamline vs. simulacin en diferencias finitas.
26
1.3.4 Ventajas de la simulacin streamline.
28
1.3.5 Desventajas de la simulacin streamline.
33
2. SIMULACIN STREAMTUBES: PRIMEROS PLANTEAMIENTOS

HACIA LA SIMULACIN STREAMLINE
34
2.1 NO LINEALIDAD
35
2.1.1 La aproximacin de Riemann en la formulacin streamtubes.
35
2.2 CLASES DE PROBLEMAS
36
2.3 MODELO MATEMTICO
36
2.3.1 Condiciones lmites.
37
2.3.2 El papel de la streamfuncin en simulacin streamtubes.
39
2.4 TRAZADO DE SOLUCIONES 1D
42
2.4.1 Mtodo A.
42
2.4.2 Mtodo B.
44
2.4.3 Trazado de las soluciones 1D sobre el grid cartesiano 2D.
44
2.5 SOLUCIONES EN DOS DIMENSIONES
46
3. CONCEPTOS Y MODELO MATEMTICO DE LA SIMULACIN

STREAMLINE.
47
3.1 DEDUCCIN DE LAS ECUACIONES
47
3.2 SOLUCIN A LA ECUACIN DE PRESIN
50
3.3 CONDICIONES LMITES
51
3.4 SOLUCIN NUMRICA DE LA MATRIZ DE PRESIN
52
3.5 DETERMINACIN DE LA VELOCIDAD DEL CAMPO
53
3.6 MTODO PARA LA DETERMINACIN DE LAS VAS Y EL

TRAZADO DE LAS Streamlines
53
3.7 EL TIEMPO DE VUELO
58
3.8 TRANSFORMACIN DE LA COORDENADA A LO LARGO

DE LAS STREAMLINES
58
3.9 SOLUCIN ANALTICA DE LA ECUACIN DE CONTINUIDAD
60
3.9.1 Clculo de las propiedades en cada bloque del grid de simulacin
61
3.9.2 Bloques que no son atravezados por streamlines.
63
3.9.3 Clculo del tiempo real.
63
3.10 SOLUCIN NUMRICA DE LA ECUACIN DE CONTINUIDAD
65
3.10.1 Planteamiento de las nuevas soluciones numricas.
66
3.10.2 Bloques que no son atravesados por streamlines.
70
3.10.3 Errores en el balance de materiales.
71
3.11 ACTUALIZACIN DEL MODELO DE LNEAS DE FLUJO
71
3.12 PASOS DE TIEMPO
73
3.13 EL CONCEPTO DE GRAVEDAD EN SIMULACIN STREAMLINE
74
3.14 FLUJO INCOMPRESIBLE Y COMPRESIBLE
76
3.15 DESCRIPCIN DEL MTODO STREAMLINE
78
4. APLICACIONES DE LA SIMULACIN STREAMLINE EN LA

EVALUACIN DE YACIMIENTOS
81
4.1 CLCULOS DEL VOLUMEN DE BARRIDO
81
4.2 CLASIFICACIN DE LOS MODELOS GEOESTADSTICOS
86
4.3 APLICACIN DE LA SIMULACIN STREAMLINE EN EL

PROCESO DE UPSCALING
89
4.4 TASA DE DISTRIBUCIN Y BALANCE DE PATRONES
94
4.4.1 Balance de patrones.
95
4.4.2 Eficiencia del inyector.
99
4.5 MODELAMIENTO DE TRAZADORES DE FLUJO Y PROCESOS

DE INYECCIN DE AGUA
101
4.6 EL PAPEL DE LA SIMULACIN STREAMLINE EN EL PROCESO

DE AJUSTE HISTRICO
102
4.6.1 Aproximaciones streamline para ajuste histrico.
104
5. SIMULADORES COMERCIALES STREAMLINE
123
5.1 SIMULADOR 3DSL STREAMSIM TECHNOLOGIES.INC.
123
5.1.1 Caractersticas bsicas del simulador 3DSL.
124
5.2 SIMULADOR FRONTSIM. ECLIPSE- SCHLUMBERGER
126
5.2.1 Caractersticas bsicas de Frontsim.
126
5.2.2 Estructura de una simulacin Frontsim.
128
6. ULTIMOS AVANCES EN SIMULACIN STREAMLINE
130
6.1 YACIMIENTOS NATURALMENTE FRACTURADOS
130
6.1.1 Formulacin matemtica modelo streamline en YNF.
130
6.1.2 Soluciones numricas a las ecuaciones de saturacin.
134
6.2 NUEVOS MTODOS DESARROLLADOS PARA MANEJAR FLUJO

COMPRESIBLE EN SIMULACIN STREAMLINE
139
6.2.1 Mtodo secuencial.
140
6.2.2 Mtodo implcito.
142
6.3 SISTEMAS DE ACEITE NEGRO CON EFECTOS CAPILARES
143
6.3.1 Modificaciones a la ecuacin de presin.
143
6.3.2 Modificaciones a la ecuacin de saturacin.
144
6.4 SISTEMAS COMPOSICIONALES EN SIMULACIN STREAMLINE
146
6.4.1 Ecuaciones bsicas.
146
6.4.2 Soluciones unidimensionales a las ecuaciones de flujo

composicional.
148
7. PLANTEAMIENTO DE DOS CASOS BASE
150
7.1 MODELO MULTIPOZO. PRIMER CASO BASE
150
7.1.1 Condiciones de desarrollo del modelo.
150
7.1.2 Resultados de la simulacin.
151
7.1.3 Simulador Frontsim vs. simulador Eclipse 100.
164
7.2 MODELO HETEROGNEO. SEGUNDO CASO BASE
169
169
8. PERSPECTIVA DE LA SIMULACIN STREAMLINE EN COLOMBIA
180
9. CONCLUSIONES Y RECOMENDACIONES
183
9.1 CONCLUSIONES
183
9.2 RECOMENDACIONES
185
BIBLIOGRAFA
186
LISTA DE TABLAS
pg.
Tabla 1. Tiempo de cmputo del simulador Frontsim para cada
paso de tiempo empleado.
163
Tabla 2. Produccin obtenida por FRONTSIM para diversos

pasos de tiempo.
164
Tabla 3. Tiempo de simulacin empleado en cada simulador.
168
LISTA DE FIGURAS
pg.
Figura 1. reas de aplicacin de la simulacin streamline vs.
simulacin en diferencias finitas.
27
Figura 2. Representacin del movimiento de los fluidos en el

modelo streamline y en el modelo de diferencias finitas.
29
Figura 3. Visualizacin de las streamlines en el modelamiento

de flujo de fluidos.
32
Figura 4. Descripcin grfica del modelo streamtube.
35
Figura 5. Condiciones lmites. Modelo matemtico streamtube.
39
Figura 6. Streamlines y la streamfuncin.
39
Figura 7. Mtodo de integracin simple a lo largo de los

streamtube para determinar el volumen poroso acumulativo
en el Punto X.
43
Figura 8. Representacin del mtodo B.
44
Figura 9. Trazado de la solucin unidimensional a lo largo de

los streamtubes sobre el grid cartesiano fundamental.
45
Figura 10. Regiones de flujo rpidas y lentas.
46
Figura 11. Mtodo de Pollock para el trazado de las streamlines

conociendo el punto de entrada.
54
Figura 12. Trazado de las lneas de flujo.
56
Figura 13. Celdas no ortogonales y el mtodo de Pollock.
57
Figura 14. Forma grfica del mtodo empleado para calcular la

saturacin en el grid de simulacin en cada celda.
63
Figura 15. Discretizacin de la ecuacin de continuidad en

direccin del tiempo de vuelo para nodos centrados.
67
Figura 16. Discretizacin de la ecuacin de continuidad en

direccin del tiempo de vuelo para nodos distribuidos.
68
Figura 17. Trazado del frente de saturacin para la

actualizacin de las lneas de flujo.
73
Figura 18. Influencia la gravedad sobre los vectores velocidad

de las fases.
75
Figura 19. Cambio de las geometras streamline debido a

cambios en la movilidad y en las condiciones de los pozos.
79
Figura 20. Diagrama de flujo para un simulador streamline.
80
Figura 21. Distribucin del tiempo de vuelo y reas de barrido

a diferentes tiempos para modelos de cinco puntos.
84
Figura 22. Clculo del rea de barrido basado en el tiempo de

vuelo streamline para patrones de cinco puntos homogneos.
85
Figura 23. Clculo del rea de barrido basado en el tiempo de

vuelo streamline para patrones de cinco puntos heterogneos.
85
Figura 24. Ranking basado en el porcentaje de recobro vs.

ranking basado en la eficiencia volumtrica de barrido a la ruptura.
88
Figura 25. Ranking basado en el porcentaje de recobro vs. ranking

basado en la eficiencia volumtrica de barrido momentos
despus de la ruptura.
88
Figura 26. Ranking basado en el porcentaje de recobro vs. ranking

basado en la eficiencia volumtrica de barrido momentos
despus de la ruptura.
89
Figura 27. Resultado de una tcnica de escalamiento.
90
Figura 28. Streamlines para dos modelos escalados y el modelo

de referencia a fina escala.
93
Figura 29. Comparacin del volumen asociado con los pozos entre
el modelo a fina escala y modelos robustos mediante streamlines.
94
Figura 30. Diagramas streamlines y factores de distribucin de

pozos en el yacimiento.
95
Figura 31. Streamlines y el proceso de balance de patrones.
96
Figura 32. Influencia del tiempo de vuelo sobre las zonas barridas.
97
Figura 33. Influencia del rea ligada a 7 inyectores.
97
Figura 34. Volumen de influencia cuando existen 6 pozos

productores verticales.
98
Figura 35. Volumen de influencia cuando existen 6 pozos

horizontales productores.
99
Figura 36. Agua inyectada vs, desplazamiento de aceite

producido.
100
Figura 37. Saturacin de aceite promedio por pozo vs aceite

Producido.
101
Figura 38. Ajuste de la curva de flujo fraccional.
110
Figura 39. Diagrama de flujo para el mtodo propuesto por

Mickaele Le Ravalec Dupin y Darryl Fenwick.
117
Figura 40: Ilustracin de la tcnica de renormalizacin

basada en la caracterizacin de Dykstra-Parsons.
119
Figura 41. Resultados del ajuste histrico para el corte de agua

empleando un mtodo convencional y el mtodo de AHA.
121
Figura 42. Bosquejo del entorno del simulador 3DSL.
124
Figura 43. Representaciones grficas construidas a partir de

los datos arrojados por el simulador 3DSL.
125
Figura 44. Bosquejo del entorno del simulador FRONTSIM.
129
Figura 45. Efecto de la gravedad en un bloque de la matriz

rodeado por fracturas parcialmente llenas de agua.
135
Figura 46. Esquema del grid para el primer caso base a simular.
150
Figura 47. Esquemas de distribucin de presin para el primer

caso base.
151
Figura 48. Grficas de presin variando el nmero de pasos de

tiempo en Frontsim. Primer caso base.
152
Figura 49. Streamlines determinadas en el primer caso base a las

condiciones lmites iniciales.
152
Figura 50. Patrn de lneas inicial en el primer caso base.
153
Figura 51. Primer ajuste a las tasas de inyeccin.
153
Figura 52. Ultimo ajuste conseguido a las tasas de inyeccin.
154
Figura 53. Factor de recobro ajustando las tasas de inyeccin.
154
Figura 54. Balance estableciendo un nuevo pozo inyector.
155
Figura 55. Aumento en el factor de recobro con la creacin de

un nuevo pozo inyector.
155
Figura 56. Distribucin de las streamlines despus de cambiar

las condiciones lmites. Primer caso base.
156
Figura 57. Factor de recobro. Primer caso base.
156
Figura 58. Tasa de produccin de aceite por pozo para

P1, P2, P3, P4 y P5.
157
Figura 59. Corte de agua por pozo para P1, P2, P3, P4 y P5.
157
Figura 60. Tasa de produccin de aceite por pozo para

P6, P7, P8, P9 y P10.
158
158
Figura 62. Presin con y sin efectos de gravedad. Primer

caso base.
159
Figura 63. OIP con y sin efectos de gravedad para el primer

caso base.
160
Figura 64. Produccin total de agua con y sin efectos de gravedad

para el caso base uno.
160
Figura 65. Produccin total de aceite con y sin efectos de gravedad.

Primer caso base.
160
Figura 66. Distribucin de saturaciones halladas por el simulador

Frontsim para el primer caso base.
161
Figura 67. Porcentaje de error en el balance de materiales.
162
Figura 68. Anlisis del porcentaje de error en el balance de

materiales con cambios en las condiciones lmites.
163
Figura 69. Grfica de presin para los simuladores ECLIPSE 100 y

FRONTSIM.
164
Figura 70. Comparacin de la distribucin de presin obtenida

por los dos simuladores para el mismo tiempo de simulacin.
165
Figura 71. Distribucin de la saturacin obtenida por los dos

simuladores para el mismo tiempo de simulacin.
166
Figura 72. Tasa de produccin de aceite obtenida para los dos

simuladores.
166
Figura 73. Tasa de inyeccin de agua para los dos simuladores

considerados.
167
Figura 74. Produccin total de agua obtenida por Eclipse 100

y Frontsim.
167
Figura 75. Produccin total de aceite obtenida por Eclipse 100

y Frontsim.
167
Figura 76. Distribucin de permeabilidad para el modelo

heterogneo considerado.
169
Figura 77. Distribucin de presin en el modelo heterogneo.
170
Figura 78. Presin para cada uno de los pasos de tiempo empleados.
Modelo heterogneo.
170
Figura 79. Distribucin de streamlines para el modelo heterogneo.
171
Figura 80. Tiempo de vuelo obtenido para el modelo establecido

al cabo de los 3650 das de inyeccin.
172
Figura 81. Distribucin de saturacin con y sin gravedad. Modelo

heterogneo.
172
Figura 82. Distribucin de presiones hallada por el simulador streamline

con y sin efectos de gravedad para el caso base dos.
173
Figura 83. FOIP para el caso dos con y sin efectos de gravedad.
173
Figura 84. Tasa de produccin de aceite para el caso dos con y sin
efectos de gravedad.
174
Figura 85. Produccin total de aceite para el caso dos con y sin
174
Figura 86. Inyeccin total de agua para el caso dos con y sin
174
Figura 87. Produccin total de agua para el caso dos con y sin
175
Figura 88. Saturacin para un solo paso de tiempo considerado.

Caso heterogneo.
175
Figura 89. Saturacin para 4 pasos de tiempo considerados.

Caso base dos.
176
Figura 90. Saturacin para 8 pasos de tiempo considerados.

Caso base dos.
176
Figura 91. Errores en el balance de materiales tomando

diferentes tamaos en los pasos de tiempo. Modelo dos.
177
Figura 92. Corte de agua tomando diferentes tamaos en los

pasos de tiempo. Modelo heterogneo.
177
Figura 93. Inyeccin total de agua tomando diferentes tamaos

en los pasos de tiempo. Modelo heterogneo.
178
Figura 94. OIP tomando diferentes tamaos en los pasos

de tiempo. Modelo heterogneo.
178
Figura 95. Produccin total de aceite tomando diferentes tamaos

178
Figura 96. Tasa de produccin de aceite tomando diferentes

tamaos en los pasos de tiempo. Modelo heterogneo.
179
Figura 97. Produccin total de agua tomando diferentes tamaos

179
RESUMEN
TITULO: SIMULACIN STREAMLINE: ESTADO DEL ARTE Y APLICACIONES

ALTAMENTE HETEROGNEOS*
AUTORES:
ADRIANA MARCELA MNDEZ BOHRQUEZ.
HERNANDO ABRIL PREZ. **
PALABRAS CLAVES: SIMULACIN STREAMLINE, UPSCALING, AJUSTE

HISTRICO, PROCESOS GEOESTADSTICOS.
El objetivo de esta tesis es suministrar una fuente de informacin para el

conocimiento y comprensin de la simulacin streamline, empleada en los ltimos
aos en el modelamiento y simulacin de yacimientos que por sus caractersticas
de extensin, complejidad y heterogeneidad, requieren de herramientas ms
fuertes, veloces y sofisticadas en el manejo de fluidos multifsicos. El
planteamiento de esta nueva tecnologa se inicia con las generalidades referentes
al tema, en las que se destacan los tipos de desplazamiento preferidos por el
mtodo, un contexto histrico, una comparacin entre la simulacin streamline y
las tcnicas convencionales, y las ventajas y desventajas del mtodo.
En este trabajo se exponen los argumentos matemticos y conceptuales que

perfilan esta nueva formulacin. Seguidamente, se presentan las aplicaciones
desarrolladas en esta tcnica de simulacin en el desarrollo y evaluacin de
procesos como son el ajuste histrico, upscaling y validacin de procesos
geoestadsticos entre otros. Se presentan los simuladores comerciales ms
conocidos y empleados sobre simulacin streamline, desarrollando mediante uno
de ellos un ejemplo prctico de validacin y anlisis del proceso de simulacin.
Finalmente se describen los ltimos avances desarrollados para el tratamiento de
sistemas ms complejos como fracturados, composicionales, compresibles y
altamente influenciados por presin capilar.
Los simuladores streamline han demostrado ser herramientas poderosas en la

caracterizacin de yacimientos, especialmente cuando se trabajan con modelos
dominados por desplazamientos convectivos. Sus ventajas son consecuencia del
mtodo aplicado y de las suposiciones inherentes a este.
____________
* Tesis.
**
Facultad de Ingenieras Fisico-Qumicas. Escuela de Ingeniera de Petrleos.
Director: Ing.Olga Patricia Ortiz Cancino.
ABSTRACT
TITLE: STREAMLINE SIMULATION: STATE OF THE ART AND FUNDAMENTAL

APPLICATIONS IN THE MODELLING AND STUDY OF HIGHLY
HETEROGENEOUS RESERVOIRS. *
AUTHORS:
ADRIANA MARCELA MNDEZ BOHRQUEZ.
HERNANDO ABRIL PREZ. **
KEYWORDS: STREAMLINE SIMULATION, UPSCALING, HISTORY MATCH,
GEOESTATISTICAL PROCESS.
The objective of this work is to give a source of information for the knowledge and
understanding of the streamline simulation, employee in the last years in the
modelling and simulation of reservoirs that for its extension characteristics,
complexity and heterogeneity, require of more strong, speedy and sophisticated
tools in the handling of multiphase fluids. The position of this new technology
begins with the relating generalities to the topic, in those that stand out the favorite
displacement types for the method, a historical context, a comparison among the
streamline simulation and the conventional techniques, and the advantages and
disadvantages of the method.
In this work the mathematical and conceptual definitions are exposed.

Subsequently, are shown the applications developed in this simulation technique in
the development and evaluation of process like the history match, upscaling and
validation of geoestatistical process among others. Next are exposed commercial
simulators well known and employees in streamline simulation, developing by
means of one of them a practical example of validation and analysis of the
simulation process. Finally the last advances developed are described for the
treatment of more complex systems as having fractured, compositionals, compress
and highly influenced by capillary pressure.
The streamline simulators have demonstrated to be powerful tools in the

characterization of reservoirs, especially when one works with models dominated
by convective displacements. Their advantages are consequence of the applied
method and of the inherent suppositions to this.
____________
* Thesis.
**
School of Physical-Chemical Engineering. School of Petroleum Engineering.
Director: Ing. Olga Patricia Ortiz Cancino.
INTRODUCCION
Los procesos de caracterizacin de yacimientos requieren en la actualidad de

la compresin, entendimiento e interpretacin de todos los fenmenos estticos
y dinmicos que pueden estar ocurriendo en un yacimiento. En el pasado estos
estudios se fundamentaban en la informacin que podan brindar solo los
parmetros geolgicos y estticos y se buscaba a partir de estos dar una
explicacin a un comportamiento que indudablemente es de carcter dinmico,
obtenindose como resultado modelos con un gran ndice de falencias e
incertidumbres.
En la actualidad el panorama de la ingeniera de yacimientos se ha

encaminado hacia el desarrollo de herramientas de simulacin que aparte de
brindar las garantas que ofrecen los simuladores convencionales, permitan
interactuar aspectos estticos y dinmicos del yacimiento, logrando avanzar
cada da hacia la bsqueda y determinacin de modelos de yacimientos ms
completos y reales que involucran la consideracin de aspectos complejos
tanto en lo que se refieren al tamao, geologa y composicin. La simulacin
streamline es justamente una de estas alternativas. Su fortaleza ha sido
demostrada a travs de los ltimos aos mediante estudios y aplicaciones de
campo, sin dejar de mencionar los esfuerzos y expectativas que se tienen en el
futuro en cuanto a la posibilidad de expansin y profundizacin en aspectos
que hasta el momento no se consideran viables.
Este proyecto contempla a travs de sus captulos los aspectos tericos,

conceptuales, campos de aplicacin y presentacin de herramientas
necesarias para llevar a cabo una simulacin streamline, esperando como
consecuencia dar una iniciativa para el desarrollo de futuras investigaciones.
22
1. GENERALIDADES
1.1 PLANTEAMIENTO DEL PROBLEMA

Los mtodos de simulacin streamtubes y streamlines han sido usados y
aceptados ampliamente en los ltimos aos, logrando con sus avances un
creciente inters por parte de toda la comunidad petrolera a nivel mundial, lo
cual se debe principalmente a:
Con el desarrollo en la caracterizacin de yacimientos, se han

generado modelos de alta resolucin constituidos por millones de
celdas. Este avance ha ocasionado inconvenientes entre el
modelamiento geolgico y el proceso de simulacin.
Con el incremento en la resolucin de los modelos, as mismo ha

aumentado el grado de incertidumbre que rodea stos modelos.
Es por eso que mediante nuevas tecnologas se busca valorar
esta incertidumbre y predecir el comportamiento de un yacimiento
a partir de modelos ms precisos y realistas.
Los simuladores numricos convencionales son inadecuados para satisfacer

estas necesidades y especialmente no son lo suficientemente rpidos y
eficaces para manejar stos modelos.
Los modelos streamline han cambiado dramticamente el trabajo de simulacin
y se han hecho mejores predicciones.
Existen varias tecnologas que pueden mejorar substancialmente los
pronsticos de simulacin y reducir los porcentajes de error. Algunas de estas
nuevas tecnologas son:
Streamline basada en la simulacin del flujo

Tcnicas Asistidas de History Matching (AHM)
Ssmica 4D
Clculo en Paralelo
Integracin de la simulacin de flujo con efectos geo-mecnicos.
1.2 DESPLAZAMIENTOS DOMINADOS POR CONVECCIN

Los desplazamientos dominados por conveccin se refieren a desplazamientos
que tienen una dependencia de primer orden de la permeabilidad y de la
movilidad total, y en un segundo orden o pequea dependencia de la presin
(compresibilidad), difusin, dispersin y presin capilar. Ejemplos de
desplazamientos dominados por conveccin son los mecanismos de recobro
secundario y terciario tales como inyeccin de agua o desplazamientos tipo
WAG.
23
1.3 SIMULACIN STREAMLINE BASADA EN EL FLUJO.
1.3.1 Contexto histrico. En la actualidad la simulacin streamline ha sido

precedida por varios mtodos para el modelamiento del flujo dominado por
conveccin en el yacimiento.
Los mtodos de lnea fuente han sido ampliamente usados por la industria del
petrleo. Estos mtodos usan soluciones analticas a la distribucin de presin
y velocidad en el yacimiento. La limitacin primaria de estos mtodos es que se
requieren propiedades homogneas y espesores del yacimiento constantes.
Los mtodos de rastreo de partculas tambin han predominado ya que

permiten modelar el comportamiento de trazadores en yacimientos de
hidrocarburos. Estos mtodos rastrean el movimiento de un grupo significante
de partculas a lo largo de caminos apropiados, trabajando de forma excelente
en frentes escarpados, pero no en perfiles uniformes.
Los mtodos streamtubes son ms generales y han sido aplicados
exitosamente en modelos escalados de campo principalmente en procesos
miscibles y de inyeccin de agua. Un estudio detallado sobre est tecnologa y
sus alcances es presentado por Thiele*. En estos mtodos, el dominio de
inters es dividido dentro de un nmero de streamtubes y los clculos de
saturacin de fluidos son realizados a lo largo de stos streamtubes.
Lamentablemente este tipo de simulacin ha sido limitada a sistemas de dos
dimensiones o a alguna forma de aproximacin hbrida para tener en cuenta
efectos 3D.
Las tecnologas streamtube y streamline, en gran magnitud, han sido

manejadas por el hecho de que la heterogeneidad controla los factores de
recobro de muchos campos. Esta evolucin caus la derivacin de modelos
geolgicos ms complejos, pero desgraciadamente tambin ocasion un gran
problema entre el detalle geolgico y la capacidad de simulacin.
La tecnologa streamtube fue originalmente desarrollada en los aos sesenta.
Inicialmente se generaron dos modelos adimensionales streamtube para
trabajar con patrones de flujo regular de permeabilidad homognea, como por
ejemplo un modelo de cinco puntos. Los modelos streamtube fueron ms tarde
generados para posiciones irregulares de pozo y para yacimientos arealmente
heterogneos. Los viejos modelos streamtube solo permitan trabajar con tasas
y posiciones constantes (es decir, no se podan incluir ni cerrar pozos y sus
tasas tanto de inyeccin como de produccin deban permanecer constantes.).
Los efectos de gravedad y por consiguiente la eficiencia de barrido vertical
nunca fue considerada.
____________
*
THIELE, Marco R. Modeling Multiphase Flow in Heterogeneous Media Using Streamtubes. Stanford
University, Dic 1994.
24
Como resultado de estas limitaciones los modelos iniciales streamtube solo

permitan evaluar la eficiencia de barrido areal de un solo patrn.
Para calcular la eficiencia del barrido vertical, Chevron desarroll los modelos
hbridos streamtube que permitan evaluar en forma eficiente el barrido areal y
vertical en el yacimiento. En esta tcnica la seccin transversal vertical era
simulada inicialmente y posteriormente combinada con los modelos areales
streamtube ya desarrollados. Sin embargo, la introduccin de nuevos pozos y
los grandes cambios en las tasas de produccin-inyeccin significaban que la
geometra streamtube deba cambiar originndose limitaciones con esta
tcnica. Es por eso que la tecnologa streamline es en el momento prctica en
muchas aplicaciones de campo debido a que esta incluye:
Gravedad
Efectos 3D
Cambio en las condiciones del pozo.
Flujo Multifsico
La tecnologa streamline incluye efectos de gravedad y permite cambios en la

tasa de los pozos. Esto les permite a ingenieros realizar un proceso sencillo
que evala tanto la eficiencia areal como la vertical y permite cambios en los
pozos. Desde entonces han surgido una explosin de estudios, resaltndose el
trabajo dirigido por Thiele*, quien es uno de los que ha resaltado la utilidad de
la simulacin streamline.
La actual popularidad de la simulacin streamline se debe ms a un

resurgimiento oportuno, dado que las streamlines y su definicin han
permanecido en la bibliografa desde la publicacin de los artculos tcnicos de
Muskat** en 1937 y ha recibido continua atencin desde entonces.
La simulacin streamline descansa sobre seis principios fundamentales:
Trazado de las streamlines tridimensionales (3D) en trminos del

tiempo de vuelo.
Replanteamiento de las ecuaciones de conservacin de masa a lo
largo de las streamlines en trminos del tiempo de vuelo.
____________
*
THIELE R, Marco. Streamline Simulation. 6 International Forum on Reservoir Simulation. Schloss

Fusch, Austria. 3-7 Sep. 2001.
**
MUSKAT, M. Flow of Homogeneous Fluids. International Human Resources Development Corporation.
Boston, MA. 1937.
25
Actualizacin peridica de las streamlines.

Soluciones numricas de transporte 1D a lo largo de las streamlines.
Clculo de los efectos de gravedad.
Extensin a flujo compresible.
1.3.2 Formulacin general de la simulacin streamline. Los simuladores

streamline aproximan clculos de flujo de fluidos en sistemas 3D mediante la
suma de soluciones 1D a lo largo de las streamlines. La seleccin de las
direcciones de las streamlines para los clculos en una dimensin hace de la
aproximacin extremadamente efectiva para modelar procesos dominados por
conveccin en el yacimiento. Este es el caso en el que la heterogeneidad es el
factor predominante que gobierna el comportamiento del flujo.
Segn Datta-Gupta*, un concepto fundamental en simulacin streamline es la
independencia de los clculos del flujo de aspectos como son los efectos de
heterogeneidad geolgica. Matemticamente se logra esta independencia por
el uso de lo que se conoce como tiempo de vuelo. El planteamiento expresa un
sistema de coordenadas donde todas las streamlines son lneas rectas y la
distancia es reemplazada por el tiempo de vuelo. El impacto de la
heterogeneidad es integrado en el tiempo de vuelo y en la trayectoria de las
streamlines. Los clculos de los procesos fsicos son reducidos a soluciones
unidimensionales a lo largo de las cuales son distribuidas en el espacio con
una alta resolucin areal y transversal.
1.3.3 Simulacin streamline vs. simulacin en diferencias finitas. En ingeniera

de yacimientos se dispone de un juego de herramientas para modelar un
yacimiento, variando por simple analoga los mtodos de simulacin. La
herramienta depende de la disponibilidad de los datos, la cantidad de tiempo
destinada para la simulacin y la exactitud de los resultados requeridos. Si en
un estudio se desea realizar una simulacin, entonces se hace necesario
evaluar si es preciso usar simulacin streamline o simulacin en diferencias
finitas, entendiendo que no solo una tcnica de simulacin puede ser aplicada
en todos los casos.
La simulacin streamline basada en el flujo es conveniente para modelar

yacimientos extensos que son dominados por procesos de desplazamientos
convectivos, es decir procesos de inyeccin de agua, flujo miscible, WAG,
donde el comportamiento PVT no depende demasiado de la presin (sistemas
quasi-incompresibles).
____________
DATTA-GUPTA, AKHIL. Streamline Simulation: A Technology Update. SPE, Texas A&M University.
Dic. 2000.
26
De otra forma, la simulacin en diferencias finitas es conveniente en modelos

pequeos, donde son importantes los detalles fsicos del flujo tales como
compresibilidad, efectos capilares y permeabilidades relativas de histresis. Por
estas razones, antes de realizar una simulacin Baker* en su artculo tcnico
sugiere tener en cuenta los siguientes aspectos:
Cun importantes son los detalles geolgicos en el proceso de

simulacin y qu tan extenso es el yacimiento a modelar?
Cul es el grado de importancia que se le da al comportamiento del
flujo en el modelamiento del yacimiento?
En simulacin composicional pueden originarse preguntas similares; es decir si

es ms importante considerar el nmero de grid-blocks o el nmero de seudo
componentes.
Figura 1. reas de aplicacin de la simulacin streamline vs. simulacin en diferencias

finitas.
DOMINIO DEL SIMULADOR
, , = f (P, x )
Dominio de
trminos
Convectivos
(heterogeneidad)
Complejidad de Fluidos
Dominio de
trminos
Compresivos
Simulador
Diferencias
finitas
Simulador
Streamline
V = f (K )
Tamao del Modelo

Tomada de StreamSim Technologies Inc. 1999.
____________
BAKER, KUPPE, CHUGH, BORA, BATYCKY, STOJANOVIC. Full-Field Modeling Using StreamlineBased Simulation: 4 Case Studies. Streamsim Technologies Inc. SPE 66405. Texas, Feb. 2001.
27
En una simulacin convencional en diferencias finitas, hay un segmento tanto

de presin como de saturaciones a resolver. En diferencias finitas se resuelve
para la presin y entonces se calcula el flujo basado en la distribucin de
presin determinada (Formulacin IMPES), pero el transporte de flujo se hace
bloque por bloque, mientras que en un modelo de simulacin streamtube o
streamline, los fluidos son transportados a lo largo de streamlines como se
muestra en la figura 2.
Debido a que el problema de transporte tiene un comportamiento altamente nolineal, el mtodo de solucin en diferencias finitas puede ser muy sensible al
tamao y a la orientacin de las celdas, afectando la simulacin. En una
simulacin streamline, la ecuacin de presin es resuelta sobre un grid de
simulacin, similar al mtodo de una simulacin convencional. Despus, las
streamlines son calculadas ortogonal mente a los contornos de presin.
Por lo tanto, una red de transporte natural es construida y el fluido es

transportado a lo largo de cada streamline, rastreando el movimiento del gas, el
aceite y el agua en el yacimiento.
Las streamlines tienen por lo tanto una ventaja inherente porque el fluido es
transportado en la direccin del gradiente de presin a lo largo de las
streamlines y no entre grid blocks como sucede en diferencias finitas y se
observa en la figura 2. Baker* denota la importancia de este mtodo y brinda
una perspectiva de este tipo de simulacin.
1.3.4 Ventajas de la simulacin streamline. Los modelos streamline poseen

algunas aplicaciones importantes sobre las simulaciones convencionales.
Entre los mltiples beneficios que trae consigo la simulacin streamline
sobresalen:
o Velocidad computacional y exactitud mejorada. Una ventaja de la simulacin

streamline sobre muchas otras aproximaciones tradicionales es su inherente
eficiencia computacional. La eficiencia es entendida como: memoria y eficiencia
computacional. Segn Thiele**, la memoria es el resultado de dos aspectos
fundamentales:
____________
*
BAKER, Richard. Streamline Technology: Reservoir History Matching and Forecasting= its Success,
Limitations and Future. JCPT. Abril. 2001.
**
THIELE R, Marco. Streamline Simulation. 6 International Forum on Reservoir Simulation. Schloss
Fusch, Austria. 3-7 Sep. 2001, Pg. 9-10.
28
Figura 2. Representacin del movimiento de los fluidos en el modelo streamline y en el

modelo de diferencias finitas
Modelo Streamline
Productor
Los fluidos en modelos
Streamline son transportados
a lo largo de Streamlines
Inyector
Modelo Diferencias finitas
Productor
Los fluidos en modelos
Diferencias finitas son
transportados a travs de los
grid blocks
Inyector
Tomada de la tesis de postgrado de VARGAS, JOSE A.
La simulacin streamline es una formulacin tipo IMPES y por lo

tanto involucra solo la solucin implcita de la presin.
El trazado de las streamlines y la solucin del problema de transporte

es hecho de una forma secuencial. Slo se guarda una streamline en
memoria a un tiempo dado.
La eficiencia computacional se logra porque:
Las ecuaciones de transporte a lo largo de las streamlines a menudo

pueden ser resueltas analticamente.
El nmero de streamlines incrementa linealmente con el nmero de

celdas
activas.
Para
desplazamientos
dominados
por
heterogeneidades, el tiempo de clculo a menudo se mide casi
linealmente con el nmero de grid-blocks, haciendo de este el
mtodo preferido para simulaciones geolgicas a pequea escala.
29
Las soluciones numricas 1D a lo largo de las streamlines no se ven

afectadas por los criterios de estabilidad geolgica del grid,
permitiendo de esta forma largos pasos de tiempo.
Las streamlines no necesitan ser actualizadas frecuentemente.
Una simulacin streamline puede ser 100 veces ms rpida que una simulacin
convencional.
o Visualizacin cuantitativa del flujo.

La caracterstica ms atractiva para
muchos ingenieros es el poder visual que tienen las streamlines para perfilar
los patrones de flujo. En lugar de tener una secuencia de tiempos de los
cambios de saturaciones, las streamlines ofrecen una visualizacin del flujo del
campo, mostrando claramente como los pozos, la geometra del yacimiento y
las heterogeneidades interactan para representar el movimiento de fluidos de
inyectores a productores.
o Modelamiento total del campo. La mejor aproximacin que se puede realizar
es modelar totalmente el campo permitiendo a los patrones evolucionar los
aspectos impuestos por las interacciones de las ubicaciones de los pozos, las
tasas, la arquitectura del yacimiento y la heterogeneidad del medio. Pero la
habilidad de optar por este tipo de modelamiento requiere de una eficiente
aproximacin de simulacin, ambas en trminos de la memoria de
almacenamiento as como del tiempo de cmputo.
Los modelos completos para un campo pueden ser notoriamente extensos,

incluso cuando se usa un nmero limitado de celdas entre pozos. Mientras la
simulacin de flujo streamline simplifica las suposiciones para lograr eficiencia,
en muchos casos el modelo streamline para todo un campo se hace preferible
a utilizar que un modelo sectorizado bajo aproximaciones tradicionales. Esto se
debe a que pueden originarse errores ms significativos en la determinacin
aproximada de los lmites del sector que aquellos introducidos por el amplio
modelo streamline.
o Descripcin de la fsica del flujo. Existen suposiciones ocultas en la
formulacin de la simulacin streamline, particularmente con respecto a la
fsica del flujo. Esto se debe a que la tcnica se desarroll dentro de una
estructura incompresible, con el principal objetivo de capturar el
comportamiento del flujo producto de la configuracin del yacimiento y de las
interacciones entre
pozos; problemas para los cuales los simuladores
tradicionales de diferencias finitas no satisfacen adecuadamente,
particularmente cuando los modelos llegan a ser muy extensos y heterogneos.
30
Debido a esto, el enfoque que se dio a la simulacin streamline fue determinar

la eficiencia de desplazamiento para los diferentes tiempos de ejecucin
aumentando progresivamente la complejidad fsica. Esto se debe a que la
complejidad fsica tiende a incrementar el nmero de streamlines que necesitan
ser actualizadas y el tiempo requerido para resolver el problema de transporte
a lo largo de cada streamline. Esta metodologa favorece la investigacin de
problemas iniciando con el modelo ms sencillo y progresivamente agregar la
fsica de flujo requerida, llegndose a considerar la compresibilidad y
comportamientos de fase complejos. Contrario a lo que se cree, el hbito en
simulacin de yacimientos por mtodos tradicionales, ha sido incluir tanta
complejidad fsica como el simulador permita, es decir, iniciando con el modelo
fsico ms complejo.
o Incompresibilidad y controles de pozo.

En sistemas realmente
incompresibles todo lo que se requiere es una diferencia de presin para
calcular la velocidad total usando la ley de Darcy. Aunque no hay sistemas
realmente incompresibles, la suposicin de incompresibilidad es
matemticamente poderosa y puede ser usada cuando sea posible. Para
sistemas con fuertes mecanismos de desplazamiento con agua, sistemas que
tienen una relacin de movilidades cercana a uno, o sistemas que permanecen
por encima del punto de burbuja, la suposicin de incompresibilidad ha sido
usada con gran xito. stos son los sistemas donde se ha probado que las
streamlines trabajan particularmente bien.
Una consecuencia muy atractiva de los sistemas incompresibles es que las

tasas histricas de los pozos pueden ser conservadas, sin tener que asegurar
previamente que los modelos de los pozos puedan dar una presin de fondo
fluyendo, es decir, P>0. Esto tiene implicaciones muy importantes en la etapa
de ajuste histrico. Los modelos incompresibles permiten al ingeniero iniciar el
ajuste con las tasas de las fases observadas sin considerar la presin.
o Generacin de nueva informacin. Los simuladores streamline van ms all

de su atractivo visual, debido a la produccin de nuevos datos no disponibles
con simuladores convencionales. Esta, es posiblemente la contribucin ms
interesante y valiosa de los simuladores streamline al rea de simulacin de
yacimientos, aunque la industria no aproveche an esta informacin.
Ya que las streamlines inician en una fuente y terminan en otro punto, es
posible determinar qu inyectores o parte de un acufero estn soportando a un
productor en particular. Un alto corte de agua en un pozo productor puede por
lo tanto ser rastreado brindando informacin de influencia de nuevos pozos y
presencia de nuevos lmites. Recprocamente, es posible determinar
exactamente, cuanto volumen de una inyeccin particular de un pozo est
contribuyendo a los productores, informacin particularmente valiosa cuando se
trata de hacer balance de patrones.
31
Las streamlines tambin permiten identificar el volumen del yacimiento

asociado a cualquier inyector o productor en los sistemas. Por primera vez, es
posible dividir el yacimiento dentro de zonas de drenaje dinmicamente
definidas atadas a pozos. Todas las propiedades normalmente asociadas con
los volmenes del yacimiento pueden ahora ser expresadas sobre una base
por pozo, tales como OIP, WIP, presin promedio, etc.
o Ajuste histrico rpido o integracin de los datos de produccin dentro de

modelos de yacimientos de alta resolucin.
Para estudios de campo, el
periodo de ajuste histrico puede ser reducido de dos a cinco veces lo que
durara una simulacin tradicional. Considerando que los tiempos de ejecucin
de una simulacin convencional incrementan con el nmero de grid blocks en
una forma exponencial (n=2 a 3), la simulacin streamline aumenta casi
linealmente con el nmero de grid blocks. De est forma, la tecnologa
streamline permite la simulacin rpida de modelos ms extensos y complejos.
o Mejor identificacin de las reas de drenaje y habilidad para mostrar modelos

geolgicos detallados. Una clara ventaja de la simulacin streamline es que
permite una fcil visualizacin de los factores de ubicacin de las reas de
drenaje y las relaciones entre inyectores y productores. Este comportamiento
se aprecia en la figura 3. Esta visualizacin es extremadamente til en la
optimizacin de desplazamientos con agua o gas debido a que los beneficios
de la inyeccin pueden ser fcilmente cuantificados. Tambin es til para
determinar como sera el comportamiento del yacimiento si se perforaran ms
pozos o se cerraran algunos de ellos.
En cuanto al modelamiento geolgico, la simulacin streamline es exitosa ya
que permite una mejor resolucin areal/vertical del yacimiento, debido al gran
nmero de celdas empleadas.
Figura 3. Visualizacin de las streamlines en el modelamiento de flujo de fluidos.
Tomada del Simulador Frontsim.
32
1.3.5 Desventajas de la simulacin streamline.

Entre las desventajas que
puede tener la simulacin streamline se encuentran:
Dificultades para incorporar procesos fsicos complejos
Dificultades en el manejo de mecanismos transversales a las

streamlines.
Son menos exactos en el balance de materiales que los simuladores

de volmenes finitos.
No tienen en cuenta la presin capilar y como consecuencia no

pueden ser usados en casos en los cuales domine este aspecto.
El clculo del tiempo puede llegar a ser similar a los simuladores

convencionales siempre que los fluidos sean altamente
compresibles.
No ofrecen todas las posibilidades y solidez de los simuladores

tradicionales para procesamiento de pozos y manejo global de todos
los factores de produccin.
33
2.
SIMULACIN STREAMTUBES: PRIMEROS PLANTEAMIENTOS

HACIA LA SIMULACIN STREAMLINE.
Los mtodos streamtube y streamline han sido usados en el clculo de flujo de

fluidos por muchos aos. Las primeras aplicaciones para simulacin de
yacimientos petroleros fueron reportadas por Higgins* et al. a principios de los
sesenta, LeBlanc** et al. y Martin*** et al. en los aos setenta. Estos mtodos
numricamente resuelven modelos de flujo de fluidos complejos y multifsicos,
en medios porosos con un balance razonable entre eficiencia computacional y
fsica modelada. Antes de hablar sobre la simulacin streamline actualmente
desarrollada, se considera lgico dar una mirada a la formulacin streamtubes,
ya que fue sobre est tcnica donde empezaron a trabajarse conceptos como
las definiciones de streamlines.
El objetivo primario de la simulacin streamtube es facilitar rpida y

exactamente soluciones numricas a los desplazamientos a travs de
sistemas fuertemente heterogneos mientras se conservan los detalles de los
modelos fsicos. La suposicin fundamental en el uso de la aproximacin
streamtube reposa sobre la creencia de que los desplazamientos de escala de
campo son dominados por la heterogeneidad del yacimiento.
Las regiones de flujo rpido y lento en el yacimiento pueden ser representadas

usando streamtubes semi-unidimensionales. Los streamtubes, como se aprecia
en la figura 4, pueden ser visualizados como un grupo de tuberas, teniendo
geometras variables y conectando los pozos inyectores con los productores.
La forma de cada streamtube es dictada por la geologa del yacimiento y lo
ms importante es que en cada tubera se asume que existe conservacin de
masa.
____________
*
HIGGINS, R.V., BOLEY, D.W., AND A.J. LEIGHTON. Aids to Forecasting The Performance of Water
Floods," JPT. Sep 1964. Pg 1076-1082.
**
LEBLANC, J.L. AND CAUDLE, B.H.A Streamline Model for Secondary Recovery," Society of Petroleum
Journal. March 1971. Pg 7-12.
***
MARTIN, J.C., WOO, P.T., AND WEGNER, R.E. Failure of Stream Tube Methods to Predict Waterflood
Performance of an Isolated Inverted Five-Spot at Favorable Mobility Ratios,". JPT .Feb 1973. Pg 151153.
34
Figura 4. Descripcin grfica del modelo streamtube.
Tomada de Internet.
2.1 NO LINEALIDAD
La dificultad fundamental en resolver las ecuaciones diferenciales parciales que
gobiernan el flujo a travs de medios porosos es su formulacin no-lineal. En
otras palabras, para explicar la fsica relevante del flujo de fluidos, los
coeficientes que aparecen en las ecuaciones que gobiernan el proceso, se
convierten en funciones de las variables independientes del problema,
generalmente saturaciones de las fases y/o composiciones totales.
Un caso especial ocurre cuando los coeficientes son asumidos constantes con
respecto a las variables independientes, asumiendo una relacin de movilidad
igual a 1.
En este caso, los streamtubes son fijados con el tiempo y el flujo es supuesto
como lineal. Para una inherente no-linealidad de otros desplazamientos, la
aproximacin streamtube peridicamente actualiza la informacin de los
streamtubes y el trazado de la solucin unidimensional, utilizando la
aproximacin de Riemann.
2.1.1 La Aproximacin de Riemann en la formulacin streamtubes. El mtodo

ms usado para modelar algunos mecanismos de desplazamiento (flujo patrn,
flujo de dos fases inmiscibles, flujo miscible, flujo composicional), se centra
sobre la idea de que un streamtube es un objeto quasi-unidimensional. Las
soluciones en dos dimensiones son entonces construidas por trazado de
soluciones unidimensionales a travs de ecuaciones de conservacin de
masa, a lo largo de cada streamtube.
35
Debido a que los streamtubes son tratados como objetos unidimensionales, las
ecuaciones de conservacin son resueltas usando condiciones lmites tipo
Riemann, y trazando o rastreando a lo largo de los streamtubes dichas
soluciones.
La aproximacin de Riemann es introducida para manejar las dificultades

asociadas con las condiciones iniciales de tipo general que se originan a lo
largo de la actualizacin peridica de los streamtubes.
2.2 CLASES DE PROBLEMAS

La aproximacin streamtube pretende resolver problemas que son dominados
por las heterogeneidades del yacimiento y las fuerzas convectivas.
2.3 MODELO MATEMTICO

Por definicin, una streamline es una lnea tangente al vector velocidad en un
instante de tiempo dado.
Dos streamlines adyacentes definen un streamtube, el cual tiene un volumen y

transporta una tasa de flujo fija. Por definicin, un Streamtube se puede ver
como una tasa de flujo volumtrica, la cul est dada por la diferencia de la
streamfuncin asociada con las streamlines limitantes.
Para Thiele*, la
fundamentales:
aproximacin
streamtube
descansa
sobre
ideas
Generacin de las streamlines y streamtubes para un dominio

particular de inters.
Trazado o seguimiento de una solucin unidimensional a lo largo de

cada streamtube.
Para problemas 2D es posible resolver para la streamfuncin usando

directamente la siguiente ecuacin:

x
+
x
y
= 0
..
(1)
Y estableciendo el valor de la streamfuncin, , a intervalos iguales.

____________
THIELE, MARCO R. Modeling Multiphase Flow in Heterogeneous Media Using Streamtubes. Dic 1994.
36
2.3.1 Condiciones lmites. La ecuacin (1) es una ecuacin diferencial parcial

elptica, la cual requiere condiciones lmite tipo Neumann o Dirichlet. Solo
dominios de seccin transversal son considerados, y por consiguiente, las
condiciones lmite a ser consideradas son no flujo en el tope y en el fondo del
dominio, y presin constante o tasa uniforme en cualquier extremo.
La tasa de flujo volumtrica en cualquier punto en el dominio puede ser escrita

en forma diferencial como:
dQ = u d A
.... (2)
Donde d A es un rea arbitraria entre dos streamlines definida como:
d A = ez d s ..... (3)
es el vector unitario perpendicular al plano XY, y ds es la longitud de
d A.
La tasa de flujo volumtrica total que atraviesa un rea arbitraria A, entre dos
streamlines A y B es simplemente dada por:

A u dA = QAB = A u ez ds
u dy
x
... (4)
u y dx
d = B
....... (5)
A
En otras palabras, la tasa de flujo total entre dos streamlines es dada por la
diferencia en el valor de la streamfuncin asociada con cada streamline.
Usando este aspecto, las condiciones lmites para el dominio de seccin
transversal es relativamente sencillo de encontrar. La diferencia entre el valor
de la streamfuncin entre el lmite de tope y fondo debe ser igual a la tasa de
flujo total. Una opcin obvia es utilizar el lmite de fondo a =0, y el lmite de
tope a =QTotal. Similarmente, una distribucin de tasa uniforme a lo largo de
la cara de la entrada o de la salida del flujo debe darse por una distribucin
lineal de desde 0 a QTotal de la siguiente forma:
entrada
=
salida
yQ
Total
0 y 1 .. (6)
37
Un lmite de presin constante se presenta cuando el gradiente en la direccin

Y sea cero, es decir, a partir de la ecuacin de Cauchy-Riemann se tiene:
P
1
....
=
Y
y X
(7)
P
= 0 .... (8)
= 0
X
Y
Las posibles condiciones lmites existentes se resumen en la figura 5.
Retomando el planteamiento matemtico propuesto por Thiele*, la solucin
numrica de la ecuacin (1) en un dominio heterogneo con condiciones lmites
especificados en la figura (5), discretizada para un arreglo de 5 puntos da como
resultado la siguiente formulacin en diferencias finitas:
i +1, j
+ B
i 1, j
+ C
i , j +1
+ D
i , j 1
(A+ B +C + D)
i, j
=0
.. (9)
Donde los coeficientes A, B, C y D estn dados por:

A=
B=
C=
D=
x i , j +1
... (12)
x i , j 1
... (13)
... (11)
i1
y
,j
2
.. (10)
i +1
y ,j
Estas ecuaciones hacen referencia a la movilidad total y nunca puede ser igual
a cero a menos que la permeabilidad absoluta del bloque sea cero. Un
promedio armnico es usado si se desea encontrar el valor en los puntos
internodales.
____________
*
Pg 19-26.
38
Figura 5. Condiciones lmites. Modelo matemtico streamtube
= yQ
total
1 1
=0
+
x y x y x y
total
= 0
= yQ
total
=0
=0
Tomada de THIELE, MARCO R. Modeling Multiphase Flow in Heterogeneous Media

Using Streamtubes. Dic 1994
2.3.2 El papel de la streamfuncin en simulacin streamtubes. La

streamfuncin es un simple escalar desconocido que puede brindar una
completa descripcin del flujo. Todas las propiedades de flujo (las velocidades,
la presin entre otras) pueden estar relacionadas con esta funcin.
Cuando la streamfuncin es constante la relacin entre las coordenadas x y y
describen una sola streamline. Por el contrario, si la streamfuncin da un nuevo
valor, la relacin entre coordenadas x y y describen una streamline diferente.
Ver figura (6).
Figura 6. Streamlines y las streamfuncin.
= Cte
v
dy
=
u
dx
ni
B
dy
vi
streamline
n
dy
A
ds
nx
-dx
ny
-dx

39
Una segunda caracterstica de la streamfuncin es que la diferencia numrica

entre dos streamlines es igual a la tasa de flujo volumtrica entre ambas
streamlines. Para analizar esto, considrense dos streamlines con valores A
y B respectivamente. Dos puntos A y B son escogidos arbitrariamente y
conectados por cualquier camino como se muestra en la figura 6.
El flujo volumtrico entre las streamlines es:
AB
. n v ds
A B
n u + n y v ds
A B
.. (14)
Donde n es un vector normal al elemento ds.
Por geometra se tienen las siguientes relaciones:
n ds = dy
n y ds = dx
...... . (15)
AB
AB
... (16)
(ud vd ) = d
y
A B
=
B
. (17)
A B
(18)
La tasa de flujo entre dos streamlines es la diferencia entre sus valores de

streamfuncin.
Posteriormente x y y son las movilidades totales dadas por:
np
kk
=
n k k
=

x
j =1
j =1
rj
. (19)
rj
... (20)
Donde kx y ky son las permeabilidades absolutas en la direccin X y Y

respectivamente, krj es la permeabilidad relativa de la fase j y j es la
viscosidad de la fase j. np es el nmero total de fases presentes.
La ecuacin (18) se resuelve para la streamfuncin con condiciones lmites de
no flujo en el tope y en el fondo, flujo total constante en la entrada y presin
constante a la salida.
40
Con la geometra de los streamtubes conocida, Thiele* expresa que es posible

definir una longitud adimensional a lo largo de cada streamtube, i, usando la
siguiente ecuacin:
Di
A ( )d ... (21)
i
Donde es una coordenada a lo largo del streamtube, es la porosidad, Ai
es el rea de seccin transversal y V p es un volumen poroso arbitrario usado

para ajuste. Similarmente el tiempo adimensional puede ser expresado como:
Di
d
=
q t ...... (22)
NV
t
Donde N es el nmero de streamtubes y qi=qt/N porque los streamtubes son

tomados a incrementos iguales de .
Aunque el V es arbitrario y es el mismo para todos los streamtubes, la mejor

p
eleccin es hacer
= V pT
donde V pT es el volumen poroso total del
sistema.
En este caso, un sistema homogneo puede tener streamtubes con una
longitud adimensional de XD=1. Con XD y tD definidos por las ecuaciones (21) y
(22), es posible trazar cualquier composicin de fluido como una funcin de XD
y tD a lo largo de un streamtube.
Si la velocidad del campo es constante en el tiempo, entonces los streamtubes
necesitan ser determinados solo una vez, y la solucin puede ser construida
para cualquier tiempo tD. De otra parte, si la velocidad del campo cambia con el
tiempo, entonces los streamtubes deben ser peridicamente actualizados.
____________
Pg 30-31.
41
Para trasladar la solucin en el tiempo para cualquier desplazamiento no lineal

se plantea el siguiente algoritmo:
1.
Resolver para la velocidad inicial del campo. Usar la distribucin inicial

de movilidad total para resolver la streamfuncin, , y determinar los
streamtubes.
2.
Trazar la solucin 1D desde 0 hasta tD. Usar la ecuacin (21) para

determinar XD y trazar la solucin 1D sobre cada streamtube para un
tiempo dado tD.
3.
Actualizar la velocidad del campo, determinar la nueva distribucin de

movilidad y volver a determinar con stos datos la streamfuncin, , y
encontrar los nuevos streamtubes.
4.
Volver a trazar la solucin 1D desde 0 hasta tD+tD. Se lleva hasta un

tiempo tD+tD para generar las verdaderas condiciones iniciales no
uniformes y generar pequeos porcentajes de error en la aplicacin del
mtodo.
Con el tiempo adimensional y la distancia definida a lo largo de cada

streamtube, tambin es posible definir una velocidad adimensional como:
Di
X
t
Di
Di
A ( )d
i
N
NV p
0
=
qt t
A ( )d
qt
i
........ (23)
2.4 TRAZADO DE SOLUCIONES 1D
2.4.1 Mtodo A. Con la geometra del streamtube determinada, encontrando

ya sea la distancia adimensional o la velocidad adimensional es sencillo
evaluar la siguiente integral (expresin del volumen poroso) a lo largo de cada
streamtube:
s
A ( )d
i
.. (24)
Una primera aproximacin se encuentra definiendo el rea A como la diferencia

en la coordenada Y de dos streamlines definiendo un streamtube para un
42
valor particular de X. Con esta definicin, la integral a resolver o el volumen

poroso como una funcin de X, puede ser aproximado a:
Ai ( )d (Y A Y B )dX
s
Donde
X I
)
( +
2 i =1 Y Ai Y Ai +1 Y Bi Y Bi +1
... (25)
Y y Y son las coordenadas Y de las streamlines lineales, I es un

nodo tal que I (N + 1) , y la porosidad, , es asumida como constante. La
B
interpretacin de esta ecuacin se muestra grficamente en la figura 7.

Idealmente, el rea de seccin transversal A debe coincidir con el isbaro del
punto particular en el cual se desea la integracin.
Por otra parte puede argumentarse que el error es pequeo desde que el flujo
est principalmente en la direccin X, y los contornos de presin resultantes
sean justamente verticales. Adems incrementando el nmero de streamtubes
se obtiene una mejor aproximacin del rea de seccin transversal,
obtenindose un error significativamente reducido.
Figura 7. Mtodo de integracin simple a lo largo de los streamtube para determinar el

volumen poroso acumulativo en el punto X.

Lo ms importante es que el error ocasionado no es acumulativo. En cambio,

este es solo funcin de la aproximacin del rea, A, al lmite superior de la
integral, tal que:
s
Ai ( )d
X
2
( Y + Y ) A( X ) Y ( X )
I
i =1
i +1
43
.. (26)
Donde
Y = Y Y
A
.. (27)
2.4.2 Mtodo B. Aunque el error en la bsqueda del volumen poroso de un

streamtube en particular usando el mtodo A descrito anteriormente es
pequeo, puede tener la desventaja de declarar implcitamente que los frentes
pueden ser trazados como lneas verticales a lo largo de los streamtubes.
Esto puede resultar en algo similar a cortes irregulares de los frentes,

especialmente si el flujo en esta zona es fuertemente vertical. Un mtodo que
puede evitar esto se muestra en la figura 8. En lugar de usar una sola lnea
vertical para aproximar el rea de seccin transversal, A, dos lneas
perpendiculares son bajadas desde la streamline central hasta las streamlines
lmites. Cuando una lnea perpendicular no pueda ser bajada sobre una de las
streamlines, se asume un segmento vertical como en el mtodo A.
Figura 8. Representacin del mtodo B. Este mtodo define de una forma ms real los
frentes.
METODO A
METODO B

2.4.3 Trazado de las soluciones 1D sobre el grid cartesiano 2D. El paso final
en la construccin de la solucin en dos dimensiones para dominios
heterogneos de inters es trazar la solucin unidimensional sobre cada
streamtube.
La idea es que cada grid block cartesiano posee N*N puntos regularmente
distribuidos en su interior, como se puede observar en la figura 9.
44
Cada punto debe por lo tanto caer dentro de un streamtube particular y puede
ser asociado con un valor adimensional de volumen poroso XD de un
streamtube. Para un tiempo particular tD, un valor de concentracin/saturacin
es dado por la solucin unidimensional.
El valor promedio de cada grid block se calcula a partir de todos los valores de
los puntos existentes en el enmallado.
Cuando se dan suficientes streamtubes y suficientes puntos dentro de cada
grid block, los errores causados por el algoritmo de trazado pueden ser
menores.
Figura 9. Trazado de la solucin unidimensional a lo largo de los streamtubes sobre

el grid cartesiano fundamental
NODOS USADOS PARA EL TRAZADO
GRID BLOCK CARTESIANO

La ventaja de usar streamtubes versus streamlines radica en que los

streamtubes ofrecen una interpretacin visual de la velocidad de flujo local,
mientras que las streamlines no lo hacen. Un streamtube por otra parte,
permite la identificacin de regiones de flujo lento y rpido: secciones de gran
espesor de un streamtube corresponden a regiones de flujo lentas, secciones
delgadas a regiones de flujo rpido.
La geometra de los streamtubes por lo tanto, captura la distribucin de la

velocidad de flujo impuesta por la permeabilidad del campo como se muestra
en la figura 10.
45
Figura 10. a: Regiones de flujo rpidas. b: Regiones de flujo lentas.
a
b
b

2.5 SOLUCIONES EN DOS DIMENSIONES

La aproximacin streamtube ya definida gira en torno a la suposicin de que un
yacimiento heterogneo en dos dimensiones que es dominado por
heterogeneidades a gran escala puede ser descompuesto dentro de unas
series de streamtubes 1D. Debido a que las soluciones 1D son bien entendidas
para una variedad de mecanismos de desplazamiento, combinando stos con
los streamtubes permite una extensin de la fsica detallada representada por
aquellas soluciones de dominios 2D y 3D.
Thiele* en varios de sus estudios dedicados a la simulacin streamtubes
expone diversos tipos de desplazamiento a modelar con esta tecnologa.
____________
*
THIELE, MARCO R, BLUNT, MARTIN J. M. FRANKLIN, JR.ORR Modeling Flow In heterogeneous
Media Using Streamtubes I. Miscible and Immiscible Displacements. Stanford University. March 3, 1995.
____
Modeling Multiphase Flow in Heterogeneous Media Using Streamtubes. Dic 1994.
46
3. CONCEPTOS Y MODELO MATEMTICO DE LA SIMULACIN

STREAMLINE.
Los simuladores streamline se caracterizan por determinar inicialmente la

distribucin de presiones en el campo basados en las caractersticas del grid
de simulacin y posteriormente centra su objetivo en el clculo de la
distribucin de saturaciones a travs del yacimiento a lo largo de lneas de flujo
conocidas como streamlines.
El mtodo empleado es una formulacin tipo IMPES: Implcita en presin y
explcita en saturacin.
3.1 DEDUCCIN DE LAS ECUACIONES.

El desarrollo del modelo matemtico presentado por Batycky* para simulacin
streamline es el siguiente.
La ecuacin que gobierna el flujo de un componente i con np fases fluyendo en

un medio poroso es definida por Lake como:
np
j =1 t
( w S )+ w u
ij
ij
=
q s j ij
j S j D ij ij
(28)
representa la tasa de flujo de la fuente o el sumidero, D
ij
Donde qs
caracteriza el componente de dispersin, ij es la fraccin msica del
componente i en la fase j.
u es la velocidad de la fase dada por la Ley de Darcy como:

j
uj
Kk
=
rj
Pj + gD
j
... .. (29)
Teniendo en cuenta que Pj es la presin de la fase, D es la profundidad y g es

la constante gravitacional.
____________
BATYCKY, R. A Three-Dimensional Two-Phase Field Scale Streamline Simulator. Stanford University.

Jun 1997. Pg 15-18.
47
Para simplificar la ecuacin (28), se asume que los fluidos son incompresibles
(j=Constante) y que no existe dispersin, dando como resultado:
np
wij S j + ij u j = qs ij
j =1
. (30)
Ahora se asume que la ecuacin (30) se cumple para todos los componentes y
adems:
= 1 (31)
i =1
ij
De tal forma que reemplazando (31) en (30) obtenemos la ecuacin de balance

para flujo incompresible:
ut =
... (32)
Asumiendo que la velocidad total es definida para todas las fases presentes.
ut = K t P + g
j
(33)
La presin capilar es despreciada de tal forma que P=Pj. La movilidad total t y

la movilidad debido a las fuerzas gravitacionales g son definidas por:
np
j =1
np
j =1
rj
.. (34)
k g
rj
. (35)
Combinando las ecuaciones (32) y (33) obtenemos la ecuacin que permite

determinar la distribucin de presiones en flujo multifsico incompresible
K t P + g
j
D = q
.... (36)
48
La ecuacin que gobierna la distribucin de saturacin esta regida por el

mtodo IMPES y se deriva de la ecuacin (30). Para simplificar dicha ecuacin
se asume que las fases son inmiscibles, de forma que:
ij
= 0 si i j
S
t
=
j
q f
s
= 1 para i = j
ij
.. (37)
j ,s
Sustituyendo la ley de Darcy ecuacin (29) dentro de la ecuacin anterior y

eliminando el trmino P mediante la ecuacin (33) se llega a:
k rj
k rj
S j + j u + gD j np k rm

t

K
t
n p k rm
n p k rm m=1
m
m=1
m=1
m
m
( )= q f
s
j ,s
. (38)
Definiendo el flujo fraccional mediante la ecuacin de Buckely-Leverett como:
k
f
rj
k
np
m =1
. (39)
j
rm
Y el trmino de flujo fraccional para la gravedad como:
k rj
j np k rm
gD

K
n p k rm m=1
m
m=1

=
m j Gj
.... (40)
Empleando las anteriores definiciones, la ecuacin (38) puede ser escrita

como:
S +
f
t
ut +
=
G q f
j
. (41)
j,s
49
Recordando que para flujo incompresible
= 0 , la ecuacin que gobierna la
saturacin para flujo incompresible se expresa como:
S +
=q f
f + G
u
t t
..
(42)
j ,s
Las ecuaciones (36) y (42) son las ecuaciones no lineales para la formulacin
IMPES empleadas en simulacin streamline. Estas ecuaciones aparte de ser
no lineales, los coeficientes de cada ecuacin son dependientes de las
variables no conocidas (presiones y saturaciones). En comparacin con el
mtodo IMPES de la simulacin convencional en diferencias finitas, el mtodo
streamline permite la transformacin de la ecuacin (42) en un grupo de
ecuaciones unidimensionales (1D).
3.2 SOLUCIN A LA ECUACIN DE PRESIN

En el mtodo streamline, el yacimiento es dividido dentro de un sistema de grid
cartesiano como se realiza en un simulador convencional.
La ecuacin que gobierna la presin (ecuacin 36) es resuelta mediante un

mtodo estndar de diferencias finitas sobre el grid cartesiano. La forma
discretizada de la ecuacin (36) en un sistema 3D usando un patrn de siete
puntos sobre un grid cartesiano de coordenadas i,j,k , segn Batycky* est
dada por:
z ,k 1
2
i , j ,k 1
+ T y, j 1
i , j 1,k
+ T x,i 1
i 1, j ,k
Pi, j ,k T z ,k 1 + T y, j 1 + T x,i 1 + T z,k + 1 + T y, j + 1 + T x,i+ 1

2
2
2
2
2
2
+ T z ,k + 1 Pi, j ,k +1 + T y, j + 1 Pi, j +1,k + T x,i+ 1 Pi+1, j ,k =

2
z ,k 1
2
i , j ,k 1
+ Gy, j 1
i , j 1,k
+ Gx,i 1
. (43)
i 1, j ,k
Di, j ,k Gz,k 1 + Gy, j 1 + Gx,i 1 + Gz,k + 1 + Gy, j + 1 + Gx,i+ 1

2
2
2
2
2
2
+ Gz,k + 1
i , j ,k +1
+ Gy , j + 1
i , j +1,k
+ Gx,i+ 1
i +1, j ,k
s ,i , j ,k
____________

Jun 1997. Pg 18-20
50
Esta notacin asume que los ndices i,j,k son en las coordenadas x,y,z
respectivamente. Para un grid de bloque centrado, la transmisibilidad est
definida por:
z ,k + 1
=
2
2 x k y
z
K
... (44)
z
+
K
k +1
z ,k
t ,k
t , k +1
z , k +1
Donde x , y y z representan las dimensiones del grid-block.

Similarmente, las transmisibilidades interbloque para la gravedad estn dadas
por:
z ,k + 1
=
2
2 x k y
z
K
g ,k
z ,k
(45)
z
K
k +1
g , k +1
z , k +1
En la ecuacin (43), el trmino qs,i,j,k , representa el flujo de la fuente o del

sumidero en el bloque i,j,k. La inyeccin es asumida como positiva, mientras la
produccin es asumida como negativa.
3.3 CONDICIONES LMITES

Las condiciones lmites para la solucin de la ecuacin (43) son definidas en
los pozos y en los lmites de no flujo de la superficie del modelo del yacimiento.
Para cualquier pozo, la presin o la tasa total pueden ser especificadas. El
simulador streamline asume que el pozo es modelado con un gradiente de
densidad variable en la cara del pozo, pero no hay prdidas por friccin. La
ecuacin empleada para un pozo con nt estratos, est dada por:
q = T [P P ]
nt
k =1
Donde
block.
.... (46)
es la presin en la cara del pozo y Pk es la presin en el grid-
La transmisibilidad para el estrato del pozo,
2 z
r +
s
r
k
Ln
est dada por:
(47)
t ,k
o ,k
w ,k
51
Donde Sk es el factor skin, ro,k es el radio de Peaceman, y rw,k es el radio

interno del pozo. La movilidad en la cara del pozo, es asumida como la
movilidad del grid-block para pozos de produccin, y t,wk la movilidad de la
fase de inyeccin para pozos de inyeccin.
Finalmente, para pozos de mltiples estratos, la presin de cada estrato del

pozo es relacionada con la presin del pozo en el tope del completamiento (k*).
Se asume un gradiente de densidad variable en la cara del pozo mediante la
siguiente relacin:
Pk
+ 0 .5
i=
+1
i 1
)(D D )
i
i 1
.... (48)
Donde la gravedad especfica en la cara del pozo, i, en el estrato ith puede ser
calculado por:
g ,i
= g
j
(49)
t ,i
Combinando las ecuaciones (46) y (48) se obtiene la ecuacin para

condiciones de pozo:
nt
q s = T kw
k =1
k
P * Pk + 0.5 + (Di Di 1)
i 1
i
*
k
i = k +1
Las incgnitas en la ecuacin anterior son
.. (50)
P k para un pozo especfico con una

w
tasa total constante, o qs para un pozo especfico con una presin constante.
3.4 SOLUCIN NUMRICA DE LA MATRIZ DE PRESIN

Agrupando las ecuaciones (43) y (50) y separando las presiones desconocidas
del grid-block, se establece el grupo de ecuaciones discretizadas para ser
resueltas. En forma matricial el grupo de ecuaciones aparece como:
T P = B ... (51)
52
Donde T contiene los trminos de transmisibilidad del pozo y el grid-block de

las ecuaciones (43) y (50). Para el mtodo IMPES, T es una matriz simtrica. El
vector P contiene los valores de presin desconocidos de todos los grid-blocks
y de los pozos. El vector B contiene los trminos debidos a la influencia de la
gravedad, las tasas de las fuentes o sumideros y los trminos de los pozos
donde P w son definidos.
k
La ecuacin (51) forma un grupo de ecuaciones lineales con todas las

saturaciones dependientes de los trminos evaluados en el paso de tiempo
previo, mientras que P es determinada en el actual paso de tiempo.
La solucin a la ecuacin (51) no requiere informacin sobre la presin inicial.
Sin embargo, para definir la solucin, la presin de por lo menos un pozo debe
ser definida (Condicin Lmite tipo Dirichlet). La ecuacin (51) representa un
gran sistema lineal de ecuaciones, las cuales son resueltas mediante un
mtodo iterativo.
3.5 DETERMINACIN DE LA VELOCIDAD DEL CAMPO

Una vez la presin del campo ha sido determinada, los vectores velocidad del
campo son definidos para trazar las streamlines. La Ley de Darcy aplicada
entre dos nodos de presin, define la velocidad total de la interfase del gridblock como:
t ,k + 1
=
2
T
A
z ,k + 1
k+ 1
(P
k +1
Pk ) +
G
A
z ,k + 1
k+ 1
(D
k +1
Dk ) ... (52)
Donde Ak+1/2 representa el rea de seccin transversal de la interfase del gridblock. Para definir un vector velocidad en una cara del grid-block, el paso final
requiere convertir la velocidad Darcy en velocidades intersticiales vt, dividiendo
por la porosidad del grid-block. La velocidad intersticial es entonces definida en
la cara del grid-block en una direccin normal a esta.
3.6 MTODO PARA LA DETERMINACIN DE LAS VAS Y EL TRAZADO DE

LAS STREAMLINES.
Una streamline es definida como la curva instantnea en el espacio a lo largo
de la cual cada punto es tangente al vector velocidad. El trazado de las
streamlines desde inyectores hasta productores est basado en el mtodo de
trazado de Pollock*.
____________
POLLOCK, D.W. Semianalytical Computation of Path Lines for Finite-Difference Models. Ground Water.
Nov-Dic 1988. 26, N 6, Pg 743-750.
53
La suposicin fundamental es que la velocidad del campo en cada direccin de

la coordenada vara linealmente y es independiente de las velocidades en las
otras direcciones.
El mtodo de Pollock es atractivo debido a que es analtico y consistente con la

ecuacin de balance de materiales.
El desarrollo del mtodo es el siguiente:
Se calcula la distribucin de presiones mediante la solucin de la

ecuacin planteada en el modelo matemtico, usando un esquema
numrico implcito de diferencias finitas, que puede ser resuelto
mediante mtodos iterativos como gradiente conjugado entre otros.
A partir de la distribucin de presiones hallada, determinar la

velocidad total en cada una de las caras de la malla de simulacin
mediante solucin de la ecuacin de Darcy (ecuacin 52).
Con el flujo conocido, el algoritmo se centra en determinar el punto

de salida de la streamline y el tiempo, dada cualquier posicin de
entrada que asuma una aproximacin lineal de la velocidad del
campo en cada direccin de la coordenada. Ver figura 11
Figura 11. Mtodo de Pollock para el trazado de las streamlines conociendo el punto
de entrada
Ruta de la
Streamline
Salida
Entrada
Origen
Tomada de Streamline Simulation. 6 Internacional Forum on Reservoir Simulation.
54
Las ecuaciones que plantea el mtodo de Pollock son las siguientes:

Si v es la velocidad intersticial
= u ,
entonces una descripcin de velocidad
lineal en la direccin x es:
=
x
x0
g =
x
+g
xx
(x x );
.. (53)
x 0
... (54)
Donde vx0 es la velocidad en x a un x=x0, y gx es el gradiente de velocidad en la

direccin x. Ahora se tiene que:
dx
dt
.. (55)
Se puede integrar la expresin de la velocidad en las direcciones (x,y,z) para

obtener el tiempo de salida fuera de cada cara dada por un punto de entrada
arbitrario (xi,yi,zi), y coordenadas de salida (xe,ye,ze):
ln
g x
1
x0
x0
g
+g
(x x )
1
ln
=
t
(x x )
g
e
y0
y0
(y y )
+ g ( y y )
+g
e
i
+ g (z z )
0
e
z0
z
ln
g z z 0 + g z (z i z 0 )
1
(56)
Si se conocen las coordenadas de entrada de una lnea de flujo en el bloque

(xi,yi,zi), se puede calcular las coordenadas de salida (xe,ye,ze), teniendo en
cuenta que la lnea de flujo debe salir a travs de la cara cuyo tiempo de
trnsito sea menor, entonces, el delta de tiempo de trnsito de la lnea de flujo
en el bloque es dado por:
= MIN
(t , t , t ) ...
x
(57)
Y definiendo el tiempo de vuelo como:

nDt
tv = tv = tm .... (58)
1
55
Donde nDt es el nmero de bloques que atraviesa la lnea de flujo.

Con esta suposicin y conociendo el tiempo de vuelo, las coordenadas de
salida, se determinan fcilmente a partir de las siguientes ecuaciones:
x = g ln[
1
xi
exp g
t )
m
x0
]+ x
..... (59)
ln [
g
1
yi
exp
(g t ) ]+ y
y
y0
.... (60)
0
ln[
g
1
zi
exp
(g t ) ]+ z
z
z0
.. (61)
Donde vxi, vyi y vzi son las velocidades en el punto de entrada. Repitiendo el
procedimiento a travs de todos los bloques del grid, siguiendo las lneas de
flujo desde un pozo productor hasta un pozo inyector, se obtiene la trayectoria
y los tiempos de vuelo de todas las lneas de flujo. Conociendo las
coordenadas de entrada y de salida de la lnea de flujo en cada bloque, la
trayectoria desde el punto de entrada hasta el punto de salida es trazada en
forma parablica. Ver Figura 12.
Figura 12. Trazado de las lneas de flujo.
ze
(x
y ,z
e
)
z
Trayectoria de
Lnea de Flujo
ye
xe
xi
(x , y , z )
i
yi
zi
Tomada de la Revista Fuentes. Volumen 3 de 2003.
Las ecuaciones de Pollock son derivadas asumiendo bloques ortogonales, pero

muy pocos modelos de yacimientos reales usan tal estructura cartesiana.
56
Usando una transformacin isoparamtrica, es posible transformar los grids de

geometra de punto distribuido dentro de unidades de cubos como se observa
en la figura 13, aplicar el mtodo de Pollock y entonces transformar la
coordenada de salida a un espacio fsico.
Figura 13. Celdas no ortogonales y el mtodo de Pollock.
Tomada de Streamline Simulation. 6 Internacional Forum on Reservoir Simulation
Usando el mtodo de Pollock y las modificaciones de este es posible trazar los

perfiles streamline a travs de cualquier grid real usado en simulacin de
yacimientos. Un estudio ms amplio y detallado sobre el trazado de las
streamlines en sistemas no-ortogonales y grids no estructurados fue
recientemente desarrollado por Mathieu Prvost* et al. en la Universidad de
Stanford. En esta tesis se presenta ampliamente las ecuaciones y los
procedimientos a seguir en el mtodo de Pollock adaptado para el tratamiento
de grids complejos como lo son en este momento las mallas PEBI y los grids
triangulares. Por su parte para un mejor entendimiento de este tipo de mallas,
la tesis Estudio Comparativo de las Tcnicas de Enmallado Empleadas en
simulacin de Yacimientos realizada por los Ingenieros Elquin Santaf** et al.
brinda una amplia descripcin y explicacin de estas estructuras y su forma de
simulacin.
En otro interesante estudio, Matringe*** et al. validaron las posibles fuentes de
error ocasionadas por el mtodo de trazado de Pollock y desarrollaron
ecuaciones y procedimientos encaminados a reducir estos inconvenientes. En
su anlisis comparan los resultados obtenidos por el mtodo de Pollock con
una solucin analtica, que en este caso es la encontrada por Morel Seytoux
para un patrn de cinco puntos.
____________
MATHIEU, MICHAEL G.EDWARDS, BLUNT, MARTIN J. Streamline Tracing On Curvilinear Structured

And Unstructured Grids. Imperial College,London .Jun 2002.
**
SANTA F, ELKIN, ENELSO, LUIS. Estudio Comparativo de las Tcnicas de Enmallado Empleadas en
simulacin de Yacimientos. Universidad Industrial de Santander. Nov 2004.
***
MATRINGE, S.F, GERRITSEN, M.G, On Accurate Tracing Of Streamlines. SPE Annual Technical
Conference. and Exhibition.Texas 26-29 Sept 2004 SPE N 89920.
57
3.7 EL TIEMPO DE VUELO

El tiempo de vuelo es el tiempo requerido para alcanzar una distancia s a lo
largo de una streamline basado en la velocidad del campo a lo largo de la
streamline. En los ms recientes estudios, los autores King* y Datta Gupta** han
usado el concepto de tiempo de vuelo para modelar el flujo en los yacimientos
de petrleo. Matemticamente el tiempo de vuelo se define como:
(s ) =
0
( )
d (62)
(
)
ut
La integral anterior es evaluada de forma que:
nblocks
t
j =1
c ,i
.. (63)
Donde tc,i es el incremento del tiempo de vuelo a travs de cada bloque i.

3.8 TRANSFORMACIN DE LA COORDENADA A LO LARGO DE LAS
STREAMLINES.
Otro procedimiento importante en esta etapa es el replanteamiento de las
ecuaciones de conservacin de masa en trminos del tiempo de vuelo. La
suposicin fundamental en esta derivacin es que las streamlines no pueden
cambiar en el tiempo.
Esta derivacin presentada por Thiele*** es simple e inicia con el conocimiento
de una ecuacin de conservacin uni-dimensional de una especie i dada por:
M + F
t x
i
=0
... (64)
____________
KING, M.J., MANSFIELD, MARK. Flow Simulation of Geologic Models. BP Exploration Ltd. Aug 1999.
Dic. 2000.
***
THIELE R, MARCO. Streamline Simulation. 6 Internacional Forum on Reservoir Simulation. Schloss
Fusch, Austria. 3-7 Sep. 2001. Pg 4.
**
58
La meta
es demostrar que por combinacin de las soluciones
unidimensionales a lo largo de las streamlines es posible reproducir la solucin
3D. En otras palabras que existe un vector v tal que:
M +
F
t
i
= 0 = u M i +
.. (65)
Definiendo una coordenada que es paralela a v (es decir, una streamline) es

posible escribir que:
.... (66)
Ahora considrese la definicin de tiempo de vuelo, que lleva a la siguiente

expresin:
d
=
=
... (67)
Y permite que la ecuacin de conservacin tridimensional sea reformulada

usando un flujo unidimensional a lo largo de una streamline como:
M
t
u
+ F i = 0
...
(68)
Varios aspectos deben tenerse en cuenta en esta derivacin. Entre stos se

encuentran que el fluido transportado a lo largo de cada streamline es
incompresible, que las streamlines no cambian con el tiempo y que las
soluciones 1D deben tener los mismos lmites y condiciones iniciales como el
problema 3D original. Pero la derivacin muestra que un problema de
transporte tridimensional puede ser reescrito en trminos de mltiples
problemas unidimensionales a lo largo de las streamlines.
59
Para el simple caso de una inyeccin de agua incompresible, es posible definir

la ecuacin de la siguiente forma:
S + f
t ut
j
=0=
S
j
Streamlines t
todos ( SL )
... (69)
El detalle ms importante acerca de esta ecuacin es que la velocidad total en

el problema 3D ha desaparecido dentro del tiempo de vuelo para cada
streamline. Esta es la integracin de sistemas heterogneos 3D dentro de una
serie de sistemas homogneos 1D en trminos del tiempo de vuelo que hacen
que el mtodo streamline sea atractivo.
3.9 SOLUCIN ANALTICA DE LA ECUACIN DE CONTINUIDAD

Recordando la ecuacin 69, cualquier desplazamiento que pueda ser modelado
por esta ecuacin puede ser trabajado mediante ecuaciones 1D a lo largo de
streamlines, no obstante, cada tipo de desplazamiento difiere en su
determinacin de la funcin de flujo fraccional y de la no linealidad del
problema. Para el evento en el que las condiciones iniciales o lmites
permanecen constantes durante el proceso de desplazamiento, el
planteamiento de soluciones analticas es una va alternativa para dar solucin
a las ecuaciones de flujo. Si las condiciones lmites son constantes, los
cambios en la mobilidad total del campo son los que determinan el nivel de
cambio en las vas de los streamlines.
Considerndose el caso especfico para un proceso de inyeccin de agua, por

definicin se sabe que fw es solamente funcin de Sw, entonces por regla de la
cadena, la ecuacin (69) puede ser escrita como:
S
t
f
S
S
S
t
'
w
=0
.. (70)
Donde
'
w
f w
S w
(71)
60
La ecuacin (70) es una ecuacin diferencial parcial hiperblica de primer

orden, cuya particularidad principal es que la informacin se propaga a lo largo
de lneas caractersticas a velocidad finita.
Usando relaciones apropiadas es posible obtener una solucin analtica para la

ecuacin (70).
El sistema de ecuaciones diferenciales que llevan a la solucin es el siguiente:
dt
= 1, w = t
dw
(72)
d
= fw' , = fw' w+ cte = fw' t + cte ...... (73)
dw
Como (0)=0 y t(0)=0. Entonces, cte=0. Por lo tanto, la solucin analtica para
una lnea de flujo definida por los tiempos de vuelo, viene dada por:
(x, y , z , t ) =
( x, y , z , t )
f
f
= w
S S w
'
w
t ... (74)
.... (75)
Con la ecuacin (75), es posible obtener la distribucin de saturacin de agua a

lo largo de una lnea de flujo, conociendo de antemano el tiempo de vuelo en el
punto de inters y el tiempo de simulacin, donde f es la pendiente de la
S
curva de flujo fraccional que se considera funcin solamente de Sw. La relacin

de fw y Sw es dada por la anterior ecuacin, y la distribucin de saturacin se
puede obtener numricamente usando la grfica de flujo fraccional. Est
relacin tambin es vlida para otros tipos de desplazamiento.
3.9.1 Clculo de las propiedades en cada bloque del grid de simulacin. Las
saturaciones deben ser trazadas desde el grupo de streamlines al grid de
simulacin para completar un paso de tiempo en el mtodo.
61
Deacuerdo con Mallison* et al., la saturacin en un bloque del grid de

simulacin por lo general es calculada como un promedio de las saturaciones
de las streamlines que atraviesan dicha celda. La figura 14 muestra
grficamente el procedimiento. Matemticamente se expresa como:
S celda = wi S i
(76)
i =1
Donde Si es la saturacin de un segmento de streamline y N es el nmero de

streamlines por celda. Por su parte, los pesos o variables wi se escogen
deacuerdo con su tiempo de vuelo en la celda:
wi =
... (77)
k =1 k
N
El trazado sobre el grid de simulacin claramente requiere que por lo menos

una streamline atraviese cada celda del grid de simulacin.
De igual forma, Batycky** sostiene que la movilidad total promedio y el tiempo
de vuelo promedio en cada celda pueden ser calculados a partir de las
siguientes frmulas:
tcelda =
i =1
N

i =1
.. (78)
i
celda =
ti
i i
i =1
i
k =1
N
... (79)
____________
MALLISON, B.T, GERRITSEN, M.G, MATRINGE, S.F. Improved Mappings for Streamline-Based
Simulation. Symposium on Improved Oil Recovery. Oklahoma, april. 2004. SPE N 89352. Pg 4.
**
Jun 1997. Pg 33.
62
En esta parte del procedimiento streamline pueden generarse errores en el

trazado de las saturaciones. Para resolver estos inconvenientes Mallison* et al.
de la Universidad de Stanford (2004) han desarrollado una nueva metodologa
basada en tcnicas de Krigging para mejorar el trazado y el clculo de las
saturaciones en esta etapa de la simulacin.
Figura 14. Forma grfica del mtodo empleado para calcular la saturacin en el grid de
simulacin en cada celda.
Tomada de Roxar Inc.
3.9.2 Bloques que no son atravesados por streamlines. En la mayora de los

casos de simulacin, no todas las celdas pueden ser atravesadas por
streamlines. Para estos bloques se debe calcular un tiempo de vuelo a fin de
definir las propiedades de los fluidos en estas celdas.
Para definir este tiempo de vuelo, simplemente se traza una streamline hacia la
celda ms cercana que ya tiene definido un tiempo de vuelo promedio. El
tiempo de vuelo en el bloque no atravesado es entonces calculado como:
noatravezado = bloque+ albloquecercano
.. (80)
3.9.3 Clculo del tiempo real. En cada tiempo en el que las posiciones de las
streamlines son nuevamente calculadas, se desconoce el tamao real de los
pasos de tiempo.
____________
MALLISON, B.T, GERRITSEN, M.G, MATRINGE, S.F. Improved Mappings for Streamline-Based
Simulation. IOR. Oklahoma, April. 2004. SPE N 89352. Pg 5.
63
Este problema surge debido a que la informacin sobre la saturacin en el grid

no se mueve explcitamente en el tiempo, sino que en su lugar est es
removida y una nueva solucin escalada 1D a un intervalo de tiempo t + t es
trazada a lo largo de las nuevas posiciones de las streamlines. La informacin
del registro anterior de saturacin solo se tiene en cuenta para los trminos de
mobilidad en la ecuacin de presin.
Para determinar el tiempo real se emplea la siguiente expresin:
T n+1 =
(W
inic
Wrn+1 + WIn WPn

Qn+1 f n+1 1
... (81)
Donde el flujo fraccional promedio de produccin en el campo es definido

como:
f n+1 =
(f
+ f n+1
2
(82)
Qn+1 es la tasa de inyeccin total durante el paso de tiempo n+1, Winic es el

volumen inicial de agua en el yacimiento, y Wrn +1 es el volumen de agua
presente despus de trazar todas las streamlines en el paso de tiempo n+1.
La produccin acumulativa de agua en el paso de tiempo anterior es calculada

por:
W = Q f T i
n
P
. (83)
i =1
Si todos los inyectores inyectan 100% agua, el volumen acumulativo de

inyeccin de agua es dado por la siguiente expresin:
n
W = Q i T i ... (84)
n
I
i =1
t es considerado un tiempo interno, el cual es requerido para el trazado de la

solucin 1D en las streamlines. T es el tiempo real que es externo y es el
tiempo de inters.
64
El tiempo real es calculado al final, despus de que una solucin ha sido

trazada sobre el grid de simulacin, y de nuevo es actualizado desde T hasta
un tiempo Tn+1=Tn+Tn+1.
3.10 SOLUCIN NUMRICA DE LA ECUACIN DE CONTINUIDAD

La solucin analtica presentada anteriormente es vlida bajo condiciones
iniciales uniformes. Sin embargo, los procesos de inyeccin de agua a menudo
involucran perforaciones de pozos, reconversin y cierres de otros o cambios
en la movilidad de las fases que hacen que las condiciones iniciales no se
mantengan uniformes durante toda la vida del proyecto. En estas
circunstancias es mejor emplear soluciones de tipo numrico.
Las soluciones numricas 1D a lo largo de las streamlines fueron primero
introducidas por Bommer* et al. en 1979 para resolver problemas de lixiviacin
de Uranio.
Irnicamente, en su caso, las streamlines fueron asumidas fijas y la solucin
numrica fue introducida debido a que no exista una solucin analtica para el
problema de inters. Batycky** combin esta teora en 1997 utilizando las
streamlines tridimensionales con una solucin general unidimensional, solucin
numrica en el espacio tiempo de vuelo. Esta fusin de ideas permiti usar la
simulacin streamline en casos reales de campo, donde las streamlines no
solo cambiaran debido a las diferencias de movilidad, sino tambin debido a
cambios en las condiciones de los pozos.
Con cada nuevo grupo de streamlines, las condiciones iniciales correctas
podan ser dibujadas sobre las streamlines, es decir, las condiciones existentes
al final del paso de tiempo anterior, y de esta forma avanzar numricamente en
el tiempo. Esto permiti mover correctamente los componentes en 3D a pesar
de cambios radicales y significantes en las geometras de las streamlines,
debido a cambios en las condiciones lmites, como en el caso de pozos
cerrados o la adicin de nuevos pozos. Usando una solucin numrica tambin
se hizo posible considerar cualquier solucin 1D a lo largo de las streamlines
incluyendo desplazamientos composicionales complejos.
En este estudio, se ha decidido considerar las ms recientes soluciones
numricas desarrolladas para este tipo de simulacin, que disminuyen las
fuentes de error en el mtodo y que ya han sido incluidas en los simuladores
comerciales. La razn radica en que las soluciones numricas convencionales
que inicialmente se plantearon para simulacin streamline llevaban implcitos
problemas de difusin numrica, originados por traslado de informacin entre
grids.
____________
*
BOMMER, M.P. AND SCHECHTER, R.S. Mathematical Modeling of In-Situ Uranium Leaching," Society
of Petroleum Engineers Journal (December 1979) Pg 393-400.
**
Jun 1997. Pg 77-85.
65
3.10.1 Planteamiento de las nuevas soluciones numricas. Teniendo en

cuenta la ecuacin (69) nuevamente para el caso tpico de inyeccin de agua,
esta se puede discretizar siguiendo la direccin del tiempo de vuelo utilizando
un esquema central de diferencias finitas, tomando como nodos los puntos
intermedios de la lnea de flujo en cada bloque segn el esquema de nodo
centrado mostrado en la figura 15. Esta ecuacin sera expresada como:
n +1
w
Sw
n
wi + 1
2
wi 1
.. (85)
=0
Donde, n se refiere al paso en el tiempo, i+1/2 e i-1/2 indican los lmites entre
los nodos.
La saturacin de agua al tiempo n+1 se obtiene por la siguiente ecuacin:
n +1
n
w
wi + 1
wi 1
. (86)
Esta aproximacin presenta dificultad para evaluar los flujos fraccionales en los
puntos medios que corresponden con las interfases entre dos bloques
adyacentes. Un mtodo prctico consiste en utilizar el promedio ponderado del
tiempo de vuelo entre los dos bloques, que matemticamente representara:
wi 1
(f
wi + 1
+
i
=
2
wi
(f
wi
i 1
. (87)
i 1
+
i
wi 1
wi +1
i +1
. .... . (88)
i +1
66
Figura 15. Discretizacin de la Ecuacin de Continuidad en Direccin del Tiempo de

Vuelo para Nodos Centrados.
i+1
i
i-1
i-1 i-1/2
i 1
i+1/2 i+1
Peddibhotla* et al. (1997) presentan una solucin para nodos distribuidos. Ver
Figura 16. En este caso los flujos fraccionales en los puntos medios
corresponden con los valores del bloque. En esta solucin emplearon un
esquema numrico de alta resolucin conocido como Disminucin de la
Variacin total (DVT) para minimizar los efectos de la difusin numrica y
prevenir oscilaciones.
El mtodo se describe a continuacin.
La variacin total se define como:
VT f
n +1
w
ND 1
f
i =1
n +1
wi +1
n +1
wi
... (89)
Donde ND es el nmero de nodos de la lnea de flujo.

____________
PEDDIBHOTLA, S., SPATH, J., BATYCKY, R. An Efficient PC Based Streamline Simulator for
Immiscible and Miscible Displacements. Petroleum Computing Conference. Texas, 8-11 Jun. 1997. SPE
N 38129.
67
Figura 16. Discretizacin de la Ecuacin de Continuidad en Direccin del Tiempo de

Vuelo para Nodos Distribuidos.
i+1
i
i-1
i-1
i+1
i 1
Para prevenir oscilaciones se emplea el criterio sugerido por Harten* (1983),

que expresa:
VT (F
n +1
w
) VT (F )
n
.... (90)
Donde F w es el flujo fraccional corregido por un trmino antidifusivo.
i+ 1
=
2
n
wi
+ (r )
n
wi +1
+
2
n
wi
.... (91)
____________
HARTEN, A. High Resolution Schemes for Hyperbolic Conservation Laws. J. Comp. Phys. Vol. 49 Pg
357. 1983.
68
Donde (r ) es una funcin que define el orden de discretizacin as:

(r ) = 0
.Un solo punto aguas arriba (Primer orden)
(r ) = r
.Dos puntos aguas arriba (Segundo orden)
(r ) = 1
(r ) =
Punto Medio (Segundo orden)
2 1
+ r
3 3
.Esquema de Leonard (Tercer orden)
r es la medida de la suavidad de los datos dada por:
r =
f
f
n
wi
wi +1
n
wi 1
n
wi
i +1 i
i i 1
.. (92)
Los esquemas convencionales brindan soluciones estables con acentuada

dispersin numrica. El uso de esquemas de alto orden minimiza esta
dispersin pero se pueden presentar oscilaciones. La forma de prevenirlas
consiste en escoger la funcin (r ) de tal manera que el trmino antidifusivo
sea maximizado en amplitud pero sujeto a la restriccin de DVT. Sweby* (1990)
deriv la siguiente condicin que garantiza el cumplimiento de lo anterior:
(r )
r
(r ) 2
... (93)
De acuerdo con la condicin de Sweby la funcin (r ) para las formulaciones

de alto orden son:
(r ) = max (0, min (r ,2 ))
. Dos puntos aguas arriba (Segundo orden)
(r ) = max (0, min (r ,2 ))
..... Punto medio (Segundo orden)
(r ) = max (0, min (2,2 r , (r )))
Esquema de Leonard (Tercer orden)
____________
*
SWEBY, P.K. High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws. SIAM J.
Numerical Anlisis Vol 21. Pg 995-1011. 1990.
69
Datta-Gupta* et al. (1997) mostr que el mejor esquema es el de Leonard y ha

sido implementado en su simulador S3D.
La saturacin de agua se puede calcular a travs de la siguiente ecuacin:
n +1
w
= Sw +
n
*
1
1
F i + 2 F i 2
.... (94)
Una vez las saturaciones han sido determinadas a lo largo de las streamlines,
estas puedes ser determinadas sobre el grid de simulacin. Tal y como ocurri
en el mtodo de soluciones analticas, la distribucin de saturacin en el grid de
simulacin puede ser determinada mediante la ecuacin 67, ya descrita con
anterioridad.
Por su parte, Batycky**, define que la movilidad total promedio en cada celda
es ahora determinada mediante la siguiente ecuacin:
t celda =

k
N fases rj S celda
j =1
....... (95)
3.10.2 Bloques que no son atravesados por streamlines. Como ya se explic

anteriormente en las soluciones analticas, existen celdas que no son
atravesadas por streamlines, pero en las cuales se requiere determinar su valor
de saturacin. La diferencia con las soluciones numricas es que se requiere
de una condicin lmite de inyeccin.
La condicin de inyeccin es definida trazando una streamline desde estos

bloques hasta un inyector cercano. De esta forma las propiedades a lo largo de
la streamline son calculadas de la forma ya descrita.
____________
* PEDIDIHOPTIA, S., DATTA-GUPTA, A., XUE, GUOPING. Multiphase Streamline Modeling In Three
Dimensions: Further Generalizations and a Field Application Reservoir Simulation Symposium. Texas 611 Jun 1997. SPE N 38003
**
Jun 1997. Pg 82.
70
3.10.3 Errores en el balance de materiales. El mtodo de solucin numrica

streamline por la naturaleza del trazado de las soluciones desde las
streamlines al grid de simulacin y despus determinar un promedio de
saturacin por celda basado en los intervalos de tiempo de vuelo, no asegura
ciertamente que el volumen sea conservado. Este problema en el balance de
volumen ocurre porque las que se mueven a lo largo de las streamlines son las
saturaciones y no volmenes explcitos. Para reducir este porcentaje los
diversos simuladores streamline brindan al usuario la opcin de corregir este
error en un orden del 0.01-0.1%.
3.11 ACTUALIZACIN DEL MODELO DE LNEAS DE FLUJO

Cuando se presentan perforaciones de pozos de relleno o reconversin y cierre
de pozos, el trazado de las lneas de flujo cambia. La actualizacin en el
trazado se puede realizar de una forma prctica utilizando el mtodo de Pollock
ya descrito con anterioridad. El problema consiste en obtener la distribucin de
la saturacin a lo largo de las nuevas lneas de flujo.
Con base en la solucin obtenida a travs de las lneas de flujo anteriores,

Peddibhotla* et al. emplearon un mtodo de interpolacin 3D con base en el
mtodo cuadrtico modificado utilizado para interpolacin de datos dispersos.
El mtodo es el siguiente:
Dada una funcin g con valores gm en el nodo (xm,ym,zm) para m=1 ,N.
La funcin interpolante G se define de la siguiente forma:
G(x, y, z) = Wk(x, y, z) * lk(x, y, z) ... (96)

k =1
Donde lk es una funcin cuadrtica que satisface la siguiente condicin:
lk(xk, yk, zk) = gk ....... (97)
____________
* PEDIDIHOPTIA, S., DATTA-GUPTA, A., XUE, GUOPING. Multiphase Streamline Modeling In Three
Dimensions: Further Generalizations and a Field Application Reservoir Simulation Symposium. Texas 611 Jun 1997. SPE N 38003. Pg 3-4.
71
Fijando los valores de g dentro de un conjunto de nodos vecinos, los pesos
W k( x, y, z) , deben satisfacer las condiciones siguientes:
Wk
(x , y , z ) = k
j
Wk
k =1
1 j = k
0 j k
(x , y , z ) = 1 ... (98)
j
Wk
(x, y, z) = Wk (x, y, z) = Wk (x, y, z) = 0 ...... (99)
x
y
z
Estos pesos son definidos con base en la funcin de distancia inversa:
Wk ( x , y , z ) =
Wk ( x , y , z )
.... (100)
Wm (x , y , z )
m =1
Donde,
(Rw dk ) *
Wk ( x , y , z ) =
Rw * dk
.. (101)
Definiendo estas condiciones:
(Rw dk )* =
Rw dk dk < Rw
0 dk Rw
Donde:
dk
Es la distancia euclidiana entre los nodos ( x , y , z ) y (xk , yk , zk ) .
Rw
Es el radio de influencia alrededor del nodo (xk , yk , zk ) .
72
Figura 17. Trazado del frente de saturacin para la actualizacin de las lneas de flujo.
Lnea de flujo a
= t i +1
Lnea de flujo a t = t i
Celda ijk
Tomada de LOLOMARI. The Use of Streamline Simulation in Reservoir Management:

Methodology and Case Studies
3.12 PASOS DE TIEMPO

La razn principal de por qu los mtodos streamline son ms rpidos que los
mtodos convencionales de diferencias finitas es que un simulador streamline
no tiene restricciones sobre el tamao del paso de tiempo debido a
consideraciones de estabilidad en el proceso. Sin embargo, para
desplazamientos no lineales, pueden todava existir consideraciones sobre los
lmites superiores del tamao del paso de tiempo, en orden a obtener
soluciones con significado fsico. Ante esta posibilidad es recomendable
realizar varias corridas, variando los pasos de tiempo, en el tiempo de
simulacin ya establecido y analizar la convergencia de las soluciones
obtenidas.
El proceso de paso de tiempos en simulacin streamline es particularmente

diferente de los simuladores convencionales. En algunos simuladores
streamline como el 3DSL desarrollado por Thiele*, se requiere un tamao de
paso de tiempo para resolver la presin tp , y un tamao del paso de tiempo
para el paso convectivo tcs. stos pasos de tiempo no tienen que ser del
mismo tamao.
____________
MANUAL SIMULADOR 3DSL. STREAMSIM TECNOLOGIES INC. Version 2.20 Agosto 2004.
73
El tamao del paso de tiempo de presin para un nuevo nivel de tiempo n+1,
es calculado usando un indicador que expresa como es de compresible el
sistema. En la mayora de los simuladores el indicador es el cambio en el
volumen total del fluido a condiciones estndar entre dos niveles de tiempo. El
tamao de un paso de tiempo de presin para un nuevo nivel de tiempo est
dado por:
n +1
p
= t p
n
DVOLMAX
.... (102)
Donde DVOLMAX es la fraccin mxima de cambio permitida en el volumen del

fluido, y V nf es el cambio adimensional en el volumen del fluido definido como:
n
f
(V
n
f ,s
V
n
n 1
f ,s
(103)
f ,s
Vf,s es el volumen de fluido a condiciones de superficie.

El tamao del paso de tiempo convectivo para el nuevo nivel de tiempo n+1 es
calculado basado en un indicador de rendimiento del volumen poroso.
El tamao de un paso de tiempo convectivo para un nuevo nivel de tiempo est

dado por:
n +1
cs
PV
n +1
DPVMAX
n +1
n +1
MAX Q
,Q
iny
prod
.... (104)
Donde DPVMAX es el mximo volumen poroso adimensional y PV es el

volumen poroso a la nueva presin y Q es la tasa inyectada o producida al
final de la solucin de presin.
3.13 EL CONCEPTO DE GRAVEDAD EN SIMULACIN STREAMLINE

La gravedad fue un problema en el planteamiento de la simulacin streamline.
74
Como ya se mencion, el vector velocidad total (el cual define una streamline)
es la suma de los vectores velocidad de las fases, pero los vectores de las
fases no son paralelos en presencia de la gravedad. Ver figura 18.
Una solucin fue presentada por Bratvedt* et al en 1996 usando el concepto de
un operador divisin, una idea que ya haba encontrado previa aplicacin en
frentes de trazadores.
El concepto del operador divisin en este caso es sencillo y resuelve las

ecuaciones de balance de materiales en dos pasos: primero, un paso
convectivo, que es tomado a lo largo de las streamlines, y segundo, un paso
gravitatorio, tomado a lo largo de las lneas de gravedad, es decir, lneas
paralelas al vector gravedad, g,. En el paso gravitatorio, los fluidos son
simplemente segregados verticalmente teniendo en cuenta solo la diferencia de
las densidades de sus fases.
Figura 18. Influencia la gravedad sobre los vectores velocidad de las fases. stos
vectores no son alineados con el vector velocidad total.
Petrleo
Velocidad Total
Streamline
Agua
La formulacin matemtica es la siguiente:
S
t
G (S
z
1
)=0
... (105)
____________
BRATVEDT, F., GIMSE, T. AND TEGNANDER, C. Streamline Computations for Porous Media Flow
Including Gravity. Transport in Porous Media, Vol. 25, N 1. Pg. 63-78. Oct 1996.
75
G (S ) = f ( )
np
i =1
u
u
j
... (106)
i
Dividiendo la ecuacin 105 en dos secciones, se pueden expresar los pasos

convectivo y gravitatorio como:
Paso Convectivo:
S
t
Paso Gravitatorio:
S
t
=0
1 G j
...
(S )
(107)
= 0 . (108)
Aunque el orden de la solucin no tiene importancia, la simulacin streamline

siempre resuelve el paso convectivo primero, seguido por el paso gravitatorio.
3.14 FLUJO INCOMPRESIBLE Y COMPRESIBLE

Muchos lquidos son difciles de comprimir (es decir, sus volmenes o
densidades no cambian mucho cuando se les aplica cierta presin), de manera
que su densidad es esencialmente una constante. Para dichos fluidos
incompresibles la ecuacin de continuidad se simplifica a:
v = 0
........ (109)
Todos los trabajos de streamtube y streamline en el pasado fueron restringidos

a la suposicin de flujo incompresible. La razn es que el flujo incompresible
introduce suposiciones simplificadas que son particularmente convenientes
para la simulacin streamline.
Cabe mencionar dos suposiciones en particular:
El origen y el final de las streamlines corresponden a pozos, significando

que todos las streamlines deben iniciar en una fuente (inyector) y
terminar en un sumidero (productor).
La tasa de flujo a lo largo de cada streamline es constante.
76
El problema es que no hay sistemas reales que sean incompresibles: todos los
casos reales de campo involucran caractersticas de flujo compresible. En flujo
compresible, las streamlines pueden iniciar o terminar en cualquier grid-block
que acte como una fuente o sumidero debido a la naturaleza compresible del
sistema, aun si el bloque no tiene ningn pozo.
Determinar y trazar streamlines en flujo compresible no es difcil. El algoritmo

de trazado de Pollock es vlido sin tener en cuenta cmo las velocidades de
flujo fueron determinadas. Una aproximacin muy conocida ha sido publicada
por Ingebrigtsen* (1999).
Otra diferente, pero no muy publicada aproximacin ha sido implementada

dentro de un cdigo comercialmente disponible conocido como 3DSL** y es
usada para modelar sistemas compresibles inmiscibles y miscibles de tres
fases ( Black oil e Inyeccin de Gas Miscible ).
Pero mientras las streamlines pueden modelar sistemas ciertamente

compresibles, la inherente ventaja de velocidad sobre mtodos de diferencias
finitas puede disminuir significativamente, dependiendo del tamao del
modelo y los mecanismos de desplazamiento que gobiernan el sistema.
Esto se debe a que si la presin absoluta del sistema debe ser resuelta para
capturar los trascientes, entonces, los lmites sobre el tamao global del paso
de tiempo son muy similares entre los mtodos streamlines y diferencias finitas.
Existen sin embargo ejemplos donde soluciones compresibles streamline son el
nico mtodo posible para brindar una respuesta cuando el modelamiento es
extenso en procesos de desplazamientos secundarios y terciarios.
Si las propiedades PVT son slo una dbil funcin de la presin, el trabajo
incompresible puede ser usado exitosamente.
____________
INGEBRIGTSEN, L., BRATVEDT, F., BERGE, J. A Streamline Based Approach to Solution of ThreePhase Flow. Reservoir Simulation Symposium Held in Houston. Texas 14-17 Feb. 1999. SPE N 51904.
**
MANUAL SIMULADOR 3DSL. STREAMSIM TECNOLOGIES INC. Version 2.20. Agosto 2004.
77
3.15 DESCRIPCIN DEL MTODO STREAMLINE

Durante una simulacin streamline se realiza el siguiente procedimiento:
1)
Dadas unas condiciones iniciales y de frontera convenientes, o
determinadas mediante la ecuacin 50, se determina la distribucin de presin
mediante la ecuacin 43 a partir del grid de simulacin. Posteriormente la
velocidad en las caras de las celdas es calculada numricamente de la
ecuacin 52.
2)
Despus que la velocidad en las caras de las celdas ha sido
determinada, se prosigue a desarrollar el paso convectivo, trazando a partir de
la velocidad ya encontrada las streamlines mediante el mtodo de Pollock
descrito anteriormente.
3)
Se resuelve ya sea en forma analtica o numrica la ecuacin 69
calculando las saturaciones a lo largo de las streamlines.
4)
Posteriormente se trazan las nuevas saturaciones streamline en el grid
empleando la ecuacin 76.
5)
Despus de que se ha desarrollado el paso convectivo se procede al
desarrollo del paso gravitatorio, trazando las lneas de gravedad.
6)
Se resuelve la ecuacin 108 calculando la saturacin debido a las
diferencias de densidad a lo largo de las lneas de gravedad.
7)
Posteriormente se trazan las nuevas saturaciones debidas a la gravedad
en el grid.
8)
Despus de desarrollado el paso convectivo y el gravitatorio para un
tiempo dado, se actualiza el tiempo y se vuelve a repetir el procedimiento
descrito anteriormente hasta el final de la simulacin.
Peridicamente se deben actualizar las streamlines para considerar efectos de

movilidad o cambios en las condiciones de los pozos. Ver figura 19. Una vez
las streamlines son nuevamente generadas, se recalcula el tiempo de vuelo a
lo largo de las nuevas streamlines. Finalmente, los clculos de saturaciones
son reasumidos con el tiempo de vuelo actualizado. Un paso crtico aqu es el
trazado de la informacin desde las streamlines viejas a las streamlines
nuevas. Esto puede ser una fuente potencial de error durante la simulacin
streamline.
78
Figura 19. Cambio de las geometras streamline debido a cambios en la movilidad y en

las condiciones de los pozos.
Finalmente, despus de realizada la simulacin puede calcularse el tiempo total

de la corrida T, mediante la siguiente ecuacin:
sol nsl sl
T t + t j
1
1
nts
... (110)
Donde:
nts : Nmero de pasos de tiempo empleados en la simulacin.
t sol : Tiempo requerido para resolver la ecuacin de presin a cada paso de
tiempo.
n sl : Nmero de streamlines en cada paso de tiempo.
t slj : Tiempo requerido para resolver la ecuacin de transporte para cada
streamline.
La figura 20 tomada de Baker* et al, muestra un esquema paso a paso del
procedimiento realizado en simulacin streamline.
____________
*
BAKER, R., KUPPE, F., CHUGH, S. BORA, R., STOJANOVIC, S., BATYCKY, R. Full-Field Modeling
Using Streamline-Based Simulation: 4 Case Studies. Streamsim Technologies Inc. Texas. Feb 2001. SPE
66405. Pg 11.
79
Figura 20. Diagrama de flujo para un simulador streamline.

Condiciones de
Pozo al tiempo
Condiciones Iniciales
Generar La Presin Del

Campo
Generar La Velocidad Del

Campo
Trazar Las Streamlines
Paso
Convectivo
sl
Resolver las Ecuaciones

1D para la Saturacin a lo
largo de los Streamlines
Trazar las nuevas

Saturaciones Streamline
sobre el grid de simulacin
Trazar las Lneas de

Gravedad
Paso
gravitatorio
gl
Actualizar el
tiempo a
i +1
= t i + dti
Resolver las Ecuaciones

1D para la Saturacin a lo
largo de las lneas de
Gravedad
Trazar las nuevas

Saturaciones de las lneas de
Gravedad sobre el grid de
simulacin
Fin de la Simulacin
Tomada de BAKER, Full-Field Modeling Using Streamline-Based Simulation.
80
4. APLICACIONES DE LA SIMULACIN STREAMLINE EN LA EVALUACIN

DE YACIMIENTOS.
Una vez se ha realizado el proceso de simulacin, entran en juego un sin

nmero de procedimientos encaminados a la evaluacin de los resultados
obtenidos y a la optimizacin del proceso, con el fin de poder hacer de la
simulacin una herramienta poderosa en la planificacin presente y futura del
campo.
La velocidad y la versatilidad de la simulacin streamline ha llevado a un

extenso nmero de aplicaciones durante los ltimos aos. Entre estas
aplicaciones se encuentran:
4.1 CLCULOS DEL VOLUMEN DE BARRIDO

La simulacin streamline provee una gran ventaja en el clculo del volumen de
barrido y las reas barridas bajo las condiciones ms generales. Esta es una
consecuencia directa de la formulacin del tiempo de vuelo. El tiempo de vuelo
streamline refleja la propagacin del frente de fluido a varios tiempos. Para un
tiempo dado, la conectividad en el tiempo de vuelo streamline provee una
medida directa del barrido volumtrico para una heterogeneidad y configuracin
de pozo arbitraria. El volumen de barrido es una cantidad fundamental que es
esperada para ser correlacionada con el recobro.
Recordando la ecuacin (62) correspondiente a la definicin de tiempo de

vuelo, la velocidad del campo para un medio general tridimensional puede ser
expresada en trminos de streamfunciones de la siguiente forma:
u =
.... (111)
Una streamline es definida por la interseccin de un valor constante para

con un valor constante para . En aplicaciones bidimensionales, se usan las
formas funcionales simplificadas, = ( x, y ) , = z , llevando de esta forma a las
siguientes expresiones conocidas:
.... (112)
x
81
Donde si se recuerda el planteamiento matemtico para el modelo streamtubes
corresponde a la streamfuncin.
Las tcnicas streamline estn basadas en la transformacin de la coordenada

del espacio fsico a la coordenada de tiempo de vuelo, donde todas las
streamlines pueden ser tratadas como lneas continuas de longitud variada.
Esta transformacin de la coordenada es facilitada por el hecho que el
Jacobiano empleado asume una forma extraordinariamente sencilla:
( , , )
= ( ) = u = .. (113)
(x, y, z )
Adems se tienen las siguientes relaciones entre el espacio fsico y las
coordenadas del tiempo de vuelo siguiendo la direccin del flujo:
dxdydz
= d d d
. (114)
Es fcil darse cuenta que la transformacin de la coordenada tambin

preserva el volumen poroso, una caracterstica esencial para calcular la
eficiencia del barrido volumtrico.
Considerando el caso bidimensional para un patrn de cinco puntos con un

pozo productor y un inyector ubicado cada uno en un extremo, el contorno del
tiempo de vuelo brinda las reas de barrido a diferentes tiempos. Ver figura 21.
El planteamiento matemtico que permite obtener el rea y el volumen de
barrido en un patrn determinado, fue desarrollado por Idrobo* et al (2000). El
planteamiento es el siguiente:
Para cualquier tiempo particular, el rea de barrido puede ser calculada

mediante la siguiente ecuacin:
barrido
(t ) = dxdy [t (x, y )]
.. (115)
____________
IDROBO, EDUARDO A., CHOUDHARY, MANOJ K., DATTA-GUPTA, A. Swept Volume Calculations
And Ranking Of Geostatistical Reservoir Models Using Streamline Simulation. Texas A&M University.
California, 1923 June 2000. SPE N 62557. Pg 2-3.
82
Cambiando el sistema de coordenadas de ( x, y ) a ( , ) , la ecuacin anterior se

reduce a:
barrido
(t ) = dd [t ( )]
. (116)
La ltima rea de barrido asociada con un pozo productor puede ser obtenida
considerando t :
A (t ) = ( )d
Final
(117)
Las ecuaciones (116) y (117) indican que las reas de barrido a cualquier
tiempo dado o la ltima rea de barrido asociada con un productor dado,
pueden ser calculadas por simple integracin del tiempo de vuelo en el
productor contra la streamfuncin normalizada para patrones de cinco puntos
homogneos y heterogneos como puede apreciarse en las figuras 21, 22 y 23.
En flujo tridimensional, se puede derivar anlogamente las anteriores
expresiones para la definicin de volmenes de barrido:
barrido
(t ) = dxdydz (t (x, y , z ))
barrido
(t ) = ddd (t ( , )) .... (118)
En lugar de evaluar las integrales complejas en 3D, Idrobo et al obtuvieron la

siguiente expresin para calcular el volumen de barrido examinando la
conectividad en el tiempo de vuelo streamline:
barrido
(
(t ) = d
i
) (t )q
( ) (119)
i
Donde es la funcin Heviside y q(i) es la tasa de flujo volumtrica asignada

al streamline i. Por otra parte, la eficiencia volumtrica de barrido puede ser
calculada simplemente dividiendo el volumen de barrido obtenido por el
volumen poroso total. Los clculos obtenidos para reas y volmenes de
barrido, as como eficiencias obtenidas son informacin que puede ser til en
el desarrollo de estrategias, tales como nuevas perforaciones y conversin de
patrones.
83
Figura 21. Distribucin del tiempo de vuelo y reas de barrido a diferentes tiempos
para modelos de cinco puntos homogneos y heterogneos
200 das
508 das
0.72 PV
Homogneos
400 das
451 das
200
das
200
das
200
das
Heterogneos
Figuras tomadas del artculo tcnico SPE 62557.
84
Figura 22. Clculo del rea de barrido basado en el tiempo de vuelo streamline para
patrones de cinco puntos homogneos.
Tiem po de Vuelo vs. Stream Funcin Norm alizada

Para un Patrn de Cinco Puntos Hom ogneo
3
2.5
1.5
0.5
AREA DE BARRIDO A 1 P V
0
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
S t r e a m Fu n c i n N o r m a l i z a da
Figura 23. Clculo del rea de barrido basado en el tiempo de vuelo streamline para
patrones de cinco puntos heterogneos.
Tiem po de Vuelo vs. Stream Funcin Norm alizada Para un Patrn de Cinco
Puntos Heterogneo
4
3.5
2.5
1.5
0.5
AREA DE BARRIDO A
0
0
0.1
0.2
0.3
0.4
0.5
0.6
S t r e a m F u n c i n N or m a l i z a d a

85
0.7
0.8
0.9
4.2 CLASIFICACIN DE LOS MODELOS GEOESTADSTICOS (RANKING)

Los simuladores streamline tienen la ventaja de ser mucho ms rpidos que
los simuladores convencionales de diferencias finitas, lo cual hace de stos una
herramienta ideal para clasificar extensos modelos geolgicos.
Lolomari* et al (2000) creen que los mtodos de simulacin streamline son
herramientas apropiadas para clasificar los modelos geolgicos por varias
razones importantes.
Una de ellas es la capacidad que tienen los simuladores streamline para

manejar modelos geolgicos sin necesidad de realizar un proceso de upscaling
previo. Esto preserva las heterogeneidades geolgicas a pequea escala en el
modelo, lo cual tiene significante impacto en el flujo del fluido. Otra razn
resulta de la combinacin de los desarrollos en los algoritmos streamlines y el
avance en los computadores. Los algoritmos de simulacin streamline son ms
exactos y al mismo tiempo los computadores llegan a ser ms poderosos.
Esto significa que extensos modelos geolgicos pueden ser simulados en

menor tiempo con mayor grado de exactitud.
Gilman** et al. (2002) han usado mtodos de simulacin streamline para
clasificar modelos geolgicos a gran escala y para dirigir estudios de
incertidumbre. Ellos declaran que los mtodos de simulacin streamline son
bien reconocidos como componentes importantes de los estudios de
caracterizacin integrada de yacimientos y creen que estos mtodos son
satisfactorios para resolver problemas a gran escala, ya que un nmero de
suposiciones simplificadas pueden ser tenidas en cuenta en la fsica del flujo.
En el estudio integrado de yacimientos, la construccin de los modelos no es

suficiente, estos modelos deben ser clasificados y su incertidumbre debe ser
evaluada basada en mtodos estticos o dinmicos.
____________
*
LOLOMARY, T., BRATVEDT, K., CRANE, M., MILLIKEN, J. The Use of Streamline Simulation in
Reservoir Management: Methodology and Case Studies. SPE 63157. Texas. Oct 2000.
**
GILMAN, J., MENG, H., ULAND, M., DZURMAN, P., COSIC, S. Statistical Ranking of Stochastic
Geomodels Using Streamline Simulation: A Field Application. Annual Technical Conference and
Exhibition. SPE 77374. Texas. 29 Sep-2 Oct 2002.
86
Uno de los mtodos estticos ms comunes es la correlacin estadstica entre

los datos de entrada y de salida. Tericamente todos los algoritmos de
modelamiento estocstico son capaces de reproducir estadsticas de los datos
de entrada (Isaaks y Srivastava 1989). Si hay una correlacin pobre entre los
datos estadsticos de entrada y de salida y esta diferencia no es justificada,
entonces significa que existen problemas con el modelo.
Otro mtodo sencillo de clasificacin es la visualizacin 3D del modelo

geolgico. Durante este proceso, los modelos pueden ser verificados para
corregir la geometra, la orientacin de los principales cuerpos geolgicos y
otras caractersticas anormales.
La clasificacin tambin puede realizarse sobre los valores de volumen poroso.

Aunque todos los modelos estocsticos son generados desde los mismos
datos, y tienen las mismas probabilidades de ocurrir, ellos tienen diferentes
volmenes porosos.
Por otra parte la clasificacin dinmica puede ser hecha usando un simulador
streamline. La conectividad en el tiempo de vuelo suministra una medida
favorable para la clasificacin de mltiples modelos geoestadsticos. Debido a
que la eficiencia volumtrica de barrido es fundamental en todos los procesos
de recobro, una clasificacin basada en la conectividad en el tiempo de vuelo
puede ser vlida sin tener en cuenta los procesos de recobro considerados. El
primer paso es seleccionar un adecuado nmero de pasos de tiempo para
optimizar el tiempo necesario para ejecutar la simulacin. Esto se puede hacer
simplificando los perfiles de produccin del pozo tomando una movilidad
promedio para las tasas de produccin e inyeccin, mientras se capturan los
principales eventos de produccin. Se deben intentar primero con varios pasos
de tiempo para obtener el nmero ptimo.
Despus de esto, diferentes modelos tales como P(10), P(50) y P(90) pueden
ser simulados a travs de un simulador streamline. Los resultados simulados
pueden ser analizados de acuerdo con el ajuste de stos con los datos
observados en el campo.
El anlisis debe ser dirigido por la generacin de grficas como las observadas
en las figuras 24, 25 y 26, entre los datos observados y simulados para los
tiempos de ruptura por pozo. Tambin los resultados obtenidos de movilidad
del fluido, eficiencia de barrido y produccin total en el yacimiento deben ser
graficados para todos los modelos simulados. Finalmente, el mejor modelo
geolgico debe ser seleccionado para el uso en el desarrollo futuro del campo.
87
Figura 24. Ranking basado en el % recobro vs. ranking basado en la eficiencia

volumtrica de barrido a la ruptura.
60
R s = 0 .9 9 7
50
40
30
20
10
0
0
10
15
20
25
30
35
40
45
50
R / %R e c

volumtrica de barrido momentos despus de la ruptura.
60
R s = 0 .9 6 5
50
40
30
20
10
0
0
10
15
20
25
R / %R e c
88
30
35
40
45
50

volumtrica de barrido momentos despus de la ruptura.
60
R s = 0 .9 0 4
50
40
30
20
10
0
0
10
15
20
25
30
35
40
45
50
R / %R e c
Tomada del artculo tcnico SPE 62557.
4.3 APLICACIN DE LA SIMULACIN STREAMLINE EN EL PROCESO DE

UPSCALING
Los modelos computacionales de yacimientos de aceite han llegado a ser cada
vez ms complejos en orden a representar la realidad geolgica y su impacto
sobre el flujo de fluidos. Las limitaciones de memoria y tiempo de CPU en
simuladores de diferencias finitas o volmenes finitos obligan a la utilizacin de
modelos ms robustos (coarse-scale model) a travs de procedimientos de
upscaling.
La simulacin de flujo de fluidos en un medio poroso sobre una escala de

campo requiere simulaciones numricas a gran escala. Si bien se han
desarrollado continuamente computadoras ms poderosas y tcnicas de
simulacin ms avanzadas, el tamao de las celdas usadas en las
simulaciones de flujo del campo es demasiado grande para tener en cuenta
explcitamente el efecto de las heterogeneidades de escala pequea. Por lo
general el nmero de celdas de un modelo geolgico es del orden de 106 o
107 celdas, mientras que los simuladores comerciales solamente manejan
modelos de simulacin del orden de 105 celdas. Por lo tanto es necesario la
aplicacin de una tcnica de escalamiento (upscaling) para reducir
(reagrupar) el nmero de celdas del modelo geolgico a una cantidad
manejable por un simulador comercial.
89
El upscaling es una tcnica cuyo objetivo es reducir el nmero de celdas de

la malla (grid) de simulacin, con el fin de ahorrar costos computacionales,
reducir el tiempo de simulacin y el manejo de datos; el modelo geolgico
detallado (fine grid) es reagrupado (coarse grid ) a un tamao razonable para
la simulacin de flujo.
El escalamiento (upcaling) es necesario para enlazar los modelos geolgicos

y de simulacin.
Dada una descripcin del yacimiento a escala muy fina y un modelo de

simulacin, un algoritmo de escalamiento asigna valores apropiados de
porosidad, permeabilidad y algunas funciones de flujo a las celdas gruesas
(coarse grid-block) de la malla de simulacin. Por lo general esta malla de
simulacin resultante se denomina malla gruesa (coarse grid), como se
muestra en la figura 27.
Figura 27. Resultado de una tcnica de escalamiento.
Fine Grid-malla Fina

6x4x6
Coarse Grid-malla Gruesa

2x2x2
Upscaling
Tomada de Internet.
El procedimiento de upscaling puede conducir a dificultades significantes en el

estudio de yacimientos:
Mientras el modelo geolgico a fina escala (fine-scale) es

construido a partir de datos petrofsicos, ssmicos y de
registros, su comportamiento dinmico nunca es verificado.
Como resultado, el estudio del yacimiento a una escala ms
robusta (coarse-scale) puede ser vinculado al modelo
geolgico a fina escala, pero los dos pueden ser
inconsistentes en su comportamiento dinmico.
90
Inversamente, el modelo escalado puede no ser propiamente

probado si el flujo y el comportamiento de produccin al nivel
de una escala muy fina no estn disponibles.
Se requiere un gran nmero de modelos por sectores en el

ptimo diseo de los patrones de pozo.
La simulacin streamline es una alternativa atractiva para superar algunas de

las desventajas ya mencionadas; adems, ofrece eficiencia computacional
mientras minimiza los efectos de difusin numrica y de orientacin del grid.
Esto permite la integracin de modelos geolgicos a fina escala dentro del
estudio del comportamiento del flujo en la ingeniera de yacimientos.
La simulacin dinmica del flujo es todava un problema en muchos estudios

integrados de yacimientos, los cuales tratan de reconciliar el modelo geolgico
con los datos ssmicos y de pozo. Herramientas tridimensionales, datos
ssmicos de alta resolucin as como aplicaciones mejoradas de modelamiento
esttico 3D, producen modelos que son cada vez ms detallados y permiten
identificar ms aspectos que la generacin previa de modelos estticos. En la
actualidad, los modelos a fina escala estn comnmente en el rango de 1 a 10
millones de celdas. De otra parte, la tecnologa de simulacin de flujo basada
en volmenes finitos o diferencias finitas ha madurado notoriamente.
Ante este desarrollo, el proceso de upscaling de los modelos geolgicos a fina

escala permanece como una realidad para muchos estudios, causando un
deterioro significativo en el modelo geolgico. Bajo condiciones reales del
yacimiento, rigurosos procesos de upscaling llegan a ser extremadamente
difciles forzando al ingeniero a realizar aproximaciones dubitativas.
Una metodologa que permita soluciones sobre el modelo geolgico original es

por lo tanto deseable y la simulacin streamline es una alternativa disponible.
El simulador streamline en la mayora de los estudios de yacimientos compara
el comportamiento dinmico de los grids a escala fina y a escala gruesa bajo
diferentes condiciones de flujo (sistemas incompresibles, compresibles y blackoil), apoyndose en comparaciones grficas de los modelos. La meta de un
proceso de upscaling es obtener un modelo de flujo representativo del grid
completo.
Samier* et al (2001) en su artculo tcnico exponen a grandes rasgos con
ejemplos prcticos, las caractersticas del proceso de upscaling y su validacin
mediante el simulador streamline.
____________
*
SAMIER, QUETTIER, THIELE, MARCO. Applications Of Streamline Simulations To Reservoir Studies

SPE Reservoir Simulation Symposium Held In Houston, Texas, 1114 Feb 2001. SPE N 66362.
91
La simulacin streamline es una herramienta til para validar el proceso de

upscaling de permeabilidades absolutas y determinar rpidamente el mejor
mtodo de upscaling. Los resultados en trminos de las curvas de produccin,
saturaciones y presiones graficadas pueden ser comparadas entre el grid
geolgico y el grid de simulacin. Ante esto surgen dos ventajas importantes:
La simulacin sobre el modelo geolgico a fina escala puede

ser mejorada en un tiempo de CPU razonable permitiendo
generar una solucin de referencia.
Las comparaciones grficas usando el tiempo de vuelo, el

tiempo de drenaje o los volmenes porosos por pozo son
tiles para verificar la calidad del proceso de upscaling. La
premisa es que para que los sistemas escalados tengan un
comportamiento dinmico similar al modelo de escala fina, los
pozos deben estar drenando volmenes similares y estos
volmenes deben estar conectados en una forma similar. Las
streamlines proveen est informacin.
Comparando los volmenes conectados sobre el tiempo, para

pozos individuales se visualiza una nueva forma de analizar el
comportamiento relativo de los algoritmos de upscaling. La
figura 28 muestra que los patrones streamline de los modelos
escalados reproducen las streamlines a fina escala solo en
ciertas partes del campo. La figura 29 cuantifica esta
discrepancia y compara este con el comportamiento del flujo
fraccional del pozo.
El pozo P2C es un buen ejemplo donde un buen ajuste del

volumen poroso sobre el tiempo tambin lleva a un buen
ajuste sobre el flujo fraccional de agua. Por otra parte, el pozo
P3 no tiene un buen ajuste del volumen poroso a fina escala,
pero el flujo fraccional es aceptable.
Un anlisis del proceso de upscaling basado en volmenes de pozo y

geometra puede al final ser ms relevante que una simple comparacin de los
perfiles de produccin. Para casos reales de campo, es bien conocido que el
upscaling puede eliminar completamente pequeas fallas, barreras de
transmisibilidad
y
otras
caractersticas
geomtricas
cambiando
significativamente los patrones de flujo y los volmenes de drenaje por pozo.
92
En este caso no hay tratamiento de los parmetros de flujo que puedan

remediar este problema, aunque puede ser posible un ajuste de los parmetros
de flujo no fsicos. El problema es geomtrico, y hasta que los patrones de flujo
en el modelo escalado no sean adecuadamente ajustados a lo que se observa
en el modelo a fina escala, se debe planear un acercamiento basado en otros
parmetros.
Thiele* (2001) en su amplio consenso dedicado a la simulacin streamline,
expresa que adems de las ventajas mencionadas anteriormente, un simulador
streamline puede ayudar al ingeniero de yacimientos a verificar las pseudo
permeabilidades relativas por comparacin directa del modelo de simulacin
escalado con el modelo geolgico y geoestadstico.
Figura 28. Streamlines para dos modelos escalados y el modelo de referencia a fina
escala. Un buen upscaling produce patrones de streamlines similares entre el modelo
a fina escala y los modelos escalados.
Escalado 1
Escalado 2
Fina Escala
____________
*

93
Figura 29. Comparacin del volumen asociado con los pozos entre el modelo a fina
escala y los modelos de escala robusta mediante streamlines.
Volumen Poroso por Pozo (0-10%)
Flujo Fraccional de Agua (0-100%)
Escala Fina
Escalado 1
Escalado 2
Tiempo
4.4 TASA DE DISTRIBUCIN Y BALANCE DE PATRONES

La aproximacin streamline puede ayudar en el manejo de yacimientos
suministrando informacin importante, como relaciones inyector/productor y
factores de distribucin para inyectores. Esta informacin es brindada
naturalmente por los modelos streamlines pero no por los simuladores
numricos convencionales. La figura 30 ilustra esta aplicacin sealada por
Datta-Gupta* (2000) en uno de sus estudios. Debido a que cada streamline
est asociada con una tasa de flujo, se puede fcilmente calcular los factores
de distribucin para los inyectores (es decir, qu fraccin del fluido inyectado
est dirigindose a los productores). Esta informacin puede ser til en el
balance de patrones y en el manejo de frentes de inundacin. Una de las
aplicaciones ms importantes de la informacin suministrada por las
streamlines se refiere a los factores de distribucin de pozos y a los volmenes
porosos asociados.
____________
*
Dic. 2000. Pg 70.
94
Figura 30. Diagramas streamlines y factores de distribucin de pozos en el yacimiento.
Factores de Ubicacin de Pozo
Tomada del artculo tcnico: Streamline Simulation. A Technology Update. 2000.
Los factores de distribucin de pozo cuantifican la cantidad de flujo en un pozo

particular debido a la presencia de otros pozos en el sistema. Los volmenes
porosos son los volmenes del yacimiento asociados con cada pozo individual.
4.4.1 Balance de patrones. Tradicionalmente los nmeros de distribucin han

sido determinados usando metodologas empricas. Para Thiele* (2001) y la
gran mayora de autores dedicados al estudio y desarrollo de la teora de lneas
de flujo, las streamlines ofrecen una solucin inmediata y rigurosa a este
problema: por cada inyector, la cantidad de flujo inyectado que soporta
cualquier productor en el campo es conocida exactamente y por consiguiente la
distribucin de fluidos entre inyectores y productores es obtenida
automticamente como parte cualquier simulacin. Esto tambin significa que
las streamlines pueden sealar inmediatamente cualquier prdida de fluidos a
los pozos fuera del patrn.
Las streamlines automticamente permiten determinar la distribucin del flujo

entre pozos mediante la suma del flujo de todas las streamlines asociadas con
un pozo en particular, pares de pozos o grupos de pozos como puede
observarse en la figura 31. Usando esta informacin y el despliegue visual de
las streamlines se obtiene un balance de patrones ms correcto y eficiente que
aquellos obtenidos con tcnicas convencionales.
____________

Fusch, Austria. 3-7 Sep. 2001. Pg 16-17.
95
Adems de los principios bsicos de un simulador streamline tres tipos de

resultado son generados tilmente para el balance de patrones:
El tiempo de vuelo indica el tiempo necesario para que una partcula

inyectada alcance un punto dado en el yacimiento, y es una
representacin instantnea del rea de barrido. Como se puede
apreciar en la figura 32 las zonas rojas son zonas que pueden ser
barridas mientras que las zonas blancas no lo son.
Figura 31. Streamlines y el proceso de balance de patrones. Las tasas son cambiadas
progresivamente para obtener un buen balance.
En este caso los inyectores se ven manteniendo la presin sobre los

costados pero no pueden barrer mucho el petrleo del rea central
hacia el productor. Un patrn de inyectores en el rea central es
aconsejable para hacer ms satisfactorio el incremento del recobro del
petrleo.
96
Figura 32. Influencia del tiempo de vuelo sobre las zonas barridas.
Otros resultados tpicos de un simulador streamline son los volmenes

porosos asociados con cada pozo como se indica en la figura 33, el cual
muestra los volmenes del yacimiento asociados con los diferentes
pozos inyectores. Las regiones son mucho ms extensas que el rea de
barrido debido a que la representacin de los volmenes porosos
asociados con cada pozo es una construccin puramente geomtrica.
Esto no dice nada acerca de cuanto tiempo se tardara en barrer el rea
completa. Este tipo de figuras ayudan a encontrar la zona de influencia
geomtrica de los inyectores.
Figura 33. Influencia del rea ligada a 7 inyectores
97
En la figura 33 se puede notar que los inyectores tienen poca influencia

sobre el rea central y que incrementando la tasa inyectada o
aumentando el tiempo de inyeccin no se obtiene un efecto significativo
en el barrido. La figura tambin indica que el inyector (color gris parte
superior izquierda) slo tiene un pequeo volumen poroso asociado a l,
sugiriendo que este inyector no es eficiente y puede ser ubicado en
cualquier otra parte. Este tipo de resultados es difcil de obtener usando
un simulador convencional de volmenes finitos.
Un tercer tipo de resultados tpicos es el volumen de influencia del

productor como se indica en la figura 34. Cada color indica el volumen
del yacimiento ligado a un productor dado. La figura 34 describe el
volumen del yacimiento asociado con cada productor para seis pozos
verticales mientras la figura 35, para seis pozos horizontales. En este
caso donde todas las fallas estn en comunicacin el esquema de pozo
horizontal tiene el potencial de drenar casi todo el petrleo del campo
pero no se indica el tiempo necesario para lograr este objetivo. Los
volmenes asociados con los dos pozos ubicados al norte son menores
usando pozos horizontales. Esto indica que la implementacin de dos
pozos verticales en esta regin del yacimiento es mucho ms eficiente.
Superponiendo la figura 35 con la grafica OIP ayuda a detectar reas no
drenadas y localizar posiciones para nuevos pozos. Otro aspecto
importante es examinar los volmenes de influencia de los productores
en el tiempo cuando nuevos pozos son abiertos. Este tipo e informacin
es extremadamente til para preparar un primer estimativo del patrn de
pozo y para generar esquema de desarrollo y optimizacin del manejo
de los pozos y la produccin.
Figura 34. Volumen de influencia cuando existen 6 pozos productores verticales.
98
Figura 35. Volumen de influencia cuando existen 6 pozos horizontales productores.
4.4.2 Eficiencia del inyector. Un componente importante en la optimizacin del

comportamiento del campo es poder comparar y clasificar la eficiencia de
inyectores. Inyectores ms eficientes pueden probablemente recibir una alta
porcin del agua de inyeccin disponible que aquellos inyectores menos
eficientes.
Existen diferentes formas para determinar la eficiencia de un inyector. Segn

Thiele* (2001), una de stas es determinar la cantidad de desplazamiento de
petrleo producido como una funcin del volumen inyectado. Esta es
exactamente el tipo de informacin que las streamlines proporcionan
naturalmente. La figura 36 presenta datos reales para una inyeccin avanzada
de agua. Graficando el volumen inyectado contra el desplazamiento de aceite
producido para cada inyector se obtiene la eficiencia de los inyectores para el
campo total. Escogiendo un porcentaje de corte de aceite producido a agua
inyectada, la grfica puede ser dividida en cuatro cuadrantes.
El cuadrante representado por puntos amarillos presenta los inyectores ms
eficientes, es decir, aquellos inyectores que producen la cantidad ms alta de
petrleo por barril de agua inyectada.
____________
*

99
El cuadrante de puntos rojos incluye los pozos menos eficientes, pozos que
inyectan una buena cantidad de agua pero recuperan poco aceite. Estos pozos
son los primeros candidatos a ser cerrados, particularmente en casos donde la
cantidad de agua es limitada o puede ser usada de forma ms eficiente en otra
rea. Finalmente los cuadrantes 2 y 3 representan casos intermedios en los
cuales se pueden plantear nuevas estrategias de desarrollo.
Si las streamlines pueden ayudar a cuantificar la eficiencia de los inyectores,

tambin puede identificar la eficiencia de los productores.
El flujo fraccional de un pozo es un indicador de eficiencia. Un pozo con alto
flujo fraccional es menos eficiente que uno con bajo flujo fraccional.
Figura 36. Agua inyectada vs. desplazamiento de aceite producido.
2500
2000
1500
1000
500
0
0
500
1000
1500
2000
2500
3000
3500
4000
A g ua Inyect ad a ( ST B / D )
Pozos produciendo a altas tasas estn relacionados con saturaciones promedio

de aceite bajas. Similarmente, pozos produciendo a bajas tasas se relacionan
con altas saturaciones de aceite promedio. La figura 37 es un ejemplo de este
tipo de anlisis. En este tipo de anlisis la grfica es dividida usando lneas
diagonales, con la lnea amarilla representando la relacin lineal ms eficiente
entre tasa de aceite y saturacin de aceite promedio. Las lneas verde y roja
son lneas paralelas a la lnea amarilla.
100
Generando un corte para un valor especfico de aceite producido se generan

en la grfica sub-reas en las que se pueden identificar los productores ms
eficientes (amarillos), los menos eficientes (rojos) y aquellos que pueden tener
un potencial para mejorar su produccin mediante procesos de workover o
side-tracking.
Figura 37. Saturacin de aceite promedio por pozo vs. aceite producido.
6000
5000
4000
3000
2000
1000
0
10
15
20
25
30
35
40
45
50
55
S a t u r a c i n d e A c e i t e P r o m e d i o p o r P o z o ( %)
4.5 MODELAMIENTO DE TRAZADORES DE FLUJO Y PROCESOS DE

INYECCIN DE AGUA
Los mtodos streamline han sido ampliamente usados en el pasado en el
modelamiento e interpretacin de pruebas de trazadores.
Tales pruebas, tpicamente involucran la inyeccin de un pequeo bache de

trazador radiactivo. De esta forma, minimizar la dispersin numrica es
esencial para inferir la heterogeneidad sobre las bases de la prueba de
trazador. Los mtodos streamline pueden ser particularmente ventajosos en
este aspecto. Una amplia investigacin sobre este y otros tipos de
desplazamiento en simulacin streamline es presentada por Batycky* (1997) en
su tesis. Pero sin lugar a dudas los desplazamientos preferidos por el mtodo
streamline son los procesos de inyeccin de agua, por las caractersticas de
incompresibilidad de los fluidos involucrados. Desde el inicio de la teora han
sido los procesos por excelencia, y su xito se aprecia en los mltiples estudios
y aplicaciones de campo que se han publicado.
____________
*

Jun 1997. Pg 15-18.
101
4.6 EL PAPEL DE LA SIMULACIN STREAMLINE EN EL PROCESO DE

AJUSTE HISTRICO
El proceso de ajuste histrico juega un papel muy importante en la prediccin
del comportamiento futuro del yacimiento. Su propsito principal es construir un
modelo numrico del yacimiento consistente con todos los datos disponibles
como el conocimiento geolgico y geoestadstico as como los datos de
produccin (cortes de agua, tasas de flujo, presiones, etc).
Esencialmente, el modelo numrico del yacimiento consta de un modelo de grid

tridimensional. Para cada bloque del grid se atribuyen valores de porosidad y
permeabilidad. Por razones prcticas, las distribuciones de estas propiedades
son usualmente aproximadas mediante mtodos Gausianos. Una vez un
modelo de yacimiento es ajustado basado en el conocimiento a priori
geoestadstico, este debe ser modificado para tener en cuenta los datos de
produccin. Este proceso conocido como ajuste histrico es tradicionalmente
dirigido como un problema de optimizacin.
Como primera medida, una funcin objetivo es definida para cuantificar el

desajuste entre los datos actuales y las respuestas paralelas calculadas por el
modelo numrico del yacimiento. Segundo, el modelo del yacimiento es
iterativamente modificado en orden a minimizar esta funcin objetivo.
Un eficiente procedimiento de optimizacin debe poseer las siguientes

cualidades:
Debe ser capaz de modificar los modelos establecidos.
Debe preservar la variabilidad espacial del modelo cualquiera

que sea la modificacin
Debe requerir un nmero reducido de simulaciones previas o

iteraciones.
Tradicionalmente el procedimiento de ajuste histrico perturba los valores de

porosidad y permeabilidad, cualquiera que sea la tcnica de parametrizacin, y
se ejecuta una corrida de simulacin para estimar el impacto sobre la funcin
objetivo.
Entre algunas aproximaciones de ajuste histrico se encuentran las tcnicas de

annealing, los coeficientes de sensibilidad y los mtodos de gradiente.
102
Las tcnicas de coeficientes de sensibilidad calculan la sensibilidad de la

funcin objetivo al cambio de permeabilidad de una celda o un grupo de celdas
y resuelve un sistema inverso que puede ser muy extenso y adems difcil de
construir. Los mtodos de coeficientes de sensibilidad pueden ser tambin
computacionalmente costosos si los coeficientes son evaluados
numricamente mediante mltiples corridas de simulacin.
Vasco* et al. (1999), combin streamlines y una aproximacin de coeficientes
de sensibilidad mientras integra datos dinmicos de produccin.
Ellos emplean streamlines para estimar analticamente los coeficientes de

sensibilidad, acelerando de este modo el procedimiento.
El anlisis streamline permite alinear los fluidos inyectados en el instante de la

ruptura en los pozos productores y de esta forma ajustar la historia de
produccin.
Los coeficientes de sensibilidad

tambin han sido empleados en la
identificacin de la geometra de las caractersticas geolgicas tales como
fallas y canales de flujo.
La tcnica annealing aplicada por Datta-Gupta et al. perturba la permeabilidad

en un grupo de bloques del grid y evala la energa de las funciones objetivo o
el grado de desajuste entre los resultados simulados y los deseados. Este
proceso es estocstico y no garantiza que la perturbacin pueda disminuir el
nivel de energa. La decisin de aceptar o no la perturbacin, est basada en el
cambio de energa causada por esta perturbacin. Perturbaciones que
incrementan el grado de desajuste son aceptadas con una frecuencia que
disminuyen con el incremento en el error. Usualmente se requieren muchas
iteraciones para obtener una solucin aceptable, aumentando con este aspecto
los costos de cmputo.
Muchos de los procedimientos de optimizacin que se usan estn basados en
los mtodos de gradiente. Estos mtodos involucran el clculo de las derivadas
de la funcin objetivo relativas a los parmetros desconocidos.
Estos mtodos pueden ser directos o indirectos. Los directos buscan el ptimo
analticamente fijando el gradiente de la funcin objetivo en cero. Los mtodos
indirectos buscan el ptimo saltando a lo largo de la funcin objetivo
movindose en la direccin del gradiente.
____________
VASCO, D.W, YOON, S., DATTA-GUPTA, A. Integrating Dynamic Data Into High-Resolution Reservoir
Models Using Streamline-Based Analytic Sensitivity Coefficients. Texas A&M U. Dic. 1999.
103
Ambos mtodos requieren de la informacin del gradiente con base en la

expansin de series de Taylor y dependiendo del orden de las derivadas que
se considere pueden ser de primero o segundo orden. Los de primer orden
consideran despreciables los trminos despus de la primera derivada usando
solamente el vector jacobiano. Los de segundo orden utilizan los trminos de
la segunda derivada y requieren el clculo de la matriz Hessiana.
Una limitacin de estos mtodos es que pueden encontrar ptimos locales y

que se requiere la existencia de las derivativas. Un aspecto crtico es el clculo
de los coeficientes de sensibilidad a los parmetros de la funcin objetivo. Es
decir, la variacin en el resultado del problema directo motivado por cambios en
los parmetros del modelo.
La tcnica streamline es un mtodo aproximado de simulacin de yacimientos

que tambin ha desarrollado tcnicas especiales para la optimizacin del ajuste
histrico. Estas tcnicas se describen a continuacin con ms detalle.
4.6.1 Aproximaciones Streamline para ajuste histrico.

Una de las
aplicaciones ms prometedoras de las streamlines est en el rea del ajuste
histrico. Debido a que las streamlines estn unidas a los pozos y al mismo
tiempo marcan reas de alta influencia en el yacimiento, se debe suponer que
debe existir informacin escondida a lo largo de las streamlines que pueden
ayudar en el proceso de ajuste histrico. Adicionalmente, el ajuste histrico del
comportamiento de un yacimiento es un proceso no-lineal que requiere de
muchas simulaciones previas, hecho que puede conseguirse rpida y
eficazmente con los simuladores streamline por las caractersticas intrnsecas
de los mismos. Entre las aproximaciones recientemente desarrolladas para el
proceso de ajuste histrico se encuentran:
o Aproximacin de Yuandong Wang* et al. Una nueva aproximacin de ajuste
histrico basada en streamlines fue propuesta por Y. Wang y A:R: Kovscek
(2000). Ellos proponen un mtodo de dos pasos fundamentales para ajustar los
datos dinmicos de produccin y determinar la heterogeneidad del yacimiento.
El primer paso es modificar la distribucin de permeabilidad al nivel streamline,

basado en la diferencia entre los resultados de la simulacin y los datos de
campo para corte de agua, cada de presin y tasa de flujo.
____________
WANG, YUANDONG, KOVSCEK, ANTHONY R. A Streamline Approach For History-Matching

Production Data . Improved Oil Recovery Symposium held in Tulsa, Oklahoma, 35 Abril 2000. SPE N
59370
104
Ajustando la curva de flujo fraccional a travs de la manipulacin de la

permeabilidad del campo, se trata de capturar la heterogeneidad del
yacimiento.
El segundo paso en esta aproximacin streamline es trasladar la modificacin

de la permeabilidad streamline hacia los grid-blocks. De esta forma, se realiza
la simulacin del flujo para verificar el ajuste. El anterior proceso es de carcter
iterativo hasta lograr la convergencia.
Este mtodo propuesto es similar al mtodo desarrollado por Vasco et. Al

(1999), sin embargo, en esta aproximacin, no se calculan los coeficientes de
sensibilidad. Cuando la curva de flujo fraccional del resultado de la simulacin
para una permeabilidad dada no ajusta con la curva de corte de agua del
campo, se deducen las streamlines responsables de la diferencia entre las
curvas de flujo fraccional. Posteriormente, basndose en la relacin entre el
tiempo de ruptura de las streamlines y la permeabilidad efectiva, se puede
calcular una modificacin de la permeabilidad efectiva a lo largo de las
streamlines para ajustar los datos de produccin.
El planteamiento matemtico de esta tcnica de ajuste es el siguiente:
La funcin objetivo como ya se conoce indica el error entre la comparacin de

los resultados de la simulacin y los datos de campo. Este error puede
expresarse como:
Np
=
n =1
(E + E
t ,n
p ,n
+ E q ,n ... (120)
Donde n representa el subndice del productor, Np es el nmero total de

productores en el estudio, Et,n,Ep,n,Eq,n son los errores en el tiempo de ruptura
adimensional de las streamlines individuales, en la presin y en la tasa de flujo
en el productor i respectivamente.
Para ajustar la presin y la tasa de flujo, se modifica la permeabilidad efectiva

de toda la regin. Sin embargo, para capturar la heterogeneidad, se necesita
resolver un sistema inverso para ajustar la curva de flujo fraccional.
Continuando con el modelo, el grado de desacuerdo entre los resultados de
referencia y los resultados del ajuste histrico se calcula mediante la siguiente
ecuacin:
105
Et =
1 N sl
Et
sl i =1
.... (121)
BT , i
Donde N sl representa el nmero de streamlines conectadas con el productor.

El error E se refiere al tiempo de ruptura de la streamline i definido de la
t
siguiente forma:
BT , i
Et
= t D , BT ,i t D , BT ,i .. (122)
C
BT , i
En donde
C
D , BT ,i
R
D , BT ,i
son los tiempos adimensionales de ruptura calculados y
de referencia respectivamente para el i-simo streamline.

Sistema inverso. Esta es la parte ms importante de este planteamiento. La
modificacin de la permeabilidad altera el tiempo de ruptura no solo del i-simo
streamline, sino tambin de otras streamlines. Por consiguiente todas las
streamlines deben ser consideradas y se debe resolver un sistema de
ecuaciones.
El sistema a resolver para cada productor es el siguiente:
a11 a12
a21 a22
a31 a32
aN1 aN2
a13
a23
a33
aN3
a1N
a2N
a3N
aNN
k1
k2
k3
..
kN
Et,BT,1
Et,BT,2
- Et,BT,3
Et,BT,N
(123)
Donde kj es la modificacin de la permeabilidad efectiva a lo largo de la

streamline j requerida para obtener un ajuste, y aij es la sensitividad del tiempo
de ruptura del i-simo streamline a la permeabilidad efectiva del j-simo
streamline.
106
Estas derivadas son definidas como:
ij
t
k
D , BT ,i
.. (124)
Debido a que las streamlines no se encuentran comunicadas, las derivadas

pueden ser aproximadas aplicando el mtodo de Dykstra-Parsons para estratos
no comunicados de diferente longitud. El mtodo relaciona el tiempo de ruptura
de diferentes estratos con la permeabilidad efectiva de cada estrato.
El tiempo de ruptura para una streamline i es calculada por:
D , BT , i
N sl
A l x D , k
k =1
k
=
N sl
A l
k =1
.. (125)
Donde l es la longitud de la streamline, xD,k es la posicin adimensional del

frente de la fase desplazada del k-simo streamline cuando el i-simo
streamline llega a la ruptura.
Los smbolos k y A k representan la porosidad promedio y el rea de seccin

transversal promedio de la streamline k respectivamente y pueden ser definidos
como:
A = A( x D )d xD
.... (126)
= ( x D )d x
..... (127)
107
Teniendo en cuenta las anteriores ecuaciones se puede expresar la ecuacin

125 de la siguiente forma simplificada:
D , BT ,i
N sl
V
k =1
D ,k
X D ,k ..... (128)
Aplicando la regla de la cadena

siguiente planteamiento:
tD,BT,i
k j
Nsl
tD,BT,i xD,k
k =1
xD,k k j
se evala la ecuacin
Nsl
xD,k
k =1
k j
= VD,k
124 mediante el
... (129)
El mtodo de Dykstra-Parsons (1950) proporciona valores de xD,k en trminos

de la permeabilidad efectiva de todas las streamline.
Las formulas para calcular xD,k para relaciones de movilidad de uno y menores
a uno son diferentes la una de la otra.
Para una relacin de movilidad de uno, la presin del campo as como la

distribucin de streamlines permanecen invariables a travs de los procesos de
desplazamiento para condiciones lmites constantes.
Cuando la ruptura ocurre en una streamline i , la posicin del frente de la

streamline k es calculada por:
x D , k = c ik
kk
.. (130)
Donde cik es una constante relacionada con la longitud de las streamlines i y k .
108
Aplicando esta definicin y llegando a la derivacin parcial se llega a la

siguiente solucin
t D ,i
=
k j
x D ,k
=
k j
N sl
V
k =1
D ,k
ijV
D , j
1
2
ki
N sl
k =1, k i
ik
VD ,k k k
Si i=j
.. (131)
si ij
Este procedimiento pude ser repetido para relaciones de movilidades diferentes

de uno obtenindose resultados estndar de Dykstra-Parsons, donde xD,k es
una funcin de la permeabilidad de los streamline i y k y del valor final de la
relacin de movilidad.
La metodologa a seguir en esta aproximacin involucra los siguientes pasos:
1. Obtener una permeabilidad inicial del campo mediante estimaciones o
prcticas geoestadsticas.
2. Ejecutar una primera simulacin sobre la permeabilidad inicial del
campo. Verificar si los resultados de la simulacin ajustan con los datos
de campo incluyendo la curva de flujo fraccional en los pozos
productores, la tasa de flujo y la presin. Si no se presenta ajuste, se
debe modificar la permeabilidad de acuerdo a los pasos siguientes.
3. Trabajar sobre las streamlines. Calcular el tiempo de vuelo o el volumen
poroso asociado para todas las streamlines. Clasificar las streamlines en
orden ascendente segn el volumen poroso.
4. Calcular la diferencia en el flujo fraccional, la tasa de flujo y la presin
entre el resultado de la simulacin y los datos de referencia. Relacionar
diferencias en la curva de flujo fraccional para el correspondiente
streamline
5. Resolver el sistema descrito en la ecuacin 123 para calcular la
modificacin de la permeabilidad efectiva de cada streamline requerida
para un buen ajuste de los datos disponibles.
6. Modificar la distribucin de permeabilidad a lo largo de cada lnea de
flujo.
7. Repetir los pasos 2 al 6 hasta obtener un ajuste satisfactorio.
109
Las aplicaciones publicadas lo muestran como un mtodo fuerte que converge

rpidamente. Sin embargo, presenta la limitacin que la variacin de
permeabilidad se hace en la misma proporcin a lo largo de la lnea de flujo,
por lo que puede no hacer honor a la informacin geolgica. Una forma de
anlisis del proceso y resultados grficos se presenta en la figura 38.
Figura 38: Ajuste de la curva de flujo fraccional

1
0.9
0.8
Flujo Fraccional
0.7
0.6
Aproximacin inicial
Segunda modificacin
0.5
Tercera modificacin
Primera modificacin
0.4
0.3
0.2
0.1
0
0
0.2
0.4
0.6
0.8
1.2
1.4
1.6
Tiem po Adim ensional
Tomada del artculo tcnico SPE 59370
o Procedimiento combinado de geoestadstica y simulacin streamline.

Aproximacin de Mickaele Le Ravalec-Dupin* et al (2002).
La caracterstica
esencial de esta aproximacin es que una simulacin del flujo no se asocia a
una sola perturbacin del modelo del yacimiento, sino a un grupo de
perturbaciones del modelo.
Esta ventaja es una consecuencia de los dos pasos estructurados del

procedimiento de ajuste.
____________
DUPIN, RAVALEC, DARRYL, MICKAELE H., FENWICK. A Combined Geostatistical and StreamlineBased History-Matching Procedure. SPE Annual Technical Conference. And Exhibition. Texas, 29 sep. 2
oct. 2002.
110
El primer paso hace referencia a la identificacin de grupos de streamlines y a

la estimacin de correcciones o perturbaciones para agregar a las
permeabilidades efectivas de estos grupos de streamlines en orden a reducir el
desajuste entre los flujos fraccionales medidos y computados. El segundo paso
es dedicado a la propagacin de la correccin estimada a las permeabilidades
absolutas del grid block. Durante el segundo paso, la distribucin de las
streamlines se asume como constante. Se define una funcin objetivo
intermedia para cuantificar la discrepancia entre las permeabilidades efectivas
de las streamlines deseadas y algunas calculadas para la permeabilidad
absoluta del campo dado. En este punto, se implementa un mtodo basado en
el gradiente para minimizar la funcin objetivo intermedia y para modificar las
permeabilidades absolutas del grid block usando el mtodo de deformacin
gradual. Desarrollando esta metodologa se obtiene la permeabilidad absoluta
del campo para la cual las permeabilidades efectivas de las streamlines son
cercanas a los valores deseados. Al mismo tiempo, la permeabilidad del campo
todava respeta la variabilidad espacial del modelo. Durante el segundo paso,
no se requieren simulaciones adicionales. El algoritmo del mtodo puede ser
planteado de la siguiente forma:
1. Dibujar un modelo de yacimiento inicial respetando el variograma

experimental.
2. Correr una simulacin de flujo para identificar las streamlines;
identificar grupos de streamlines y calcular las correcciones para
agregar a las permeabilidades efectivas de estos grupos en orden a
reducir el desajuste entre los flujos fraccionales medidos y
computados.
3. Propagar la correccin estimada al grid block de permeabilidades
absolutas. Este paso involucra una perturbacin iterativa de las
permeabilidades absolutas sobre la base del mtodo de deformacin
gradual. Este paso no requiere ninguna simulacin de flujo.
4. Volver al paso 2, si el ajuste de los flujos fraccionales no es
satisfactorio.
El planteamiento matemtico del mtodo es el siguiente:
Correccin de las permeabilidades efectivas de las streamlines.
Esta
seccin se enfoca en la determinacin de la correccin para agregar a las
permeabilidades efectivas a lo largo de las streamlines. El principal propsito
de este primer paso es ajustar las tasas de fluido inyectado medidas a la
produccin de los pozos. Este dato es especialmente sensible a la
heterogeneidad del yacimiento y puede ser aprovechado para determinar vas
de flujo preferenciales y barreras de flujo.
111
Como el error en las tasas de flujo est relacionado con el error en las
permeabilidades efectivas de las streamlines, el ajuste debe ser mejorado
variando las permeabilidades efectivas de las streamlines. La consideracin de
grupos de streamlines es la base fundamental para aplicar este procedimiento
a modelos de yacimientos con un extenso nmero bloques.
La permeabilidad efectiva de un grupo de streamlines es un promedio armnico

de las permeabilidades absolutas de los grid blocks interceptados por las
streamlines de este grupo. Est correccin es expresada como:
N sli N k
eff
i
q
N N q
k
k =1 j =1
sli
k, j
k, j
k, j
k, j
k =1 j =1
..... (132)
Donde i expresa el i-simo grupo de streamlines. Nsli es el nmero de

streamlines en el i-simo grupo i y Nk es el nmero de grid blocks interceptados
por el k-simo streamline de este grupo.
qkj es la tasa de flujo para el k-simo streamline y el j-simo grid block

interceptado. kj es el tiempo de vuelo para el k-simo streamline a travs del
j-simo grid block interceptado.
Est ecuacin es una ecuacin emprica basada en el promedio volumtrico a

lo largo de las streamlines.
Para mejorar el ajuste del flujo fraccional, la permeabilidad efectiva para este
grupo debe ser igual a:
eff
i , deseada
eff
i, deseada
t
t
i ,calculada
eff
i , computada
.. (133)
i , medida
eff
i, computada
son las permeabilidades efectivas deseadas y computadas
para el i-simo grupo de streamlines.

modelo considerado mientras
eff
i, deseada
eff
i,computada
depende de la permeabilidad del
es el valor deseado para mejorar el
112
ajuste.
i, medido
i,computado
son los tiempos de ruptura asociados con el incremento
q para las tasas de flujo medidas y computadas.

i
Propagacin de la correccin deseada a las permeabilidades del grid block.

Una vez las permeabilidades efectivas deseadas son conocidas, est
correccin debe ser propagada a las permeabilidades absolutas del grid block.
Para hacer esto se plantea un problema de optimizacin intermedio basado en

el gradiente, dirigido a minimizar la siguiente funcin objetivo intermedia:
1 N
FOI ( ) =
2 i =1
(K
eff
i , deseada
K i , computada ( ) .. (134)
eff
N es el nmero de grupos de streamlines. Las incgnitas de este problema de

optimizacin no lineal son todas las permeabilidades absolutas del grid block
reunidas en el vector . Introduciendo la tcnica de parametrizacin gradual
dentro del procedimiento de optimizacin se puede reducir el nmero de
parmetros desconocidos mientras se garantiza la preservacin de la
variabilidad espacial del modelo. La funcin objetivo intermedia se define
nuevamente como:
FOI ( ) =
1 N
2 i =1
(K
K i , computada ( ) ... (135)

i , deseada
eff
eff
El vector est constituido por los parmetros de la deformacin gradual. Para

disear un eficiente proceso de minimizacin, se necesitan las derivadas de la
funcin objetivo intermedia relativas a los parmetros de deformacin. Para
esto se hace uso del clculo del gradiente.
Clculo del gradiente. En el planteamiento del clculo del gradiente se
consideran dos grupos de desviaciones normales independientes llamadas
interferencias Gaussianas: interferencia Gaussiana inicial z1 e Interferencia
Gaussiana complementaria z2. Un nuevo grupo de desviaciones normales
independientes z es obtenido de la siguiente relacin:
z ( ) = z1 cos( ) + z 2 sen( ) ..... (136)

113
Variando los parmetros de deformacin se puede describir una cadena de

interferencia Gaussianas.
Cuando se consideran dos interferencias Gaussianas la derivada de z relativa

al parmetro de deformacin, se expresa como:
z
= z1 sen ( ) + z 2 cos ( ) ..... (137)
Entonces el grupo de desviaciones z es convertido en una realizacin

Gaussiana y con una funcin de covarianza C del siguiente producto de
convolucin:
y = f * z .. (138)
La funcin f resulta de la descomposicin de la funcin de covarianza como
C=f f, donde f * es la transpuesta de f . La derivada de y esta dada por:
= f*
z
...... (139)
La siguiente relacin permite convertir y en una realizacin normal w :
w = m + y .. .. (140)
Y las derivadas son:
w
y
=
...... (141)
Cuando los datos fuertes (valores de permeabilidades medidos a las

ubicaciones dadas) estn disponibles, stos pueden integrarse dentro la
realizacin w empleando tcnicas de Kriging.
114
Una realizacin wc es obtenida de:
wc = wdK + (w wK )
... (142)
El estimativo para un Kriging dual es formulado como:
wK (x ) = pi C ( x xi ) + m ..... (143)
i =1
En donde C es la funcin de covarianza, n es el nmero de observaciones

puntuales y x es la ubicacin.
La derivada de wc es obtenida de:
w
wc
(x ) =
(x ) wK (x )
n
wK
(x ) = pi C (x xi )
i =1
pi
C (x
n
j =1
xi ) =
w
(x ) ....... (144)
wc es transformada en un registro normal de permeabilidades de campo k

mediante :
k = Exp(wc )
.. (145)
El gradiente de permeabilidad est dado por:
k wc
=
k ... (146)
115
Finalmente los gradientes de permeabilidad efectiva son:

N sli N k q k , j k , j k j
K ieff
= K ieff
k =1 j =1
..... (147)
k 2j
N sli N k qk , j k , j
k =1 j =1
El procedimiento para realizar el problema intermedio de optimizacin es el

siguiente:
1. Asignar un parmetro de deformacin a cada grupo de streamlines.

2. Trazar aleatoriamente una interferencia gaussiana inicial.
3. Trazar una interferencia gaussiana complementaria independiente.
4. Combinar las anteriores interferencias gaussianas usando el mtodo de
deformacin gradual y calcular los gradientes asociados.
5. Convertir las interferencias gaussianas resultantes a la permeabilidad
del campo y calcular los gradientes asociados.
6. Determinar las permeabilidades efectivas de los grupos de streamlines
asumiendo que sus geometras son invariables. Calcular los gradientes
asociados.
7. Estimar la perturbacin para aplicar los parmetros de deformacin a fin
de reducir la funcin objetivo intermedia.
8. En este punto hay varias alternativas:
Si la FOI no es lo suficientemente pequea y no converge,

actualizar los parmetros de deformacin gradual e ir al paso 4.
Si la FOI no es muy pequea, pero converge, actualizar la

interferencia gaussiana inicial e ir al paso 3.
Si no se presentan ninguna de las anteriores alternativas, se debe

detener el procedimiento.
El diagrama de flujo para esta aproximacin de ajuste histrico se presenta a

continuacin.
116
Figura 39. Diagrama de flujo para el mtodo propuesto por Mickaele Le Ravalec Dupin
y Darryl Fenwick
Modelo
Permeabilidad
Correccin para las

Permeabilidades
Efectivas a lo largo
de los Streamlines
Simulacin del
Flujo de
Fluidos
Paso 1: Estimacin de la Correccin

Interferencia
Gaussiana
Inicial
Actualizacin Interferencia Gaussiana y trazado aleatorio de IGC
Interferencia
Gaussiana
Complementaria
Deformacin
Gradual
Modelo
Permeabilidad
Funcin
Objetivo
Intermedia
Optimizacin PDG
Paso 2: Optimizacin Intermedia-Propagacin de la Correccin
Tomada de DUPIN, A Combined Geostatistical and Streamline-Based HistoryMatching Procedure.
o Procedimiento de ajuste histrico asistido (AHA). Aproximacin de Emanuel*

et al (1998). El ajuste histrico asistido (AHA) es una metodologa por medio
de la cual datos 3D de las streamlines son usados para ayudar en el proceso
de ajuste histrico.
En general, el modelo de simulacin en consideracin es un modelo

convencional de diferencias finitas aunque el mtodo tambin aplica modelos
de desplazamiento streamline.
El proceso AHA est basado en la suposicin de que el ajuste histrico es

alcanzado por alteracin de las propiedades geolgicas a lo largo de las vas
de flujo que conectan pozos productores con sus fuentes de flujo.
Las fuentes pueden ser pozos inyectores de agua, gas, acuferos o capas de
gas, sin embargo, el mecanismo de conduccin debe ser un desplazamiento
que se d a lo largo de una va definible.
____________
EMANUEL, A.S, AND MILLIKEN, W.J. Application of 3D Streamline Simulation to Assist History
Matching. SPE 49000. New Orlens. Oct 1998.
117
En grids muy robustos (coarse grids), un solo grid block puede interceptar el
flujo de varios pozos y en este caso la aplicacin del mtodo de AHA no es
muy aconsejable.
Una vez se determinan las vas de las streamlines, se identifican los grid blocks
a travs de los cuales dichas streamlines se mueven y entonces se asignan al
productor en el cual termina la streamline. Se debe usar un nmero suficiente
de streamlines para obtener un buen alcance de los grid blocks del modelo.
Algunos grid blocks pueden ser interceptados por ms de una streamline y en
algunos casos mltiples streamlines a travs de un grid block particular,
pueden terminar en diferentes pozos. En estos casos se debe hacer una
asignacin de acuerdo con reglas ya determinadas. Las dos reglas aplicadas
en estos casos son:
1. Asignacin fundamentada en una base de prioridad predeterminada

(ejemplo, tasa de produccin).
2. Asignacin al pozo con el mayor nmero de streamlines o la streamline
ms extensa en el grid block.
La tcnica no trabaja bien para mecanismos de expansin de fluidos,

expansin de roca o aumento de gas en solucin debido a que estos
mecanismos no definen muy bien las vas de flujo.
Ajuste del modelo geolgico. El simulador streamline es usado para
determinar las vas de flujo y asignar a cada productor los grid blocks en el
modelo que principalmente afectan el pozo. Una vez estos grid blocks son
identificados, se hace necesario determinar e implementar los cambios para
lograr el ajuste histrico deseado.
Permeabilidad. La permeabilidad del bloque es probablemente el parmetro
que con mayor frecuencia se cambia durante el ajuste histrico. Esto es lgico
ya que la permeabilidad generalmente tiene el efecto ms fuerte sobre el flujo
de fluidos y al mismo tiempo es el parmetro menos conocido. Los
multiplicadores convencionales para las permeabilidades horizontales y
verticales estn implementados en el mtodo. Los multiplicadores horizontales
y verticales de permeabilidad pueden ser diferentes, consecuentemente, la
relacin de permeabilidades K v / K h puede ser asignada sobre una base pozo
por pozo.
En adicin a los multiplicadores de permeabilidad, el mtodo incorpora cambios

sistemticos a la heterogeneidad de la distribucin de permeabilidad.
118
Construido sobre una aproximacin 2D de Emanuel y Milliken, el mtodo usa

parmetros de Dykstra-Parsons como una medida de la heterogeneidad. Para
cada columna vertical a lo largo de la va de flujo del pozo, se analizan los grid
blocks conectados al pozo y se determina un coeficiente caracterstico de
Dykstra-Parsons. Las permeabilidades son entonces ajustadas para hacer
coincidir el nuevo coeficiente de Dykstra-Parsons. La figura 40 presenta una
ilustracin de la tcnica de re-normalizacin empleada.
En lugar de multiplicar la permeabilidad por un factor especfico, est puede ser

alterada volviendo a derivar la transformada porosidad-permeabilidad.
Especficamente el mtodo calcula la distribucin de permeabilidad-porosidad

desde los datos de entrada y permite al usuario reconstruir la distribucin de
permeabilidad de entrada usando un estimativo diferente (alto o bajo) de la
transformada de permeabilidad-porosidad. La tcnica minimiza el peligro
inherente en la aplicacin de un factor de multiplicacin que genere valores de
permeabilidad inconsistentes con la porosidad del yacimiento.
Figura 40: Ilustracin de la tcnica de renormalizacin basada en la caracterizacin de

Dykstra-Parsons.
120
100
Rest ablecimient o de la
variacin de permeabilidad
80
60
Int erpolacin Geoest ad st ica
40
20
0
0
10
20
30
40
50
60
70
80
90
Fr e c u e nc i a A c u l u l a t i v a
Tomada de EMANUEL, A.S, AND MILLIKEN, W.J. Application of 3D Streamline

Simulation to Assist History Matching. SPE 49000.
119
Porosidad. El mtodo tambin ha sido diseado para realizar cambios a la

porosidad o al volumen poroso del modelo. Dos tipos de cambios pueden ser
hechos:
Cambios que preservan el volumen

Cambios que no preservan el volumen.
En el primer caso, se asignan los multiplicadores pozo por pozo para

incrementar o disminuir la porosidad. Dichos cambios pueden afectar el
volumen poroso total del modelo.
En el segundo caso, se implementa una tcnica de re-normalizacin. En este

caso, el volumen poroso se incrementa a lo largo de las streamlines
conectadas a ciertos productores, pero el volumen poroso total es preservado
por re-normalizacin.
Asignaciones regionales.
Finalmente, la tcnica AHA es idealmente
conveniente para asignar regiones de propiedades al modelo teniendo en
cuenta cada pozo. En este mtodo, todos los grid blocks asociados con un
pozo en particular son simplemente asignados a la regin designada por ese
pozo.
La decisin de determinar el tipo y grado de cambios que deben ser hechos al

modelo todava dependen de la experiencia del ingeniero de simulacin. Sin
embargo, incorporando la metodologa streamline, el proceso del ajuste
histrico se ve ampliamente simplificado. La experiencia ha demostrado que se
pueden alcanzar ajustes satisfactorios en menor tiempo de lo que tardara
realizar el proceso por mtodos convencionales. Un ejemplo de esto puede
apreciarse en la figura 41 en donde la lnea blanca representa el corte histrico
de agua para el periodo de inters, la lnea roja representa el ajuste histrico
realizado por mtodos convencionales y la lnea naranja el nuevo mtodo de
AHA. En ella puede observarse que se obtienen excelentes resultados
mediante la aplicacin del AHA, an sobre los resultados del mtodo
tradicional.
Guas para la seleccin de los parmetros.
Multiplicadores de grid block. Los multiplicadores convencionales para

permeabilidad horizontal incrementan o disminuyen la conectividad entre
la fuente y el productor. Un incremento en la permeabilidad normalmente
resulta en tiempos tempranos de ruptura del fluido desplazado y un gran
mantenimiento de la presin.
120
Figura 41. Resultados del ajuste histrico para el corte de agua empleando un mtodo
convencional y el mtodo de AHA.
Multiplicadores para la porosidad incrementan o disminuyen la cantidad

relativa de volumen de yacimiento entre el productor y la fuente, y
ocasionan un efecto opuesto sobre el tiempo de ruptura.
En contraste al efecto directo de la permeabilidad horizontal, la

permeabilidad vertical influye fuertemente en la importancia de las
fuerzas viscosas y gravitatorias en el proceso de desplazamiento.
Incrementos en la permeabilidad vertical dan como resultado un gran
efecto de segregacin gravitacional. Esto generalmente mejora la
eficiencia de barrido si el yacimiento tiene una pendiente considerable.
Disminuciones en la permeabilidad vertical mejora el flujo estratificado y
ste tiende a incrementar los efectos de heterogeneidad y digitacin.
Transformada porosidad-permeabilidad.
Grandes cambios a la
permeabilidad pueden generar pares de valores de porosidadpermeabilidad que tienden al rango mximo de los valores medidos. La
asignacin de una alta permeabilidad a un grid block de baja porosidad
puede tambin inducir a inestabilidades numricas en la solucin de
saturacin. Un cambio a la transformada de porosidad-permeabilidad
reduce la posibilidad de que ocurran este tipo de problemas.
121
ndices de heterogeneidad.
Para modelos de yacimientos
caracterizados por permeabilidades correlacionadas, incrementos en la
heterogeneidad (incrementos en el coeficiente de Dykstra-Parsons)
inducen a rupturas tempranas y reducen la eficiencia de barrido. Para
distribuciones de permeabilidades aleatorias, incrementos en la
heterogeneidad demoran el tiempo de ruptura ya que se crean vas
tortuosas para el desplazamiento del fluido.
122
5. SIMULADORES COMERCIALES STREAMLINE.
Los principales
disponibles son:
programas
de
simulacin
streamline
comercialmente
3DSL de StreamSim Technologies.Inc.

Frontsim de Schlumberger
El modelo Streamsim fue desarrollado por la Universidad de Stanford y es

conocido como 3DSL. El modelo de streamline de Schlumberger hace parte del
conocido programa de simulacin Eclipse.
A continuacin se presentaran algunas generalidades y caractersticas

particulares de cada uno de stos simuladores. No obstante, se profundizara
un poco ms en el simulador Frontsim de Schlumberger por ser el simulador
streamline con el que se cuenta en la facultad.
5.1 SIMULADOR 3DSL STREAMSIM TECHNOLOGIES.INC.

3DSL* es el simulador streamline ms avanzado de la industria y representa
una herramienta novedosa en la simulacin de yacimientos altamente
heterogeneos. 3DSL permite la simulacin rpida de modelos ampliamente
extensos y completos en computadores personales.
Es particularmente eficiente en el modelamiento y estudio de diversos tipos de

desplazamientos y est diseado para simular modelos constituidos por
millones de celdas con miles de pozos y muchos aos de historia,
desarrollando la simulacin en tiempos de ejecucin aceptables y reduciendo
los costos en PC. 3DSL es un simulador de aceite negro que trabaja con
sistemas trifsicos en tres dimensiones y se extiende a sistemas compresibles
de aceite negro.
____________
*
MANUAL SIMULADOR 3DSL. STREAMSIM TECNOLOGIES INC. Version 2.20 Agosto 2004.
123
5.1.1 Caractersticas bsicas del simulador 3DSL. Como ya se mencion, el

modelo implementado por 3DSL es un modelo de tres fases, tridimensional que
tiene en cuenta la gravedad y permite cualquier combinacin de
compresibilidad entre las tres fases presentes. El gas puede estar presente en
la fase aceite, y por otra parte el gas libre y el agua son consideradas fases
independientes. Un sin nmero de desplazamientos pueden ser modelados
mediante la formulacin 3DSL, entre los que se pueden mencionar:
Flujo Trazador
Flujo Miscible
Flujo Inmiscible
Flujo Trifsico con Aceite Vivo
Sistemas de Porosidad Dual
Figura 42. Bosquejo del entorno del simulador 3DSL.
Tomada del simulador 3DSL.
Dependiendo de la complejidad geolgica del yacimiento, el simulador 3DSL

trabaja sistemas de punto centrado, punto distribuido, grids estructurales o
configuraciones estrictamente cartesianas con bloques activos e inactivos.
124
Trabaja sobre curvas de permeabilidades relativas en tres fases y permite

definir por regiones el yacimiento. Tiene en cuenta aspectos intrnsecos a la
configuracin de cada uno de los pozos como la presencia de pozos verticales,
horizontales, desviados y completamientos o trabajos de work-over practicados
a los pozos pertenecientes al dominio de inters.
Una caracterstica importante del simulador 3DSL es su produccin de nuevos

datos adicionales asociados con las streamlines que pueden ser fundamentales
en el estudio y prediccin del comportamiento del campo. Entre stos datos se
encuentran los factores de ubicacin de pozos y los volmenes porosos
asociados a cada pozo. Los factores de ubicacin de pozos identifican el
desplazamiento de los fluidos producidos debido a los volmenes inyectados y
son usados por el simulador para construir grficas eficientes y mapas de
patrones de flujo como puede observarse en la figura 43. Los volmenes
porosos relacionados con cada pozo cuantifican la presencia de fluidos sobre
cada dominio de los pozos y la extensin areal de las regiones de inyeccin y
produccin.
En cuanto a las caractersticas operativas del simulador 3DSL es sencillo de

operar y su programacin est basada en la utilizacin de keywords. Es
compatible con otros formatos de simuladores como Gocad y Eclipse.
Para un anlisis e informacin ms detallada sobre este simulador se puede
consultar en su pgina oficial www.streamsim.com en donde se dispone del
manual para usuario general, demos gratuitos y de una amplia variedad de
artculos tcnicos asociados con el tema.
Figura 43. Representaciones grficas construidas a partir de los datos arrojados por el
simulador 3DSL.
Mapa de
Patrones de
Flujo
Volumen Poroso Por
Pozo
Tomada del simulador 3DSL.
125
5.2 SIMULADOR FRONTSIM. ECLIPSE- SCHLUMBERGER

Este simulador streamline es el ms conocido y el ms usado en el
modelamiento de yacimientos por el gran dominio que se tiene sobre la
plataforma Eclipse en el medio ingenieril. Su caracterstica fundamental como
ocurre con todos los simuladores streamline es su habilidad para poder
visualizar las vas de flujo en el yacimiento y las relaciones entre pozos
productores e inyectores. Su programacin y ejecucin guarda el ncleo bsico
de los simuladores Eclipse, es decir el planteamiento de la simulacin se
desarrolla mediante keywords especficos similares a los empleados en Eclipse
100.
FRONTSIM* es un simulador de yacimientos basado en una formulacin tipo
IMPES (implcita la presin y explcita la saturacin) y en un concepto de
seguimiento del frente para la determinacin de las saturaciones.
Como ya se mencion, la ecuacin de presin es resuelta implcitamente a

travs de un mtodo de diferencias finitas y dependiendo de la estructura del
grid se pueden emplear volmenes de control o formulaciones estndar sobre
el mtodo.
Por otra parte una vez definida la distribucin de presiones, las saturaciones
para cada una de las streamlines a determinados tiempos son encontradas
mediante un procedimiento de seguimiento del frente sobre una escala macro
en el yacimiento.
FRONTSIM conserva la estructura matemtica y procedimental ya planteada

en el capitulo 3 en donde cada streamline trazado representa una tasa
volumtrica constante y acta como un espacio unidimensional para hallar
saturaciones.
5.2.1 Caractersticas bsicas de Frontsim.
Trabaja especialmente sistemas bifsicos, aunque existe una versin

disponible para trabajar sistemas de tres fases.
Es ampliamente aplicado en modelamiento de procesos de trazadores y

flujo inmiscible.
Trabaja sistemas compresibles e incompresibles.
____________
*
MANUAL ECLIPSE FRONTSIM GEOQUEST.SCHLUMBERGER. 2001.
126
Tiene en cuenta efectos de gravedad.
Ampliamente flexible en el manejo de geometras de grid de punto

distribuido.
Por su naturaleza no presenta dispersin numrica ni efectos del grid.

Trabaja un ptimo refinamiento local del grid.
Tiene en cuenta aspectos altamente heterogneos y de anisotropa

geolgica, as como datos petrofsicos que ayudan a una mejor
descripcin del yacimiento.
Incluye tensores completos de permeabilidades y multiplicadores de

transmisibilidad, lo cual es bsico en el trabajo de grids no ortogonales.
La ecuacin de presin es resuelta mediante mtodos iterativos

precondicionados.
Ecuacin de saturacin resuelta mediante un mtodo de seguimiento del

frente sobre streamlines tridimensionales.
Trabaja tasas de inyeccin o produccin variables controladas por pozo,

grupo o a nivel de campo.
Permite el reclculo automtico o definido por el usuario de las

streamlines o de la distribucin de presin.
Brinda la opcin de especificaciones regionales de la presin y de los

fluidos presentes en el sitio.
Trabaja como herramienta de clasificacin de posibles modelos

geolgicos, reduciendo la incertidumbre en trminos de la valoracin del
comportamiento dinmico del yacimiento.
Maneja adecuadamente procesos de upscaling, reduciendo

cuidadosamente el nmero de celdas y permitiendo simular yacimientos
completos.
Define claramente relaciones entre inyectores y productores y cuales no

estn contribuyendo en el desarrollo del campo. As mismo, permite la
identificacin de reas de posible inters.
Trabaja con procesos de Ajuste Histrico Asistido.
Su programacin es mediante keywords y posee un formato libre para la

entrada de archivos de datos.
127
5.2.2 Estructura de una simulacin Frontsim. Como todos los simuladores

Eclipse, FRONTSIM posee la siguiente secuencia estructural para su
programacin:
RUNSPEC: Es una seccin requerida, en donde se especifican el

nombre del proyecto, las dimensiones del problema, cambios a realizar
y nmero de fases presentes entre muchas otras caractersticas.
GRID: Seccin requerida donde se especifica la geometra del grid

computacional y su esquema de distribucin; as mismo se definen
ciertas propiedades de la roca como porosidades y permeabilidades
absolutas entre otras.
EDIT: Seccin en la que se realizan modificaciones a los valores de

volmenes porosos calculados, profundidades de los nodos centrales y
transmisibilidades. Es de carcter opcional en la simulacin. Sus
keywords son algunas de las empleadas en la seccin GRID y
desempean la misma funcin.
PROPS: Seccin de carcter requerido en FRONTSIM. En ella se

especifican tablas de propiedades de la roca y de los fluidos como
funcin de la presin y de la saturacin (densidad, viscosidad,
permeabilidades relativas, etc.).
REGIONS: Su funcin es dividir el grid computacional en regiones para

calcular propiedades PVT (densidades y viscosidades), propiedades de
saturacin (Permeabilidades relativas), condiciones iniciales y
determinar por zonas el comportamiento del flujo y el OOIP. Si el
yacimiento presenta variaciones en sus propiedades por zonas es
aconsejable emplear esta seccin. En su definicin emplea los mismos
keywords de la seccin EDIT y algunos adicionales.
SOLUTION: Seccin de carcter requerido en donde se especifican las

condiciones iniciales del yacimiento y de los fluidos presentes.
SUMMARY: Esta seccin no es empleada por FRONTSIM.
SCHEDULE: Especifica las operaciones a ser simuladas, condiciones y

los tiempos a los cuales se exigen reportes de los resultados arrojados
por el simulador.
128
Para un mejor desarrollo de la simulacin es aconsejable seguir la secuencia

de secciones anteriormente descrita. La programacin de la estructura del
problema de simulacin puede desarrollarse en el block de notas o
directamente en ECLIPSE OFFICE y posteriormente correrse en el entorno del
simulador.
Figura 44. Bosquejo del entorno del simulador FRONTSIM.
Tomada del simulador Frontsim.
La visualizacin grfica de los resultados arrojados por FRONTSIM pueden

observarse mediante la herramienta GRIDSIM de ECLIPSE.
Para una informacin ms detallada sobre el funcionamiento y planteamiento

de estas secciones y keywords el lector puede dirigirse directamente al manual
del usuario ECLIPSE-FRONTSIM 2001A que se encuentra en el paquete de
Simulacin de ECLIPSE.
129
6. ULTIMOS AVANCES EN SIMULACIN STREAMLINE.
6.1 YACIMIENTOS NATURALMENTE FRACTURADOS

La simulacin streamline ha demostrado ser exitosa en el modelamiento de
yacimientos que presentan un nico sistema de porosidad. Sin embargo,
debido a sus caractersticas, la simulacin streamline puede ser fcilmente
adaptable a modelos de porosidad dual mediante el uso e implementacin de
funciones de transferencia entre la matriz y la fractura en las ecuaciones
bsicas que constituyen el modelo matemtico streamline. De esta forma, los
sistemas matriz y fractura son tratados en forma independiente, trazando en
cada uno de ellos grupos de streamlines. Entre los mltiples autores que han
estudiado y desarrollado est metodologa de aplicacin a sistemas fracturados
se encuentran Ahmed H Al-Huthali* et al (2004) quienes a su vez han centrado
sus investigaciones en la aplicacin de funciones de transferencia estudiadas
por autores como Kazemi**, Sonier*** y Litvat**** entre otros.
6.1.1 Formulacin matemtica modelo streamline en YNF. En yacimientos

naturalmente fracturados pueden presentarse los siguientes sistemas:
9 DPDP (Sistema de porosidad dual y doble permeabilidad): se
presenta cuando los sistemas fractura y matriz interactan entre
s y en cada uno de ellos se presenta almacenamiento y flujo de
fluidos.
9 DPSP (Sistema de porosidad dual y nica permeabilidad): Se
presenta cuando la fractura adems de almacenar fluidos es la
nica va principal de flujo. La matriz solo alberga fluidos.
9 SPSP (Sistema de nica porosidad y nica permeabilidad): Se
presenta cuando el almacenamiento y el flujo de fluidos solo
ocurre a travs de la fractura.
Con base en estas consideraciones Ahmed Al-Huthali y Datta-Gupta (2004)
formularon las siguientes ecuaciones para la determinacin de presiones y
saturaciones en YNF empleando simulacin streamline:
____________
*
AL-HUTHALI, ARAMCO, AHMED, SAUDI, DATTA-GUPTA, AKHIL. Streamline Simulation of Water
Injection in Naturally Fractured Reservoirs. Symposium on Improved Oil Recovery. Oklahoma. 17-21 abril.
2004. SPE N 89443.
**
KAZEMI, H. Analytical and Numerical Solution of Oil Recovery from Fractured Reservoirs with Empirical
Transfer Functions. SPEJ. May 1992. Pg 219.
***
SONIER, F. Numerical Simulation of Naturally Fractured Reservoirs. SPE 15627. New Orleans. Oct
1986.
****
LITVAK, B.L. Simulation and Characterization of Naturally Fractured Reservoirs. Dallas. April 29-May
1 1985.
130
o Presiones.
9 DPDP: Las streamlines deben ser trazadas tanto para el sistema
matriz como para el sistema fractura.
Fractura: k f
Matriz:
tf
r
Pf + gf Z ) + t = q sf (148)
r
k m ( tm Pm + gm Z ) t = q sm .... (149)
Donde la movilidad total t y la movilidad debida a la gravedad

g se expresan como:
t = o + w
g = og + wg
... (150)
Y la funcin de transferencia total t est dada por:
t = o + w .. (151)
9 DPSP: En esta ecuacin el trmino de transferencia no afecta la
trayectoria de las streamlines, ya que este se da en igual
magnitud en la matriz. Las streamlines solo se trazan para el
sistema fractura.
Sistema Fractura:
r
k f (tf Pf + gf Z ) = q sf
. (152)
o Saturaciones en funcin del Tiempo de vuelo.

9 DPDP
S wf
Fractura:
f wf
f
131
r
Gf
= 0 .... (153)
Matriz:
r
S wm f wm G m w
+
+
+
= 0 .. (154)
t
m
m
m
9 DPSP
Fractura:
Matriz:
9 SPSP
Fractura:
Ecuacin 152
w = m
S wm
t
.... (155)
r
S w f w G
+
=0
+
..... (156)
Como puede observarse, las ecuaciones para distribucin de saturacin llevan

implcita la determinacin de funciones de transferencia matriz-fractura. Estas
funciones de transferencia pueden ser convencionales o empricas.
Funciones de Transferencia Convencionales.
9 DPDP:
Fase Agua:
w = Fs k m wmf (Pwf Pwm )
..... (157)
Fase aceite:
o = Fs k m omf (Pof Pom )
........ (158)
Donde
Pwf = Pof Pcf

Pwm = Pom Pcm
.. (159)
132
Fs es el factor de forma que define la conectividad entre la

fractura y los bloques de la matriz que rodean dicha fractura.
Matemticamente se expresa como:
1
1
1
Fs = 4 2 + 2 + 2 ... (160)
I
x Iy Iz
Y wmf y omf representan las relaciones de movilidad entre los
sistemas matriz-fractura para cada una de las fases.
9 DPSP:
wmf omf
(Pcm Pcf
wmf + omf
w = Fs k m
.... (161)
Funciones de Transferencia Empricas. Este tipo de funciones de

transferencia son usadas cuando la imbibicin contra corriente es la fuerza que
domina el desplazamiento de aceite desde la matriz; no obstante solo pueden
aplicarse en sistemas DPSP. En este caso la funcin se define como:
9 DPSP:
w @ Sw =1 = Q e t .... (162)
Donde es una tasa constante determinada en laboratorio y
definida como el recproco del tiempo requerido para que la matriz
libere el 63% del aceite recuperable.
DeSwaan propuso de igual forma la siguiente funcin de

transferencia emprica:
t
w = Q e
(t )
S wf ( )
.... (163)
133
Otro aspecto importante a considerar en este tipo de yacimientos es la

influencia de la gravedad en el sistema matriz-fractura, como puede observarse
en la figura 45, hecho que puede ocasionar movimiento de fluidos.
En estos casos, los autores Ahmed Al-Huthali y Datta-Gupta (2004) han

planteado las siguientes ecuaciones para el clculo de funciones de
transferencia:
9 DPSP:
w = Fs k m
wmf omf
omf + wmf
(P
cm
Pcf + Pgh ) .... (164)
Donde la influencia de la gravedad Pgh , y la saturacin de agua

normalizada en el sistema fractura y en la matriz, S wnf y S wnm
respectivamente son determinadas de la siguiente forma:
Pgh = L z (S wnf S wnm )( w o )g
S wnf =
S wf S wcf
1 S orf S wcf
S wnm =
S wm S wcm
1 S orm S wcm
.. (165)
6.1.2 Soluciones Numricas a las Ecuaciones de Saturacin.

9 DPDP:
Para resolver las ecuaciones 153 y 154 para los sistemas matrizfractura, Ahmed Al-Huthali y Datta-Gupta (2004) plantearon dividir
dichas ecuaciones en dos secuencias: primero resolver un paso
convectivo y despus un paso correctivo. El paso convectivo
incluye la evolucin de la saturacin a lo largo de las streamlines
debido a las fuerzas viscosas. El paso correctivo incorpora el
trmino de transferencia entre los sistemas matriz-fractura. Estos
pasos matemticamente son definidos como:
134
Figura 45. Efecto de la gravedad en un bloque de la matriz rodeado por fracturas

parcialmente llenas de agua.
Agua en la Fractura
Aceite en la Fractura
Agua en la Matriz
Aceite en la Matriz
Paso Convectivo:
S wf
t1
f wf
f
=0
(166)
S wm f wm
+
= 0 ..... (167)
t1
m
Paso Correctivo:
S wf
t 2
S wf
t 2
f
w
=0
(168)
=0
. (169)
Para resolver las anteriores ecuaciones se emplea un esquema

explcito de diferencias finitas para discretizar los trminos
convectivos en los sistemas fractura y matriz.
135
Sistema Fractura:
n +1
wf ,i
= t
n
wf ,i
f wfn ,i f wfn ,i 1
f
.. (170)
Sistema Matriz:
n +1
wm,i
n
wm,i
= t
n
n
f wm
,i f wm,i 1
.. (171)
Donde i representa los nodos a lo largo de las streamlines.
De igual forma se emplea el mismo esquema numrico para

resolver el paso correctivo en los sistemas matriz y fractura.
n +1
wf ,i
n
wf ,i
F k
= t s m
wmf omf
+
wmf
i omf
(Pcm Pcf
.. (172)
Donde i representa el nmero de grid blocks para el paso

correctivo.
Las saturaciones de la fractura y de la matriz pertenecientes al

paso convectivo son trazadas sobre el grid y usadas como
condiciones iniciales para el paso correctivo. Para trazar la
saturacin de los sistemas fractura y matriz en el grid se han
propuesto las siguientes ecuaciones:
Saturacin Promedio en la Fractura:

nsl
S wf , grid =
S
i =1
wf ,i
f ,i
. (173)
nsl

i =1
f ,i
136
Saturacin Promedio en la Matriz:

nsl
S wm, grid =
S
i =1
wm ,i
m,i
. (174)
nsl
m,i
i =1
En esta ecuacin nsl es el nmero de streamlines que pasan a

travs de un grid block y es el tiempo de vuelo a travs de dicho
grid block.
9 DPSP:
En este tipo de sistemas las streamlines son generadas y trazadas
nicamente en el sistema fractura. Dependiendo de si en el
planteamiento
se
emplean
funciones
de
transferencia
convencionales o empricas las ecuaciones para determinar
saturacin en los sistemas matriz y fractura son las que se
presentan a continuacin:
Funciones de transferencia convencionales (CTF):
Saturacin en la Fractura:
n +1
wf ,i
n
wf ,i
fn fn
F k
wf ,i
wf ,i 1
= t
+ s m

f
wmf omf
+
wmf
i omf
(Pcm Pcf
) (175)
n
Saturacin en la Matriz:
n +1
wm ,i
n
wm ,i
F k
= t s m
wmf omf
+
wmf
i omf
137
(Pcm Pcf
.. (176)
Funciones de transferencia empricas (ETF):

Saturacin en la Fractura:
n +1
wf ,i
n
wf ,i
1
1 Q
f n f wfn ,i 1 Q

t wf ,i
e
SUM n 1e t (177)
+
=
+

f
f
f
t

)]
SUM n1 = SUM n2 + S wfn ,i S wfn1,i e t para
n 1
SUM n1 = 0
Saturacin en la Matriz:
n +1
n
n 1 t
S wm
S
t
SUM
e
,i
wm,i
.. (178)
f
De igual forma como ocurre en el sistema DPDP la saturacin de

agua en el sistema fractura a lo largo de las streamlines puede ser
trazada sobre el grid empleando la ecuacin 173. Por su parte la
saturacin en la matriz a lo largo de las streamlines puede ser
trazada sobre los grid blocks usando el siguiente promedio
aritmtico:
S wm, grid =
1 nsl
S wm,i
nsl i =1
. (179)
Aunque muchos autores han desarrollado mtodos alternos para

tratar este tipo de yacimientos altamente complejos, este es el
planteamiento bsico que se realiza en el estudio de yacimientos
naturalmente fracturados empleando la simulacin streamline. Un
estudio bastante detallado en cuanto al uso de las Funciones de
Transferencia en Yacimientos Naturalmente fracturados puede
observarse en los artculos tcnicos escritos por Di-Donato* et al.
____________
DONATO, GINEBRA, BLUNT, MARTIN. A Streamline-Based Dual Porosity Simulator for Fracture Flow
Simulation. ITF Project. Jun. 9 2003.
DONATO, GINEBRA, HUANG, WENFEN, BLUNT, MARTIN. Streamline-Based Dual Porosity Simulation
of Fractured Reservoirs. Annual Technical Conference and Exhibition. Colorado, 5 8 Oct 2003. SPE N
84036.
138
6.2 NUEVOS MTODOS DESARROLLADOS PARA MANEJAR FLUJO

COMPRESIBLE EN SIMULACIN STREAMLINE
Con miras a ampliar la aplicabilidad de la simulacin streamline, Ingebrigtsen*
et al (1999) desarrollaron dos novedosos mtodos para adaptar sistemas
compresibles en la formulacin streamline. Los mtodos planteados se
describen a continuacin.
Considrese flujo de tres fases compresibles en un medio poroso. Por

simplicidad los efectos de gravedad y las fuerzas capilares son despreciables.
Asumiendo una fase agua, una fase aceite, y fase gas de tal forma que una
parte del gas est disuelta en la fase aceite, mientras que la otra se presenta
como una fase libre. Las ecuaciones de conservacin para cada componente
son:
( s wb w) +
b f v
t
( s ob o) +
b f v
t
= qw
( (s g b g + R s s o b o ))
t
= qo
...... (180)
.. (181)
+ bg
f v
g
R s f
s
vt =
o
(182)
g
Sumando las ecuaciones de conservacin para cada componente se obtiene la

ecuacin de presin:
P
+
v
=
Q
P .. (183)
ct t t
Donde vt es la velocidad de Darcy y esta dada por:
vt = t P .. (184)
____________
INGEBRIGTSEN, L., BRATVEDT, F., BERGE, J. A Streamline Based Approach to Solution of ThreePhase Flow. Reservoir Simulation Symposium Held in Houston. Texas 14-17 Feb. 1999. SPE N 51904.
139
6.2.1 Mtodo Secuencial. Con este mtodo se desea obtener una expresin
explcita para cada una de las saturaciones de las fases. La ecuacin para el
componente del gas contiene un tiempo derivado de la saturacin de aceite.
Para eliminar este trmino se asume que la cantidad de gas en el aceite debe
ser constante en el tiempo. Si se asume que la porosidad es constante en el
tiempo se obtiene una expresin para la saturacin del gas de la siguiente
forma:
bg
sg
+ v t b g
t
(b f ) v + b
g
vt R s = q
.... (185)
Despreciando el trmino v t a lo largo de las streamlines en la ecuacin

anterior y aplicando un operador igualdad finalmente se obtienen las
ecuaciones de saturacin a lo largo de las streamlines:
bo
(s o)
t
bw
bg
(s w)
t
(s g )
t
(bo
(bw
= q .... (186)
o
(b g
= q .... (187)
w
+ bo
(R s )
o
... (188)
g
Las ecuaciones obtenidas son similares a las ecuaciones para flujo

incompresible y se pueden tener en cuenta para efectos compresibles.
El procedimiento de solucin del mtodo es dividido en cinco pasos:

1. Se resuelve la ecuacin de presin (ecuacin 183), conociendo
(Po,sgo,swo,soo) como condiciones iniciales. La presin (P1(so)) es
obtenida usando un esquema implcito de diferencias finitas.
2. la velocidad total es calculada mediante la ecuacin 184 y las
streamlines son trazadas.
140
3. Se ajustan las saturaciones para obtener conservacin de masa

despus de que la presin ha sido actualizada. Suponga las
saturaciones iniciales ya establecidas tal que la fraccin msica para
cada componente a una presin inicial esta definida por:
n = s b (P , P ) .... (189)
o
o,0
n = s b (P , P ) + R (P , P )n
g ,0
. (190)
Las nuevas saturaciones que corresponden a los valores de

Rs(P1(so),Pb) y bj(P1(so),Pb) pueden ser obtenidas, calculando primero
los volmenes de las fases a partir de las siguientes ecuaciones:
n R (P (s ), P )n
~
(
)
=
s P
b (P (s ), P )
1
g ,0
~
sw,0 = sw,0
.... (191)
n
~
s (P ) = ( ( ), )
b Ps P
o,0
.. (192)
..... (193)
Seguidamente, la presin del punto de ebullicin Pb

debe ser
nuevamente calculada. Despus de determinar los volmenes de las
fases (ecuaciones 191, 192 y 193), stos se suman y las nuevas
saturaciones son calculadas mediante la siguiente ecuacin:
~
s
s =~
s
j
.... (194)
4. Las saturaciones son propagadas a lo largo de las streamlines

resolviendo las ecuaciones de saturacin y obteniendo las saturaciones
sj . Las saturaciones son trazadas en la malla fundamental.
141
5. Como las saturaciones son propagadas a lo largo de las streamlines, la

presin debe haber cambiado. Debido a que se estn resolviendo las
ecuaciones en una forma diferente a la tradicional y los mtodos de
solucin para la saturacin producen resultados que no satisfacen la
ecuacin de volumen constante que indica
s
j
= 1 .. (195)
Se debe sustentar esta discrepancia en el volumen en la ecuacin de

presin dada inicialmente (183)
6.2.2 Mtodo implcito. El procedimiento de solucin de este mtodo para

trabajar sistemas de flujo compresible es el siguiente:
1. Resolver la ecuacin de presin sobre la malla fundamental, calculando

la velocidad total y trazando las streamlines.
2. Ajustar las saturaciones a la presin actualizada como se hace en el
paso 3 del mtodo descrito anteriormente.
3. Resolver para la presin y la saturacin usando un esquema implcito
para cada streamline.
4. Trazar las saturaciones en la malla fundamental.
Recordando las ecuaciones (180), (181) y (182) para saturacin, se pueden

resolver estas ecuaciones sustituyendo para la ecuacin de volumen constante
(ecuacin 195), donde sg, sw y P son las incgnitas primarias. El nuevo
sistema de ecuaciones para este mtodo, despreciando y acomodando
trminos, es:
((
1 sg sw
t
)b ) + (b f ) q
o
= 0 . (196)
(s wb )+ (b
(s b + R (1 s s )b ) + (b f + R b f ) q
f
g
) q
=0
(197)
142
= 0 .... (198)
En este mtodo, se asume que las saturaciones llenan el volumen poroso

exactamente. Para hacer este mtodo completo se tendra que resolver
implcitamente Rs en el caso de flujo sub-saturado. Para est implementacin,
Rs es recalculada al final de cada paso de tiempo local.
Este mtodo es algo inmaduro, pero los resultados obtenidos son bastante
prometedores.
6.3 SISTEMAS DE ACEITE NEGRO CON EFECTOS CAPILARES

Aunque la mayora de los casos aplicados a la simulacin streamline
desprecian por simplicidad del mtodo los efectos capilares, una nueva
formulacin que permite incluir tanto fuerzas capilares como gravedad ha sido
desarrollada por Roman A. Berenblyum* et al (2003), pertenecientes a la
Universidad Tcnica de Denmark y la Universidad de Stanford. Como ya se ha
presentado en el desarrollo del planteamiento bsico del modelo matemtico
streamline, anteriormente no se tenan en cuenta los efectos de flujo en
direccin transversal a las streamlines, hecho generado por diferencias de
presin capilar entre las fases presentes. Estas circunstancias cambiaron con
un nuevo planteamiento de ecuaciones tanto de presin como de saturacin en
el mtodo, que involucra dos conceptos reales e influyentes en el flujo de
fluidos en el yacimiento: presin capilar y gravedad.
Esta nueva aproximacin suministra una herramienta rpida y confiable para

simular procesos de inyeccin de agua en yacimientos heterogneos de baja
permeabilidad.
La introduccin de las fuerzas capilares requiere la modificacin de las

ecuaciones de presin y de saturacin.
6.3.1 Modificaciones a la Ecuacin de Presin. La ecuacin de continuidad

para flujo incompresible considerada en esta nueva aproximacin es escrita
como:
u t = (u w + u o ) = 0
...... (199)
____________
BERENBLYUM, R.A., SHAPIRO, A., JESSEN, K., STENBY, E., ORR, F. Black Oil Streamline Simulator
with Capillary Effects. Technical University of Denmark and Stanford University. SPE 84037. Oct. 2003.
143
La velocidad de la fase j es obtenida de:
u j = k j Pj
.... (200)
La presin capilar es definida como la diferencia de la presin entre la fase no

mojante (aceite) y la fase mojante (agua), de forma que:
Pc = Po Pw
... (201)
Redefiniendo la ecuacin de continuidad (ecuacin 199) para cada una de las

fases, teniendo en cuenta la ecuacin de velocidad (200) y la definicin de
presin capilar (201), se llega a:
u t = (kt Pw + ko Pc ) = 0
. (202)
Si en la ecuacin (202) se tiene en cuenta el efecto de la gravedad, finalmente

se llega a la nueva ecuacin para distribucin de presiones streamline con
efectos capilares y gravitacionales:
ut = (kt Pw + ko Pc + k g D) = 0
... (203)
6.3.2 Modificaciones a la Ecuacin de Saturacin. El planteamiento tiene en

cuenta la velocidad total del flujo y parte de la definicin de gradiente de
presin para el agua, introduciendo en ella el concepto de presin capilar:
Pw =
1
u t + o Pc
k t
t
... (204)
De igual forma, la ecuacin para el balance de masa para la fase agua es

formulada como:
S w
+ u w = 0
t
.... (205)
144
Introduciendo los conceptos de velocidad del agua dado por la ecuacin (200) y
gradiente de presin para el agua (ecuacin 204), en la ecuacin de balance de
masa se llega a la siguiente expresin:
S w
+ w ut + (1 f )k w Po = 0 .... (206)
t
t
Donde:
o t w
=
= (1 f )
t
t
. (207)
Reagrupando y haciendo uso de la ecuacin de continuidad, se llega a la

expresin para la distribucin de saturacin con efectos capilares:
S w
+ u t f + [(1 f )k w Pc ] = 0 ... (208)
t
La ecuacin de saturacin con fuerzas gravitacionales es obtenida de la misma

manera:
S w
+ u t f + Tg D = 0
t
. (209)
Donde:
Tg = kg
wo
(o w )
t
.... (210)
La ecuacin de saturacin es finalmente resuelta usando una tcnica de

operador divisin:
S
t1
+ ut f = qi
. (211)
145
S
t 2
+ G = 0 . (212)
En estas ecuaciones qi hace referencia a fuentes o sumideros y G a la accin

de la gravedad y las fuerzas capilares, de forma que:
G = (1 f )kwPc + Tg D
..... (213)
La solucin de la ecuacin 211 es la forma tradicional del mtodo streamline

ya descrita en el captulo 4. Por su parte, la ecuacin 212 es resuelta
explcitamente sobre el grid de simulacin de diferencias finitas, empleando
una aproximacin de primer orden donde los gradientes de gravedad y
capilaridad son discretizados usando un esquema Stencil de 7 puntos.
Una explicacin ms detallada del mtodo es presentada por Roman, Shapiro,

Jessen, Erling. y Orr en su artculo tcnico.
6.4 SISTEMAS COMPOSICIONALES EN SIMULACIN STREAMLINE

Cuando los mtodos tradicionales de inyeccin de agua no alcanzan un
porcentaje de recobro ptimo, la implementacin de otros procesos de
inyeccin, requieren de un correcto modelamiento de los efectos de
transferencia de masa.
El simulador streamline en su formulacin composicional tiene en cuenta estos

efectos, y su operatividad y aplicabilidad en hechos reales ha sido demostrada
por Thiele* et al (1997) y Wei Yan** et al (2004).
6.4.1 Ecuaciones bsicas.

streamline es el siguiente:
El planteamiento composicional en simulacin
____________
THIELE, MARCO R., BATYCKY, ROD P, BLUNT, MARTIN J. A Streamline-Based 3D Field-Scale

Compositional Reservoir Simulator. Annual Technical Conference and Exhibition. Texas, Oct 5-8, 1997.
SPE N 38889.
**
WEI, Y., MICHELSEN, M., STENBY, E., BERENBLYUM, R., SHAPIRO, A. Three-phase Compositional
Streamline Simulation and Its Application to WAG. University of Denmark. SPE 89440. April 2004.
146
Las ecuaciones de conservacin de masa multidimiensional, para una especie

qumica, i , fluyendo bajo fuerzas convectivas puede expresarse como:
Ci
+ u t Fi = 0
t
i = 1.............nc .. (214)
Donde Ci, son las moles locales del componente i , definidas como:
np
Ci = xij j S j
j =1
j = 1............n p
(215)
Para esta y las dems ecuaciones del planteamiento composicional, np es el

nmero de fases presentes, nc es el nmero de componentes, xij es la fraccin
molar del componente i en la fase j, y j y Sj son la densidad y la saturacin de
la fase j respectivamente. Por otra parte, el flujo molar del componente i, se
define matemticamente como:
np
Fi = u D xij j f j
j =1
j = 1............n p
.. (216)
Donde uD es la velocidad volumtrica total adimensional y fj es el flujo fraccional

de la fase j.
Usando la definicin de tiempo de vuelo (ecuacin 62), es posible definir la

ecuacin 214, de la siguiente forma:
Ci Fi
+
=0
t
i = 1.............nc
(217)
La velocidad total en el campo, ut, necesaria para determinar el tiempo de

vuelo, puede ser determinada resolviendo la ecuacin de balance de masa
para la presin y entonces, usar la Ley de Darcy para determinar la velocidad
local del campo.
147
La aproximacin formulada para este evento es la siguiente:
K (t P ) = 0
.... (218)
Donde K es el tensor de permeabilidad local, y t es la movilidad total.
6.4.2 Soluciones unidimensionales a las ecuaciones de flujo composicional.

La ecuacin 217 es resuelta para las moles totales del componente i, usando
un esquema explcito con limitaciones de primer orden en el tiempo y de
segundo orden en el espacio. El esquema puede ser mejorado usando una
simple modificacin, de forma que:
Fi = u D Fi
.. (219)
Con esta relacin, la ecuacin 217 puede ser escrita en forma de diferencias
finitas como:
n +1
i ,k
n
n
t

u
F
u
F
=C +
D i
D i

k 1
k + 1
2
2
n
i ,k
.... (220)
Donde los flujos inter.-celdas, Fi , son aproximados por la siguiente ecuacin:

n
i ,k +
1
2
i ,k
n
+
i ,k +
1
2
n
n
F
1 t .. (221)
F
i
k
i
k
+
,
1
,
es el lmite de Van Leer, definido como:
r+ r
1+ r
.. (222)
Definiendo a r como la funcin del radio de diferencias de flujo adyacentes,

expresada mediante la siguiente ecuacin:
148
i ,k +
1
2
F ni + Fi n1
= n
Fi +1 Fi n
.... (223)
Donde k, es el nodo opuesto de la streamline discretizada y n, es el nivel de

tiempo local a lo largo de cada streamline. Por otra parte, las fracciones
molares totales para cada componente, son obtenidas por:
zl =
Cl
nc
i =1
Ci
.. (224)
Finalmente, la velocidad total adimensional uD, se expresa mediante la

siguiente relacin:
n
n +1
i =1 t Ci ,k Ci ,k + u D Fi 1
k
2
=
nc
un
D,k +
1
2
F
nc
i =1
i ,k +
.... (225)
1
2
Para la metodologa planteada, las composiciones de las fases, densidades y

saturaciones son determinadas mediante clculos de equilibrio usando la
ecuacin de estado de Peng-Robinson y clculos flash. Las viscosidades de
las fases son calculadas usando la correlacin de Lohrenz-Bray-Clark.
149
7. PLANTEAMIENTO DE DOS CASOS BASE
Para demostrar tanto las ventajas como desventajas de la simulacin

streamline y corroborar lo dicho en teora, se planteo la simulacin de dos
modelos, uno homogneo multipozo, y otro heterogneo compuesto por un
pozo inyector y uno productor.
7.1 MODELO MULTIPOZO. PRIMER CASO BASE

Este modelo, es un yacimiento areal y homogneo, compuesto por 10 pozos
productores (verticales) y 5 pozos inyectores de agua, con 400 celdas activas
(20x20x1), porosidad de 0.25, permeabilidad en todas las direcciones igual a
100 md y saturacin inicial de agua de 0.25. El objetivo de este modelo, aparte
de validar el procedimiento de simulacin streamline, es poder presentar la
utilidad de esta tecnologa en el estudio de patrones de inyeccin-produccin
para sistemas de mltiples pozos. El esquema del yacimiento modelado es el
siguiente:
Figura 46. Esquema del grid para el primer caso base a simular.
7.1.1 Condiciones de desarrollo del modelo.

Inicialmente, se planteo el
desarrollo del modelo con todos los pozos tanto productores como inyectores
abiertos a produccin e inyeccin respectivamente, durante un periodo de
aproximadamente cinco aos.
150
Despus de este tiempo se cambiaron las condiciones de operacin del campo

cerrando 7 de los pozos productores (P1,P2,P3,P7,P8,P9 y P10). Se contino
con el proceso de inyeccin de agua con el mismo nmero de inyectores por un
periodo de 5 aos ms. El modelo fue simulado en FRONTSIM y en ECLIPSE
100.

Anlisis de la veracidad del mtodo
streamline. Despus de realizada la simulacin en FRONTSIM se observaron
los siguientes resultados.
o Determinacin de la distribucin de presin.
Definiendo diferentes tamaos

en los pasos de tiempo, se observo el mismo comportamiento en el clculo y
determinacin de la distribucin de presin en el campo.
Grficamente se observaron los mismos esquemas de presin para los dos

eventos modelados, y los datos arrojados por el simulador para cada tamao
de tiempo empleado difieren aproximadamente en 0.2 PSI. Estos
comportamientos constantes a travs del tiempo se observan en las figuras 47
y 48.
Figura 47. Esquemas de distribucin de presin para el primer caso base. Izquierda
evento con todos los pozos productores abiertos. Derecha evento con 3 pozos
productores.
151
Figura 48. Grficas de presin variando el nmero de pasos de tiempo en Frontsim.

Primer caso base.
Tomada de Eclipse Office.
o Determinacin de las streamlines.

Teniendo en cuenta que durante los
primeros 5 aos todos los pozos tanto productores como inyectores se
encontraban en operacin a ciertas condiciones, el simulador pudo determinar
el siguiente patrn de lneas de flujo para esas condiciones lmites constantes
(ver figura 49). De igual forma se puede apreciar el avance del proceso de
inyeccin de agua a travs de todo el campo e identificar las reas de
influencia de cada inyector, los pozos productores que soporta, las reas que
no pueden ser barridas y los volmenes de yacimiento asociados con cada
pozo productor.
Figura 49. Streamlines determinadas en el caso base uno a las condiciones lmites
iniciales. Izquierda: patrn de lneas establecidas. Derecha: progreso de la inyeccin

152
Cambiando las tasas tanto de produccin como de inyeccin de los pozos

basados en la perspectiva visual que brindan las streamlines, se puede lograr
un mejor balance de los patrones tanto de inyeccin y produccin, mejorando el
barrido en zonas inalcanzables, y mejorando el porcentaje de recobro de
aceite. En este caso, el parmetro fundamental que permiti una mejor
optimizacin del esquema inicial de produccin fue la disminucin en las tasas
de inyeccin de algunos pozos, especialmente del inyector I3, de forma que
disminuyendo esta tasa, se lograr una mejor cobertura areal. El progreso del
balance de patrones establecido y el aumento en el factor de recobro se
observa en la figuras 50, 51, 52 y 53.
Figura 50. Patrn de lneas inicial. Primer caso base. Produccin = 9.32543 e4 STB
Figura 51. Primer ajuste a las tasas de inyeccin. Produccin = 1.01110e5 STB
153
Figura 52. Ultimo ajuste conseguido a las tasas de inyeccin. Primer caso base.
Produccin = 1.03417e5 STB.
Figura 53. Factor de recobro ajustando las tasas de inyeccin y produccin. Primer
caso base.
Otra forma de lograr un balance de patrones ptimo para el modelo multipozo

establecido y aumentar el porcentaje de recobro, es abrir un pozo inyector en el
rea de no influencia. De esta forma se definen nuevas zonas de drenaje antes
no consideradas. Este comportamiento puede apreciarse en las figuras 54 y 55.
154
Figura 54. Balance estableciendo un nuevo pozo inyector. Produccin = 1.10736 e5

STB
Figura 55. Aumento en el factor de recobro con la creacin de un nuevo pozo inyector.
Verificando el comportamiento del simulador, se corrobor visualmente la

aplicacin del mtodo de solucin numrica streamline (ver figura 56), ya que
despus de transcurridos los primeros cinco aos, se observ un cambio en la
distribucin de las lneas de flujo en el campo, esto como consecuencia del
cierre de algunos pozos productores que ocasion un nuevo planteamiento en
las condiciones lmites, y por consiguiente un nuevo esquema de barrido del
modelo.
155
No obstante, si se hubieran conservado las condiciones de produccin inicial,

se habra obtenido un mayor porcentaje de recobro de aceite, como puede
observarse en la figura 57.
Figura 56. Distribucin de las streamlines despus de cambiar las condiciones lmites.
Izquierda: nueva distribucin de lneas de flujo. Derecha: progreso de la nueva
inyeccin. Primer caso base.
Figura 57. Factor de recobro. Lnea roja: continuando con el patrn de pozos inicial.
Lnea verde: cambio en las condiciones de los pozos.
Todos los pozos productores abiertos
Solo 3 productores abiertos
156
Otra ventaja que presentan los simuladores streamline es el planteamiento de

los resultados a la escala de pozo, caso no observado en el simulador
convencional, lo cual permite un anlisis ms detallado del comportamiento
tanto de pozos productores como de inyectores. Para el patrn de pozos
productores establecidos, se obtuvo las siguientes grficas:
Figura 58. Tasa de produccin de aceite por pozo para P1, P2, P3, P4 y P5.
157
Figura 60. Tasa de produccin de aceite por pozo para P6, P7, P8, P9 y P10.
Mediante estas figuras es posible establecer que los mejores productores son
los pozos P5 y P6, ya que no experimentan un marcado descenso en su tasa
de produccin de aceite, y al mismo tiempo no presenta un alto corte de agua.
De igual forma puede decirse inversamente que los pozos P1, P2, P7 y P10 no
son muy productivos, declinan aceleradamente su tasa de produccin, y
experimentan un alto corte de agua debido a sus proximidades a los pozos
inyectores.
158
o Influencia de la gravedad.
En el caso base planteado, el efecto de la
gravedad no se evidencia fuertemente, debido a las condiciones del modelo
(areal) y a las densidades de los fluidos presentes que eran muy similares. Los
resultados obtenidos con y sin efecto de la gravedad para distribucin de
saturacin y comportamiento de produccin mostraron resultados similares,
con una disminucin poco importante en la distribucin de presin de
aproximadamente 9 PSI para el caso de no gravedad. Esta diferencia es
consecuencia de los trminos gravitacionales tenidos en cuenta en la ecuacin
de presin, para el caso de considerarse este efecto. Recordar ecuaciones 43 y
45. Este comportamiento se aprecia en las figuras 62, 63, 64 y 65.
Se debe resaltar que la gravedad siempre debe ser considerada en la

simulacin a fin de obtener resultados ms exactos y satisfactorios,
especialmente cuando se trabajan yacimientos con marcados efectos de
segregacin gravitacional o marcada diferencia en la densidad de las fases
presentes.
Figura 62. Presin con y sin efectos de gravedad. Primer caso base.
159
Figura 63. OIP con y sin efectos de gravedad. Primer caso base.
Figura 64. Produccin total de agua con y sin efectos de gravedad. Primer caso base.
Figura 65. Produccin total de aceite con y sin efectos de gravedad. Primer caso base.
Figuras tomadas de Eclipse Office.
160
o Determinacin de la distribucin de saturacin en el modelo. Para el caso

base establecido, se comparo la eficiencia del simulador en cuanto a la
determinacin de las saturaciones. Se tomaron diferentes tamaos en el paso
de tiempo y se analizaron los eventos de produccin y cierre. Para todos los
tamaos de pasos de tiempo la distribucin de saturaciones obtenida fue
similar, aunque en ellas puede evidenciarse un mnimo efecto de difusin
numrica (grafica arriba) cuando se aumenta el nmero de pasos de tiempo
empleados. Este efecto de difusin ocasiona a su vez errores mnimos en el
balance de materiales, aunque al algoritmo de control empleado por el
simulador que busca disminuir este porcentaje. Estos efectos se observan en
las figuras 66 y 67.
Figura 66. Distribucin de saturaciones halladas por el simulador Frontsim para el
primer caso base.
Arriba: 8 pasos de tiempo. Abajo: 2 pasos de tiempo.
Tomadas de Eclipse Office.
161
Figura 67. Porcentaje de error en el balance de materiales. Primer caso base. Lnea
roja: 8 pasos de tiempo. Lnea verde: 2 pasos de tiempo.
8 pasos de tiempo
2 pasos de tiempo
Segn los resultados obtenidos, en estas figuras claramente se observa que

tomando pasos de tiempo relativamente largos, el simulador realiza los clculos
en una forma sencilla y los datos obtenidos van a ser ampliamente confiables y
libres de efectos de difusin numrica, ocasionados por la naturaleza del
mtodo. Esta es justamente la caracterstica o ventaja que brinda la simulacin
streamline: robustez, eficiencia y velocidad a partir de pasos de tiempo
extensos. Esta ventaja solo se consigue con este tipo de simuladores, ya que
como se explico anteriormente, los clculos de saturaciones no son funcin
directa de la geometra ni de las dimensiones de las celdas.
o Anlisis de los errores en el balance de materiales introducidos por cambios

en las condiciones de los pozos. Para analizar este efecto, se consideraron
dos casos: el primero el definido inicialmente en el modelo con todos los pozos
productores en operacin y luego un periodo de cierre de siete de estos pozos;
el segundo, el planteamiento de las mismas condiciones anteriores y otros
periodos de apertura y cierre de algunos de estos pozos por el mismo periodo
de tiempo establecido. El resultado grfico de este anlisis se observa en la
figura 68. Este aumento en el porcentaje de error (aunque se mueve en un
rango aceptable) a medida que se varan las condiciones de operacin del
modelo ocurre debido al traslado de la informacin de las condiciones de
saturacin antiguas a las nuevas lneas de flujo establecidas por el simulador,
que ocasionan errores de este tipo. Cuando las streamlines se actualizan
frecuentemente por cambios en las condiciones de operacin, los errores en el
trazado limitan la exactitud total del mtodo streamline.
162
Figura 68. Anlisis del porcentaje de error en el balance de materiales con cambios en
las condiciones lmites. Primer caso base.
Condiciones lmites muy variables
Condiciones lmites poco variables
o Anlisis del tiempo de simulacin empleado. De igual forma se analiz los

datos en cuanto al tiempo de CPU empleado en cada caso. La tabla 1
presenta los resultados obtenidos. Tal y como se esperaba, a medida que se
incrementa el tamao de los pasos de tiempo, se disminuye notoriamente el
tiempo de simulacin y los resultados determinados por el simulador son tan
confiables como aquellos simulados mediante pasos de tiempo cortos. Esta
evidencia en los resultados se aprecia en la tabla 2 y en las grficas de
produccin. No obstante cabe resaltar que para desplazamientos que exhiben
un comportamiento altamente no-lineal, previamente se debe establecer un
tamao de pasos de tiempo ptimo, a fin de reducir los porcentajes de error en
el mtodo.
Tabla 1. Tiempo de cmputo del simulador Frontsim para cada paso de tiempo
empleado.
Pasos de tiempo empleados
Tiempo de CPU (seg.)
1.28
46
2.55
92
4.66
Tomada de datos de Frontsim.
163
Tabla 2. Produccin obtenida por FRONTSIM para diversos pasos de tiempo.

Pasos de tiempo empleados
Produccin (STB)
9.32931 e4
46
9.32672 e4
92
9.32768 e4
Tomada de datos de Frontsim.
7.1.3 Simulador Frontsim vs. simulador Eclipse 100.

Para determinar las
ventajas del simulador streamline, el modelo del caso base inicialmente
planteado se construyo en el simulador Eclipse 100 (Black-Oil). Para ambos
casos, se compararon los resultados en cuanto a comportamientos de
produccin obtenidos, tiempo de simulacin, memoria empleada, distribucin
de presiones y saturaciones.
o Distribucin de presiones y saturaciones.

Para el comportamiento de la
presin arrojado por los dos simuladores y presentado grficamente en la figura
69, se observan que los resultados obtenidos son bastante consistentes,
presentando una diferencia de presin a lo largo del yacimiento de cerca de 8
PSI entre los dos simuladores. Grficamente esta distribucin de presiones
puede observarse en la figura 70 para el mismo tiempo de simulacin.
Figura 69. Grfica de presin para los simuladores ECLIPSE 100 y FRONTSIM.
Primer caso base.
164
Figura 70. Comparacin de la distribucin de presin obtenida por los dos simuladores
para el mismo tiempo de simulacin. Izquierda simulador Eclipse 100. Derecha
simulador Frontsim.
Tomada del simulador Eclipse 100 y Frontsim.
Para la distribucin de saturaciones obtenida por los dos simuladores, se

observo una diferencia de cerca de 0.04 para los valores de saturacin en las
celdas obtenidos por Eclipse 100 y Frontsim, adems de observar en la
distribucin hallada por Eclipse 100 un marcado efecto de difusin numrica
(izquierda). Esta diferencia es el resultado de las formulaciones matemticas
aplicadas en cada simulador, no obstante los resultados son bastante
satisfactorios, indicando un buen planteamiento del modelo por parte de los dos
simuladores. Este comportamiento es apreciado en la figura 71.
o Anlisis de los datos de produccin obtenidos. Para los dos simuladores, el

comportamiento en cuanto a los datos de produccin fue bastante similar y
consistente con los eventos de simulacin considerados (tiempos de
produccin y cierres).
165
Figura 71. Distribucin de la saturacin obtenida por los dos simuladores para el
mismo tiempo de simulacin. Izquierda Eclipse 100. Derecha Frontsim.
Tomada del simulador Eclipse 100 y Frontsim .
La ligera discrepancia puede ser consecuencia del esquema empleado en cada

simulador (Eclipse 100: esquema de punto centrado. Frontsim: esquema de
punto distribuido.) y de la solucin a las ecuaciones de flujo empleada en cada
mtodo. Los resultados obtenidos se observan en las figuras 72, 73, 74 y 75.
Figura 72. Tasa de produccin de aceite obtenida para los dos simuladores.
166
Figura 73. Tasa de inyeccin de agua para los dos simuladores considerados.
Figura 74. Produccin total de agua obtenida por Eclipse 100 y Frontsim.
Figura 75. Produccin total de aceite Obtenida por Eclipse 100 y Frontsim.
167
o Tiempo de simulacin y memoria empleada. Finalmente, para mostrar las

ventajas de los simuladores streamline, se compararon los tiempos de
simulacin y la memoria requerida por cada software en cada corrida. Se vario
el numero de pasos de tiempo para verificar la eficiencia streamline en cuanto a
velocidad sobre el simulador tradicional.
Los resultados se observan en la tabla 3.
Tabla 3. Tiempo de simulacin empleado en cada simulador.

Nmero de
Pasos de
Tiempo
TCPU
Frontsim
seg.
TCPU
Eclipse
seg.
% Velocidad de
Frontsim sobre
Eclipse
Memoria
Usada por
Frontsim
MB
Memoria
Usada por
Eclipse MB
92
4.78
13.125
0.675
0.477
46
2.83
7.563
0.675
0.477
0.98
1.609
2.75 veces ms
rpido
2.68 veces ms
rpido
1.6 veces ms
rpido
0.675
0.477
Tomada de datos de los simuladores Eclipse 100 y Frontism..
En estos datos, claramente se puede apreciar la ventaja en cuanto a velocidad

del simulador streamline sobre el tradicional, ya que se evidencia que an
tomando el mismo nmero de pasos de tiempo, el simulador streamline, por
sus algoritmos y su formulacin 1D es ms rpido y a medida que se aumentan
los tamaos en los pasos de tiempo, se reduce el efecto de difusin numrica,
efecto contrario al comportamiento de Eclipse. Si esto se traduce a ejemplos
reales de campo, donde los modelos considerados son bastante extensos, el
simulador streamline podramos decir que sobrepasara en velocidad a un
simulador convencional, y esto es aprovechable especialmente cuando se
realizan procesos de ajuste histrico, validacin de procesos de upscaling,
geoestadstica entre otros, que reduciran en trminos econmicos los gastos
en tiempo de simulacin.
En cuanto a la memoria requerida por Frontsim, es obvio que por las

caractersticas del mtodo streamline se requiera un poco ms de memoria
disponible en el equipo, no obstante el simulador posee comandos
especializados para no almacenar en forma permanente estos datos, hecho
que no podra considerarse como una desventaja del mtodo.
168
7.2 MODELO HETEROGNEO. SEGUNDO CASO BASE

Para apreciar la especialidad del simulador streamline en cuanto al manejo de
yacimientos altamente heterogneos, se planteo un modelo de yacimiento de
este tipo, compuesto por un pozo inyector de agua y un pozo productor. El
modelo est compuesto por 500.000 celdas activas (100x100x50) con
porosidad de 0.2 en todas las celdas y saturacin de agua inicial igual a 0.2. La
distribucin de permeabilidad en el modelo establecido puede apreciarse en la
figura 76.
Figura 76. Distribucin de permeabilidad para el modelo heterogneo considerado.
Tomada de GridSim.
El objetivo de esta simulacin, aparte de validar el mtodo streamline como se

hizo con el caso anterior, es demostrar la ventaja de este simulador en cuanto
al tratamiento de este tipo de yacimientos altamente heterogneos.
Una primera ventaja obtenida al inicio de la simulacin fue el mismo

planteamiento del modelo, ya que por la extensin y complejidad de este
(500.000 celdas activas), no fue posible correrlo en Eclipse 100, ya que el
simulador arrojaba problemas de convergencia. Para poder realizar esta
simulacin en un simulador convencional se requerira de establecer un
sistema de simulacin en paralelo.
7.2.1 Resultados de la simulacin. Anlisis de la veracidad del mtodo

streamline. Despus de realizada la simulacin en FRONTSIM se observaron
los siguientes resultados.
169
o Determinacin de la distribucin de Presin. Para los diferentes pasos de

tiempo establecidos, se obtuvo la siguiente distribucin de presin. Ver figura
77
Figura 77. Distribucin de presin en el modelo heterogneo.
Figura 78. Presin para cada uno de los pasos de tiempoeEmpleados. Modelo
heterogneo.
170
Segn la figura 78, puede observarse una ligera discrepancia entre los datos
obtenidos por el simulador streamline para cada tamao de paso de tiempo, lo
cual verifica la veracidad del mtodo en cuanto a la determinacin de la presin
en yacimientos altamente heterogneos.
o Determinacin de las streamlines. Para cada set de pasos de tiempo

empleado, se observ la misma distribucin de lneas streamline graficada por
el simulador, identificando claramente las reas de flujo en el modelo. En estas
figuras adems de apreciarse mejor el mtodo del trazado de Pollock, se
observa con detalle el avance del proceso de inyeccin y la saturacin a travs
de las lneas de flujo definidas. Esto se observa en la figura 79.
Figura 79. Izquierda: streamlines definidas para el modelo simulado. Derecha: avance
de la inyeccin en las streamlines. Modelo heterogneo.
Otra propiedad apreciable en el simulador es el tiempo de vuelo que

presentaran los fluidos en cada instante de tiempo, como se observa en la
figura 80. Mediante ella se puede concluir que las reas laterales superior del
modelo (lneas amarillas y verdes) son zonas difciles de barrer y que requieren
de aproximadamente 7000 a 8000 das de inyeccin para lograr que se
experimente un efecto notorio sobre esta zona.
171
Figura 80. Tiempo de vuelo obtenido para el modelo establecido al cabo de los 3650
das de inyeccin.
o Influencia de la gravedad. En este caso planteado, por las caractersticas

estructurales y de heterogeneidad del modelo, el efecto de la gravedad juega
un papel importante en la prediccin del comportamiento del modelo,
especialmente en lo que se refiere a la determinacin de la presin y de la
distribucin de saturaciones en el modelo. Este efecto se observa en las figuras
81 y 82.
Figura 81. Izquierda: distribucin de saturacin sin efectos de gravedad.
Derecha: distribucin de saturacin con efectos de gravedad. Modelo heterogneo.

172
Figura 82. Distribucin de presiones hallada por el simulador streamline con y sin
efectos de gravedad para el caso base dos.
Para las figuras 81 y 82 el efecto de la gravedad se hace un poco notorio en la

forma del frente de saturacin y en la diferencia de presin de 60 PSIA, hallada
en la distribucin de presin. Para la prediccin del comportamiento de
produccin los resultados obtenidos son bastante confiables y divergen en un
pequeo porcentaje, evidenciando la validez del modelo construido (ver figuras
83, 84, 85, 86 y 87).
Figura 83. FOIP para el caso dos con y sin efectos de gravedad.
173
Figura 84. Tasa de produccin de aceite para el caso dos con y sin efectos de
gravedad.
Figura 85. Produccin total de aceite para el caso dos con y sin efectos de gravedad.
Figura 86. Inyeccin total de agua para el caso dos con y sin efectos de gravedad.
174
Figura 87. Produccin total de agua para el caso dos con y sin efectos de gravedad.
o Determinacin de la distribucin de saturacin en el modelo. Para los tres

tamaos en los pasos de tiempo considerados para las corridas, se obtuvieron
las siguientes distribuciones de saturacin, evidencindose nuevamente, pero
de una forma ms marcada el efecto de difusin numrica en la forma externa
del frente de invasin y en los valores de saturacin de algunas celdas, a
medida que se aumentan los pasos de tiempo en la simulacin.
Figura 88. Saturacin para un solo Paso de tiempo considerado. Tp = 3650 das.
Modelo heterogneo.
175
Figura 89. Saturacin para 4 pasos de tiempo considerados. Tp=912.5 das. Modelo
heterogneo.
Figura 90. Saturacin para 8 pasos de tiempo considerados. Tp= 456.25 das. Modelo
heterogneo.
176
o Prediccin del comportamiento de la produccin. Aunque al aumentar el

nmero de pasos de tiempo se originan efectos de difusin numrica, estos
errores son relativamente pequeos y el simulador busca reducirlos mediante
algoritmos internos. No obstante para casos heterogneos como el
considerado, se debe hacer un anlisis ms detallado a fin de determinar el
nmero de pasos de tiempo ptimo, sin dejar a un lado la ventaja de
seleccionar pasos de tiempo algo extensos. Los comportamientos obtenidos
para cada uno de los pasos de tiempo seleccionados son bastante consistentes
y aproximados y se observan en las siguientes figuras:
Figura 91. Errores en el balance de materiales tomando diferentes tamaos en los

pasos de tiempo. Modelo heterogneo
Figura 92. Corte de agua tomando diferentes tamaos en los pasos de tiempo. Modelo
heterogneo.
177
Figura 93. Inyeccin total de agua tomando diferentes tamaos en los pasos de
tiempo. Modelo heterogneo.
Figura 94. OIP tomando diferentes tamaos en los pasos de tiempo. Modelo dos.
Figura 95. Produccin total de aceite tomando diferentes tamaos en los pasos de
178
Figura 96. Tasa de produccin de aceite tomando diferentes tamaos en los pasos de
Figura 97. Produccin total de agua tomando diferentes tamaos en los pasos de
tiempo. Caso base dos.
179
8. PERSPECTIVA DE LA SIMULACIN STREAMLINE EN COLOMBIA
En la actualidad, teniendo en cuenta la gran necesidad que requiere Colombia

en cuanto a mejorar el rendimiento y produccin de sus yacimientos, una mejor
comprensin y definicin de cada uno de ellos a partir de tecnologas como la
simulacin streamline, puede ayudar al pas a aprovechar al mximo sus
campos, y as mismo garantizar a aquellos que estn por venir, un mejor
desarrollo y operacin.
Las amplias ventajas que brinda la simulacin streamline de la mano de

nuevas tecnologas como son procesos geoestadsticos, algoritmos genticos,
mtodos mejorados de integracin de datos entre otros, en los ltimos aos
han sido aplicadas con xito en algunos campos de nuestro territorio,
especialmente yacimientos sometidos a inyeccin de agua.
En recientes estudios realizados en conjunto con la UIS y ECOPETROL, se

han generado investigaciones que interactan los procesos anteriormente
mencionados, como es el caso del estudio realizado en un piloto de inyeccin
del campo La Cira Infantas. En este estudio diseado para cuantificar la
incertidumbre asociada en la prediccin de parmetros de produccin, se
integraron el modelamiento geoestadstico con la simulacin streamline,
generando un procedimiento bastante robusto que permiti predecir
eficientemente cortes de agua y tasas de produccin. Este estudio realizado
por el Ingeniero Jorge Mantilla V* (2003), fue tan exitoso que en la actualidad
toda el rea 7 del campo la Cira Infantas ya ha sido simulado mediante est
tcnica, arrojando buenas recomendaciones para mejorar el proceso de
inyeccin-produccin.
Otro interesante proyecto fue el desarrollado por el Ingeniero Jos Arnobio
Vargas Medina** (2003), quin combin simulacin streamline, modelamiento
Geoestadstico y tcnicas de inversin dinmica de datos en un modelo
sinttico, a fin de obtener una novedosa metodologa para optimizar el modelo
de permeabilidad de yacimientos altamente heterogneos, obtenindose
nuevamente resultados satisfactorios. Segn el autor, esta metodologa permite
la apropiacin para la industria petrolera nacional de una tecnologa de punta
til en la caracterizacin de nuestros yacimientos.
____________
MANTILLA, V., JORGE. Cuantificacin de la Incertidumbre Asociada en la Prediccin del

Comportamiento de Produccin de un Yacimiento altamente Heterogneo Sometido a Inyeccin de
Agua. Universidad Industrial de Santander. Bucaramanga. 2003.
**
VARGAS MEDINA, JOSE ARNOBIO, Optimizacin del Modelo de Permeabilidad de un Yacimiento
Heterogneo Mediante Inversin Dinmica de Datos Basada en Simulacin Streamline. Universidad
Industrial de Santander. Bucaramanga. 2003.
180
Concluyendo con este grupo de estudios por su parte, el Ingeniero Santiago

Gonzlez
Fernndez*
(2003)
desarroll
un
programa
llamado
INVERDINMICA, generado en plataforma PC y bajo ambiente Windows que
integra est vez simulacin streamline con formulaciones de algoritmos
genticos aplicables a yacimientos en produccin primaria como en
yacimientos sometidos a inyeccin de agua. La metodologa desarrollada
mostr ser una tcnica apropiada para encontrar descripciones de yacimientos
que hacen honor a los datos de produccin, informacin de los pozos y a los
rasgos geolgicos, adems de manejar gran cantidad de parmetros, lo cual es
un problema cuando se trabaja en modelos de simulacin de yacimientos.
Segn el autor est metodologa ayudar a los ingenieros de yacimientos a

realizar el ajuste histrico del modelo propuesto, disminuyendo las horashombre requeridas para realizar este proceso.
En cuanto a empresas del sector privado, la OXY, est aplicando actualmente

la simulacin streamline en Cao-Limn y usa el simulador de Schlumberger
Frontsim.
Muchos de nuestros campos, sobre todo los ms antiguos han experimentado
un marcado efecto de disminucin de su produccin. Si a esto se le agrega que
no se cuenta con un modelo de yacimiento bien establecido, que describa tanto
esttica como dinmicamente el comportamiento real de los fluidos, no se
puede esperar un buen aprovechamiento de estos recursos.
La simulacin streamline sera clave en estos casos, permitiendo una buena

definicin y comprensin de nuevos modelos que realmente representen estos
fenmenos y que adems brinden la posibilidad de mejorar o innovar los
procesos de recobro con el fin ltimo: aumentar en forma eficiente la
produccin.
Como es de apreciarse, esta tcnica tiene un excelente futuro en nuestro pas

debido a su aprovechamiento, novedad y su rapidez de solucin. Actualmente
se utiliza para procesos de inyeccin-produccin, pero ya se han hecho
desarrollos para aplicar a procesos relacionados con deplecin por gas en
solucin.
____________
*
GONZLEZ FERNANDEZ, SANTIAGO. Inversin Dinmica de Datos para Caracterizar un Yacimiento
Fluvial Usando Geoestadstica y Algoritmos Genticos. Universidad Industrial de Santander.
Bucaramanga. 2003.
181
Dadas las amplias caractersticas de heterogeneidad de nuestros yacimientos,

sus extensiones, complejidad geolgica, y la inmensa mayora de campos
sometidos a procesos inmiscibles, esta rea de la ingeniera se perfila como
una herramienta clave en el buen manejo y prediccin del comportamiento de
nuestros campos.
182
9. CONCLUSIONES Y RECOMENDACIONES.
9.1 CONCLUSIONES
Los simuladores streamline como herramientas en la caracterizacin de

yacimientos, especialmente modelos altamente heterogneos, son
ampliamente ventajosos, ya que permiten simular procesos basados
directamente en el modelo geolgico original, empleando sencillamente un
procesador convencional.
La simulacin streamline sobrepasa los lmites establecidos por los

simuladores convencionales: permite trabajar con mallas grandes y a pasos de
tiempo extensos, haciendo eficiente uso de la memoria disponible.
La diferencia en tiempo de CPU entre una simulacin streamline y una

simulacin de diferencias finitas es debida a la habilidad del mtodo streamline
para desacoplar un problema 3D en una serie de problemas unidimensionales,
que solo son funcin de las propiedades geolgicas de las celdas, el tiempo de
vuelo y los flujos fraccionales de los fluidos presentes. Los efectos de
orientacin del grid son substancialmente reducidos en el mtodo streamline.
La tecnologa streamline brinda nuevas caractersticas grficas como la

visualizacin del tiempo de vuelo y las zonas de influencia de inyectores y
productores. Estos nuevos tipos de resultados pueden ser ampliamente tiles a
la hora de verificar procesos de upscaling, consistencia de modelos geolgicos
y dinmicos, clasificacin de modelos geoestadsticos y optimizacin de los
patrones de inyeccin-produccin establecidos, entre otros.
Los simuladores streamline han demostrado ser eficientes para pasos de

tiempo largos y extensos. Aunque tomando pasos de tiempo cortos se
evidencia en la distribucin de saturacin hallada un leve efecto de difusin
numrica, los resultados en cuanto a prediccin del comportamiento de la
produccin arrojados por el simulador son ampliamente consistentes y
congruentes con el modelo establecido. Si en el modelo establecido no se
varan en forma significativa las condiciones lmites o de operacin, y adems
el simulador brinda la opcin de solucin analtica, est posibilidad origina
soluciones ampliamente satisfactorias y libres de efectos de difusin.
183
El efecto de la gravedad es un factor que de una u otra forma puede influir en el

desplazamiento de los fluidos, especialmente, si estos son totalmente
diferentes en trminos de densidad y composicin y si el modelo a simular es
ampliamente estratificado.
Aunque la metodologa streamline en los ltimos aos se ha extendido y se ha
desarrollado para trabajar procesos compresibles y composicionales, su amplio
campo de aplicacin se ha mantenido en el rea de procesos de inyeccin de
agua, debido a que son sistemas de desplazamiento altamente estables y muy
acordes con el mtodo streamline.
Los cambios constantes en las condiciones lmites o de operacin del modelo

establecido, llevan implcito errores en el balance de materiales, reduciendo de
sta forma la eficiencia del mtodo streamline. No obstante, cada simulador
streamline posee algoritmos especiales que buscan reducir este porcentaje, o
dado el caso detener el proceso de simulacin cuando el error al final de algn
paso de tiempo es muy alto.
En sistemas altamente compresibles, la ventaja inherente de velocidad del

simulador streamline se ve reducida debido a la necesidad de actualizacin
peridica de la presin y de las lneas de flujo. Lo anterior no quiere expresar
error en los resultados obtenidos por el simulador streamline, solo que su
comportamiento en cuanto a velocidad sera semejante al esperado en un
simulador convencional.
Mediante los casos base establecidos, se verific la efectividad, veracidad y

aplicabilidad del mtodo de simulacin streamline. Se confront la teora con la
prctica, obtenindose los resultados esperados. Sin embargo, se debe resaltar
que cada tipo de simulador streamline comercialmente disponible, posee su
forma de tratamiento de informacin y anlisis de datos. 3DSL es muy amplio y
en los ltimos aos le han sido implementadas muchas aplicaciones no
encontradas en Frontsim, que aumentan las ventajas de este tipo de
simulacin.
Mecanismos dominados por fuertes efectos de presin capilar originan flujo

perpendicular a las streamlines. En estos casos, el uso de un simulador de
diferencias finitas es la opcin ms recomendable.
184
9.2 RECOMENDACIONES.
En la actualidad, se han desarrollado modelos matemticos streamline

adaptados para tratar sistemas capilares y gravitacionales. El desarrollo de un
software que tenga en cuenta estos aspectos sera una alternativa interesante.
Cuando se requieran trabajar procesos altamente no-lineales, es aconsejable

seleccionar el tamao y el nmero de pasos de tiempo ptimo, basados en los
datos arrojados por el simulador y el comportamiento histrico-real del campo.
Ninguno de los simuladores streamline comercialmente reconocidos brinda la
posibilidad de trabajar mallas no estructuradas ni tipo PEBI. El planteamiento
matemtico ya ha sido desarrollado y sera excelente adaptarlo a algn
simulador.
La simulacin streamline es solo una de las herramientas aplicadas en la

actualidad para la caracterizacin de yacimientos. La sistematizacin de
algoritmos que integren en un solo programa que incluya simulacin
geoestadstica, simulacin streamline y esquemas de inversin dinmica puede
ser ventajoso y se evitara la manipulacin y transferencia manual de
informacin entre las subrutinas.
Difundir en el medio estudiantil este mtodo de simulacin, presentando las

amplias aplicaciones y ventajas que posee, pero sobre todo el amplio campo
de aplicacin que tiene en nuestro pas por ser una tecnologa nueva, que
responde a las necesidades de optimizacin de los procesos de recobro y por
ende al mejoramiento de la produccin de los campos.
185
BIBLIOGRAFIA
AGARWAL, BIJAN, BLUNT, MARTIN. A Streamline-Based Method for

Assisted History-Matching Applied to an Arabian Gulf Field. SPE Annual
Technical Conference. And Exhibition. Denver, Colorado. 5-8 oct. 2003. SPE
N 84462
AL-HUTHALI, ARAMCO, AHMED, SAUDI, DATTA-GUPTA, AKHIL. Streamline

Simulation of Water Injection in Naturally Fractured Reservoirs. Symposium on
Improved Oil Recovery. Oklahoma. 17-21 abril. 2004. SPE N 89443.
ATES, HARUN, BAHAR, ASNUL, CHARFEDDINE, MOHSEN MOHAN,

KELKAR, DATTA-GUPTA, AKHIL. Ranking And Upscaling Of Geostatistical
Reservoir Models Using Streamline Simulation: A Field Case Study. Texas
A&M University 9-12 Jun 2003. SPE N 81497
BAKER, RICHARD. Streamline Techology: Reservoir History Matching and

Forecasting= its Success, Limitations and Future. JCPT. Abril. 2001.
-----------Full-Field Modeling Using Streamline-Based Simulation: 4 Case
Studies. Streamsim Technologies Inc. SPE 66405. Texas, Feb. 2001.
BARTHELEMY, MAUD, BLUNT, MARTIN, FENWICK, DARRYL,. StreamlineBased History Matching Using Sensitivity Computations and the Gradual
Deformation Method. Institut Francais Du Petrole, Imperial Collage. Jun. 2003.
BATENBURG, VAN, BRUINING, J., BAKKEN, H. The Effect of Capillary

Forces in Heterogeneous Flow-Units: A Streamline Approach. SPE Annual
Technical Conference. And Exhibition.Texas, oct. 5-9 1991. SPE N 22588.
BATYCKY, ROD P., THIELE R, MARCO, BLUNT, MARTIN J.. A StreamlineBased

Reservoir
Simulation
of
the
House
Mountain
Waterflood. Stanford University - SUPRIC Research Group. SCRF, 1997.
---------------A Three-Dimensional Two-Phase Field Scale Streamline
Simulator. Stanford University. Jun 1997. Pg 15-18.
BERENBLYUM, R.A., SHAPIRO, A., JESSEN, K., STENBY, E., ORR, F. Black
Oil Streamline Simulator with Capillary Effects. Technical University of
Denmark and Stanford University. SPE 84037. Oct. 2003.
186
BRATVEDT, F., GIMSE, T. AND TEGNANDER, C. Streamline Computations

for Porous Media Flow Including Gravity. Transport in Porous Media, Vol. 25,
N 1. Pg. 63-78. Oct 1996.
CHAWATH, ADWAIT, TAGGART, LAN. Insights Into Upscaling Using 3D

Streamline. Reservoir Simulation Symposium. Texas, 11-14 feb. 2001. SPE N
66379.
CHRISTIE, MIKE, SUBBEY, SAM, SAMBRIDGE, MALCOLM, THIELE R,

MARCO. Quantifying Prediction Uncertainty in Reservoir Modelling Using
Streamline Simulation. 15 ASCE Engineering Mechanics Conference.
Columbia University, New York. 2-5 Jun. 2002.
DATTA-GUPTA, AKHIL. Streamline Simulation: A Technology Update. SPE,

Texas A&M University. Dic. 2000.
DONATO, GINEBRA, BLUNT, MARTIN. A Streamline-Based Dual Porosity

Simulator for Fracture Flow Simulation. ITF Project. Jun. 9 2003.
-------------Streamline-Based Dual Porosity Simulation of Fractured Reservoirs.
Annual Technical Conference and Exhibition. Colorado, 5 8 Oct 2003. SPE N
84036.
DUPIN, RAVALEC, DARRYL, MICKAELE H., FENWICK. A Combined

Geostatistical and Streamline-Based History-Matching Procedure. SPE Annual
Technical Conference. And Exhibition. Texas, 29 sep. 2 oct. 2002.
EMANUEL, A.S, AND MILLIKEN, W.J. Application of 3D Streamline Simulation

to Assist History Matching. SPE 49000. New Orlens. Oct 1998.
GEOQUEST, FRONTSIM Case Study: FrontSim Streamline Simulation. Sept

2001.
GILMAN, J., MENG, H., ULAND, M., DZURMAN, P., COSIC, S. Statistical
Ranking of Stochastic Geomodels Using Streamline Simulation: A Field
Application. Annual Technical Conference and Exhibition. SPE 77374. Texas.
29 Sep-2 Oct 2002.
187
GONZLEZ FERNANDEZ, SANTIAGO. Inversin Dinmica de Datos para

Caracterizar un Yacimiento Fluvial Usando Geoestadstica y Algoritmos
Genticos. Universidad Industrial de Santander. Bucaramanga. 2003.
HARTEN, A. High Resolution Schemes for Hyperbolic Conservation Laws. J.

Comp. Phys. Vol. 49 Pg 357. 1983.
HE, ZHONG, DATTA-GUPTA, AKHIL, YOON, SEONGSIK. Streamline-Based

Production Data Integration Under Changing Field Conditions Annual
Technical Conference and Exhibition held in New Orleans, Louisiana, 30 Sep.3
Oct. 2001. SPE N 71333
IDROBO, EDUARDO A., CHOUDHARY, MANOJ K., DATTA-GUPTA, A.

Swept Volume Calculations And Ranking Of Geostatistical Reservoir Models
Using Streamline Simulation. Texas A&M University. California, 1923 June
2000. SPE N 62557.
INGEBRIGTSEN, L., BRATVEDT, F., BERGE, J. A Streamline Based

Approach to Solution of Three-Phase Flow. Reservoir Simulation Symposium
Held in Houston. Texas 14-17 Feb. 1999. SPE N 51904.
KAZEMI, H. Analytical and Numerical Solution of Oil Recovery from Fractured

Reservoirs with Empirical Transfer Functions. SPEJ. May 1992. Pg 219.
KING, M.J., MANSFIELD, MARK. Flow Simulation of Geologic Models .BP

Exploration Ltd. 4 de Agos 1999.
LITVAK, B.L. Simulation and Characterization of Naturally Fractured

Reservoirs. Dallas. April 29-May 1 1985.
LOLOMARI, TONY, BRATVEDT, KYRRE, CRANE, MARTHA, MILLIKEN,

WILLIAM J., TYRIE, J.J. The Use of Streamline Simulation in Reservoir
Management: Methodology and Case Studies. SPE Annual Technical
Conference. And Exhibition. Texas 1-4 oct. 2000. SPE N 63157.
---------------The Use of Streamline Simulation in Reservoir Management:
Methodology and Case Studies. SPE 63157. Texas. Oct 2000.
M.A. CHRISTIE, HERIOT-WATT U, BLUNT, M.J. Tenth SPE Comparative

Solution Project: A Comparison of Upscaling Techniques. SPE. Agos. 2001.
188
MALLISON, B.T, GERRITSEN, M.G, MATRINGE, S.F. Improved Mappings for

Streamline-Based Simulation. Symposium on Improved Oil Recovery.
Oklahoma, 17-21 abril. 2004. SPE N 89352.
MANTILLA, V., JORGE. Cuantificacin de la Incertidumbre Asociada en la

Prediccin del Comportamiento de Produccin de un Yacimiento altamente
Heterogneo Sometido a Inyeccin de Agua. Universidad Industrial de
Santander. Bucaramanga. 2003.
MANUAL ECLIPSE FRONTSIM GEOQUEST.SCHLUMBERGER. 2001.
MANUAL SIMULADOR 3DSL. STREAMSIM TECNOLOGIES INC. Version 2.20

Agost 2004.
MASCHIO, CLIO, SCHIOZER, DENIS JOS. Analysis of the Performance of

Streamline Simulation. UNICAMP. Buenos Aires - Argentina, Nov 5-6, 2002.
MARATHE, R.V. A Novel Streamtube Model for Tracer Flow In An Areally

Heterogeneous Reservoir. Oil and Gas conference and Exhibition. India. 17-19
Feb 1998. SPE N 39550.
MATHIEU, MICHAEL G.EDWARDS, BLUNT, MARTIN J. Streamline Tracing

On Curvilinear Structured And Unstructured Grids. Imperial College,London
.Jun 2002.
MATRINGE, S.F, GERRITSEN, M.G, On Accurate Tracing Of Streamlines.

SPE Annual Technical Conference. And Exhibition.Texas 26-29 Sept 2004 SPE
N 89920.
PEDIDIHOPTIA, S., DATTA-GUPTA, A., XUE, GUOPING. Multiphase

Streamline Modeling In Three Dimensions: Further Generalizations and a Field
Application Reservoir Simulation Symposium. Texas 6-11 Jun 1997. SPE N
38003
----------------------An Efficient PC Based Streamline Simulator for Immiscible
and Miscible Displacements. Petroleum Computing Conference. Texas, 8-11
Jun. 1997. SPE N 38129.
POLLOCK, D.W. Semianalytical Computation of Path Lines for FiniteDifference Models. Ground Water. Nov-Dic 1988. 26, N 6, Pg 743-750.
189
PONTING, D. K. Hybrid Streamline Methods. Asia Pacific Conference on

Integrated Modelling for Asset Management. Malasya, 23-24 marzo. 1998. SPE
N 39756.
PORTELLA, R.C.M., HEWETT, T.A. Upscaling, Gridding, and Simulating Using
Streamtubes .Stanford U. Sept 2000.
RUIJIAN, LI, REYNOLDS, A.C., OLIVER, D.S. History-Matching of ThreePhase Flow Production Data. Reservoir Simulation Symposium Held in
Houston. 11-14 feb. 2003.
SAMIER, QUETTIER, THIELE, MARCO. Applications Of Streamline

Simulations To Reservoir Studies SPE Reservoir Simulation Symposium Held
In Houston, Texas, 1114 Feb 2001. SPE N 66362.
SANTA F, ELKIN, ENELSO, LUIS. Estudio Comparativo de las Tcnicas de

Enmallado Empleadas en simulacin de Yacimientos. Universidad Industrial de
Santander. Nov 2004.
SARKAR, AMIT K., FARUQI, SOHAIL A., BARMAN, I., WANG, M. K. Efficient
Use of Streamline Simulation and Inversion Techniques to Improve Reservoir
Modeling. SPE Annual Technical Conference. And Exhibition. Texas 1-4 oct.
2000. SPE N 63291.
SONIER, F. Numerical Simulation of Naturally Fractured Reservoirs. SPE

15627. New Orleans. Oct 1986.
SWEBY, P.K. High Resolution Schemes Using Flux Limiters for Hyperbolic
Conservation Laws. SIAM J. Numerical Anlisis Vol 21. Pg 995-1011. 1990.
THIELE R, MARCO. Streamline Simulation. 6 Internacional Forum on

Reservoir Simulation. Schloss Fusch, Austria. 3-7 Sep. 2001.
---------------------------Modeling Flow In heterogeneous Media Using
Streamtubes. Miscible And Immiscble Displacements Department of
Petroleum Engineering Stanford University, Stanford California 1995
--------------------------- Simulating Flow In Heterogeneous Systems Using
Streamtubes And Streamlines . SPE, Stanford University. Feb 1996.
-----------------------A Streamline-Based 3D Field-Scale Compositional Reservoir
Simulator. Annual Technical Conference and Exhibition. Texas, Oct 5-8, 1997.
SPE N 38889.
190
-----------------------Modeling Flow In heterogeneous Media Using Streamtubes I.

Miscible and Immiscible Displacements. Stanford University. Marzo 3, 1995.
----------------------- Modeling Multiphase Flow in Heterogeneous Media Using
Streamtubes. Dic 1994.
TRAN, THOMAS T , WEN, XIAN-HUAN, BEHRENS, RONALD A. Efficient

Conditioning of 3D Fine-Scale Reservoir Model To Multiphase Production Data
Using Streamline-Based Coarse-Scale Inversion and Geostatistical
Downscaling. Annual Technical Conference and Exhibition. Texas, 36 Octr
1999. SPE N 56518.
VARGAS, JOSE A. Optimizacin Del Modelo De Permeabilidad De Un

Yacimiento Heterogneo Mediante Inversin Dinmica De Datos Basada En
Simulacin Streamline. Universidad Industrial De Santander, Bucaramanga
2003.
VASCO, D.W, YOON, S., DATTA-GUPTA, A. Integrating Dynamic Data Into

High-Resolution Reservoir Models Using Streamline-Based Analytic Sensitivity
Coefficients. Texas A&M U. Dic. 1999.
WANG, YUANDONG, KOVSCEK, ANTHONY R. A Streamline Approach For

History-Matching Production Data . Improved Oil Recovery Symposium held in
Tulsa, Oklahoma, 35 Abril 2000. SPE N 59370
-----------------------------A Streamline Approach For Ranking Reservoir Models
That Incorporates Production History .. SPE Annual Technical Conference.
And Exhibition held in San Antonio. Texas 29 Sept-2 Oct 2002. SPE N 77377.
-----------------------------A Streamline Approach For History-Matching Production
Data. Improved Oil Recovery Symposium held in Tulsa, Oklahoma, 3 5 April
2000. SPE N 59370.
WEI, Y., MICHELSEN, M., STENBY, E., BERENBLYUM, R., SHAPIRO, A.

Three-phase Compositional Streamline Simulation and Its Application to
WAG. University of Denmark. SPE 89440. April 2004.
WU, ZHAN, DATTA-GUPTA, A. Rapid History Matching Using a Generalized

Travel-Time Inversion Method .Texas A&M U. SPE. Jun 2002
191

Simulacion Streamline

Cargado por

Copyright:

Formatos disponibles

Simulacion Streamline

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 Streamline

Cargado por

Copyright:

Formatos disponibles

SIMULACIN STREAMLINE: ESTADO DEL ARTE Y APLICACIONES

FUNDAMENTALES EN EL MODELAMIENTO Y ESTUDIO DE YACIMIENTOS

ADRIANA MARCELA MNDEZ BOHRQUEZ

UNIVERSIDAD INDUSTRIAL DE SANTANDER

SIMULACIN STREAMLINE: ESTADO DEL ARTE Y APLICACIONES

ADRIANA MARCELA MNDEZ BOHRQUEZ

Tesis para optar al ttulo de

UNIVERSIDAD INDUSTRIAL DE SANTANDER

A mi madre, por el don de la vida, por su incansable entrega, amor y esfuerzo,

A mi hermano Vctor por su compaa y cario y por ensearme a su corta

A Evangelista, por su apoyo y ayuda en los duros y gratos momentos y por la

A toda mi familia, por los ms cercanos y los ms distantes, porque siempre

A mis luces celestiales, quienes desde el ms all me aportaron su ayuda

A mis amigos, compaeros de estudios y aventuras, a ellos gracias por su

A mi madre porque me dio el regalo ms preciado que es la vida.

A mis tos Hector Melndez y Samuel Gmez por su ayuda y colaboracin ya

A Marlene por su amor y apoyo incondicional sin esperar nada a cambio.

A mis primos en especial a Joan Sebastin y Eduardo, por aceptarme como un

A mi compaera, por su fuerza de voluntad, ya que sin ella este trabajo no se

Los ingenieros Nicols Santos y Fernando Calvete, por la atencin y

A la Escuela de Ingeniera de Petrleos, en cabeza de su directora, Dra. Zuly

A los ingenieros Jorge Mantilla V. y Jos Arnobio Vargas Medina, quienes

1.1 PLANTEAMIENTO DEL PROBLEMA

1.2 DESPLAZAMIENTOS DOMINADOS POR CONVECCIN

1.3 SIMULACIN STREAMLINE BASADA EN EL FLUJO

1.3.1 Contexto histrico.

1.3.2 Formulacin general de la simulacin streamline.

1.3.3 Simulacin streamline vs. simulacin en diferencias finitas.

1.3.4 Ventajas de la simulacin streamline.

1.3.5 Desventajas de la simulacin streamline.

2. SIMULACIN STREAMTUBES: PRIMEROS PLANTEAMIENTOS

2.1.1 La aproximacin de Riemann en la formulacin streamtubes.

2.2 CLASES DE PROBLEMAS

2.3 MODELO MATEMTICO

2.3.1 Condiciones lmites.

2.3.2 El papel de la streamfuncin en simulacin streamtubes.

2.4 TRAZADO DE SOLUCIONES 1D

2.4.3 Trazado de las soluciones 1D sobre el grid cartesiano 2D.

2.5 SOLUCIONES EN DOS DIMENSIONES

3. CONCEPTOS Y MODELO MATEMTICO DE LA SIMULACIN

3.1 DEDUCCIN DE LAS ECUACIONES

3.2 SOLUCIN A LA ECUACIN DE PRESIN

3.3 CONDICIONES LMITES

3.4 SOLUCIN NUMRICA DE LA MATRIZ DE PRESIN

3.5 DETERMINACIN DE LA VELOCIDAD DEL CAMPO

3.6 MTODO PARA LA DETERMINACIN DE LAS VAS Y EL

3.7 EL TIEMPO DE VUELO

3.8 TRANSFORMACIN DE LA COORDENADA A LO LARGO

3.9 SOLUCIN ANALTICA DE LA ECUACIN DE CONTINUIDAD

3.9.1 Clculo de las propiedades en cada bloque del grid de simulacin

3.9.2 Bloques que no son atravezados por streamlines.

3.9.3 Clculo del tiempo real.

3.10 SOLUCIN NUMRICA DE LA ECUACIN DE CONTINUIDAD

3.10.1 Planteamiento de las nuevas soluciones numricas.

3.10.2 Bloques que no son atravesados por streamlines.

3.10.3 Errores en el balance de materiales.

3.11 ACTUALIZACIN DEL MODELO DE LNEAS DE FLUJO

3.12 PASOS DE TIEMPO

3.13 EL CONCEPTO DE GRAVEDAD EN SIMULACIN STREAMLINE