2019 1 Ip
2019 1 Ip
2019 1 Ip
Director
JUAN EDUARDO RIVERA DE LA OSSA
Ingeniero de Petróleos
________________________________
________________________________
________________________________
________________________________
________________________________
________________________________
________________________________
________________________________
________________________________
3
DIRECTIVAS DE LA UNIVERSIDAD
4
Las directivas de la Universidad de América, los jurados calificadores y los docentes
no son responsables por las ideas y conceptos emitidos en el presente documento.
Estos corresponden únicamente al autor.
5
DEDICATORIAS
6
AGRADECIMIENTOS
7
CONTENIDO
pág.
RESUMEN 25
INTRODUCCIÓN 26
OBJETIVOS 28
2. MARCO TÉORICO 44
2.1 GENERALIDADES DE LA INYECCIÓN WATER ALTERNATING GAS 44
2.2 HISTORIA DE LA INYECCIÓN CÍCLICA DE AGUA Y GAS (WAG) 46
2.3 FACTORES QUE AFECTAN LA INYECCIÓN WAG 48
2.3.1 Características geológicas y físicas del yacimiento 48
2.3.2 Propiedades de los fluidos 50
2.2.3 Relación entre las movilidades 51
2.3 GENERALIDADES DE LA SIMULACIÓN NÚMERICA DE YACIMIENTOS 53
8
2.3.1 Clasificación de los simuladores de yacimientos 55
2.3.2 Ecuaciones que usa el simulador 59
9
4.3.6 Permeabilidades relativas 90
4.4.1 Array properties 97
4.4 CONFIGURACIÓN DEL MODELO 99
4.4.2 Generalidades del modelo 99
4.4.4 Rock Fluid – Relación roca fluido 109
4.4.5 Carga del modelo dinámico 111
4.4.6 Sección condiciones iniciales 117
4.4.7 Petróleo original in situ 118
4.5 COTEJO HISTÓRICO 119
4.5.1 Proceso de sensibilidades 119
7. CONCLUSIONES 146
8. RECOMENDACIONES 148
10
BIBLIOGRAFÍA 149
ANEXOS 152
11
LISTA DE FIGURAS
pág.
Figura 1. Localización del área de interés 32
Figura 2. Columna estratigráfica generalizada para la cuenca del Valle Medio
del Magdalena 35
Figura 3. Proceso de inyección alternada de gas y agua 45
Figura 4. Distribución del frente de invasión para la inyección WAG 49
Figura 5. Espesor desfavorable para la inyección WAG 50
Figura 6. Diseño planta de inyección de agua del campo 68
Figura 7. Área piloto de inyección de agua 70
Figura 8. Mapa estructural del campo objeto de estudio 77
Figura 9. Sector Model escogido por la compañía 78
Figura 10. Ubicación de pozos 79
Figura 11. Líneas base de arenas y arcillas para las Arenas B de la
Formación Mugrosa 80
Figura 12. Espesores del modelo 98
Figura 13. Net to Gross ratio 100
Figura 14. Porosidad del modelo 101
Figura 15. Porosidad promedio arrojada por Builder 102
Figura 16. Permeabilidad del modelo 104
Figura 17. Promedio de permeabilidad del modelo 105
Figura 18. Permeabilidad vertical 106
Figura 19. Permeabilidad vertical promedio 106
Figura 20. Saturación de agua del modelo 109
Figura 21. Permeabilidades relativas para las 3 fases en el tipo de roca 1 111
Figura 22. Composición del fluido en WinProp 112
Figura 23. Presión crítica, Temperatura crítica y factores acéntricos 113
Figura 24. Fracción molar de componentes en WinProp 114
Figura 25. Composition, WinProp 114
Figura 26. Prueba de liberación diferencial en WinProp 115
Figura 27. Sección condiciones iniciales para Builder 118
Figura 28. Factor de recobro de porcentajes de inyección WAG 126
Figura 29. Respuesta de la inyección WAG comparada con la inyección de
agua en el arreglo I-57 129
Figura 30. Producción de gas asociado al pronóstico de producción de los
próximos 10 años 131
Figura 31. Corte de agua para los próximos 10 años 132
Figura 32. Flujo de caja neto del proyecto 144
Figura 33. Permeabilidad relativa gas-petróleo 155
Figura 34. Curva de permeabilidad relativa gas-petróleo 156
12
LISTA DE TABLAS
pág.
13
Tabla 37. Inversión 137
Tabla 38. Secciones requeridas para transporte del Flue Gas desde la
Refineria de Barrancabermeja hasta disposición. 138
Tabla 39. Mano de obra de instalación de tubería 138
Tabla 40. Costos de inyección de agua 139
Tabla 41. Costo de Flue Gas proveniente de la Refinería de
Barrancabermeja. 140
Tabla 42. Costos operacionales 141
Tabla 43. Flujo de caja operativo del proyecto para 10 años 143
Tabla 44. Datos de porosidad y propiedades de la roca. 153
Tabla 45. Datos de porosidad para las arenas b de la formación mugrosa. 154
Tabla 46. Cronograma de inyección WAG. 157
14
LISTA DE ECUACIONES
pág.
15
LISTA DE GRÁFICAS
pág.
Gráfica 1. Historial de producción acumulada 42
Gráfica 2. Uso de la inyección Water Alternating Gas en el mundo 47
Gráfica 3. Producción anual histórica del campo objeto de estudio 64
Gráfica 4. Historial de inyección del patrón I-35 71
Gráfica 5. Historial de inyección del patrón I-57 73
Gráfica 6. Historial de inyección arreglo I-111 74
Gráfica 7. Comportamiento del punto de burbuja para cada % molar de
inyección 86
Gráfica 8. Comportamiento del factor de hinchamiento en el análisis PVT. 88
Gráfica 9. Viscosidad del fluido con respecto a la presión. 89
Gráfica 10. Rs fluido 90
Gráfica 11. Curvas de permeabilidades relativas para los diferentes tipos de
roca. 97
Gráfica 12. Permeabilidades relativas del aceite y el agua para el tipo de roca
1. 110
Gráfica 13. Viscosidad de los fluidos vs presión en WinProp 116
Gráfica 14. Bo y GOR vs presión en WinProp 117
Gráfica 15. Sensibilidades de BHP para ajuste histórico 121
Gráfica 16. Producción acumulada para ajuste histórico. 122
Gráfica 17. Tasa liquida para ajuste histórico. 123
Gráfica 18. Tasa de petróleo para ajuste histórico. 124
Gráfica 19. Pronóstico de producción a 10 años 130
16
LISTA DE ANEXOS
Anexo A. Datos de porosidad y propiedades de la roca. 153
Anexo B. Datos de porosidad para las arenas b de la formación mugrosa. 154
Anexo C. Permeabilidad relativa gas-petróleo. 155
Anexo D. Curva de permeabilidad relativa gas-petróleo. 156
Anexo E. Cronograma de inyección WAG. 157
17
LISTA DE ABREVIATURAS
18
Sw Saturación del agua
Swi Saturación inicial de agua
WAG Water Alternating Gas (Inyección alternada de gas y agua)
19
GLOSARIO
1
SCHLUMBERGER. Arcilla. Oil Field Glossary en linea], 2018 [revisado 19 de noviembre 2018].
Disponible en Internet: https://www.glossary.oilfield.slb.com/es/Terms/s/arcilla.aspx.
2
Op. Cit. MANRIQUE, Jesus. p. 36
3
SENER. GLOSARIO DE TERMINOS PETROLEROS. Secretaria de Energia de Mexico[en linea],
2015 [revisado 19 de noviembre 2018]. Disponible en Internet:
20
COMPRESIBILIDAD: es una propiedad que tienen todos los cuerpos, la cual
genera que estos disminuyan su volumen al someterlos a distintas condiciones de
presión o compresión manteniendo constante otros parámetros.
CONCORDANTE: serie de estratos que representan un período de sedimentación
continúo.
CONTINUIDAD DE LAS PROPIEDADES DE LA ROCA: en todo caso de
simulación es necesario tener en cuenta en relación con la permeabilidad y la
continuidad vertical.
CORRELACIÓN: es toda relación generada a partir de ensayos y por métodos
empíricos que permiten identificar y generar variables para calcular sus valores.
CUENCA SEDIMENTARIA: una depresión de la corteza terrestre, formada por la
actividad tectónica de las placas, en la que se acumulan sedimentos. La persistencia
de la depositación puede producir un grado adicional de depresión o subsidencia.
Las cuencas sedimentarias, o simplemente cuencas, pueden tener forma de cubeta
o de fosa alargada. Si las rocas ricas generadoras de hidrocarburos se combinan
con condiciones adecuadas de profundidad y duración de sepultamiento, pueden
generarse hidrocarburos en la cuenca.4
DAÑO (SKIN): reducción de la capacidad que tiene naturalmente un yacimiento
para producir los fluidos presentes en él, puede generarse por una disminución de
la porosidad o de la permeabilidad, o de ambos. Puede presentarse cerca de la cara
del pozo o ser profundo en la roca.
DISCORDANCIA: una superficie geológica que separa las rocas más antiguas de
las rocas más modernas, es cualquier periodo de tiempo en el cual no se generó
una depositación de sedimentos.5
EMPUJE POR AGUA: los yacimientos con empuje de agua, son aquellos que
tienen una conexión hidráulica entre el yacimiento y una roca porosa saturada con
agua, el cual se denomina acuífero.
ESPESOR NETO: resultado de restar al espesor total de todas aquellas porciones
de la roca que no tienen posibilidades de producir hidrocarburos
https://www.gob.mx/cms/uploads/attachment/file/8317/GLOSARIO_DE_TERMINOS_PETROLERO
S_2015.pdf.
4
SCHLUMBERGER. Cuenca Sedimentaria.Oil Field Glossary en linea], 2018 [revisado 19 de
noviembre 2018]. Disponible en Internet:
https://www.glossary.oilfield.slb.com/es/Terms/s/sedimentary_basin.aspx
5
Ibid. https://www.glossary.oilfield.slb.com/es/Terms/u/unconformity.aspx
21
ESPESOR: es la medida de una capa o de un estrato sedimentario en forma
perpendicular a su extensión lateral, presumiendo la existencia de depositación
sobre una superficie horizontal.6
ESTRATIGRAFIA: es la rama de la geología que trata del estudio e interpretación
de las rocas sedimentarias, metamórficas y volcánicas estratificadas, y de la
identificación, descripción, secuencia, tanto vertical como horizontal, cartografía y
correlación de las unidades estratificadas de rocas.7
EXPLORACION: actividad o conjunto de actividades que se valen de métodos
directos, incluyendo la perforación de pozos, encaminadas a la identificación,
descubrimiento y evaluación de Hidrocarburos en el Subsuelo, en un área definida. 8
FACTOR DE RECOBRO: es la cantidad recuperable de petróleo existente en un
yacimiento, normalmente es expresado como un porcentaje.
FALLA NORMAL: es una falla con desplazamiento vertical en el cual el bloque de
techo se desplaza hacia abajo en relación con el bloque de muro 9.
FALLA: una falla es una grieta en la corteza terrestre. Generalmente, las fallas
están asociadas con, o forman, los límites entre las placas tectónicas de la Tierra. 10
GEM CMG: según la compañía CMG (Computer Modelling Group) el software
preprocesador GEM es un simulador de yacimientos que se rige mediante la
ecuación de estado y es líder en el mundo para el modelado de yacimientos
composicionales, químicos y no convencionales.
GEOLOGÍA DEL PETRÓLEO: es una aplicación especializada de la geología que
estudia todos los aspectos relacionados con la formación de yacimientos
petrolíferos y su prospección.
GRAVEDAD API: es la clasificación que se le da a los hidrocarburos, según su
densidad relativa y su calidad. Esta clasificación fue creada por el American
Petroleum Institute.
HIDROCARBURO: es un compuesto orgánico el cual está constituido por la mezcla
de carbonos e hidrógenos, contando también, con distintas sustancias que los
acompañan o que se derivan de estos.
6
Ibid. https://www.glossary.oilfield.slb.com/es/Terms/b/bed_thickness.aspx?p=1
7
Ibid. https://www.glossary.oilfield.slb.com/es/Terms/s/stratigraphy.aspx
8
RUSSELL, RANDY. Exploracion Petrolera. Ventanas al universo [en linea], 19 de Enero de 2010
[revisado 19 de noviembre 2018]. Disponible en Internet:
https://www.windows2universe.org/earth/geology/fault.html&lang=sp
9
Op. Cit. SENER. Disponible en:
https://www.gob.mx/cms/uploads/attachment/file/8317/GLOSARIO_DE_TERMINOS_PETROLERO
S_2015.pdf
10
Ibid. https://www.windows2universe.org/earth/geology/fault.html&lang=sp
22
HUMECTABILIDAD: es una propiedad importante en un yacimiento, debido a que
afecta el comportamiento capilar y el del desplazamiento de los fluidos en la roca.
En otras palabras, la humectabilidad es la tendencia que tiene un fluido miscible en
presencia de otro fluido miscible a adherirse a la roca.
INYECCIÓN ALTERNADA DE AGUA Y GAS (WAG): es un proceso de
recuperación mejorada de petróleo, el cual consiste en inyectar agua y gas durante
distintos periodos de tiempo para lograr aumentar el factor de recobro de un
reservorio.
KERÓGENO: es el constituyente orgánico de las rocas sedimentarias, no soluble
en solventes orgánicos comunes, es sólido y resulta de la alteración diagénetica de
la materia orgánica. Suele dividirse en Kerógeno tipo l, Kerógeno tipo ll y Kerógeno
tipo lll, los cuales se caracterizan dependiendo de su contenido de hidrocarburos y
oxígeno. Siendo el tipo l el de más alto contenido de hidrocarburos y más bajo
contenido de oxígeno.
LIMOLITA: roca sedimentaria, su textura es de grano muy fino de diversos
minerales. El tamaño del grano esta entre 0.0086-0.0036 mm. Se compone
principalmente de óxido de hierro, cuarzo, calcita entre otros. Este tipo de rocas son
el producto de la compactación de sedimentos compuestos por partículas, que por
su tamaño, se encuentra en un rango intermedio entre las arcillas y las arenas. 11
LITOLOGÍA: es toda la naturaleza macroscópica del contenido mineral, la
granulometría, la textura y el color de las rocas
LUTITA: es una roca sedimentaria de origen detrítico compuesta por partículas del
tamaño de la arcilla tienden a ser plásticas y el tamaño de grano es menor a 1/256
mm .Su compone de feldespatos y mica principalmente.12
MALLA DE SIMULACIÓN: es una aproximación numérica que consiste en dividir
en un número finito de unidades de volumen discretas (celdas) un yacimiento, lo
cual permite modelar las estructuras y propiedades presentes en el mismo.
MARCO GEOLOGÍCO: es un resumen de la geología presente en un lugar en
específico en el cual se incluye la estratigrafía, las rocas presentes en el área de
estudio. Se puede hacer de manera regional o puntual dependiendo de la escala en
la que se esté trabajando.
MIGRACIÓN: término utilizado en exploración para denominar al movimiento que
sufre el aceite, del lugar donde se formó hasta la roca almacenadora o trampa. 13
11
Op. Cit. MANRIQUE, Jesus p. 36
12
Op. Cit. MANRIQUE, Jesus. p. 36
13
RACEFN. Glosario de Geología. Real Academia de Ciencias exactas [en linea], 26 de enero
2010 [revisado 19 de noviembre 2018]. Disponible en Internet:
https://www.ugr.es/~agcasco/personal/rac_geologia/rac.htm
23
MISCIBILIDAD: propiedad que tienen los fluidos, trata de la habilidad que tienen
dos o más fluidos para mezclarse entre sí.
MOVILIDAD: es la facilidad con la cual un fluido se mueve en el yacimiento. Se
calcula como la relación entre la permeabilidad efectiva de la roca a un fluido y la
viscosidad de éste.
POZO: perforación para el proceso de búsqueda o producción de petróleo crudo
gas natural o para proporcionar servicios relacionados con los mismos. Los pozos
se clasifican de acuerdo a su objetivo y resultado como: pozos de aceite y gas
asociado, pozos de gas seco y pozos inyectores. 14
PROYECTO PILOTO: es aquél proyecto que se lleva a cabo en un pequeño sector
representativo de un yacimiento, en donde se efectúan pruebas similares a las que
se llevarían a cabo en toda el área del yacimiento. El objetivo es recabar información
y/o obtener resultados que puedan ser utilizados como base de estudios
convencionales o de simulación matemática de todo el yacimiento.15
REFÍNERIA: planta industrial que tiene como objetivo la refinación y separación de
distintos componentes y derivados de los hidrocarburos.
SISTEMA PETROLÍFERO: son todos componentes geológicos y los procesos
necesarios para generar, migrar y almacenar hidrocarburos; esto incluye una roca
generadora madura, un trayecto de migración, una roca yacimiento, una trampa y
un sello.16
SUPRAYACENTE: se extiende por encima de una superficie geológica.
TRAMPA ESTRATIGRAFICA: una variedad de contenedor geológico sellado con
capacidad para retener hidrocarburos, formado por los cambios producidos en el
tipo de roca o por acuñamientos, discordancias, o rasgos sedimentarios, tales como
los arrecifes.17
14
Ibid. https://www.ugr.es/~agcasco/personal/rac_geologia/rac.htm
15
Ibid. https://www.gob.mx/cms/uploads/attachment/file/8317/GLOSARIO_DE_TERMINOS
PETROLEROS _2015.pdf
16
Ibid. https://www.glossary.oilfield.slb.com/es/Terms/p/petroleum_system.aspx
17
SCHLUMBERGER.Trampa estratigrafica.oil Field Glossary en linea], 2018 [revisado 19 de
noviembre 2018]. Disponible en Internet:
https://www.glossary.oilfield.slb.com/es/Terms/s/stratigraphic_trap.aspx.
24
RESUMEN
Palabras clave:
Inyección WAG
Recobro mejorado
Flue Gas
Refinería de Barrancabermeja
Valle Medio del Magdalena
25
INTRODUCCIÓN
26
ingresos que se tendrán durante este periodo de tiempo, los cuales se generaran
mediante la predicción establecida y los costos aproximados de una implementación
de un proceso de inyección de gas y agua.
27
OBJETIVOS
OBJETIVO GENERAL
Evaluar técnica y financieramente de la inyección tipo WAG con Flue gas
proveniente de la refinería de Barrancabermeja en las arenas B de la Formación
Mugrosa en un campo mediante una simulación numérica.
OBJETIVOS ESPECÍFICOS
5. Evaluar la viabilidad de una inyección tipo WAG mediante el uso del simulador
GEM para fluidos composicionales de la empresa CMG.
6. Evaluar la viabilidad financiera del uso del uso de la Inyección tipo WAG con
Flue Gas proveniente de la refinería de Barrancabermeja en las arenas B de la
formación mugrosa en un campo, mediante el uso del Valor Presente Neto
(VPN).
28
1. GENERALIDADES GEOLÓGICAS DE LA CUENCA Y ÁREA DEL
PROYECTO
La Cuenca sedimentaria del Valle Medio del Magdalena es una de las más
importantes en cuanto a producción de hidrocarburos en Colombia, esta involucra
el área en la cual se localiza el Campo objeto de estudio de este proyecto, cercana
a la refinería de Barrancabermeja.
En este capítulo se presentan aspectos asociados a la localización geográfica,
marco geológico e historia de producción del área énfasis del proyecto.
18
ECOPETROL S.A. Contrato de la Concesión de Mares. Ecopetrol. [en linea], 23 de septiembre de
2014 [revisado 10 agosto de 2018]. Disponible en Internet:
https://www.ecopetrol.com.co/wps/portal/es/ecopetrol-web/nuestra-empresa/sala-de-
prensa/publicaciones/otras-publicaciones/cronica-de-la-concesion-de-mares/contrato-de-la-
concesion-de-mares
29
En el año 1922 inicia la historia de la Refinería de Barrancabermeja, construida en
un punto geográfico muy importante para la industria petrolera durante sus primeros
años de existencia, esta, inicia sus operaciones con equipos diseñados
especialmente para destilar petróleo mediante el uso de leña.19
En 1928 entra en operación la Refinería de Barrancabermeja, construida
inicialmente por la compañía Tropical Oil, después de haberse construido un
oleoducto entre el Campo Infantas y la ciudad de Barrancabermeja.
La Cuenca del Valle Medio del Magdalena ha tenido una amplia historia exploratoria
y de producción que ha dado lugar a la operación de 64 campos petroleros, los
cuales han generado producción de hidrocarburos en las formaciones Tablazo, La
Luna, Lisama, La Paz, Esmeraldas, Mugrosa y Colorado.
Al mes de octubre de 2018, la Cuenca del Valle Medio del Magdalena cuenta con
64 campos en operación a lo largo de la Cuenca del Valle Medio del Magdalena se
encuentran 21 son campos productores de gas, 10 producen crudo liviano, 22
crudos medianos y 31 se encargan de producir crudos pesados.
En consideración al sector de énfasis de este proyecto, la Tabla 1 relaciona los
principales campos descubiertos cercanos a la Refinería de Barrancabermeja.
CAMPO FORMACIÓN PRODUCTORA AÑO DE PERFORACIÓN PROFUNDIDAD ALCANZADA POZO DESCUBRIDOR POZOS PERFORADOS
La Cira Mugrosa y Colorado 1926 1926 Ft La Cira-58 1367
Infantas Colorado 1918 2285 Ft Infantas-1 876
Gala Mugrosa y Colorado 1985 5600 Ft Cardenales-1, Yuma-1 423
Galan Mugrosa y Colorado 1945 5885 Ft Galan-1 286
Llanito Mugrosa y Colorado 1955 5561 Ft Llanito-1 223
Casabe Mugrosa, Colorado y La Paz 1941 6500 Ft CSBE-1 456
1.2 LOCALIZACIÓN
La Cuenca del Valle Medio, está ubicada en Colombia, con un área aproximada de
24.000 km2, cubre parte de los departamentos de Santander, Cundinamarca,
Boyacá, Antioquia, Bolívar, Cesar, Norte de Santander y Casanare.
19
VECINO, Erminda. Refinería de Barrancabermeja celebra 93 años de refinar crudo para el país.
Ecopetrol. [en linea], 18 de febrero de 2015 [revisado 10 agosto de 2018]. Disponible en Internet:
https://www.ecopetrol.com.co/wps/portal/es/ecopetrol-web/nuestra-empresa/sala-de-
prensa/boletines-de-prensa/Boletines%202015/Boletines%202015/refineria-de-barrancabermeja
30
El área de estudio de este proyecto, la cual no se detalla por confidencialidad, queda
en cercanía de la Refinería de Barrancabermeja, hacia la zona occidental de la
cuenca sedimentaria del Valle Medio del Magdalena, al sur del rio Sogamoso y en
la parte oriental del Rio Magdalena.
El acceso a dicha área puede darse por vía área, tomando un avión desde Bogotá
hasta Barrancabermeja, que en tiempo requiere 30 minutos, y de allí desplazarse
por vía terrestre al sector de estudio aproximadamente 14 km. Por vía terrestre
partiendo desde Bogotá D.C, se debe realizar el siguiente recorrido: salir hacia el
Noroccidente de la ciudad por la Calle 80, en dirección hacia el Municipio de la Vega,
por la carretera nacional #45 por un trayecto de alrededor de 120 kilómetros hacia
el Noroccidente hasta llegar al municipio de Villeta. En dirección al Norte recorridos
unos 85 km, después de sobrepasar el municipio La concepción, aproximarse al
municipio de Guaduas. Al sobrepasar este municipio por la Vía principal, tomar la
carretera nacional #56 hacia el Nororiente. Al recorrer unos 36 km, y al llegar al
municipio de Puerto Vargas. Se debe tomar la Ruta del Sol, tomando dirección al
Nororiente, recorrer alrededor de 100 km. En el municipio de Puerto Libre, se debe
tomar salida hacia el Oriente en dirección hacia Puerto Boyacá, sobrepasado este
municipio recorrer unos 130 km hasta encontrar un desvió 8 km antes del municipio
de Cimitarra en el cual se debe girar hacia el Norte en dirección a la Vía Cimitarra –
Puerto Araujo. Recorrer alrededor de 80 km y tomar la Troncal Barrancabermeja –
Magdalena hacia el Norte. Pasados unos 30 km dirección al Norte, se habrá llegado
al destino de interés en la ciudad de Barrancabermeja.
En total el camino por vía terrestre desde Bogotá D.C. y Barrancabermeja es de
unos 476 km con una duración aproximada de entre NUEVE (9) y DIÉZ (10) horas
de recorrido.
En la Figura 1 se puede ver una aproximación geográfica de la zona de estudio,
también se pueden evidenciar las rutas de interés cercanas al área del proyecto,
junto con un punto de referencia tomado en Barrancabermeja, ya que allí se
encuentra la refinería de la cual se dispondrá el Flue Gas para la inyección WAG.
31
Figura 1. Localización del área de interés.
Fuente: Mapa físico político de Colombia. [En línea]. [Citado 3 de agosto de 2018]: Disponible en
Instituto Agustín Codazzi. Disponible en: http://ssiglwps.igac.gov.co/mapas_de_colombia/IGAC.
Modificado por el autor.
32
1.3 MARCO GEOLÓGICO
En el año de 1983 David Kingston denominó a la Cuenca del Valle Medio del
Magdalena como una cuenca poli-histórica20, debido a que su formación se dio a
partir de distintas etapas en distintas épocas de la prehistoria.
Su formación tuvo lugar a través de distintas etapas que tuvieron una relación
directa con distintos eventos tectónicos de la región noroeste de Suramérica, las
cuales ocurrieron durante el periodo Triásico Tardío, el Cretácico Medio, el
Paleógeno Temprano y el Neógeno Medio.21
Esta cuenca se extiende a través del Rio Magdalena y está delimitada al norte por
el sistema de fallas Espíritu Santo y al sur por el Cinturón de Girardot. Al noreste
está limitada por las fallas de Bucaramanga-Santa Marta y al sureste por el sistema
de fallas Bituima y La salina. El límite occidental está marcado por la superposición
más occidental del relleno de la cuenca del Neógeno en la Serranía de San Lucas
y el sótano de la Cordillera Central.22
El Valle Medio del Magdalena muestra una sucesión de depósitos continentales del
Jurásico superpuestos por sedimentos cretácicos, tanto calcáreos como
siliclásticos, de origen transicional a marino. La secuencia del Paleógeno está
formada por rocas siliclásticas depositadas principalmente bajo condiciones
continentales con alguna influencia marina.23
La Cuenca del Valle Medio del Magdalena es una de las cuencas más exploradas
en el país, en ella se han descubierto más de 60 campos en distintos depósitos de
hidrocarburos formados durante el Paleógeno.
Las áreas de exploración potencial se relacionan principalmente con fallas normales
invertidas y anticlinales de subdesarrollo. Las sutiles trampas estratigráficas
asociadas con las incrustaciones del Mioceno-Eoceno superior, los canales incisos
y los truncamientos también son objetivos principales para futuras exploraciones. 24
En cuanto al área de énfasis de este proyecto, la Formación Mugrosa remonta su
depositación hacia el periodo terciario, más exactamente hacia inicios del
Oligoceno. Esta se encuentra definida por areniscas de grano fino con altas
intercalaciones de Lutitas y Limolitas.
20
KINGSTON, D., DISHROON, C., WILLIAMS, P. Global basin classification system. AAPG
Bulletin. V.67. p. 2175
21
ANH. Op.cit., p. 45
22
Ibid., p. 45
23
Ibid., p. 80
24
Ibid., p. 81
33
El nombre de la Formación Mugrosa fue dado por A. K. Gill y referenciado por
Morales y otros en 1958, refiriéndose a la parte inferior del grupo Chuspas. Su
nombre fue tomado de la quebrada Mugrosa.
La Formación Mugrosa es una de las unidades litológicas más importantes de la
Cuenca del Valle Medio del Magdalena, además de ser reconocida como uno de los
reservorios más importantes en Colombia.
1.3.1 Columna estratigráfica. En la Figura 2 se muestra la columna estratigráfica
generalizada para la cuenca del Valle Medio del Magdalena, resaltando la
Formación Mugrosa, énfasis de estudio en este proyecto. La sección perforada por
los pozos del campo objeto de estudio la conforman rocas del Cretácico superior al
reciente, no obstante, la Formación Lisama no fue perforada en dicho campo debido
a una discordancia de esta formación presente en el campo objeto de estudio de
este trabajo de investigación.
34
Figura 2. Columna estratigráfica generalizada para la cuenca del Valle Medio del
Magdalena.
35
1.3.2 Estratigrafía del área de estudio y geología del petróleo. A continuación,
describen las principales características de las formaciones estratigráficas
involucradas en el sistema petrolífero del área de estudio, en orden de la más
antigua a la más reciente, con base en información de JULIVERT & DAPORTA
principalmente.
Se realiza un énfasis en la Formación Mugrosa, objetivo de este proyecto.
Adicionalmente, se relacionan los principales parámetros asociados con su sistema
petrolífero.
1.3.2.1 Formación Tablazo. Su edad según M. Julivert25 esta formación fue
depositada en el Cretácico Inferior, en el Albiano. Está compuesta principalmente
por shale negra compacta, y gran cantidad de calizas duras con un tono azul. Su
espesor varía entre los 400 y 450 metros. Tiene un ambiente de depositación
marino. La Formación Tablazo se encuentra en contacto concordante con la
Formación Paja que la infrayace.
La Fm. Tablazo tal como se muestra en la Figura 2 presenta características de
generadora potencial para la cuenca y el campo, ya que en esta formación se
presentan valores de TOC entre el 2% y 4%.
1.3.2.2 Formación Simití. Según M, Julivert26, esta formación tuvo su depositación
en el Cretáceo Inferior entre el Cenomaniano y el Turoniano. Cuenta con shales con
colores grisáceos y negros. Se presentan capas delgadas de conglomerados en el
tope de la formación, hacia el Este de la misma se presenta una predominancia de
arenas. Su espesor varía entre 250 m y 650m. El ambiente de depositación de esta
formación es marino de somero a profundo. La Formación Simití infrayace
concordantemente con la Formación La Luna.
Junto con las Formaciones La Luna y Tablazo, esta formación se destaca en el
sistema petrolífero de la Cuenca del Valle Medio del Magdalena al tener parte de
las rocas potencialmente generadoras de hidrocarburos, aunque cuenta con
distintos tipos de calizas y shales potenciales, estas no se consideran como unas
de las principales generadoras de la cuenca, ya que cuentan con valores de TOC
menores al 2%.
1.3.2.3 Formación La Luna. Según M. Julivert27, la Formación La Luna remonta su
edad al Cretácico superior, entre el Turoniano y el Coniaciano. Esta unidad se
encuentra constituida por calizas de tono gris oscuro, con pequeñas intercalaciones
de arcillas, lutitas de color gris en capas delgadas con delgadas intercalaciones de
calizas arcillosas. Se caracteriza por tener una estructura subangular en la cual se
25
JULIVERT. Manuel. Léxico estratigráfico de Colombia – Precámbrico, Paleozoico, Mesozoico.
Vol. 5, Paris. Centre National de la Reserche Scientifique. 1968. p. 260-269.
26
Ibid. p. 442-443.
27
Ibid. p. 307-309.
36
presentan grandes intercalaciones de arcillas, las cuales tienen tonalidades de color
gris verdoso, rojizo y violeta. Esta formación cuenta con un espesor aproximado de
2650 pies aproximadamente a lo largo de la zona de estudio, aunque se presentan
varias discontinuidades estratigráficas a lo largo del área de estudio. El ambiente de
depositación es marino de aguas relativamente poco profundas. La Formación La
Luna infrayace concordantemente con la Formación Umir.
La Formación La Luna juega un papel en el sistema petrolífero de la Cuenca del
Valle Medio del Magdalena, ya que es en la que se cuenta con un mayor potencial
de generación de hidrocarburos. El TOC de esta formación varía entre 2% y 6%,
Kerógeno tipo ll.
1.3.2.4 Formación Umir. Su edad es del Cretácico Superior, entre el Turoniano y
el Coniaciano. En esta formación se encuentran distintos tipos de shales, aunque
en el tope de la misma se presentan intercalaciones de areniscas de grano fino a
medio, el espesor promedio de esta formación es de unos 1000 metros. Su
ambiente de depositación es marino somero. La Formación Umir infrayace
concordantemente a la Formación Lisama hacia el Este, aunque en la zona de
estudio de este trabajo, se presenta una discordancia de la Formación Lisama, por
lo que en las secciones perforadas en el campo objeto de estudio la Formación Umir
infrayace discordantemente a la Formación Esmeraldas-La Paz.
La Formación Umir tiene un papel muy importante dentro del sistema petrolífero de
la Cuenca del Valle Medio del Magdalena, ya que sus shales, al ser rocas de grano
muy fino tienen una porosidad y permeabilidad muy bajas, por lo que no permiten el
paso del hidrocarburo, actuando como sello de las Formaciones Esmeraldas-La
Paz, Mugrosa y Colorado.
1.3.2.5 Formación Esmeraldas-La Paz. Según J. De porta28, fue depositada en el
Eoceno Inferior. En esta formación predominan sedimentos limo-arcillosos con
coloraciones grises y grises claras, en su tope cuenta con delgadas intercalaciones
de areniscas blancas, de grano medio a fino con distintas inclusiones de materia
orgánica. Su espesor promedio aproximado es de 500 pies y su ambiente de
depositación es continental-fluvial. Infrayace concordantemente a la Formación
Mugrosa.
Se consideran distintas acumulaciones de arenas presentes en esta formación
como rocas potenciales almacenadoras de hidrocarburos. Sus arenas cuentan con
porosidades entre 8% y 19% y con permeabilidades entre 20md y 450 md.
1.3.2.6 Formación Mugrosa. La Formación Mugrosa, según J. De Porta29, data su
edad de depositación en el Eoceno Superior y el Oligoceno inferior, su espesor varía
28
DE PORTA, Jaime. Léxico estratigráfico de Colombia – Terciario a Cuaternario. Vol. 5, Paris.
Centre National de la Reserche Scientifique. 1974. p. 327-330.
29
Ibid. p. 373-376.
37
entre 500 y 800 m y se formó entre hace unos 30 y 25 millones de años en un
ambiente fluvial en su mayoría, aunque conto con diferentes tipos de facies, entre
las cuales se destacan facies de lago, canales y llanuras.
La Formación Mugrosa pertenece a la base del grupo sedimentológico denominado
¨Chuspas¨. Esta es una de las unidades estratigráficas de distribución más amplia
y constante en la Cuenca del Valle Medio del Magdalena. La Formación Mugrosa
infrayace concordantemente a la Formación Colorado.
Es de gran importancia para el sistema petrolífero de la Cuenca del Valle Medio del
Magdalena debido a que se considera una de las rocas con mayor potencial de
almacenamiento de hidrocarburos, sus arenas cuentan con porosidades promedio
entre 9% y 17% y su permeabilidad varía entre 20 md y 600 md.
En consideración a los objetivos de este trabajo, el cual se basa en las Arenas B de
la Formación Mugrosa, a continuación se hace una breve descripción de esta
sección de dicha formación.
Arenas B Formación Mugrosa. La Zona B de la Formación Mugrosa se
caracteriza principalmente por contar con una secuencia limo-arcillosa con una
coloración amarillenta tornando a gris. Presenta intercalaciones de areniscas con
un grano fino.
Contiene fracciones de areniscas de grano fino, estas tienden a aumentar hacia
la base de la zona, justo antes de llegar a tope de la Zona C de la Formación
Mugrosa. El espesor promedio de la zona es de unos 1400 pies.
La Zona B de la Formación Mugrosa, según la empresa ECOPETROL S.A. tiende
a subdividirse en tres secciones de areniscas principales debido a la alta intrusión
de arcillas que se presenta en esta región, las características de estas zonas se
presentan en orden estratigráfico a continuación.
Zona B3. Cuenta con un espesor aproximado de 350 pies, cuenta con una
continuidad lateral de areniscas y depósitos de llanura. En su base se presenta un
poco continuidad.30
38
compone de arcillas de color gris claro, negro, carbonáceas con rasgos rojos
intercaladas con areniscas de grano medio poco potentes. Su espesor varía
alrededor de los 2500 pies. Cuenta con un ambiente de depositación fluvial. Esta
formación infrayace concordantemente al Grupo Real.
La Formación Colorado es considerada como roca almacenadora probada ya que
distintos campos de la zona como Campo Gala, Galán, La Cira, entre otros, cuentan
con producción de hidrocarburos en esta formación. Cuenta con una porosidad
promedio entre 10% y 15% y una permeabilidad entre 150 md y 550 md.
1.3.2.8 Grupo Real. Según J. De Porta33, el Grupo Real fue depositado a lo largo
del Mioceno Medio y Superior. El Grupo Real se subdivide en 5 (CINCO)
formaciones, las cuales se establecen de la siguiente forma: Formación Lluvia, la
cual se compone de conglomerados de cuarzo, fragmentos de carbón y mineral
arcilloso de hierro; Chontorales Shale, que se compone de areniscas
conglomeráticas con estratificación cruzada intercalada con shales con rasgos
grises y azules; la Formación Hiel, se compone de shales con rasgos grises y rojos
intercalados con areniscas en capas delgadas; Enrejado Shale, se compone de
troncos carbonizados; y Formación Bagre, que se compone de shales con rasgos
rojos, grises, morados y pardos, intercalados con areniscas en capas delgadas. La
sucesión termina con areniscas conglomeráticas con estratificación cruzada,
madera carbonizada y hojas. Su espesor varía entre 500m y 3900m. Su ambiente
de depositación es continental-fluvial. Le infrayace discordantemente la Formación
Mesa.
1.3.2.9 Formación Mesa. Según J. De Porta34, la Formación Mesa fue depositada
en el Cuaternario. Está compuesta por una serie de conglomerados, granodiorita,
cuarzo, chert, andesita y distintas rocas metamórficas con presencia de lentes de
arenas gruesas. Su espesor promedio es de unos 250m. Tiene un ambiente de
depositación continental. Le infrayace discordantemente el Grupo Real.
1.3.2.10 Migración. Durante la época del Eoceno se presenta una discordancia que
fue sumamente importante, ya que, a raíz de ella, se formó un sistema ideal, en el
cual se facilitó la migración de los hidrocarburos. En la Cuenca del Valle Medio del
Magdalena presento varias rutas de migración en las cuales se destacan las
siguientes:
Migración vertical directa de hidrocarburos, ya que la Formación La Luna atraviesa
la inconformidad del Eoceno.
Migración lateral a lo largo de las areniscas del Eoceno.
33
Ibid. p. 142-144
34
Ibid. p. 356-362
39
Migración vertical a través de las fallas en áreas donde la Formación La Luna no
atraviese la inconformidad del Eoceno.35
35
BARRERO. Darío. Colombian Sedmentary Basins. Edición 1. Bogotá, Colombia. ANH Agencia
Nacional de Hidrocarburos [En línea] 2007. [Consultado 19 de noviembre de 2018]. Disponible en:
http://www.anh.gov.co/Informacion-Geologica-y-Geofisica/Cuencas-
sedimentarias/Documents/colombian_sedimentary_basins.pdf
36
Análisis de yacimientos sometidos a procesos de inyección alternada de gas y agua (WAG)
mediante simulación numérica. Universidad Industrial de Santander. 2009, p. 134.
40
gas en solución y que posterior a este se produce un empuje leve de algún acuífero
cercano al reservorio.
En el año 2007 se iniciaron 3 (TRES) pilotos de inyección de agua por patrones
irregulares directos en distintos sectores, el cual aumento considerablemente la
producción de hidrocarburos en el campo.
Actualmente algunos pozos del campo presentan producción por medio de distintos
métodos de levantamiento artificial, como levantamiento por cavidades progresivas.
1.4.2 Tiempo de producción. El campo objeto de estudio remonta sus inicios en
producción de hidrocarburos hacia el año de 1960, es decir tiene una producción
acumulada de 58 años.
1.4.3 Número de pozos. El campo tiene un área aproximada de unos 70 Km 2 de
extensión, en los cuales se han perforado alrededor de 223 pozos. En la Tabla 2 se
observa el estado de los pozos del campo objeto de estudio.
NÚMERO
ESTADO
DE POZOS
Productores 75
Abandonados parcialmente 31
Abandonados 24
Inactivos 93
TOTAL 223
41
Gráfica 1. Historial de producción acumulada
ZONA
PARAMETRO
A B C D
Profundidad (Ft) 4000 5200 6200 7000
Espesor (Ft) 26 50 40 18
Porosidad promedio (%) 24 21 20 18
Permeabilidad promedio (%) 480 340 240 160
Swi (%) 42 51 46 47
Tipo de producción Gas Solución Empuje Agua + Gas Solución
Tipo de crudo Naftenico
42
El énfasis de estudio del presente trabajo de investigación se centra en las arenas
de la zona B de la Formación Mugrosa de la Cuenca del Valle Medio del Magdalena.
Esta sección de la formación se subdivide principalmente en tres unidades
principales, las cuales se denominan como B1, B2 y B3.
La Tabla 4 representa las propiedades petrofísicas de cada una de las unidades de
arenas presentes en la Zona B de la Formación Mugrosa.
Tabla 4. Propiedades petrofísicas promedio del área de estudio.
ESPESOR ESPESOR POROSIDAD PERMEABILIDAD
FORMACIÓN ZONA UNIDAD
BRUTO (Ft) NETO (Ft) PROMEDIO (%) (Md)
Mugrosa B 1 100 16.68 21.68 94.26
Mugrosa B 2 300 24.25 21.92 90.17
Mugrosa B 3 400 64.33 21.47 96.16
ZONA
PARAMETRO
A B C D
Gravedad API 18.4 20 22 24
Viscosidad (cP) 50 19.4 12.1 5
Factor volumetrico 1.058 1.094 1.105 1.2
Factor volumetrico @Pb 1.061 1.094 1.09 1.206
GOR (SCF/STB) 150 500 400 600
Punto de burbuja (psi) 1600 2175 2240 2900
43
2. MARCO TÉORICO
37
PARIS DE FERRER, Magdalena. Inyección de agua y gas en yacimientos petrolíferos. Segunda
edición. 2007, p.327
38
Larsen, J.A. and Skauge, A. Simulation of the Immiscible WAG Process Using Cycle-Dependent
Three-Phase Relative Permeabilities, citado por MARCHESIN, Dan. Wave Structure in WAG
Recovery. En: Society of petroleum engineering. Junio, 2011.
44
gas (alta eficiencia microscópica) y también se logra complementar las desventajas
con las que cuenta cada uno de los procesos de inyección. Y esta complementación
entre la inyección de ambos fluidos se ve reflejada en el aumento del factor de
recobro.39
El uso del agua para la inyección, se justifica de manera que está, ayuda a controlar
la movilidad y logra estabilizar el frente de invasión. Además, que la gran mayoría
de los procesos de inyección WAG parten a raíz de un proceso de inyección de
agua que se lleva por un largo tiempo atrás en los distintos campos de estudio.
Se establece una inyección por baches, debido a que de esta forma se puede
controlar más fácilmente la cantidad de fluido que se necesita inyectar, además de
tener un mayor control sobre el desplazamiento del frente de invasión para cada
fluido.
Figura 3. Proceso de inyección alternada de gas y agua
39
MENEZ DIAZ,Robinson; MUÑOZ NAVARRO,Samuel Fernando y Rivera De La Ossa, Juan
Eduardo. Inyección Alternada De Agua Y Gas (Wag): Generalidades Teóricas, Casos
Históricos,ProblemasOperacionales Comunes Y ScreeningEstadístico Para Su Aplicación.
Universidad Industrial de Santander, 2007. p. 1-55
45
inyectada, que es aquel que controla principalmente la movilidad de los fluidos
dentro del yacimiento.
40
MONROY ZAPATA, Diana; CORONADO NAVARRO, Jorge Luis. Análisis de yacimientos
sometidos a procesos de inyección alternada de gas y agua (WAG) mediante simulación numérica.
Universidad Industrial de Santander. 2009, p.4.
41
L.M. Surguchev ; RAGNHILD Korbol ; SIGURD Haugen ; O.S. Krakstad. Screening of WAG
Injection Strategies for Heterogeneous Reservoirs. SPE-25075-MS. 1992.
46
En la Gráfica 2 se establece una proporción de los distintos proyectos de
recuperación de petróleo mediante la inyección alternada de gas y agua.
42
Ibid. L.M. Surguchev ; RAGNHILD Korbol ; SIGURD Haugen ; O.S. Krakstad.
43
D. H. Tehrani, A. Danesh, M. Sohrabi And G. Henderson. Enhanced Oil Recovery Y Water
Alternating Gas (Wag) Injection. Department Of Petroleum Ingineering, Heriot-Watt University
Edinburgh, Uk. Paper 01,04 y 06. 1999-2001
47
2.3 FACTORES QUE AFECTAN LA INYECCIÓN WAG
Los procesos de inyección alternada de gas y agua deben contar con un largo
estudio técnico, el cual logre identificar todos los posibles factores que puedan
afectar o mejorar el proceso de recuperación dentro de un yacimiento.
Este tipo de proceso se ve afectado por una gran cantidad de variables, siendo las
condiciones geológicas del reservorio y las propiedades de los fluidos las más
importantes. Dentro de todos estos factores podemos encontrar aquellos que se ven
afectados por las propiedades geológicas y físicas del yacimiento, y los que se
presentan por las características de los fluidos.
2.3.1 Características geológicas y físicas del yacimiento. El uso de la inyección
tipo Water Alternating Gas tiene una gran dependencia de las características
geológicas que pueda tener un yacimiento. A continuación se establecen las
características que puedan beneficiar o complicar el uso de este sistema de recobro
mejorado en un reservorio petrolero.
2.3.1.1 Heterogeneidad del yacimiento. La estructura y los estratos de rocas y
minerales que se puedan presentar en el yacimiento son el factor más importante
que se deba tener en cuenta. Debido a que cada una de las capas o estratos con
los que se cuenta en el reservorio tiene propiedades distintas a los demás.
A lo largo de un yacimiento, aunque como ya se mencionó anteriormente sus
estratos sedimentológicos tienen distintas propiedades, se tiene la ventaja que
estos están comunicados entre sí, además que cuentan con ambientes y
condiciones de depositación muy similares.
Es indispensable determinar todas aquellas zonas favorables para la inyección de
agua y gas, ya que es necesario generar una caracterización eficiente del
yacimiento.
La principal característica para tener en cuenta debe ser la permeabilidad vertical y
horizontal del yacimiento, ya que, si se cuenta con una alta permeabilidad vertical
se pueden formar frentes de empuje muy inestables, debido a que los fluidos con
mayor densidad como el caso del petróleo y el agua, se moverán hacia las zonas
medias y bajas del yacimiento. Así como también, el gas tendera a desplazarse
hacia la zona superior del reservorio.
Un yacimiento muy heterogéneo es aquel que cuenta con grandes variaciones de
sus propiedades petrofísicas, tales como; la permeabilidad, la porosidad y la
cementación. Estas influyen importantemente en el desplazamiento que puedan
tener los fluidos a través de su paso por el reservorio.
48
2.3.1.2 Segregación gravitacional. Este factor, de segregación gravitacional,
afecta principalmente la eficiencia vertical de desplazamiento, debido a que si esta
se junta con una permeabilidad vertical alta, el gas que pueda ser inyectado va a
tender a desplazarse hacia la parte superior del yacimiento, así como el agua va
desplazarse hacia la parte inferior.
Los objetivos principales del proceso de inyección WAG son aprovechar las
propiedades macroscópicas con las que cuenta el agua, tales como, desplazar el
banco de petróleo hacia los pozos productores, estabilizar las caídas de presión
naturales del yacimiento.
Al realizar una inyección WAG, con el tiempo, el frente de desplazamiento tiende a
formar tres zonas diferentes: el gas por un lado tiende a situarse en la zona superior
de la zona de inyección, por su parte el agua tiende a bajar hacia la base y en la
mitad de estas dos fases se presenta una mezcla de agua y gas que se genera a
raíz de la inyección de estos fluidos.
49
formación, no generando eficientemente su función, y luego retirarlo de estas zonas
presenta una gran dificultad. En conclusión, se estaría desperdiciando el gas del
cual se va a proveer.
2.3.1.3 Espesor de la formación. El espesor de la formación a producir para un
proceso de inyección WAG es determinante en gran medida, ya que de esta
depende la estabilidad que pueda llegar a tener el frente de invasión, debido a que
entre más espesa sea la zona productora, las fuerzas gravitacionales van a ejercer
mayor acción e influencia sobre los fluidos inyectados.
Al tener las fuerzas gravitacionales, mayor influencia sobre los fluidos inyectados,
la estabilidad del frente de invasión se vería altamente reducida debido a que el gas
tendera a situarse en la parte superior del reservorio y el agua en la parte inferior
del mismo, como se evidencia en la Figura 5.
50
propiedades de los fluidos que se puedan ver afectadas dentro de la inyección tipo
WAG se encuentran las permeabilidades relativas de los fluidos, la viscosidad del
crudo en el yacimiento, expansión o hinchamiento del petróleo y disminución de la
saturación residual del petróleo (Sor).
2.3.2.1 Permeabilidades relativas de los fluidos. La permeabilidad relativa de
cada una de las fases de los fluidos está directamente relacionada con la saturación
con la que cuenta cada uno dentro del reservorio, esto quiere decir que la
permeabilidad relativa para cada fase depende de sí misma en cuanto a su
saturación, este fenómeno no ocurre con la permeabilidad relativa al aceite, ya que
esta depende directamente de las saturaciones de agua y gas.44
2.3.2.2 Viscosidad del crudo en el yacimiento. Debido a que durante la inyección
tipo WAG se inyectan gases de combustión al yacimiento, estos gases actúan de
forma en que este gas empieza a absorber distintos componentes livianos del crudo,
generando que este disminuya su viscosidad al ser más liviano y poderse movilizar
más fácil. En sí, el gas se disuelve dentro del petróleo.
Esto genera una ventaja ya que relaciona la razón de movilidad y aumentar también
la eficiencia del barrido. Este efecto genera más notabilidad en fluidos en crudos no
muy livianos, en los cuales las caídas de viscosidad son más fuertes.
2.3.2.3 Expansión o hinchamiento del petróleo. Como ocurre con cualquier gas
que se inyecta en un yacimiento, el gas de combustión que se disuelve produce una
expansión en el petróleo, como consecuencia de esto, el volumen poroso que ocupa
el crudo es mayor y en consecuencia, será mayor la cantidad del petróleo
desplazado.
Este efecto se puede medir en un laboratorio a través de distintos estudios de
comportamiento de las fases, basándose como indicativo en el factor volumétrico
del petróleo. Aunque también, es posible estimar esto mediante correlaciones y por
modelamiento de fluidos.
2.3.2.4 Disminución de la saturación residual de petróleo. Cuando se generan
barridos de desplazamiento por un frente miscible, la saturación residual del
petróleo tiende a tener un valor muy bajo, sin embargo es mucho menos al que se
tendría solamente con una inyección de agua. Esto se debe al hinchamiento ya
mencionado por variaciones en las características del crudo por efecto del gas que
se inyecta.
2.2.3 Relación entre las movilidades. La movilidad es la relación existente entre
la permeabilidad efectiva de un fluido y su viscosidad. Al realizarse una inyección
de fluidos distintos a los presentes en el yacimiento, y al tener estos una interacción
con los componentes in situ en el reservorio, se generara una relación de
movilidades entre estos, que se expresara como el cociente entre las movilidades
44
MONROY ZAPATA., Op. Cit., p. 36
51
de cada uno de los fluidos, como se expresa en la siguiente Ecuación 1, la cual se
refiere a la relación de movilidades entre el agua inyectada y el petróleo in situ en la
formación.
Ecuación 1. Relación de movilidad entre el petróleo y el agua.
(𝐾𝑤 𝜇𝑤 )
𝑀=
(𝐾𝑜 𝜇𝑜 )
45
Schlumberger. Oilfield Glossary. Yacimiento de doble permeabilidad, [En línea], 2018, [revisado
el 28 de agosto de 2018]. Disponible en:
https://www.glossary.oilfield.slb.com/es/Terms/f/fingering.aspx
52
La digitación viscosa como tal, puede generar que los tiempos de ruptura de los
fluidos sean mucho menores a los esperados y que el barrido del fluido no sea
eficiente.
46
ESCOBAR,Freddy H. Fundamentos De Ingeniería De Yacimientos. 2011. p. 8.
47
RAMOS,Jorge Luis. 2018.
53
La simulación de un yacimiento se realiza mediante distintas técnicas aplicadas a la
ingeniería de yacimientos, tales como:
Curvas de declinación. Según Magdalena Paris de Ferrer, existen tres tipos de
curvas de declinación, hiperbólica, exponencial y armónica.
48
PARIS DE FERRER,Magdalena. Fundamentos De Ingeniería De Yacimientos. 2009. p. 505.
49
ESCOBAR H, Freddy. Fundamentos de Ingeniería de Yacimientos. 2011.p. 110
54
Además de los dos métodos mencionados anteriormente la simulación también se
debe contar con distintos métodos analíticos, los cuales se obtienen de la literatura,
acerca del comportamiento de un reservorio. Estos, permiten identificar el porqué
de los distintos comportamientos de un yacimiento a lo largo de su vida productiva.
Por último, se debe tener en cuenta los distintos modelos numéricos, estos son
importantes a la hora de trabajar en cualquier simulación, ya que mediante ellos se
puede dar uso de las distintas correlaciones planteadas por los autores para cada
tipo de yacimiento, estas correlaciones nos permiten de algún modo identificar
distintas propiedades de un reservorio, así como predecir de algún modo su
comportamiento a corto y largo plazo. Estos modelos numéricos, también son
necesarios para dar una solución a la ecuación de difusividad que es aquella que
nos permitirá reconocer como será el comportamiento de los fluidos presentes en el
yacimiento.
50
RAMOS, Jorge Luis. 2018.
51
LOPEZ, Lorena. Introducción a la simulación numérica de Yacimientos. Ing de Yacimientos v.
[En línea], 3 de octubre de 2009 [revisado 14 de agosto de 2018]. Disponible en: http://lmla-
ingyacv.blogspot.com/2009/10/clase-n1-introduccion-la-simulacion.html
55
La clasificación de un simulador de yacimientos según su modelo de flujo se da
básicamente en dos partes, un modelo Black Oil o uno composicional.
El modelo de flujo Black Oil, se refiere a que en el existen solo 3 fases distintas
dentro del yacimiento, las cuales son petróleo, gas y agua. Dentro de este, se debe
asumir que tanto el crudo como el agua son inmiscibles entre sí, es decir que son
fases independientes la una de la otra, a pesar de ser líquidos. Además, se debe
asumir también que el gas es soluble en el petróleo, mas no en el agua.52
Estos tipos de simuladores tienen la capacidad de realizar la simulación de cualquier
yacimiento en el que estén presentes estas tres fases en cualquier proporción. Sin
usados principalmente en yacimientos con producciones únicamente de petróleo y
donde la producción se gas o agua no generan grandes depleciones de presión en
el yacimiento.
La simulación de yacimientos de petróleo negro o Black Oil es especialmente usada
para procesos de inyección de agua o de gas en los cuales estos fluidos son
inmiscibles con el petróleo, y donde no se espera que estos generen grandes
cambios en la composición de los fluidos iniciales en el reservorio. En este se puede
modelar el flujo de agua, petróleo o de gas en donde se toman variaciones de la
solubilidad del gas en el petróleo en función de la presión.53
También, es importante resaltar que los simuladores que se rigen mediante un
modelo de flujo Black Oil incluyen procesos de recobro primario, es decir, aquellos
en los que el yacimiento genera una producción de hidrocarburos por medio de su
naturaleza, dentro de estos podemos encontrar la recuperación por expansión de la
roca y los fluidos, que aunque esta no influya en gran porcentaje en el recobro final
de un yacimiento y su efecto no sea muy prolongado existe en todos los yacimientos
existentes.
En cuanto al modelo de flujo composicional, se considera esencialmente como una
generalización del modelo Black Oil, esto debido a que en este tipo de modelo los
efectos de transferencia de masa son bastante significativos durante el proceso de
desplazamiento del fluido dentro del yacimiento. 54
En el modelo composicional por lo general las tres fases de fluidos en el yacimiento
están bien definidas, la gran diferencia con el modelo Black Oil es que en este se
definen las composiciones reales de las fases del petróleo y gas, y el
comportamiento de estos dentro del análisis PVT lleva un mayor trabajo y sus
especificaciones son más exactas.
52
CARRILLO, Lucio. Simulación de petróleo negro (BLACK OIL). PETROBLOGGER.COM. [En
línea], 30 de julio de 2010, [revisado 14 de agosto de 2018]. Disponible en:
http://www.ingenieriadepetroleo.com/simulacion-de-petroleo-negro-black-oil/
53
CARRILLO, Lucio. Ibid.
54
La comunidad petrolera. Tipos de modelos de simulación de yacimiento. La comunidad petrolera.
[En línea], 27 de abril de 2009, [revisado 15 de agosto de 2018]. Disponible en:
https://www.lacomunidadpetrolera.com/2009/04/tipos-de-modelos-de-simulacion-de.html
56
Los gases, por su parte llevan cada uno de sus componentes al realizar un análisis
PVT, es decir se especifica cada componente presente tanto en los gases presentes
en el yacimiento como los que se desean inyectar si el proceso que se va a plantear
en el simulador es una inyección de algún gas.
También es importante rastrear cada uno de los componentes del petróleo, es decir
especificar si este cuenta con algún gas en solución y los componentes del mismo.
También especificar el tipo de crudo con el que se cuenta y sus distintas
propiedades.
Dentro de este modelo de flujo se pueden realizar distintos procesos de producción
tales como la inyección de cualquier gas en el cual el petróleo sea móvil dentro del
yacimiento, la inyección de un gas en yacimientos que están cercanos a los puntos
críticos del reservorio, procesos especiales como recirculación de gas en cualquier
yacimiento, entre otros.
2.3.1.2 Según régimen de flujo en la roca. Los simuladores de yacimientos
también son clasificados según el régimen de flujo en la roca del yacimiento que se
requiera estudiar, es decir los tipos de canales por los que se van a movilizar los
fluidos dentro del yacimiento, estos canales representados por la porosidad y la
permeabilidad.
Primero, existen los simuladores especializados para un régimen de flujo
convencional en la roca, es decir aquellos en los que el yacimiento tiene una
porosidad simple, que no es más que la porosidad natural con la que cuenta el
reservorio. Aquella que tiene inicialmente la roca después de toda su historia de
compresión y que queda definida luego de la liberación de la presión inicial del
yacimiento en el momento de este ser perforado.
Cabe resaltar que la porosidad simple, al ser generada en los sedimentos
resultantes después de la depositación y la compactación inicial esta depende
principalmente de las distintas características con las que pueda contar los
sedimentos, tales como el tamaño de sus granos, su forma y su clasificación. Esta
tiende a tener un decrecimiento, debido a que la roca cada vez se comprime más y
tiene a sepultarse con el tiempo.55
Los simuladores de yacimientos también son utilizados para otro régimen de flujo
en la roca, el cual se refiere a aquellos yacimientos que tienen una porosidad dual
o secundaria, estos están clasificados como yacimientos fracturados, esta se
caracteriza por tener una relación con los esfuerzos tectónicos que tienen los
reservorios después de la etapa de sepultación y compactación de la roca, también
se puede generar por las fuerzas generadas por los distintos acuíferos o fuerzas de
empuje presentes en el yacimiento.
55
MUÑOZ, Cristofer. Porosidad. Ingeniería Petrolera. [En línea], 4 de febrero de 2017, [revisado el
15 de agosto de 2018]. Disponible en: http://ingpetol.blogspot.com/2017/02/porosidad.html
57
Los yacimientos fracturados también se pueden clasificar en cuanto estos presentan
una permeabilidad dual, es decir, que en el reservorio los fluidos son capaces de
desplazarse a través de la porosidad primaria y la secundaria.56
2.3.1.3 Según modelado de métodos EOR. Otro tipo de clasificación de los
simuladores de yacimiento, también se refiere al modelado de los métodos EOR
(enhaced oil recovery) o procesos de recobro mejorado. Dentro de ellos la primera
categorización se expresa a través de los procesos de recobro químico en los cuales
se involucran la inyección de distintos materiales que no están presentes en un
yacimiento. Dentro de estos tipos de materiales se pueden encontrar polímeros,
emulsiones o combinaciones entre estos.57
Dentro de las invasiones químicas para procesos de recobro mejorado, también,
históricamente se han utilizado distintos materiales como el amonio líquido,
alcoholes, entre otros.
Dentro de los métodos de recobro mejorado, también se pueden encontrar procesos
como los térmicos, que son usados principalmente para calentar el reservorio y de
esta manera disminuir la viscosidad del petróleo para que este pueda fluir más
fácilmente a través del yacimiento. Para realizar este proceso se utiliza
principalmente la inyección de agua caliente al yacimiento o de vapor, que bien,
puede ser de manera continua o alternada.
También, dentro de los métodos térmicos de recobro mejorado está el proceso de
combustión del yacimiento, el cual es calentado mediante resistencias eléctricas las
cuales ayudan a calentarlo para que el petróleo empiece a perder componentes
para aligerarlo o bien para que este pierda su viscosidad, y así, lograr que este fluya
más fácilmente.
Por último, dentro de los modelos EOR, también están los procesos miscibles, es
decir, en los cuales se genera una inyección de fluidos en el yacimiento, los cuales
actúan como un agente desplazante al petróleo. El desplazamiento miscible puede
darse en el primer contacto, esto cuando los dos fluidos son miscibles en todas sus
proporciones. O bien, se puede generar con distintos tipos de gases, los cuales
toman un tiempo en mezclarse dentro del yacimiento.
Dentro de los procesos de recobro por inyección de fluidos miscibles se encuentra
el de inyectar tapones miscibles, es decir inyectar algún tipo de solvente líquido para
que contacte en primer instancia al banco de petróleo y luego enviar por baches
distintas cantidades de gas y agua para que estos empujen el banco mezclado entre
el petróleo y el solvente liquido hacia el pozo productor.
56
Schlumberger. Oilfield Glossary. Yacimiento de doble permeabilidad, [En línea], 2018, [revisado
el 16 de agosto de 2018]. Disponible en: https://www.glossary.oilfield.slb.com/es/Terms/d/dual-
permeability_reservoir.aspx
57
PARIS DE FERRER, Magdalena. Inyección de agua y gas en yacimientos petrolíferos. Segunda
edición. 2007. p. 313.
58
También dentro de los procesos miscibles de recobro de petróleo, clasifica la
inyección alternada de agua y gas o modelo Water Alternating Gas. Este proceso
permite controlar la inestabilidad del frente de desplazamiento y con esto mejorar la
eficiencia de barrido vertical, al disminuir la razón de movilidad, y por lo tanto,
aumentar la eficiencia de barrido volumétrico.58
2.3.1.4 Procesos especiales. Dentro de la clasificación de los simuladores también
se encuentran aquellos que son usados para procesos especiales, es decir aquellos
que se refieren a procesos geomecánicos, API TRACKING, entre otros.
Los procesos geomecánicos se refieren a aquellos en los que se trata de deformar
las rocas y los minerales presentes en el yacimiento, por medio de fracturas
generadas a presión, entre otros. Con el fin de generar nuevos canales de
desplazamiento para los fluidos del yacimiento. Estos se clasifican como procesos
especiales, debido a su alta complejidad geológica.
2.3.2 Ecuaciones que usa el simulador. Los simuladores de yacimientos manejan
distintos métodos de modelos numéricos para cumplir con su objetivo, que en sí, es
predecir y representar el flujo de los fluidos presentes en el reservorio a través de
este.
La simulación numérica de un yacimiento, se rige principalmente por principios de
balance de materiales, donde se toma en cuenta la heterogeneidad presente en el
yacimiento y la dirección del flujo de los fluidos. 59
2.3.2.1 Ecuación de balance de materiales. La ecuación de balance de materiales,
se define como el balance del volumen que se presenta al igualar la producción de
fluidos de un yacimiento, que se representa como el vaciamiento del mismo y la
expansión que presentan los fluidos con las caídas de presión que se generan
durante la producción. Fue desarrollada en el año 1941 por Schilthuis; en esta, se
establece una diferencia entre el volumen total de los fluidos iniciales en el
yacimiento, y la cantidad de fluidos remanentes deben ser iguales a la cantidad de
fluidos producidos.60
Se debe tener en cuenta las siguientes consideraciones para aplicar el balance de
materia:
Volumen poroso constante.
No existe compactación ni subsidencia.
El PVT es representativo del yacimiento
Proceso isotérmico
Las compresibilidades del agua y la formación son despreciables.
58
PARIS DE FERRER, Magdalena. Ibid. p.327.
59
RAMOS, Jorge Luis. Op. Cit.
60
La comunidad petrolera. Ecuación de balance de materiales. La comunidad petrolera. [En línea],
18 de enero de 2009, [revisado 17 de agosto de 2018]. Disponible en:
https://www.lacomunidadpetrolera.com/2009/01/la-ecuacin-de-balance-de-materiales.html
59
Se considera equilibrio termodinámico entre el gas y el petróleo a presión y
temperatura de yacimiento
Dimensión cero61
2.3.2.2 Ley de Darcy. La reconocida Ley de Darcy, descubierta por el francés Henry
Darcy en 1856. En está se describe que la velocidad a la que avanza un fluido en
61
La comunidad petrolera. Ibid.
60
un medio poroso, la Ley de Darcy es proporcional a la permeabilidad y al gradiente
de presión, e inversamente proporcional a la viscosidad del fluido.
La ecuación de Darcy para flujos monofásicos se describe en la Ecuación 4.
Ecuación 4. Ley de Darcy.
𝒒 𝒌 𝒅𝒑
𝒖=
=
𝑨 𝝁 𝒅𝒙
Fuente: EcuRed. Ley de Darcy. [En línea]. [Revisado 22 de enero de 2019].
Disponible en: https://www.ecured.cu/Ley_de_Darcy.
Donde:
μ= Viscosidad del fluido, cp k= permeabilidad de la roca, darcy
A= Sección transversal de la roca, cm2 dp= Densidad del fluido + (Densidad
del fluido)*(gravedad
q= Tasa volumétrica del flujo en especifica)*(distancia recorrida).62
movimiento, cm3/seg
62
PARIS DE FERRER, Magdalena. Ibid. .49.
61
Ecuación 5. Ecuación de estado para gases ideales.
63
Schlumberger. Ecuación de difusividad. Schlumberger Oilfield Glossary. [En línea] 2018,
[revisado 29 de octubre de 2018]. Disponible en:
https://www.glossary.oilfield.slb.com/es/Terms/d/diffusion_equation.aspx
62
3. DESARROLLO E HISTORIA DE PRODUCCIÓN DEL CAMPO
63
establece la producción anual de cada uno de los 242 pozos presentes en el campo
de estudio, desde el inicio de la producción en 1961.
Gráfica 3. Producción anual histórica del campo objeto de estudio
PROMEDIOS DE PRODUCCIÓN
Producción anual (BLS) 18.000.000
16.000.000
14.000.000
12.000.000
10.000.000
8.000.000
6.000.000
4.000.000
2.000.000 Total
0
<31/07/1961
1980
2000
1964
1968
1972
1976
1984
1988
1992
1996
2004
2008
2012
2016
Fecha
64
3.1.1 Problemas asociados al campo objeto de estudio. A continuación se
describen los distintos problemas operacionales que se han presentado en el campo
objeto de estudio a lo largo de su historia productiva.
Salinidad del agua de formación. La alta salinidad del agua de formación
procedente del campo de estudio ha sido un problema para la empresa
ECOPETROL S.A., debido a que su incompatibilidad con la formación genera que
no se pueda disponer este recurso para una eventual expansión del campo,
debido a esto, ECOPETROL S.A se encuentra realizando esfuerzos para definir
algún recurso hidrológico compatible con la salmuera de la Formación Mugrosa,
esto con el fin de poder utilizar este recurso en futuras inyecciones en el mismo
campo o en los campos más cercanos.
65
3.2 INYECCIÓN DE AGUA
66
3.3 HISTORIA DE INYECCIÓN DE AGUA EN EL CAMPO
67
Figura 6. Diseño planta de inyección de agua del campo
68
2250 Psi por pozo por lo cual se esperaba llegar a producir 0.98 MBLS de petróleo
incremental.
3.3.1.1 Definición del área piloto. El piloto de inyección de agua del Campo de
estudio comprende tres pozos inyectores I-57, I-35 y el I-111 formando arreglos
irregulares invertidos, como se evidencia en la Tabla 6.
69
Figura 7. Área piloto de inyección de agua
El patrón de inyección I-35 está configurado por un pozo inyector y tres pozos
productores, los cuales se describen a continuación:
Pozo Inyector I-35. El pozo I-35 esta completado en las subunidades de la
Formación Mugrosa B3-1, B4-0, B4-2 y C.
Pozo Productor P-128. Se encuentra completado en las subunidades de la
Formación Mugrosa B, B3-0, C, C2, C3 y Esmeralda. Comparte con el inyector las
arenas B4-2 y C. Se ha identificado de intervención de cañoneo adicional en
arenas Fósiles, B2-1, B3-1, C, C2, C3 con evaluación selectiva para verificar
afectación en Mugrosa C.
70
Pozo Productor P-25. Completado en las arenas B, B3-0 y B3-1. Comparte con
el inyector la arena B3-1. Se abre intervalo B3-1 en noviembre de 2014.
Pozo Productor P-86. Completado en Fósiles, B2-1, B3-0, B4-2, C, C1, C2, C3 y
Esmeralda. Comparte con el inyector la arena C. Se cañonearon las arenas B2-1,
B3-0 y C con evaluación selectiva para verificar afectación en Mugrosa C para
lograr un balanceo de la inyección del patrón I-35.
I-35
2000
1500
1000
500
0
28-feb-13
31-may-13
31-may-15
29-feb-16
31-may-16
28-feb-17
31-may-17
28-feb-18
31-may-18
31-Aug-13
31-Aug-15
31-Aug-16
31-Aug-17
31-Aug-18
30-sep-12
30-nov-15
30-nov-16
30-nov-17
Fuente: Elaboración propia.
En la Gráfica 4 puede observar que desde que empezó la inyección en el pozo I-35
desde el 2012 hasta el 2015 tuvo un aumento el caudal de inyección ya que en 2012
empezó la inyección con un caudal de aproximadamente 300 BWPD y en el 2015
llego a tener un caudal de inyección de alrededor de los 1800 BWPD esto debido a
la compleja geología que presente la formación mugrosa. A partir del segundo
semestre del 2015 se puede evidenciar en la gráfica la reducción del caudal de
inyección debido a la irrupción temprana del agua inyectada en los pozos
productores por lo cual la empresa Ecopetrol S.A decidió reducir el caudal de
inyección hasta el 2017 alrededor de los 500 BWPD para así no seguir presentado
los mismos inconvenientes. Entre mediados del 2017 y principios del 2018 se
mantiene un caudal más o menos constante entre los 1000-1200 BWPD hasta que
a mitad del 2018, a inyección cae prácticamente a 200 BWPD debido a un daño en
la bomba de inyección que fue solucionado de manera oportuna para así mantener
un caudal de aproximadamente 1200 BWPD que se está manejando en este
momento en el campo de estudio.
3.3.3.2 Patrón de inyección I-57. Es el segundo patrón de inyección construido en
el campo objeto de estudio, la composición de sus pozos se especifica en Tabla 8.
71
Tabla 8. Patrón de inyección I-57
PATRÓN DE POZO INYECTOR I-57
POZOS INYECTORES POZOS PRODUCTORES
I-57 P-112
P-84
P-109
72
Gráfica 5. Historial de inyección del patrón I-57
I-57
2500
2000
1500
1000
500
0
31-Jan-13
28-feb-14
31-may-14
28-feb-15
31-Aug-15
29-feb-16
28-feb-17
28-feb-18
31-may-15
31-may-16
31-may-17
31-may-18
31-Aug-12
31-Aug-14
31-Aug-16
31-Aug-17
31-Aug-18
31-oct-13
31-jul-13
30-nov-14
30-nov-15
30-nov-16
30-nov-17
30-Apr-13
Para efectos de este trabajo, este patrón de inyección juega un papel fundamental,
debido a que en él se realizara la evaluación de una posible inyección alternada de
gas y agua en las Arenas B de la Formación Mugrosa.
3.3.3.3 Patrón de inyección I-111. Este patrón inicio inyección en agosto del 2015,
considero el pozo P-117 en el patrón, Pero este nunca presento indicios de ser
afectado por la inyección, luego se determinó que entre el pozo inyector I-111 y el
pozo P-117 hay una falla.
En la se especifican los pozos que hacen parte de este último patrón de inyección
de agua implementado en el campo objeto de estudio.
73
Pozo Productor P-117. Completado en Fósiles, Fósiles-1, B, B2-1, B3-1, B4-1,
B4-2, C, C3 y Esmeralda comparte con el pozo inyector I-111 las arenas B y B2-
1.
En junio del año 2015 se empezó con el piloto de inyección del pozo I-111 el cual
inicio con una tasa de aproximadamente 500 BWPD y que a mediados del mes de
octubre del mismo año reporto la tasa más alta de inyección que se haya tenido en
ese arreglo con una tasa de inyección de aproximadamente 1200 BWPD. En febrero
de 2016 se presenta una caída en la tasa de inyección debido a mantenimiento
preventivo de una de las bombas de inyección del campo de estudio por lo cual la
tasa de inyección para ese mes estuvo en 400 BWPD la cual a comparación de al
inicio de la inyección pues fueron casi 800 BWPD menos. De ahí en adelante hubo
un pico de inyección de 800 BWPD la cual se presentó en el mes de abril pero se
tuvo que reducir la tasa de inyección en los siguientes meses debido a que el barrido
no era eficiente por lo cual se volvieron a inyectar a una tasa promedio de más o
menos de 400 BWPD. En abril del 2018 la tasa de inyección cayo nuevamente
debido a problemas operacionales que fueron solucionados y a junio de 2018 ya
se tenía una tasa de inyección aproximadamente de 1100 BWPD una tasa muy
parecida a la que se venía manejando al principio de la inyección, como se aprecia
en la Gráfica 6.
31-Dec-16
31-Dec-17
29-feb-16
28-feb-17
28-feb-18
31-Aug-15
31-Aug-16
31-Aug-17
31-Aug-18
30-jun-15
31-oct-15
31-oct-16
31-oct-17
30-Apr-17
30-jun-16
30-jun-17
30-jun-18
30-Apr-16
30-Apr-18
74
4. CONSTRUCCIÓN MODELO DE SIMULACIÓN
75
4.1.3 Descripción del modelo. El diseño de un proceso de simulación está
influenciado por distintos factores, tales como el tipo de proceso que se desea
realizar, los objetivos principales del estudio, es decir que información es la que se
desea obtener, la calidad de los datos obtenidos o realizados acerca del yacimiento
y la descripción del mismo.
Un factor importante a tener en cuenta es que el modelo que se desea trabajar tenga
concordancia con el modelo real del yacimiento.
4.1.4 Desarrollo de la simulación. El desarrollo de una simulación se realiza
dependiendo de los objetivos planteados al inicio de la simulación. El desarrollo de
una simulación está conformado por el ajuste del cotejo histórico, el cual está
basado en una gran cantidad de parámetros del yacimiento, que generan un modelo
que permita determinar un comportamiento que se asemeje al comportamiento real
del yacimiento.
Por otro lado, se realiza un análisis de sensibilidad de distintos parámetros que
puedan afectar el proceso de simulación. Posteriormente se puede generar un
pronóstico de del comportamiento del yacimiento en el futuro, esto se usa
principalmente para predecir producciones a futuro o identificar oportunidades de
distintos proyectos o procedimientos que se deseen realizar en el reservorio.
Como paso final del desarrollo de la simulación se deben analizar los resultados
obtenidos en la simulación, y a partir de ellos realizar un reporte general de los
mismos, además de identificar si se cumplió o no con los objetivos propuestos al
inicio del proyecto.
76
las Arenas B de la Formación Mugrosa en el Patrón inyector de agua I-57,
especificado con anterioridad.
4.2.1 Mapa estructural. En el modelo estático de este proyecto de investigación,
se cuenta con un mallado suministrado por la compañía ECOPETROL S.A., este
modelo se encuentra constituido por un total de 1.250.700 celdas, las cuales se
encuentran distribuidas en el modelo estructural que se evidencia en la Figura 8.
77
Figura 9. Sector Model escogido por la compañía
78
Figura 10. Ubicación de pozos
4.2.4 Perfiles direccionales. La compañía suministro una base de datos con los
survey direccionales realizados en cada uno de los pozos construidos en el campo,
con base a esta información se realizó un ajuste al formato .WBD requerido por el
simulador, en el cual se establece un orden de columnas en la cual se debe registrar
en orden el nombre del pozo, sus coordenadas en el eje x y eje y, su profundidad
MD, el TVD y el TVDSS.
4.2.5 Topes y bases. Para establecer las características de este ítem, fue necesario
realizar un análisis de los topes y bases, obtenidos en los registros eléctricos de los
pozos presentes en el campo. Posteriormente se realizó un ajuste de esta
información para el simulador, en donde se estableció un tope y una base a cada
pozo con su respectiva profundidad MD.
4.2.6 Volumen Shale (Vsh). La compañía ECOPETROL S.A. suministro los datos
del volumen de arcillas presentes en el modelo.
Se puede corroborar la información suministrada por la compañía en el modelo
estático, mediante la lectura de registros Gamma Ray, con los cuales se pueden
determinar las formaciones ricas en arcillas (aquellas que cuentan con altos valores
de los registros Gamma Ray), con esta información se pueden determinar las líneas
79
de arcilla y arenas. La muestra un resumen de las líneas de arcilla y arenas basadas
en un estudio de registros eléctricos realizado en uno de los pozos del campo objeto
de estudio.
Figura 11. Líneas base de arenas y arcillas para
las Arenas B de la Formación Mugrosa
LÍNEAS BASE DE ARCILLA Y ARENAS
ARENAS ARCILLAS
ZONA IC1 IC2 IC1 IC2
Mugrosa B1 46 40 98 89
Mugrosa B2 46 43 80 81
Mugrosa B3 38 39 90 70
Mugrosa B4 45 42 90 75
Al haber realizado el análisis de los registros Gamma Ray, se debe realizar el cálculo
del volumen de arcilla para la sección de las Arenas B de la Formación Mugrosa.
Como primer paso para realizar este cálculo, se debe hallar el índice Gamma Ray,
el cual se encuentra especificado en la Ecuación 7 para una profundidad
determinada de 5940 pies de profundidad
Ecuación 7. Índice Gamma Ray para una profundidad determinada
88 − 45
𝐼𝐺𝑅 @5940 𝐹𝑡 = = 0.956
90 − 45
Fuente: ASQUITH.George, Basic Well Log Analysis for Geologists, Primera Edición,
1982, p. 87.
80
4.2.7 Factor de Formación. El Factor de formación se define como una relación en
la proporcionalidad entre la resistividad que tiene una roca limpia contenida de agua
(Ro) y la resistividad del agua presente en la formación (Rw). El factor de formación
se relaciona también con la porosidad y el factor de cementación, como se evidencia
en la Ecuación 9, la cual fue propuesta por Archie.
Ecuación 9. Factor de Formación de Archie
𝑎
𝐹= 𝑚
∅
Fuente: SENERGY Ltda, Interactive Petrophysics User Manual, Volumen 3. 2008.
p. 744.
Donde:
a: Coeficiente de litología
m: Factor de cementación
81
El muestreo del fluido del yacimiento fue tomado en las arenas B de la Formación
Mugrosa, más exactamente a 6761 Ft de profundidad y a condiciones de presión de
2350 psi y 144.5°F de temperatura.
4.3.1 Análisis PVT. Las propiedades del fluido importadas al simulador fueron
obtenidas de un análisis PVT desarrollado en el campo a una muestra de fluido por
parte de la Unidad de servicios técnicos y laboratorios del Instituto Colombiano del
Petróleo.
En el análisis de datos del estudio PVT realizado en el campo objeto de estudio, se
estableció una densidad para la muestra de fondo de 0.8643 g/cc @ 5000 psi y
144.5 °F. Una prueba de liberación instantánea del gas permitió identificar que el
fluido cuenta con 288.9 pies cúbicos de gas @ 14.7 psi y 60 °F para cada barril de
petróleo residual a 60 °F y que el aceite residual tiene una gravedad API de 22,65.
Mediante un estudio de las relaciones Presión – Volumen a una temperatura de
144,5 °F realizado por la compañía se estableció que el fluido cuenta con una
presión de burbuja de 2175 psia.
La viscosidad del aceite se midió sobre un extenso rango de presiones a una
temperatura de yacimiento, llegando a la conclusión que la viscosidad a la presión
de burbuja del fluido es de 10,160 cP.
En la Tabla 11 se muestra un resumen de los datos del fluido del yacimiento
después de realizado el estudio análisis PVT.
82
Tabla 12. Resumen de datos de prueba de liberación
diferencial.
RESUMEN DATOS PRUEBA DE LIBERACIÓN DIFERENCIAL
G.O.R. 288.9 SCF/STB
Bo 1.0904 bl/STB
Densidad del aceite 0.8674 g/cc
Gravedad especifica del gas 0.7079
83
Tabla 13. Composición total del fluido de yacimiento
Gas Liquido Recombinado MW
Componentes
Mol % Wt % Wt % Mol % g/mol
Nitrógeno 20,082 0,0000 0,1287 0,8937 28,01
Dióxido de Carbono 0,1580 0,0000 0,0159 0,0703 44,01
Metano 836,869 0,0000 30,712 372,433 16,04
Etano 63,933 0,0000 0,4398 28,452 30,07
Propano 31,142 0,0143 0,3278 14,460 44,1
Isobutano 14,725 0,0291 0,2235 0,7481 58,12
Normal Butano 11,188 0,0323 0,1796 0,6009 58,12
Isopentano 0,7695 0,0894 0,2122 0,5722 72,15
Normal Pentano 0,3867 0,0661 0,1268 0,3419 72,15
Hexanos 0,5086 0,2625 0,3479 0,8057 84
Heptanos 0,2327 0,7533 0,7691 15,583 96
Octanos 0,0979 10,801 10,534 19,150 107
Nonanos 0,0229 13,496 12,926 20,780 121
Decanos 0,0029 13,805 13,166 19,112 134
Undecanos 0,0004 15,869 15,126 20,015 147
Dodecanos 0,0003 15,813 15,072 18,210 161
Tridecanos 0,0001 19,213 18,312 20,354 175
Tetradecanos 0,0000 20,559 19,595 20,060 190
Pentadecanos 0,0000 21,397 20,393 19,256 206
Hexadecanos 0,0000 21,113 20,123 17,631 222
Heptadecanos 0,0000 22,211 21,169 17,374 237
Octadecanos 0,0000 25,021 23,847 18,481 251
Nonadecanos 0,0000 26,819 25,561 18,905 263
Eicosanos 0,0000 25,787 24,577 17,384 275
Henicosanos 0,0000 24,399 23,255 15,544 291
Docosanos 0,0000 24,851 23,685 15,105 305
Tricosanos 0,0000 24,619 23,464 14,353 318
Tetracosanos 0,0000 23,918 22,796 13,396 331
Pentacosanos 0,0000 22,878 21,805 12,294 345
Hexacosanos 0,0000 17,725 16,894 0,9153 359
Heptacosanos 0,0000 23,617 22,509 11,707 374
Octacosanos 0,0000 20,598 19,632 0,9842 388
Nonacosanos 0,0000 21,666 20,650 0,9992 402
Triacontanos plus 0,0000 530,991 506,085 169,724 580
Benceno 0,0141 0,0213 0,0228 0,0568 78,1
Tolueno 0,0121 0,0149 0,0168 0,0353 92,2
84
4.3.3 Prueba de relación presión-volumen. A continuación, en la Tabla 14 se
muestra un resumen de los datos de la prueba de relaciones presión-volumen
realizada a temperatura de yacimiento.
85
En la Tabla 16 se muestran los resultados de las presiones de burbuja obtenidas
mediante una prueba de expansión a composición constante realizada al fluido
original y a los fluidos obtenidos después de cada inyección de gas.
Fuente: INSTITUTO
COLOMBIANO DEL
PETROLEO. Unidad de
servicios técnicos y
laboratorios, área de
optimización de la producción.
2009. Modificado por el autor.
86
En la Tabla 17 se evidencian los factores de hinchamiento del fluido obtenidos en
esta prueba. Dentro de los resultados se evidencia que el incremento de esta
propiedad no presenta un comportamiento típico, lo cual podría indicar una posible
formación de asfaltenos.
Tabla 17. Factor de
hinchamiento (SWF)
Pb (psia) SWF
2160 1
2606 1.0207
2961 1.0291
3415 1.0332
3915 1.0436
Fuente: INSTITUTO
COLOMBIANO DEL
PETROLEO. Unidad de
servicios técnicos y
laboratorios, área de
optimización de la
producción. 2009.
Modificado por el autor.
En la Gráfica 8 se describe el comportamiento del hinchamiento del fluido para cada
uno de los puntos de burbuja descritos anteriormente con un % molar de inyección
distinto.
87
Gráfica 8. Comportamiento del factor de hinchamiento en el
análisis PVT.
4.3.5 Propiedades del fluido. Luego de importar cada uno de los datos obtenidos
en las pruebas PVT descritas anteriormente al programa WinProp de la compañía
CMG, se obtuvo información detallada acerca del fluido del yacimiento y un
modelamiento de las posibles propiedades del mismo con una inyección del gas
preparado en laboratorio y el cual tiene una gran semejanza con el flue gas
proveniente de la Refinería de Barrancabermeja.
4.3.5.1 Viscosidad. Las viscosidades de los fluidos juegan un papel fundamental a
la hora de la elaboración del modelo dinámico en un proceso de simulación. Para
ello se obtienen valores de viscosidad para cada uno de los fluidos que hacen parte
del yacimiento.
En la Tabla 18. Viscosidad del aceite con respecto a la presiónse muestra el
comportamiento de la viscosidad del petróleo con respecto a la presión, obtenidos
de la prueba de separador realizada por la compañía, la cual se realizó a una
temperatura de 85 °F y 60 psig.
En la Gráfica 9 se plantea el comportamiento de la viscosidad con respecto a la
presión. Esta se presenta de forma coherente, ya que la viscosidad decrece al llegar
al punto de burbuja debido a que hay una expansión del gas en solución presente;
mediante que a presiones menores al punto de burbuja la viscosidad tiende a
incrementar debido a la liberación del gas.
88
Tabla 18. Viscosidad del aceite con respecto a
la presión.
VISCOSIDAD DEL ACEITE SATURADO
PRESIÓN (PSIA) VISCOSIDAD DEL ACEITE (cP)
5011 13.54
4510 12.826
4010 12.314
3510 11.651
2513 10.959
2311 10.486
Pb 2175 10.16
2009 10.196
1509 11.693
1009 13.396
509 15.191
108 16.959
89
4.3.5.2 Factor volumétrico. El factor volumétrico está definido como la relación que
existe entre el volumen del aceite a condiciones de yacimiento y el que puede
presentar a condiciones de superficie. El factor volumétrico estipulado para las
Arenas B de la Formación Mugrosa es de 1.094 RB/STB.
4.3.5.3 Gas en solución. El gas en solución está representado por la variable Rs,
esta se ve representada en la Gráfica 10 para valores de presión por debajo del
punto de burbuja. Este ítem se presenta una variación con respecto a la presión,
por valores por encima del punto de burbuja el gas en solución es constante, ya que
a partir de este no se presenta liberación del gas, pero por debajo del punto de
burbuja empieza a disminuir debido a que el gas empieza a liberarse del líquido.
Gráfica 10. Rs fluido
90
Tabla 19. Propiedades petrofísicas @4915 Ft TVDSS.
PROPIEDAD VALOR UNIDADES
Ko a Swi 85.1 mD
Kw a Sor 15.1 mD
K aire 416 mD
Swi 15,0 %
Sor 42.57 %
Porosidad 22.2 % .
Fuente: ECOPETROL S.A. 2012. Modificado por el
autor.
Los datos de permeabilidades relativas para este primer núcleo tomado a una
profundidad de 5469 pies se establecen en Tabla 20.
Tabla 20. Permeabilidades relativas @4915 Ft TVDSS
PERMEABILIDAD PERMEABILIDAD
SATURACION
RELATIVA AL RELATIVA AL Krw/Kro
DE AGUA (fr.)
ACEITE AGUA
0,150 1,0000 0,0000 0,0000
0,170 0,8528 0,0014 0,0016
0,190 0,6703 0,0061 0,0090
0,210 0,5268 0,0104 0,0196
0,230 0,4141 0,0144 0,0347
0,250 0,3254 0,0182 0,0561
0,270 0,2558 0,0221 0,0863
0,290 0,2010 0,0260 0,1293
0,310 0,1580 0,0301 0,1905
0,330 0,1242 0,0345 0,2779
0,350 0,0976 0,0394 0,4033
0,370 0,0767 0,0448 0,5835
0,390 0,0603 0,0508 0,8430
0,410 0,0474 0,0577 1,2170
0,430 0,0372 0,0654 1,7563
0,450 0,0293 0,0742 2,5338
0,470 0,0230 0,0841 3,6538
0,490 0,0181 0,0952 5,2650
0,510 0,0142 0,1077 7,5789
0,530 0,0112 0,1217 10,8956
0,550 0,0088 0,1373 15,6389
0,570 0,0069 0,1546 22,4067
0,574 0,0066 0,1583 24,0720
91
La segunda muestra de núcleos fue tomada en un intervalo entre 5023-5069 pies
TVDSS, en la se establecen las propiedades de esta muestra.
92
Tabla 22. Permeabilidades relativas @5023-5069 Ft TVDSS.
PERMEABILIDAD PERMEABILIDAD
SATURACION
RELATIVA AL RELATIVA AL Krw/Kro
DE AGUA (fr.)
ACEITE AGUA
0,150 1,0000 0,0000 0,0000
0,170 0,8528 0,0014 0,0016
0,190 0,6703 0,0061 0,0090
0,210 0,5268 0,0104 0,0196
0,230 0,4141 0,0144 0,0347
0,250 0,3254 0,0182 0,0561
0,270 0,2558 0,0221 0,0863
0,290 0,2010 0,0260 0,1293
0,310 0,1580 0,0301 0,1905
0,330 0,1242 0,0345 0,2779
0,350 0,0976 0,0394 0,4033
0,370 0,0767 0,0448 0,5835
0,390 0,0603 0,0508 0,8430
0,410 0,0474 0,0577 1,2170
0,430 0,0372 0,0654 1,7563
0,450 0,0293 0,0742 2,5338
0,470 0,0230 0,0841 3,6538
0,490 0,0181 0,0952 5,2650
0,510 0,0142 0,1077 7,5789
0,530 0,0112 0,1217 10,8956
0,550 0,0088 0,1373 15,6389
0,570 0,0069 0,1546 22,4067
0,574 0,0066 0,1583 24,0720
93
En la Tabla 24 en la se representan las permeabilidades relativas para esta muestra
de núcleo.
94
Tabla 26. Permeabilidades relativas @5166 Ft TVDSS.
PERMEABILIDAD PERMEABILIDAD
SATURACION
RELATIVA AL RELATIVA AL Krw/Kro
DE AGUA (fr.)
ACEITE AGUA
95
Tabla 27. Propiedades petrofísicas @5692 Ft TVDSS.
PROPIEDAD VALOR UNIDADES
Ko a Swi 8,6 mD
Kw a Sor 1,3 mD
K aire 34 mD
Swi 16,1 %
Sor 45,4 %
Porosidad 19,2 %
96
En la Gráfica 11 se establece el comportamiento de las permeabilidades relativas
al agua y al aceite, según las muestras de núcleos obtenidas del reservorio, a raíz
de esta gráfica se pueden identificar los distintos tipos de rocas presentes en el
yacimiento, los cuales se encuentran cargados dentro del Sector Model
suministrado por la compañía, aunque falta suministrarle a cada uno de ellos sus
respectivos datos de permeabilidades relativas, presiones capilares y sus
respectivos exponentes de tipo de roca. Esto se realiza en la sección Rock/Fluid,
presente en el software Buider de la compañía CMG.
Gráfica 11. Curvas de permeabilidades relativas para los diferentes tipos de roca.
El suministro de los datos de los tipos de roca al simulador, se realizan para cada
uno de los distintos (4) CUATRO tipos de roca presentes.
4.4.1 Array properties. Para una correcta generación del modelo de simulación, es
necesario establecer una validación de las distintas propiedades del yacimiento
establecidas en el modelo estructural proporcionado por la compañía.
97
Debido a que para efectos del cálculo del petróleo original in situ se ve afectado en
gran medida por la data obtenida por el modelo estático y el área estimada del
modelo. Se realizó un filtro de la información obtenida en la malla de simulación en
cuanto al espesor de los conglomerados de arenas más representativos en el área
de estudio. Para esto, se tomaron como base distintos registros eléctricos realizados
por la compañía, los cuales debido a cuestiones de confidencialidad no es posible
exponerlos en el documento.
Con los datos obtenidos en los registros eléctricos realizados por la compañía en
uno de los pozos del campo objeto de estudio se verifico que los espesores
planteados en el Sector Model debían contar con una actualización la cual se
especificó en la sección 4.2.1.1 Sector Mode. Los espesores presentes en el modelo
se especifican en la Figura 12.
98
4.4 CONFIGURACIÓN DEL MODELO
99
Figura 13. Net to Gross ratio
Mediante la opción de elección de las propiedades estáticas con las que cuenta el
simulador, se puede seleccionar la propiedad que se desea modelar mediante un
mapa, el programa inmediatamente seleccionada la función cambiara a una
distribución de dicha propiedad a lo largo del modelo.
4.4.3.2 Porosidad. La porosidad es una de las propiedades más importantes con
las que cuenta un yacimiento, ya que esta, indica los espacios libres dentro de
cualquier roca. Estos espacios pueden llegar a ser ocupados en su mayoría por
distintos fluidos.
Para el modelo estático del cual se va a partir la elaboración de este trabajo de
investigación no se le realizo ninguna modificación en cuanto a la porosidad, ya que
los valores arrojados por este en el modelo tienen una gran similitud con los datos
petrofísicos de los estudios realizados con anterioridad por la empresa
ECOPETROL S.A.
En la Figura 14 se presenta el mapa de distribución de porosidades en fracciones
como unidad, los valores gravitan entre 0%y 39%.
100
Figura 14. Porosidad del modelo
101
Figura 15. Porosidad promedio arrojada por
Builder.
102
4.4.3.3 Compresibilidad de la roca. Para realizar una estimación de la
compresibilidad de la roca, es necesario evaluar en primera instancia la
compresibilidad de las areniscas y de las lutitas, ya que como se mencionó
anteriormente en la descripción de la formación, se cuenta con areniscas con
intrusiones de lutitas.
Compresibilidad de la arenisca: Se calculó mediante el uso de la ecuación de
Newman Ecuación 11 para areniscas consolidadas y calizas.
103
Figura 16. Permeabilidad del modelo
104
Figura 17. Promedio de permeabilidad del
modelo.
105
Figura 18. Permeabilidad vertical
106
4.4.3.5 Saturación de agua. La saturación de agua está definida como la cantidad
de agua que se puede almacenar en los poros de una roca, en fracción o porcentaje.
Las arenas B de la Formación Mugrosa cuenta también con altos contenidos de
arcilla en su interior, por ende no se puede realizar un cálculo de la saturación de
agua mediante la ecuación de Archie, ya que dentro de sus consideraciones, se
plantea que solo aplica para formaciones limpias.
Para realizar entonces un modelo de saturación de agua para este tipo de
reservorio, el cual presenta altos valores de arcillosidad. Para lograr este objetivo
se trabajara mediante la ecuación de Juhasz, la cual se usa regularmente en
formaciones arcillosas y es una modificación del modelo de saturación propuesto
por Waxman. Este modelo, basa su principio en el efecto de la saturación con la
resistividad, en él se tiene en cuenta la resistividad de las lutitas, mediante su
capacidad de intercambio de cationes.
El intercambio catiónico solo se puede llegar a obtener mediante pruebas de
laboratorio con los núcleos obtenidos de la formación, las cuales fueron realizadas
por la compañía. La resume el modelo propuesto por Juhasz para la saturación de
agua en un medio altamente arcilloso.
Ecuación 13. Ecuación de Juhasz para saturación de agua
𝟏 ∅𝑻𝒎 ∗ 𝑺𝒘𝑻𝒏 𝑹𝒘
− = ∗ (𝟏 − 𝑩𝒏 ∗ 𝑸𝒗𝒏 )
𝑹𝒕 𝒂 ∗ 𝑹𝒘 𝑺𝒘𝑻
Fuente: SENERGY Ltda, Interactive Petrophysics User Manual, Version 3.5, 2008,
p. 909.
Dónde:
Rt: Resistividad de la zona virgen (ohm-m)
n: Exponente de Saturación
m: Factor de Cementación
ɸT: Porosidad Total (Fracción)
SwT: Saturación de Agua Total (Fracción)
Rw: Resistividad del Agua (ohn-m)
Bn: Conductancia de los Cationes de Arcilla
Qvn: Capacidad de Intercambio catiónico Normalizada
107
Ecuación 14. Capacidad de intercambio catiónico
Fuente: SENERGY Ltda, Interactive Petrophysics User Manual, Version 3.5, 2008,
p. 910.
Donde:
Vcl: Volumen de Arcilla (Fracción)
ɸTclay: Porosidad leída del registro en lutitas adyacentes (Fracción)
Donde:
T: Temperatura de Formación (°F)
Rw: Resistividad del Agua (ohm-m)
Para las muestras tomadas por la compañía, para una profundidad de 5541 pies de
profundidad se cuenta con una resistividad de agua de 0,157 ohm-m.
108
Ecuación 16. Conductancia de cationes de arcillas para muestra de núcleo.
109
Fluid del simulador, posteriormente se genera la opción de crear un tipo de roca
nuevo, dentro del cual se le generaron las tablas de permeabilidades relativas de
aceite y agua vs la saturación de agua, y adicionalmente se genera la tabla de
permeabilidades relativas de petróleo y gas vs la saturación del líquido, esta última
se especifica en el ANEXO C, además de su respectiva representación en el
ANEXO D.
En la Gráfica 12 se establece la curva de permeabilidad relativa del petróleo y el
agua vs la saturación de agua para el tipo de roca 1.
Gráfica 12. Permeabilidades relativas del aceite y el agua
para el tipo de roca 1.
110
Figura 21. Permeabilidades relativas para las 3 fases en el
tipo de roca 1.
111
Figura 22. Composición del fluido en WinProp
112
Figura 23. Presión crítica, Temperatura crítica y factores acéntricos
113
Figura 24. Fracción molar de
componentes en WinProp.
Se generara entonces en el software una gráfica de fracción para cada uno de los
componentes como se evidencia en la Figura 25.
114
Prueba de liberación diferencial en WinProp. Se deben seleccionar las
opciones Lab > Differential Liberation, el programa generara una nueva ventana
en la cual se especifican todas las propiedades referenciadas con anterioridad en
la sección 4.3.2 Liberación diferencial, este procedimiento se describe en la Figura
26.
115
aceite incrementa a medida que la presión disminuye después de haber alcanzado
el punto de burbuja.
Gráfica 13. Viscosidad de los fluidos vs presión en WinProp
116
Gráfica 14. Bo y GOR vs presión en WinProp
117
Dentro de la herramienta de Initial Conditions se deben especificar distintos datos
propios del yacimiento, tales como, la presión de referencia, la cual se estableció
por parte del departamento geológico de ECOPETROL S.A., además de la
profundidad, que para el caso de las Arenas B de la Formación Mugrosa es de 5200
pies de profundidad. También se estableció el contacto agua-petróleo, el cual se
sitúa a 5630 Ft, y por último se establece la saturación de agua, la cual establece
el patrón de la misma a lo largo del modelo.
118
en el simulador por medio de la corrida del modelo como se describe a
continuación.
El valor del Original Oil in Place obtenido del modelo suministrado por la compañía,
y el valor calculado con anterioridad para las Arenas B de la Formación Mugrosa,
presentan una diferencia del 2.758% de margen de error. Cabe resaltar que el OOIP
obtenido en el modelo de simulación está por encima del dato que suministro la
compañía. Los valores de petróleo original in situ no se expresan en este trabajo
por motivos de confidencialidad con la compañía.
119
los tipos de roca presentes en el reservorio, debido a que estos se toman en pozos
distantes de la zona de interés, razón por la cual las propiedades de interacción
roca-fluido varían razonablemente, los factores de daño para cada pozo presente
en el modelo es uno de los factores, que también influye dentro del ajuste histórico
del modelo.
120
El margen de error de esta primera sensibilización arrojo un valor del 11.8754%
comparado con los datos creados en el Field History File.
El valor arrojado por la curva de BHP (Bottom Hole Pressure) para el pozo I-57
asumiendo un valor de Skin o daño de formación, presento valores por encima del
registrado por el archivo Field History File, razón por la cual se debe realizar una
segunda sensibilidad en este parámetro para ajustarlo.
4.5.1.1 Sensibilidad 2 (S=3). La segunda sensibilidad del modelo para realizar un
ajuste en el cotejo histórico del modelo, se realizó aumentando en una unidad el
factor de daño de formación en el pozo I-57, esto con el fin de verificar si la curva
de presión de fondo de esta sensibilización logra ajustarse con un mínimo margen
de error a la curva generada para el mismo parámetro por el archivo Field History
Field.
El margen de error para esta segunda sensibilización arrojo un valor de error del
6.7851% con respecto a la curva de historia de presión del pozo, ya que arrojo una
presión de fondo de 1506,6 PSIA. Por esta razón se define esta segunda
sensibilidad, como el cotejo histórico debido a que cuenta con un margen de error
menor al 10%, el cual es aprobado por la compañía.
La gráfica del ajuste, arrojada por la herramienta Results Graph de CMG, se
describe en la Gráfica 15.
121
A raíz de las sensibilidades realizadas, en la ¡Error! No se encuentra el origen de
a referencia. se evidencia la producción acumulada de petróleo, en la cual se refleja
que el segundo caso de ajuste del Skin representa de manera similar la producción
histórica acumulada para el pozo I-57.
Gráfica 16. Producción acumulada para ajuste histórico.
122
Gráfica 17. Tasa liquida para ajuste histórico.
123
Gráfica 18. Tasa de petróleo para ajuste histórico.
124
5. VIABILIDAD DE LA INYECCIÓN WAG MEDIANTE EL USO DEL SIMULADOR
GEM
64
LUKOL; Valeev. Design of WAG parameters. Society of Petroleum Engineers. Moscú, Rusia. 17
de octubre de 2017. p. 6.
125
Figura 28. Factor de recobro de porcentajes de inyección WAG
126
Cronograma de inyección. En el ANEXO E se presenta un cronograma detallado
del proceso de inyección Water Alternating Gas para los próximos DIEZ (10) años
en el campo objeto de estudio.
5.1.1 Flue Gas inyectado. La composición del Flue Gas proveniente de la Refinería
de Barrancabermeja se especifica en la Tabla 32.
127
Ecuación 18 representa la presión mínima de miscibilidad para este proceso. Sin
embargo se disminuye este valor de presión en un 10% para asegurar que no se
genere la miscibilidad en el proceso.
𝟏𝟓𝟏𝟗
( 𝟐 .𝟕𝟕𝟐− 𝐓 +𝟒𝟔𝟎 )
𝐏𝐌𝐌 : 𝟏𝟎𝟕𝟏 . 𝟖𝟐𝟖𝟗𝟑 ∗ 𝟏𝟎 𝐨
128
Figura 29. Respuesta de la inyección WAG comparada con la inyección de
agua en el arreglo I-57.
129
Gráfica 19. Pronóstico de producción a 10 años
130
Así mismo se generan la gráfica de producción de gas entre los años 2018 y 2028,
partiendo de una base de producción de (0) CERO en el arreglo seleccionado para
realizar la inyección WAG en el simulador, la cual se estipula en la Figura 30.
El corte de agua arrojado por el simulador GEM para los próximos (10) DIEZ años
se modela en la Figura 31. La cual presenta una gran cantidad de variaciones con
respecto al tiempo, debido a los ciclos de baches de gas inyectados. Se puede
evidenciar que las caídas en el corte de agua producida son ocasionadas por los
periodos en que se inyecta el Flue Gas
131
Figura 31. Corte de agua para los próximos 10 años
132
6. EVALUACIÓN FINANCIERA DE LA INYECCIÓN WAG EN LAS ARENAS B
DE LA FORMACIÓN MUGROSA
6.1 INGRESOS
133
Precio referencia para los hidrocarburos líquidos. El precio de referencia para
el petróleo es utilizado para categorizar el crudo de acuerdo a su lugar de
proveniencia y sus características. Se usan distintos estándares como referencia,
por ejemplo, el crudo WTI se usa para petróleo proveniente del continente
americano, el BRENT para Europa, el DUBAI para el crudo proveniente de oriente
medio, TAPIS para crudo ligero del lejano oriente y MINAS para crudo pesado del
lejano oriente.
El crudo de referencia WTI se registra para el año 2018 en un promedio de $70.44
USD y este se caracteriza principalmente por contar con una gravedad API de
más de 30° y tener un contenido de azufre de alrededor del 24%. Cuando un
hidrocarburo no presenta estas características previamente establecidas, su
precio puede aumentar o disminuir.
Basic Sediment and Water (BSW). Representa el contenido de agua y
sedimentos que contiene el petróleo al momento de su producción, este factor
tiene alta importancia a la hora de comercializar el crudo, debido a que entre mayor
sea la cantidad de BSW presente en el crudo, su proceso de refinanción se vuelve
más complejo. Los costos del petróleo se reducen o incrementan a la hora de su
venta de acuerdo al porcentaje de presencia de agua y sedimentos.
134
Tabla 35. Ingresos estimados para 10 años
AÑO 1 2 3 4 5
BARRILES PRODUCIDOS PRONOSTICADOS 1,461,496 515,860 223,070 94,230 60,990
INGRESOS $ 80,382,280 $28,372,300 $12,268,850 $ 5,182,650 $ 3,354,450
AÑO 6 7 8 9 10
BARRILES PRODUCIDOS PRONOSTICADOS 44,390 32,070 27,530 28,410 21,210
INGRESOS $ 2,441,450 $ 1,763,850 $ 1,514,150 $ 1,562,550 $ 1,166,550
6.2 IMPUESTOS
6.3 REGALÍAS
En Colombia, las regalías son todo aquel pago que debe realizar cualquier empresa
extractora de recursos no renovables, como lo son la minería o la explotación de
petróleo y gas natural. Cada municipio y departamento en el cual se efectué
cualquier extracción de estos recursos debe estar incluido en la repartición de estas
contraprestaciones, así como también la recibirán todos los puertos marítimos o
fluviales, o estaciones de vías férreas por donde se transitan todos estos recursos.
El campo objeto de estudio comenzó su comercialidad hacia el año de 1960, con
una participación del 100% de la compañía ECOPETROL S.A., la compañía como
tal no debe realizar pagos de ningún tipo por realizar cualquier exploración dentro
del campo, aunque si debe realizar el pago de un porcentaje del crudo explotado en
el área.
65
COLOMBIA, CONGRESO DE LA REPUBLICA. Ley 1111 (27 de diciembre de 2006). Bogotá
D.C. 2006. p 5.
135
El presente trabajo de investigación se basara en el Decreto 1760 de 2003, el cual
establece que la ANH (Agencia Nacional de Hidrocarburos) es la encargada de la
administración y recaudación de la regalías a nivel nacional.
Regalías directas. Las regalías directas son todas aquellas compensaciones que
benefician directamente a las entidades territoriales a las cuales pertenecen, es
decir, a los municipios o departamentos en donde se adelantan trabajos de
explotación de hidrocarburos, a los municipios con puertos marítimos o fluviales
por donde se transportan estos recursos y a las áreas de influencia de estos
puertos.
Regalías indirectas. Son recursos no asignados directamente a los
departamentos o municipios productores de recursos no renovables, estos son
distribuidos por el fondo nacional de regalías, el cual es aprobado por el Congreso
de la República de Colombia.
La ley 756 de 2002, determina el valor de las regalías que deben pagar las
compañías explotadoras de recursos no renovables, el porcentaje de aporte al
Sistema Nacional de Regalías varía según la producción de petróleo, como se
establece en la Tabla 36.
6.4 INVERSIÓN
Las inversiones se refieren a todo aquel capital necesario para poner en marcha
cualquier proyecto, y darle un mantenimiento a lo largo de su desarrollo. Dentro de
las inversiones necesarias para que un proyecto se lleve a cabo, se encuentra todo
aquel capital que garantice el implementar toda la infraestructura necesaria para
que el proyecto se mantenga, que garantice una producción a lo largo de la vida útil
del mismo, y si es necesario, cualquier tipo de intervención que se tenga que
136
realizar. En la Tabla 37 se establecen las inversiones necesarias para realizar este
trabajo de estudio.
Tabla 37. Inversión
INVERSIÓN
COSTO UNIDAD INVERSIÓN TOTAL
DESCRIPCIÓN CANTIDAD
(USD) (USD)
Lineas de distribución $ 20.880 1433 $ 29.913.386
Compresor $ 450.000 2 $ 900.000
Mano de obra $ 30.000 N/A $ 30.000
Facilidades de superficie -
N/A N/A N/A
Tanque de almacenamiento
TOTAL INVERSIÓN $ 30.843.386
137
desplazar gases a través de algún conducto. Para efectos de este trabajo de
investigación se deben implementar 2 (DOS) compresores, uno a la salida de los
gases desde la Refinería de Barrancabermeja, a lo largo del desplazamiento hasta
el campo objeto de estudio y finalmente en superficie para que el gas pueda ser
inyectado hacia el reservorio.
Línea de tubería. Está compuesta por una gran cantidad de tubos, con distintas
especificaciones técnicas, que son necesarios para desplazar el Flue Gas desde
la refinería de Barrancabermeja hasta el campo objeto de estudio.
6.5 COSTOS
En esta sección, los costos operativos se representan como todo el capital que la
compañía desembolsa con el fin de desarrollar satisfactoriamente el proyecto en
cuestión, además de garantizar la economía del proyecto. Dentro del grupo de
costos, se pueden encontrar el costo de llevar un barril de volumen del fluido de
138
yacimiento hasta superficie, el mantenimiento de los equipos de superficie y
subsuelo, la energía eléctrica necesaria para que los equipos funcionen y los costos
del personal requerido para los procesos operativos.
Costos de producción de los fluidos. Es el desembolso de dinero necesario
para extraer un barril de fluido de formación, llevarlo a superficie y separarlo en
sus tres fases, consiste en su gran mayoría en el costo de separar el líquido del
gas y posteriormente separar el agua del petróleo. El proceso de separación es
indispensable para poder comercializar los hidrocarburos líquidos y gaseosos.
Costos de inyección de agua. Para el caso de estudio, se deben evaluar los
costos que requiere la inyección de todos aquellos volúmenes de agua necesarios
para llevar a cabo la inyección alternada de gas y agua en el campo objeto de
estudio. Todos estos costos se especifican en la Tabla 40.
139
transporte e inyección del mismo gas hasta el yacimiento. Se especifica que la
cantidad de CO2 a inyectar cada año durante el proceso es de 509.689 toneladas,
el costo de cada tonelada de Flue Gas es de $13.15 USD para el primer año,
aunque se especifica un aumento de acuerdo a la inflación del 4% para cada año.
140
Tabla 42. Costos operacionales
COSTOS DE OPERACIÓN
Consumo de energía
Costos de producción Mantenimiento
AÑO eléctrica de los TOTAL (USD)
(USD) (USD)
equipos (USD)
1 $ 35.000 $ 100.684 $ 130.000 $ 265.684
2 $ 36.400 $ 97.153 $ 130.000 $ 263.553
3 $ 37.856 $ 91.666 $ 130.000 $ 259.522
4 $ 39.370 $ 82.753 $ 130.000 $ 252.123
5 $ 40.945 $ 79.982 $ 130.000 $ 250.927
6 $ 42.583 $ 81.513 $ 130.000 $ 254.096
7 $ 44.286 $ 79.564 $ 130.000 $ 253.850
8 $ 46.058 $ 81.610 $ 130.000 $ 257.667
9 $ 47.900 $ 68.237 $ 130.000 $ 246.136
10 $ 49.816 $ 56.005 $ 130.000 $ 235.821
COSTOS DE OPERACIÓN TOTALES $ 2.539.379
6.6 AHORROS
141
S.A. se ve directamente relacionada al pago de este impuesto, debido a que, esta
misma es la que se encarga de administrar la Refinería de Barrancabermeja.
Para realizar el flujo de caja necesario para evaluar este proyecto de investigación,
es necesario identificar el indicador financiero que se va a usar para la evaluación
del mismo, e identificar si este es viable o no para la compañía.
142
● Flujo de caja para los próximos 10 años
Tabla 43. Flujo de caja operativo del proyecto para 10 años
/ AÑO 0 AÑO 1 AÑO 2 AÑO 3 AÑO 4 AÑO 5
INGRESOS $ 80,382,280 $ 28,372,300 $ 12,268,850 $ 5,182,650 $ 3,354,450
COSTOS $ 10,574,991 $ 5,922,540 $ 4,537,348 $ 3,975,057 $ 3,894,042
UTILIDAD BRUTA $ 69,807,289 $ 22,449,760 $ 7,731,502 $ 1,207,593 -$ 539,592
REGALÍAS $ 6,430,582 $ 2,269,784 $ 981,508 $ 414,612 $ 268,356
UTILIDAD ANTES DE IMPUESTOS $ 63,376,707 $ 20,179,976 $ 6,749,994 $ 792,981 -$ 807,948
IMPUESTOS $ 20,914,313 $ 6,659,392 $ 2,227,498 $ 261,684 -$ 266,623
UTILIDAD NETA $ 42,462,393 $ 13,520,584 $ 4,522,496 $ 531,297 -$ 541,325
INVERSIÓN $30,843,385.83
FLUJO DE CAJA NETO -$30,843,385.83 $ 42,462,393 $ 13,520,584 $ 4,522,496 $ 531,297 -$ 541,325
/ AÑO 6 AÑO 7 AÑO 8 AÑO 9 AÑO 10
INGRESOS $ 2,441,450 $ 1,763,850 $ 1,514,150 $ 1,562,550 $ 1,166,550
COSTOS $ 3,838,034 $ 3,907,470 $ 4,170,591 $ 4,107,346 $ 4,174,986
UTILIDAD BRUTA -$ 1,396,584 -$ 2,143,620 -$ 2,656,441 -$ 2,544,796 -$ 3,008,436
REGALÍAS $ 195,316 $ 141,108 $ 121,132 $ 125,004 $ 93,324
UTILIDAD ANTES DE IMPUESTOS -$ 1,591,900 -$ 2,284,728 -$ 2,777,573 -$ 2,669,800 -$ 3,101,760
IMPUESTOS -$ 525,327 -$ 753,960 -$ 916,599 -$ 881,034 -$ 1,023,581
UTILIDAD NETA -$ 1,066,573 -$ 1,530,768 -$ 1,860,974 -$ 1,788,766 -$ 2,078,179
INVERSIÓN
FLUJO DE CAJA NETO -$ 1,066,573 -$ 1,530,768 -$ 1,860,974 -$ 1,788,766 -$ 2,078,179
143
Flujo de caja. Es una representación de los flujos de entrada y de salida de
recursos en un proyecto, los cuales son denominados ingresos y egresos. En la
Figura 32 se representan los ingresos por conceptos de la venta de la producción
de petróleo estimada para los próximos 10 (DIEZ) años y los egresos por
conceptos de inversión en líneas de tubería, instalación de la misma, facilidades
de producción e inyección; además de costos operacionales como costos de
personal, mantenimiento, entre otros.
Valor presente neto (VPN). El valor presente neto está definido como el valor
equivalente de todos aquellos ingresos y egresos presentes en un proyecto. En la
Ecuación 19 se ve expresada la formula necesaria para realizar un correcto
cálculo del valor presente neto.
144
además de obtener ganancias. Si es igual a 0 (CERO), se recupera la inversión
inicial pero no se obtienen ganancias del proyecto y por último, si el valor arrojado
por la ecuación es menor a 0 (CERO) expresa que los ingresos obtenidos por el
proyecto no son los suficientes para cubrir la inversión generada inicialmente.
Tasa de descuento. La TIO es la tasa de interés mínima de ganancia que
establece cualquier inversionista para invertir en cualquier proyecto. Para efectos
de este proyecto, se establece una TIO del 10% efectivo anual, se toma este valor
por indicaciones de la compañía ECOPETROL S.A.
145
7. CONCLUSIONES
El método de producción primaria del campo fue de empuje por gas en solución.
Posterior a este se produjo un empuje leve de algún acuífero cercano al
reservorio debido a un decremento en las presiones por encima del punto de
burbuja y finalmente una estabilización de la presión del reservorio.
Se realizó una actualización de las permeabilidades relativas para las tres fases,
con el fin de especificar el comportamiento de cada una de ellas para cada uno
de los (4) CUATRO tipos de roca establecidos por el departamento geológico de
ECOPETROL S.A.
146
La presión de inyección del Flue Gas fue establecida a partir de un promedio
entre la presión máxima de inyección referida en el estudio de hinchamiento del
Flue Gas a condiciones de yacimiento (3915 psi) y la presión mínima de
miscibilidad del modelo (1388,947 psi). Con el fin de evitar una precipitación de
asfáltenos derivada del hinchamiento del gas en el reservorio y de garantizar
que el desplazamiento del fluido a raíz de la inyección, se presente mediante un
método miscible.
Se estableció una tasa de inyección de Flue Gas del 25% en relación con el agua
inyectada para el proceso de inyección Water Alternating Gas, instaurado a raíz
del éxito de los distintos procesos de este tipo de inyección en el mundo.
147
8. RECOMENDACIONES
Agregar al modelo estático todos los pozos del bloque junto con los sets de
registros eléctricos correspondientes a cada uno de ellos, con el fin de obtener
mayor confiabilidad en la data.
148
BIBLIOGRAFÍA
CARRILLO, Lucio. Simulación de petróleo negro (BLACK OIL). 2010. [En línea]
[Consultado el 14 de agosto de 2018]. Disponible en internet:
http://www.ingenieriadepetroleo.com/simulacion-de-petroleo-negro-black-oil/
149
LARSEN, Johne Alex y SKAUGE, Arne. Simulation of the immiscible WAG process
using cycle-Dependent Three-Phase Relative Permeabilities. Wave Structure in
WAG Recovery. 2011.
150
SCHLUMBERGER. Oilfield Glossary. Ecuación de difusividad. [En línea]
[Consultado el 29 de octubre de 2018]. Disponible en:
https://www.glossary.oilfield.slb.com/es/Terms/d/diffusion_equation.aspx
151
ANEXOS
152
ANEXO A.
DATOS DE POROSIDAD Y PROPIEDADES DE LA ROCA.
153
ANEXO B.
DATOS DE POROSIDAD PARA LAS ARENAS B DE LA FORMACIÓN
MUGROSA.
154
ANEXO C.
PERMEABILIDAD RELATIVA GAS-PETRÓLEO.
155
ANEXO D.
CURVA DE PERMEABILIDAD RELATIVA GAS-PETRÓLEO.
156
ANEXO E.
CRONOGRAMA DE INYECCIÓN WAG.
157