173526
173526
173526
Director
Manuel Enrique Cabarcas Simancas
M. Sc. Ingeniería de Petróleos
Nota Escaneada
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 3
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 4
DEDICATORIA
A Dios por darme la sabiduría y paciencia para terminar con éxito este proyecto y mi carrera
Ingeniería de petróleos.
A mis padres Henry Gomez y María del Rosario Rios los cuales se han esforzado toda su vida
por brindarnos lo mejor a mis hermanos y a mí. Por todas las enseñanzas, los regaños y el amor
hacia nosotros que han hecho de mí una excelente persona.
A mis hermanos Henry Felipe y Daniela que de alguna u otra manera siempre estuvieron
acompañándome y ayudándome en todo el transcurso de mi carrera.
A mi familia en general.
A Paula y Leidy mis amigas del colegio las cuales siempre han estado para apoyarme y darme
ánimo.
A Silvia por ser una excelente amiga, apoyarme siempre durante la carrera y contarme sus
experiencias y sufrimientos con la tesis.
A Sergio Andrey por el apoyo y palabras positivas durante el desarrollo del proyecto.
DEDICATORIA
por permitirme realizar este trabajo de grado con el cual culmina esta etapa.
A quienes han sido mi mayor motivación para culminar mis estudios y quienes han estado
presentes durante todo este proceso, mi padre Serafín Rojas y mi madre Sandra Molina, ellos que
con gran esfuerzo me han brindado su apoyo incondicional a lo largo de todo este proceso de
crecimiento personal y profesional, gracias por confiar en mí y creer que este sueño un día podría
ser realidad.
Tabla de Contenido
Pág.
Introducción .................................................................................................................................. 16
1. Objetivos ................................................................................................................................... 19
3.1 Medidas para la protección del medio marino por labores de la industria mar adentro. ........ 55
3.1.1 Convenio sobre la protección del medio marino de la zona del Mar Báltico (1974) .......... 56
tuberías. ......................................................................................................................................... 63
Lista de Figuras
Pág.
Figura 18. Presión de cabeza de pozo de diferentes pozos a nivel mundial. ................................ 87
Figura 20. Diseño por milla de la línea de pozo Kronos a plataforma. ........................................ 90
Figura 21. Diseño por milla de temperatura ambiente y coeficiente de transferencia de calor
Figura 24. Velocidad del gas vs longitud para el gas pobre. ........................................................ 93
Figura 25. Presión vs longitud para el gas pobre a diferentes caudales. ....................................... 94
Figura 28. Velocidad del gas vs longitud para el gas intermedio. ................................................ 97
Figura 29. Presión vs longitud para el gas pobre a diferentes caudales. ....................................... 97
Figura 32. Velocidad del gas vs longitud para el gas rico. ......................................................... 100
Figura 33. Presión vs longitud para el gas rico a diferentes caudales......................................... 100
Figura 34. Presión vs longitud para diferentes composiciones a diferentes caudales. ............... 101
Figura 35. Hold up vs longitud para diferentes composiciones a diferentes caudales ............... 102
Figura 37. Temperatura de formación de hidratos y del fluido vs longitud ............................... 103
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 12
Figura 38. Envolvente de las diferentes composiciones de gases. Perfil del fluido y de formación
Figura 39. Velocidad del gas vs longitud para diferentes composiciones a diferentes caudales. 105
Figura 40. Velocidad del líquido vs longitud para diferentes composiciones a diferentes caudales.
..................................................................................................................................................... 105
Figura 41. Mapa de patrón de flujo propuesto por Xiao et al. .................................................... 106
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 13
Lista de Tablas
Pág.
Tabla 6. Lineamientos generales para la selección del sitio y la evaluación de rutas para la
Tabla 7. Medidas de mitigación de los efectos ambientales negativos que posiblemente pueden
....................................................................................................................................................... 68
Resumen
TRANSPORTE SUBMARINO1
PALABRAS CLAVES: Offshore, flujo multifásico, modelamiento, campo de gas, aguas ultra
profundas.
ABSTRACT
DESCRIPTION: Recent discoveries of gas in the Colombian Caribbean Sea invite the
hydrocarbons industry to thoroughly investigate everything related to offshore, so this work
focuses on studying the behavior of the multiphase flow during the transport of gas through
submarine lines, proposing a flow line of the Kronos field located in the “Bloque fuerte sur”.
For the development of this project an extensive bibliographic review was carried out and the
Aspen Hysys simulation tool was used for the design and analysis of the gas line from the well
head to the platform, varying the gas compositions (poor, intermediate and rich), the diameters
and flow rates.
To analyzing the results of the software it is concluded that the pressure drop decreases as the
diameter of the pipe increases, however, for diameters greater than 16 inches, no significant
changes are observed. It is observed that the holdup increases as the diameter of the pipeline is
greater, the gas velocity decreases as the diameter of the pipeline is greater and it is concluded
that a larger diameter in the pipeline means higher investment costs therefore, it is recommended
to use an internal diameter of 16 inches, likewise the gas and liquid velocities were analyzed and
it was concluded that the flow regime found in the flow line was smooth stratigraphic.
The main contribution of this project is the proposal and simulation of a submarine flow line in
the Colombian Caribbean Sea for the Kronos field, which can be of help for the development of
nearby wells such as Purple Angel and Gorgon. In addition, it was designed including the
protected areas of the Caribbean Sea.
3
Bachelor Thesis.
4
Physicochemical Engineering Faculty. Petroleum Engineering School. Director: M. Sc. Manuel Enrique Cabarcas
Simancas.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 16
Introducción
aumento, en el 2015 fue de 3468,6 mientras que en el 2016 fue de 3542,9 billones de pies
cúbicos. Por otro lado, la producción mundial en el 2015 fue de 3538,6 y en el 2016 de 3551,6
billones de pies cúbicos (BP Energy, 2017). Lo anterior evidencia tanto el aumento de consumo,
explotación de yacimientos de gas costa afuera. La industria del offshore es una de las líneas que
En Colombia el gas es una fuente de energía que cada vez toma más fuerza siendo
actualmente la segunda fuente de energía más usada, sin embargo y, a pesar de su importancia, el
potencial gasífero del país no es muy amplio ya que apenas cuenta con 4,4 trillones de pies
panorama no es muy alentador teniendo en cuenta que se consume 10,6 billones de metros
cúbicos anuales (BP Energy 2017), por lo que la industria nacional se ve en la obligación de
importar gas para satisfacer la demanda energética y además, ampliar el panorama exploratorio
desabastecimiento energético que sufre en la actualidad, es por esto por lo que las diferentes
empresas operadoras han tomado iniciativa en el caso. Por ejemplo, Ecopetrol para el año 2020
600 millones de barriles en reservas probadas, para lograrlo, se están concentrando esfuerzos de
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 17
alto valor económico en proyectos de la costa afuera colombiana, que según estudios es una zona
atractiva y con potencial gasífero para el país; en los últimos meses se ha tenido éxito con el
pozo Purple Angel-1 (Colombia Energía, 2017); de este último se tienen buenas expectativas por
antecedentes como el de Orca y Kronos de donde se han confirmado recursos de gas natural
(UPME, 2017).
Sin lugar a duda hay un futuro promisorio en el offshore colombiano y para su desarrollo
El flujo de gas natural en aguas profundas y ultra profundas del caribe colombiano se ve
afectado por las bajas temperaturas del lecho oceánico por lo que se pueden presentar problemas
conocer la naturaleza del flujo multifásico y las condiciones en las que se encuentra.
Un flujo multifásico es definido como el flujo simultaneo de varias fases dentro de una
comportamiento de algunas de sus variables, entre las más importantes están las propiedades
físicas de las fases, diámetro e inclinación de la tubería, presión de operación y régimen de flujo
(Beggs y Brill, 1991); cada una de sus variables describen un comportamiento particular por lo
que también es importante reconocer el ambiente marino en el que va a estar expuesta la línea de
flujo.
El ambiente marino localizado en el bloque fuerte sur del caribe colombiano es un entorno
un tanto hostil para una línea de flujo, teniendo en cuenta pozos como Kronos-1, el cual se
lámina de agua de 1585 metros catalogándose como un pozo ultra profundo (Colombia Energía,
2017), con pozos de estas características es importante tener en cuenta aspectos como las altas
presiones y las bajas temperaturas a las que va a estar sometida la línea de flujo, por lo que es de
vital importancia analizar el comportamiento del flujo multifásico en este tipo de ambientes.
ecosistema marítimo que tiene Colombia en esta zona del caribe colombiano en el cual se
encuentran plantas acuáticas que en muchas partes forman pastos marinos, así como manglares,
y otros tipos de plantas comunes, también destacar que un 9% de los arrecifes de coral del
mundo se encuentran en el caribe (INVEMAR, 2013). Es importante delimitar las zonas viables
para el trazado de las líneas de flujo submarinas, líneas que se encargan de llevar el fluido desde
firme. Así mismo, el trazado de las líneas de flujo debe ser amigable y conveniente tanto con el
medio ambiente, al punto de no atentar contra la preservación de ecosistemas marinos, como con
1. Objetivos
marinos.
Proponer un esquema del trazado de las líneas de flujo submarinas para el desarrollo de
para el bloque Fuerte Sur, bajo diferentes condiciones de operación y composición del
gas.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 20
2. Marco Teórico
El flujo multifásico se presenta cuando dos o más fases viajan a través de un mismo recipiente,
multifásico consiste en una fase continua, o un medio fluido, y una fase dispersa, cualquiera de
las dos puede ser un medio líquido o gaseoso. Por esto el flujo multifásico es un fenómeno difícil
El flujo multifásico describe sistemas multicomponente en los que la interacción entre los
diferentes componentes tiene una gran influencia en la estructura en general del flujo. Para la
industria del petróleo y gas, el flujo multifásico es el flujo combinado de gas, petróleo y agua en
una tubería (Jacobsen, 2012). Hasta el momento hay muy pocos casos en los que el problema de
un flujo multifásico pueda ser simplificado y asumirse como un flujo monofásico menos
complejo. Teniendo en cuenta que aspectos como el perfil de velocidad, turbulencia y capa
limite son características difíciles para describir la naturaleza de estos flujos. El flujo multifásico
se puede presentar un rango de presiones y temperaturas muy amplio, estas presiones pueden
variar desde 15000 psia hasta casi condiciones atmosféricas y sus temperaturas pueden oscilar
desde 400°F hasta temperaturas por debajo del punto de congelación del agua (Brill, 1987).
mayoría de los primeros investigadores usaron datos obtenidos de pruebas de laboratorio, aunque
algunos usaron datos de campo. Estos datos generalmente incluían caudales volumétricos de gas
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 21
patrón de flujo que radicaba en que las fases de líquido y gas viajaban por la tubería a diferentes
velocidades, de esta forma se dieron a conocer diversos patrones de flujo que con el tiempo
dieron cabida a las correlaciones empíricas que se conocen hoy en día. (Brill, 1992). Además de
las consideraciones sobre las caídas de presión asociadas con el tamaño de tubería, se deben
tener en cuenta otros factores que se presentan en el flujo de dos fases, tales como: hold up de
líquido, densidad, velocidad superficial, deslizamiento entre fases y los patrones de flujo.
2.2.1 Gradiente de presión. La base teórica para las principales ecuaciones de flujo de fluidos
es la ecuación general de energía, una expresión para el balance o conservación de energía entre
dos puntos de un sistema. Brown (1977) citado por Abd Elmoniem y El-Banbi (2015), Como se
Donde:
= Energía interna,
= Energía potencial,
de ecuación de gradiente de presión (Allen Jr., T. y Ditsworth, 1975), finalmente quedando así:
Donde:
Donde:
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 23
De la ecuación anterior, se puede inferir que para calcular el gradiente de presión en flujo de
tuberías es necesario estudiar las caídas de presión debido a la fricción, la elevación (cambios en
Torres y Trauwvitz (2008) plantean que la componente de elevación es tomada sólo sobre la
figura 1. El componente de elevación para flujo vertical o inclinado es por mucho el más
importante de los tres componentes, ya que, para flujo vertical, contribuye generalmente en más
cualquier ángulo de inclinación. La ecuación 1.3 es usualmente adaptada para flujo multifásico
asumiendo que la mezcla gas – líquidos puede ser considerada como homogénea sobre un
Figura 1. Diagrama esquemático general del flujo. Fuente: Torres y Trauwvtitz. (2008). Flujo
2.2.2 Hold Up de líquido (Colgamiento). Puede ser definido como la relación entre el volumen
de líquido existente en una sección específica de tubería a las condiciones de flujo, y el volumen
de la sección aludida. El hold up permite conocer las propiedades físicas de una mezcla de
gas/liquido junto con la velocidad a la que fluye cada fase, lo cual es fundamental para calcular
gradientes de presión en el flujo multifásico. El Hold up ocurre cuando la fase liquida dentro de
la tubería viaja a una menor velocidad que la fase gaseosa y así provocando resbalamiento entre
Donde:
Hold up de líquido,
Figura 2. Ilustración hold up. Fuente: Cabarcas y Rojas (2011). Estudio comparativo entre
natural del flujo a mayor velocidad de una de las fases. El resbalamiento entre fases en el flujo
• La resistencia al flujo por fricción es mucho menor en la fase gaseosa que en la fase líquida.
• La diferencia de compresibilidades entre el gas y el líquido, hace que el gas en expansión viaje
a mayor velocidad que el líquido cuando la presión decrece en dirección del flujo.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 26
que el líquido viaje a menor velocidad que el gas cuando es flujo ascendente, y a mayor
2.2.4 Densidad. La densidad de una mezcla de fluido multifásico fluyendo es muy difícil de
evaluar, debido a la separación gravitacional de las fases y el deslizamiento entre las fases. Esta
puede ser calculada de las densidades del agua y el aceite y de las tasas de flujo si se asume que
Donde:
Donde:
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 27
T=Temperatura, R,
P= Presión, psi.
La densidad del agua in situ depende de los sólidos disueltos en el agua, la temperatura y
Donde:
2.2.5 Velocidad Superficial. Se puede definir como la velocidad de una de fase fluyendo en la
Donde:
Otros autores estudian la velocidad real la cual definen como la relación entre el caudal de
cada fase y el área trasversal que ocupa esa fase y la velocidad superficial bifásica lo cual se
define como la sumatoria de las velocidades superficiales de las dos fases. Rojas y Chacón
(2011).
Velocidad Real. Es la relación entre el caudal de cada fase y el área transversal que
Donde:
2.2.6 Deslizamiento entre fases (Slip). El Slip es un fenómeno físico que describe el
movimiento entre dos fases, producto de la diferencia de velocidades de cada fase. Figura 3.
Figura 3. Deslizamiento entre fases. Fuente: Rojas y Chacón (2011). Estudio comparativo entre
Los patrones de flujo son un fenómeno hidrodinámico muy complejo. Cuando el gas y el líquido
fluyen simultáneamente en una tubería, las dos fases pueden distribuirse en una variedad de
variables geométricas y las propiedades físicas de las dos fases (Yamaha y Abdallah 2010).
se puede evaluar qué fase está en contacto con la pared de la tubería, aspectos de trasferencia de
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 30
calor entre otros datos esenciales para determinar las características del flujo multifásico (Torres,
Trauwvits. 2008).
Los patrones de flujo en tuberías horizontales descritos por Beggs (1973) citado por Torres y
2.3.1 Flujo Laminar Estratificado. Este patrón de flujo se presenta relativamente a bajos flujos
de gas y líquido, para el cual las dos fases son separadas por efecto de la gravedad, donde el
2.3.2 Flujo Laminar ondulado. Básicamente presenta características similares el flujo laminar
2.3.3 Flujo Laminar anular. El flujo anular se presenta a muy altos flujos de gas. La fase
gaseosa fluye como un núcleo a alta velocidad, el cual puede llevar gotas de líquido atrapada. La
fase liquida fluye como una película delgada adherida a la pared interna de la tubería,
generalmente, esta película es más gruesa en el fondo que en la pared superior de la tubería,
2.3.4 Flujo intermitente. Este flujo es caracterizado por el flujo alternado de gas, fluyendo
tubería, superados por burbujas de gas, el cual tiene una capa estratificada de líquido que a su
2.3.5 Flujo intermitente tapón. Es considerado como el caso límite del flujo bache, cuando el
bache de líquido está libre de burbujas, este ocurre a bajos flujos de gas, cuando el flujo es
menos turbulento.
2.3.6 Flujo intermitente bache. Se produce a altos flujos de gas, cuando el flujo del bache toma
2.3.7 Flujo burbuja o burbujas dispersas. Este tipo de flujo se presenta cuando la tubería se
encuentra casi llena de líquido y la fase de gas libre es pequeña, el gas está presente en pequeñas
burbujas distribuidas aleatoriamente, al igual que sus diámetros; las burbujas se mueven a
observa sobre un rango completo de inclinación de tubería, mientras que el patrón de flujo
2.3.8 Flujo niebla. En este caso la fase continua es el gas, el cual arrastra y transporta al líquido
Figura 4. Diagrama esquemático de patrones de flujo en tuberías horizontales. Fuente: Brill, J.P,
Arirachakaran,S. J. (1992). State of The Art In Multiphase Flow. Journal Petroleum Technology.
2.4.1 Flujo burbuja. La tubería está casi completamente llena de líquido y tiene una pequeña
fase de gas libre. El gas está presente como pequeñas burbujas, distribuidas al azar, cuyos
diámetros también varían al azar. Las burbujas se mueven a diferentes velocidades dependiendo
de sus respectivos diámetros. El líquido asciende por la tubería a una velocidad básicamente
uniforme y, a excepción de su densidad, la fase gaseosa tiene poco efecto sobre el gradiente de
presión.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 33
2.4.2 Flujo tapón (Slug). El flujo tapón es un flujo muy complejo e inestable, la predicción de la
caída de presión, la fracción de vacío y otros parámetros son una tarea difícil. (Yahaya y
Abdallah, 2010). En este régimen la fase de gas es más pronunciada, aunque, la fase liquida es la
fase continua, las burbujas de gas se fusionan y forman burbujas estables de aproximadamente el
mismo tamaño y forma que son casi el diámetro de la tubería y están separados por baches de
líquido: La velocidad de la burbuja es mayor que la del líquido y puede predecirse en relación
con la velocidad del bache de líquido, A velocidades de flujo más altas, el líquido puede incluso
ser arrastrado en las burbujas de gas. Tanto las fases de gas como las líquidas tienen efectos
significativos sobre el gradiente de presión. Este flujo está compuesto por una serie de slugs de
gas, siendo la cabeza de los mismos generalmente despuntada, mientras que su final es plano con
una estela de burbujas. Dado una tasa de flujo de líquido constante, un incremento de la tasa de
flujo de gas conduce a un alargamiento y un rompimiento del tapón de gas (Sánchez, Montbrun,
El flujo más común en regímenes de flujo múltiple es el Slug o bache, el flujo tiene algunas
cambios bruscos de presión y fluctuaciones de flujo violentas. Este fenómeno puede causar
productividad, acortar el acumulado de vida operacional y generar daño por corrosión, por lo que
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 34
es clave saber cuándo y dónde se puede producir una inestabilidad de flujo como esta. (Li, Guo y
Xie, 2017)
subcategorías de patrones de flujo que se generan en una tubería tipo S o tipo Riser, slugging
severo normal (SS1), Slugging severo transitorio (SST) y slugging en flujo estable (STB)
liquidas. El ciclo completo para esta fase de slugging severo se describe en 4 pasos. Partiendo de
un crecimiento inicial del bache posteriormente una producción de líquido para que después
Figura 5. Etapas de slugging severo tipo SS1. Fuente: Xie et al. (2017). The influence of
longitud de este fenómeno es más corta que la altura del riser. Las etapas de producción de este
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 35
modelo parten de una generación de baches luego se presenta la penetración de gas y por último
el reventón de gas y recuperación de líquido. La principal diferencia entre el SST y SS1 es que la
Figura 6. Etapas de slugging severo de transición. Fuente: Xie et al. (2017). The influence of
El flujo estable (STB) se observa bajo dos condiciones de estudio: a una velocidad de gas alta
y a una velocidad de líquido relativamente alta. El flujo estable tiene una amplitud de presión
baja, pero una frecuencia más alta en relación con los dos modelos anteriormente nombrados.
Tiene un flujo continuo de gas en la base del riser y la velocidad del líquido se caracteriza por
El slugging severo puede ocurrir bajo las siguientes condiciones (Xie et al, 2017):
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 36
Tabla 1.
Parámetro Condición
Tipo de flujo Debe ser estratificado en la tubería aguas arriba de la base del riser.
Fuente: Torres y Trauwvtitz (2008). Flujo multifásico en tuberías. Trabajo de grado ingeniero de
2.4.3 Flujo Transición o “Churn”. Se describe como un flujo agitado, este patrón existe solo
cuando el flujo es ascendente y es de una naturaleza caótica y se producen cambios desde una
fase liquida continua a una fase gaseosa (Yahaya y Abdallah, 2010). Las burbujas de gas pueden
unirse y el líquido puede ser arrastrado por las burbujas, aunque, los efectos del líquido son
2.4.4 Flujo Anular. En el flujo anular, el gas fluye a través del núcleo central de la tubería,
mientras que el líquido fluye a lo largo de las paredes de la tubería como una película. El sistema
se puede considerar como flujo de gas monofásico a través de un diámetro ligeramente reducido
1967. Predicting two-phase pressure drops in vertical pipe. Esso production research Co.
Houston Texas.
desarrollado dos maneras que se describen a continuación con el fin de dar una noción de como
de interpretar el patrón de flujo que lleva un fluido a través de una sección de tubería.
empíricas o modelos mecanisticos. Existe un amplio conjunto de autores que se inclinan por el
(1952), Baxandell y Thomas (1961), Griffit, et al. (1961), Duns y Ros (1963), Fancher y Brown
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 38
(1963), Hagedorn y Brown (1964) , Orkiszewski (1967), Eaton et al (1967), Dukler et al (1969),
Aziz, Govier y Fogarasi (1972), Beggs y Brill (1973), Oliemans (1976), Gregory, Aziz y
Mandhane (1977), Mukherjee y Brill (1985) , Olgas 2P/3P (2006), HTFS citado por piping
Hysys y Baker citado por Torres y Trauwvitz (2008). A su vez, en la literatura se encuentran un
extenso grupo de investigadores que han desarrollado modelos mecanisticos para determinar el
cálculo de la caída de presión como Ozon et al (1987), Barnea (1987), Hasan y Kabir (1988),
Xiao et al. (1990), Tulsa, Chokshi (1994), Ansari at al (1994), Petalas y Aziz (2000) y TUFFP
(2003).
Autores como Yahaya y Al Gahtani (2010) plantean que hay una constante competencia y
desacuerdo entre estas dos escuelas de pensamiento: los primeros señalan que los modelos más
empíricos basados en datos tienen un mejor valor práctico debido a la complejidad infinita de los
predictivo de los modelos basados en mecanismos, que son capaces de generar nuevos
conocimientos.
Autores como Ruiz, Márquez y Brito (2014) en su trabajo no diferencian entre correlaciones
consideran las correlaciones empíricas como modelos, lo cual puede prestarse para confusiones
al momento de estudiar este tema. Hablar de correlaciones empíricas como modelos no esta del
todo erróneo debido a que estos también son modelos, con la diferencia de que se basan en datos
A continuación, se presenta una definición más amplia tanto de las correlaciones empíricas
como de los modelos mecanisticos, después se describirán aquellos que han sido más empleados
en la industria petrolera, los cuales serán empleados para las simulaciones del presente proyecto.
Las correlaciones empíricas son modelos que se formulan estableciendo una relación
matemática basada en datos experimentales con el fin de predecir las pérdidas de presión en una
tubería. La aplicación de correlaciones está limitada al rango de datos utilizados para generar la
correlación, en este caso, al tratarse de un fluido multifásico, hace que las correlaciones tengan
2.5.1 Poettman y Carpenter (1952). Abd El Moniem y El-Banbi,(2015) señalan que esta
correlación fue el primer intento de derivar una correlación de flujo multifásico. El desarrollo de
este método empírico se llevó a cabo utilizando la ecuación general de energía y datos de 34
pozos fluyentes de aceite y 15 con bombeo mecánico, el aceite y el gas fueron considerados para
estar en una sola fase, y no se trató de determinar la correlación de colgamiento. Toda la perdida
factor de perdida de energía, la cual se tomó como constante en toda la longitud de la tubería.
Para el desarrollo de este método se asumió un flujo vertical de aceite, gas y agua como
turbulento, además de no tomar en cuenta los patrones de flujo. Algunos de los parámetros
tenidos en cuenta para el cálculo del gradiente de presión por este método son (Tabla 2):
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 40
Tabla 2.
Diámetro 2, 2 ½ y 3 Pulgadas
Viscosidad <5 cp
RGL <1500
Fuente: Torres y Trauwvtitz (2008). Flujo multifásico en tuberías. Trabajo de grado ingeniero de
Donde:
Flujo másico
2.5.2 Duns y Ros (1963). El modelo de Duns y Ros fue desarrollado para flujo de dos fases en
pozos verticales teniendo en cuenta los patrones de flujo en el pozo. Consideraron tres patrones
dimensionales: el número de velocidades del líquido y del gas (NLV y NGV, por sus siglas en
ingles), el número de diámetro (ND) y el número de viscosidad del líquido (NL), ver las
ecuaciones (1.16) a (1.19). Validaron el modelo con una base de datos de 4000 pruebas con
20,000 puntos de datos, con un diámetro de tubería entre 1.25 y 5.6 pulgadas y velocidades
superficiales de líquido y gas de hasta 3.2 m/s y 100 m/s, respectivamente. Los fluidos utilizados
fueron gas y varios tipos de aceite y agua. (Ruiz, Brito y Márquez. 2014)
Donde
Una vez que se conoce el patrón de flujo en la tubería, las velocidades de los fluidos
adecuadas. Para el patrón de flujo slug y burbuja, la ecuación de caída de presión de fricción
Donde:
= Deslizamiento de liquido
2.5.3 Hagedorn y Brown (1964). Este método trató de determinar una correlación en la cual se
incluyera prácticamente todos los rangos de flujo, un amplio rango de gas líquido, todos los
tamaños de las tuberías usadas ordinariamente y los efectos de las características de los líquidos.
Esta correlación empírica se desarrolló utilizando datos de pozos verticales de 1500 pies, se
tuvieron en cuenta cinco tipos de fluidos diferentes, entre los que estaban el agua y cuatro tipos
de aceites con viscosidades que estaban en un rango de 10 a 110 cp con una temperatura de 80
Abd El Moniem y El-Banbi ,(2015) concluye que los efectos viscosos son insignificantes para
las viscosidades líquidas inferiores a 12 cp, pero deben tenerse en cuenta cuando la viscosidad
del líquido es mayor que este valor. La precisión informada de esta correlación fue de alrededor
2.5.4 Eaton Et Al (1966). La correlación propuesta por estos autores fue desarrollada a partir de
2016).
Donde:
=Flujo Másico
Tabla 3.
2.5.5 Orkiszewski (1967). Esta correlación hizo énfasis en que el colgamiento de líquido fuera
derivado del fenómeno físico observado y que el gradiente de presión fuera relacionado a la
distribución geométrica de la fase de líquido y gas, por lo que determinó que la densidad de la
Por otro lado, el factor de fricción se correlacionó con las propiedades del fluido en la fase
continua por lo que se define que esta correlación se limita a evaluar la caída de presión en dos
Orkiszewski reconoce cuatro patrones de flujo y para cada uno de ellos diseño una
correlación, con ello pudo determinar la velocidad de resbalamiento y la fricción para cada uno
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 45
de los patrones; los patrones son: Burbuja, Bache, Transición Bache-Niebla y Niebla.
Orkiszewski consideró varios conceptos teóricos, pero obtuvo los mejores resultados apoyándose
en los métodos de Duns y Ros, y Griffith y Wallis, por lo que desarrollo su método en base a
Abd El Moniem y El-Banbi (2015), señalan que la precisión informada fue de alrededor del
0,8% con una desviación estándar del 10,3%. La precisión del método se verificó cuando se
2.5.6 Beggs & Brill (1973). Estudio el efecto del ángulo de inclinación de la tubería sobre la
tenía un mecanismo que podía inclinar la tubería de horizontal a vertical; los fluidos utilizados
fueron aire y agua, entre los parámetros estudiados y los rangos de variación se encuentran los
Tabla 4.
Para cada diámetro de la tubería variaban los flujos de gas y de líquido, por lo que se lograron
observar todos los patrones de flujo, las variables que afectaban el gradiente de presión del fluido
multifásico para esta prueba experimental fueron el colgamiento y el ángulo de inclinación, por
lo que se evaluó cuidadosamente cada parámetro; el diseño experimental para poder determinar
el gradiente de presión está dado por la siguiente ecuación. (Torres y Trauwvitz, 2008):
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 47
Donde:
Abd El Moniem y El-Banbi ,2015 concluyen que la precisión informada de la correlación fue
Se han descrito algunas de las correlaciones empíricas más importantes las cuales pueden ser
Tabla 5.
Vertical
Eaton et al. 1966 Si Si
Horizontal
Horizontal e
Beggs & Brill. 1973 Si Si
inclinadas
Fuente: Modificada Piping Equipment- Hysys.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 48
La tecnología del flujo multifásico fue tomando importancia y se llegó a la conclusión de que
dando así la pauta al desarrollo de los modelos mecanísticos. (Torres y Trauwvtitz, 2008).
Los modelos mecanísticos se basan en un enfoque fenomenológico que tiene en cuenta los
desarrollaron modelos separados para predecir las características hidráulicas verticales del pozo,
2010).
2.6.1 Modelo Mecanístico de Govier, Aziz y Fogarasi (1976). Este modelo fue diseñado para
calcular la perdida de presión, retención y régimen de flujo. Fue desarrollado luego de un estudio
de caída de presión en pozos productores de gas y condensado, se analizaron a detalle 102 pozos
con relaciones de gas líquido de 3900 a 1170000 pies cúbicos por barril. (Ruiz, Brito y Márquez,
2014)
Los autores desarrollaron su propio mapa de patrones de flujo para identificar cuatro
regímenes de flujo: burbuja, slug, niebla y zona de transición. Govier y Aziz dan la ecuación de
energía mecánica para el flujo vertical de una mezcla de una o varias fases escrita para un
Donde:
Termino de corrección
Después de la estimación del régimen de flujo, el modelo correspondiente se usa para calcular
el retraso del líquido y las pérdidas de presión por fricción. Govier y Aziz desarrollaron
2.6.2 Modelo Mecanístico de Xiao et al. (1990). Es un modelo mecanístico para el flujo
predecir primero el patrón de flujo existente y luego detectar las características de flujo,
principalmente retención de líquido y caída de presión, para los patrones de flujo burbuja
grandes diámetros de campo, seleccionado de la base de datos A.G.A. para flujo multifásico y
datos de laboratorio publicados en la literatura. Los datos incluyen tanto aceite negro como
fluidos composicionales.
el rendimiento de algunas correlaciones utilizadas para el flujo bifásico en tuberías como Beggs
y Brill, Mukherjee y Brill , Dukler y Dukler- Eaton.. La evaluación de este método basado en las
caídas de presión, demuestra que el rendimiento general del modelo Xiao et al propuesto es
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 50
mejor que el de cualquiera de las correlaciones, con el mínimo porcentaje de error absoluto y la
patrones de flujo y características de flujo tales como hold up y caída de presión. El modelo se
evalúa utilizando una base de datos compuesta por 1.712 casos de pozos que cubren una amplia
variedad de datos de campo, fue desarrollado como parte del programa de investigación de Fluid
Flow Projects (TUFFP) de la universidad de Tulsa. Este modelo mecanicista se desarrolló para
flujo ascendente de dos fases y predice la existencia de cuatro patrones de flujo, a saber; flujo
Hagedorn y Brown, Duns y Ros, Orkiszewski con la corrección de Triggia, Beggs y Bri1l con
la corrección de Palmer, Mukherjee y Brill, Aziz et al y Hasan and Kabir como modelo
mecanístico.
Ansari et al (1994), concluye en su investigación que el modelo desarrollado por él, presenta
literatura. Así mismo el trabajo de Ruiz, Brito y Márquez (2014) soportan lo dicho por Ansari et
al ya que en su estudio la mejor predicción del patrón de flujo fue hecha por Beggs y Brill y
Abd El Moniem y El-Banbi (2015), señalan que la precisión informada por su estudio fue de
alrededor del 12.1% de error absoluto promedio con un 17.1% de desviación estándar.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 51
considera tres patrones de flujo: burbuja, slug y flujo anular, utiliza un enfoque de modelado de
flujo de deriva para la transición de burbuja a slug. Los datos medidos se obtuvieron de 324
pruebas de caudales muy variables. Las predicciones de caída de presión del modelo se
compararon con ocho correlaciones / modelos mecanísticos usando datos medidos y un banco de
mecanístico para predecir patrones de flujo en transición, gradientes de presión, “hold up” de
líquido y características de flujo “slug” en tubería de gas- liquido. El modelo fue validado para
todos los ángulos de inclinación desde 290 a 90° desde la horizontal. Este modelo unificado fue
El concepto del principio fundamental es la condición de las acciones del flujo “slug”, para
los límites de transición con los otros patrones de flujo. Los patrones de transición desde flujo
“slug” a flujo estratificado y/o anular son predicho utilizando la solución de la ecuación de
momento para flujo “slug”. La zona total de la película es tratada como el volumen de control y
tradicionales al usar por separado los modelos mecanísticos y correlaciones para cada transición.
relacionados.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 52
actualizado y validado.
2.6.6 Modelo Mecanístico OLGAS (2006). Es un modelo mecanístico para flujo estacionario.
Este modelo es capaz de predecir la caída de presión en las tuberías, retención de líquido y
patrón de flujo para cada sección de tubería. El modelo considera ecuaciones de continuidad e
impulso para cada fase. Se supone que la composición de la mezcla total es constante con el
tiempo a lo largo de la tubería, mientras que la composición del gas y el líquido cambia con la
correlaciones utilizadas en el modelo dependen del patrón de flujo en la tubería, esto debe
determinarse como primer paso. OLGAS considera cuatro patrones de flujo: estratificado, anular,
formulado para predecir patrones de flujo, caída de presión y “hold up” de líquido en flujo
vertical ascendente para flujo de dos fases. El modelo identifica cinco patrones de flujo: burbuja,
burbujas dispersas, tapón, “churn” y anular. Los modelos usados para la predicción de patrón de
flujo son: Ansari et al. (1994) para flujos de burbujas dispersas y anular; Chokschi (1994) para
flujo burbuja y un nuevo modelo para flujo “churn”. (Piping Equiment- Hysys).
Modelos hidrodinámicos separados para cada patrón de flujo son usados. Un nuevo modelo
hidrodinámico es propuesto para flujo churn y una versión modificada del modelo de Chokshi
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 53
para flujo tapón. Los modelos de Chokshi y Ansari et al. Son adoptados para flujo burbuja y
anular respectivamente.
El modelo ha sido evaluado utilizando una base de datos de 2052 pozos de Tulsa University
Fluid Flow Project, convergiendo en un amplio rango de datos. El modelo ha sido comparado
con los métodos de Ansari et al. (1994), Chokshi (1994), Hasan & Kabir (1994), Aziz et al.
(1972) y Hagerdon & Brown (1964), y se afirma que dieron resultados superiores.
tenido que extenderse y adaptarse a ambientes con algunas condiciones especiales. Con el inicio
de la industria offshore en el año 1891 se empiezan a gestar las primeras ideas para poder llevar a
cuenta que todas las operaciones para el desarrollo de esta industria se llevan a cabo sobre el
lecho marino es importante destacar las características que se encuentran en tal ambiente.
Para el desarrollo de operaciones relacionadas a los hidrocarburos en el lecho marino hay que
tener en cuenta las características ecológicas y de relieve que éste pueda presentar. Es importante
tener en cuenta que para llevar a cabo cada una de las operaciones, se deben diseñar estrategias
que sean amigables con el medio ambiente, esto con el fin de preservar la fauna y la flora
presentes en él.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 54
entre otras, relacionadas con hidrocarburos costa afuera eran llevadas a cabo sin ningún tipo de
protección para los ecosistemas marinos por lo que se causaban algunos daños, generando gran
de los hidrocarburos. Fuente. Isabelle M Cote, Jhon D. Reynolds. (2006). Cambridge University
involuntarias. Por ejemplo, si una tubería se coloca sobre el lecho marino en un lugar en donde
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 55
sedimentos y podría dar lugar a la formación rápida de dunas o bancos de arena que causarán
Authority, 2018)
Por ende, algunas naciones se han reunido con el propósito de mitigar los daños y
3.1 Medidas para la protección del medio marino por labores de la industria mar adentro.
Las regulaciones en materia de prevención y protección del medio marino fueron ampliamente
desarrolladas en la Convención de las Naciones Unidas sobre el Derecho del Mar, conocida
recomendados para la protección del medio marino, teniendo en cuenta las características propias
resultante de las actividades relativas a los fondos marinos, los cuales no pueden ser menos
instalaciones y estructuras para realizar actividades en el lecho y el subsuelo del mar (art. 56)
(Lares, 2017).
Así mismo, existen convenios regionales que contienen normas sobre la responsabilidad civil
3.1.1 Convenio sobre la protección del medio marino de la zona del Mar Báltico (1974). La
comisión encargada de liderar este convenio se encuentra en Helsinki y vela por la preservación
de los ecosistemas marinos que se encuentran en la zona concerniente al mar Báltico. (Barberan,
1974).
Artículo 10. Exploración y explotación de los fondos marinos y su subsuelo. Se define que,
cada parte contratante deberá adoptar todas las medidas oportunas a fin de prevenir la
contaminación del medio marino de la Zona del mar Báltico que resulte de la exploración o la
explotación de los fondos marinos. Velará también por que se tenga a mano el equipo apropiado
Artículo 11. Cooperación en la lucha contra la contaminación marina. Se establece que las
partes contratantes deberán adoptar medidas de cooperación entre sí, con miras a eliminar o
reducir al mínimo la contaminación de la Zona del mar Báltico por el petróleo u otras sustancias
del Mediterráneo (PAM), el primer acuerdo regional bajo los auspicios del Programa de
emergencia causadas por hidrocarburos y otras sustancias perjudiciales. Así mismo como el
denominado (Protocolo Offshore) aprobado en 1994 y, en 2011, luego del accidente del
Deepwater Horizon, la Unión Europea decidió adherirse al Protocolo Offshore, haciendo posible
Artículo 15. Protocolo offshore. En el punto 3 se expone que la autoridad competente exigirá
un certificado de seguridad e idoneidad para el objetivo. Dicho certificado deberá ser expedido
sistemas de carga, tuberías y demás instalaciones que pueda especificar la parte contratante.
Offshore, 2011).
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 58
Mediterráneo y proteger y mejorar el medio ambiente marino, así como contribuir a su desarrollo
3.1.3 Convenio de Cartagena. Se trata de un acuerdo que busca prevenir, reducir y controlar la
El convenio de Cartagena tiene efecto sobre la Región del Gran Caribe, que comprende el
medio marino del Golfo de México, el Mar Caribe y las zonas adyacentes del océano atlántico.
El comité encargado de velar por el cumplimiento del convenio está en la facultad de coordinar
Cartagena, 1983). Este convenio encargó al PNUMA las funciones de secretaría; y para ello fue
Artículo 8. Contaminación resultante de actividades relativas a los fondos marinos. Las partes
contratantes adoptarán todas las medidas adecuadas para prevenir, reducir y controlar la
1983).
Artículo 10. Zonas especialmente protegidas. Las partes contratantes adoptarán, individual o
conjuntamente, todas las medidas adecuadas para proteger y preservar en la zona de aplicación
del convenio los ecosistemas raros y vulnerables, así como el hábitat de las especies diezmadas,
amenazadas o en peligro de extinción. Con este objeto las partes contratantes procurarán
establecer zonas protegidas, el establecimiento de estas zonas no afectará a los derechos de otras
Artículo 12. Evaluación del impacto ambiental. Primero: En el marco de sus políticas de
ordenación del medio, las partes contratantes se comprometen a elaborar directrices técnicas y de
otra índole que sirvan de ayuda en la planificación de sus proyectos de desarrollo. De manera
que se prevenga o minimice su impacto nocivo en la zona de aplicación del convenio (Convenio
de Cartagena, 1983).
Segundo: Cuando una parte contratante tenga conocimiento de casos en que la zona de
aplicación del convenio se halle en peligro inminente de sufrir daños por contaminación, o los
haya sufrido ya, lo notificará inmediatamente a los otros estados que, a su juicio, puedan salir
afectados por dicha contaminación, así como a las organizaciones internacionales competentes.
Con el fin de proteger de forma más especifica el medio marino en el que se encuentra instaurada
afuera, algunos países han designado organismos internos, encargados de diseñar directrices de
ecosistemas marinos más vulnerables que se encuentran en el lecho marino dentro de sus límites
nacionales.
Jamaica, uno de los países que más demuestra interés por defender cada uno de sus
Conservation Authority) la facultad de administrar y proteger los recursos naturales del país e
muchos otros países en desarrollo, la calidad de los recursos naturales de algunas zonas costeras
se está degradando cada vez más debido a la alta presión demográfica y los efectos de la
Por otra parte, La NEPA (National Environmental and Planning Agency) postula en su
plataforma virtual las directrices dadas por la NRCA sobre los permisos y aspectos de ingeniería
de proyectos que involucran la presencia de tuberías en la zona costera, una guía que pretende
eliminar o mitigar los impactos potencialmente dañinos que la instalación y operación de dichos
identificar posibles rutas para la tubería. El propósito fundamental es desarrollar bases claras
preferida que sea técnica, económica y ambientalmente factible (Natural Resources Conservation
Authority, 2018).
La selección de una ruta para un oleoducto o gasoducto a menudo viene dictada, en cierta
medida, por factores como topografía costera y a la configuración del litoral. Sin embargo, el
objetivo principal en la selección de la ruta es evitar obstáculos y encontrar la ruta más suave a la
profundidad de agua deseada, manteniendo un recorrido lo más recto posible (Natural Resources
Tabla 6.
Elegir rutas que tengan el menor impacto negativo sobre hábitats críticos o
áreas de desove; rutas de migración; vegetación acuática; especies
amenazadas, raras o en peligro de extinción.
Alineación de Evitar pasar a través de los arrecifes de coral.
ruta Evitar pasar por parques marinos existentes o propuestos, sitios de buceo
o playas de baño.
Evitar pasar por lugares de interés arqueológico, histórico, cultural o
paleontológico.
Evitar lugares donde la costa este sufriendo erosión, acreción,
Estabilidad de progradación. (donde quiera que las líneas de la costa sean inestables, las
la costa tuberías enterradas podrán quedar expuestos o sobrecargarse más de lo
previsto, lo que provocaría interrupciones).
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 62
Sismicidad Evitar sitios que estén cerca de fallas geológicas activas conocidas.
Fuente. Adaptado de Guidelines for the construction, maintenance and monitoring of underwater
pipelines & cables in the coastal zone, Natural Resorces Conservation Authority, Nepa, 2018,
recuperado de
http://nepa.gov.jm/publications/guidelines/Underwater%20Cables%20&%20Pipelines/Undersea
%20Cables%20and%20Pipelines%20-%20part%201.html
Después de evaluar las posibles rutas para el tendido de la tubería, se deben comparar los
durante la operación. También es posible evaluar rutas alternativas sobre la base de minimizar o
eliminar los impactos potenciales mediante el uso de técnicas de mitigación aceptables; el nivel
Tipo y tamaño de tubería, ya sea enterrada o colocada sobre el lecho marino, así como
información especializadas en cada paso que conlleva realizar este proceso. Entre algunas de las
requiere acudir a mapas topográficos, cartas hidrográficas, mapas geológicos, mapas de riesgo de
de tuberías.
Tabla 7.
Medidas de mitigación de los efectos ambientales negativos que posiblemente pueden ocurrir
Desde 1973 con el descubrimiento de Ballena y 1976 con Chuchupa la industria colombiana se
oceánico, teniendo algunos resultados favorables como lo son Kronos, Purple Angel, Orca y
Gorgón.
Para que el progreso de este tipo de proyectos sea notable, es necesario un gran despliegue en
impacto que puede tener el desarrollo de este tipo de proyectos sobre los ecosistemas que se
Colombia cuenta con una posición privilegiada con jurisdicción y soberanía sobre dos
océanos que constituyen un área marina que ocupa aproximadamente el 50% de su territorio,
2008).
Tabla 8.
La Política Nacional para el Desarrollo Sostenible de los Espacios Oceánicos y las Zonas
anteriormente nombrados en la Tabla 8. Los bienes y servicios que ofrecen estos ecosistemas
representan activos estratégicos de gran importancia, no solo para los países ribereños sino para
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 66
el mundo ya que hoy en día, más de la mitad de la población y las más grandes e importantes
exploración de crudo. Y es que los bloques costa afuera se han convertido en una de las grandes
hidrocarburos (Direccion General Marítima, 2014). Se busca que la exploración costa fuera sea
identificada en las áreas marinas de los actuales bloques de exploración, 14 bloques costa afuera
Pacifico, 37 bloques costa afuera Caribe. Tres empresas internacionales han realizado oferta de
exploración que implica la realización de 20 mil kilómetros de sísmica 3D, que equivale a 32 mil
En el país hay contratos de evaluación técnica en cinco bloques ubicados en las cuencas
así como en la de Sinú y Tumaco offshore. Pero en la Ronda 2014 se sumaron cinco nuevas
cuencas: cuatro en aguas ultra profundas del mar Caribe y una en la cuenca Sinú offshore
XXI Cumbre de Rio en su capítulo 17 introduce por primera vez en la agenda internacional la
convención del Convenio de Diversidad Biológica (CDB, 1992) emite el mandato de Yakarta,
Costera, identificando en Manejo Integrado de Zonas Costeras (MIZC) como principal estrategia
de acción para lograr el uso sostenible de los recursos costeros y marinos y en 2001 se formula
en Colombia la PNAOCI.
Tabla 9.
costeros en Colombia.
primera regulación sistémica del sector ambiental en Colombia, se revisan las normas, que, sin
agotar el total existente, se cree, son las más relevantes para la temática abordada (Díaz, 2015).
Tabla 10.
Estatuto Síntesis
ordenación de áreas protegidas o áreas donde haya que tomar medidas especiales para conservar
El Plan de Acción de la PNAOCI 2002-2004, emitido a través del Documento CONPES 3164
Subsistema de Áreas Marinas Protegidas (SAMP), como parte del Sistema Nacional de Áreas
socioeconómica.
Existen áreas protegidas continentales, marinas insulares y marino costeras; las Áreas Marinas
(UICN) como “Cualquier área del territorio intermareal o submareal, cuyos fondos, aguas, flora
y fauna asociadas, así como rasgos históricos y culturales, han sido preservados por las leyes o
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 70
cualquier otra medida efectiva para proteger todo o parte del medio ambiente comprendido en
tuberías costa afuera. Según el Artículo 34, sección 1.3 de la reglamentación emitida por el
instalación costa afuera teniendo en cuenta las condiciones climáticas anticipadas y reportes
compresores, calderas y líneas de flujo deberán diseñarse, instalarse y mantenerse de manera tal
que se propenda por la eficiencia, seguridad de operación y protección del medio ambiente. Los
operadores deberán identificar los materiales apropiados para las líneas de flujo y justificar su
disposiciones e indicar el uso y construcción de las líneas de flujo; deberán diseñar las rutas de
las líneas de flujo con el fin de minimizar su longitud, reduciendo las obstrucciones y
condiciones peligrosas así como también presentar esquemas que detallen el flujo , la presión,
incluir disposiciones para una adecuada identificación que se ajuste a los estándares de la
industria respecto del uso de líneas de flujo (Ministerio de Minas y Energía, 2016).
2016).
costeros. La reserva del Archipiélago de Corales del Rosario y San Bernardo, se encuentra
protegida por la resolución 679 de 2005 del Ministerio de Ambiente Vivienda y Desarrollo
Territorial y es uno de los parques marinos más importantes y más cercanos a los recientes
aguas y parques marinos y, a su vez plantea que el operador costa afuera deberá demostrar una
responsabilidad financiera en el monto de $250.000.000 dólares de los Estados Unidos para cada
una de sus instalaciones, para compensar al Gobierno Nacional y a las comunidades afectadas
por los daños que pudieran ocurrir en el caso de derrame de petróleo de dichas instalaciones
Por otro lado, en las secciones 5.1, 5.3 y 5.5 del Artículo 17, se establece la creación de un
incidente y que el Fondo únicamente cubrirá daños relacionados con actividades de exploración
y producción de hidrocarburos y gas. Esto podrá incluir buques tanque, instalaciones costa afuera
el departamento de la Guajira y la masificación del uso del gas natural en el año 2005
(Belalcazar, 2017). El campo Chuchupa produce gas libre y abastece un 75% de la demanda y
Después de más de dos décadas, Colombia vuelve a reiniciar sus actividades costa afuera a
finales del año 2014 con el descubrimiento de el pozo exploratorio Orca-1 perteneciente al
bloque Tayrona ubicado a 40 kilómetros de la costa del departamento de La Guajira con una
profundidad de 4,420 metros y una lámina de agua de 674 metros (Energía Colombia,2017).
1 en donde se verifico la presencia de hidrocarburos en aguas ultra profundas del sector sur del
Caribe Colombiano, a 3,720 metros de profundidad bajo el nivel medio del mar, donde se
encontró un intervalo de arena neta gasífera de entre 40 a 70 metros (130 a 230 pies). Kronos-1
está ubicado a 53 Kilómetros costa fuera en el bloque de exploración y producción Fuerte Sur
Pozo Purple Angel -1. Este se encuentra a 4.7 kilómetros de distancia de kronos-1. El pozo
Purple Angel-1 alcanzó una profundidad total de 4.795 metros, incluida una lámina de agua de
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 73
1.835 metros, y registró intervalos con presencia de gas cuya suma total se estima entre 21 y 34
pozo exploratorio Gorgon-1 con posible presencia de gas en aguas ultra profundas en el sur del
Caribe colombiano. Este hallazgo que cuenta con una profundidad entre los 3.675 y los 4.415
metros bajo el nivel del mar ratifica la alta prospectividad hidrocarburíferos costa afuera con la
que cuenta el país (ANH, 2017). Gorgon-1 está ubicado en el bloque Purple Angel a 27
kilómetros al norte del pozo Purple Angel-1, que recientemente confirmó la extensión del
encontraron intervalos de arena neta gasífera entre 80 y 110 metros (260 a 360 pies) (Boletín
Ecopetrol, 2017).
En abril del 2017 se descubrió el pozo Siluro, localizado en el bloque RC-11 en el caribe a 45
kilómetros del cabo de la vela y 40 kilómetros de Manaure se estima que el pozo debe explorarse
a una profundidad de 2.090 metros (protafolio,2017), Este pozo fue perforado y resulto seco. En
el segundo semestre del año 2017 se descubrió el pozo Brahma y Ecopetrol comenzó la
perforación del pozo costa afuera Molusco, este está ubicado en unas de las cuencas más
El mar caribe se trata de una cuenca suboceánica del océano Atlántico, justo en la región
tropical. El mar Caribe es un gran cuerpo de agua salada ubicado al sureste del Golfo de México
y al oeste del océano Atlántico, entre las latitudes 9 º y 22 º Norte y las longitudes 89 º y 60 º
Oeste. Limita al sur con Colombia, Venezuela y Panamá, al oeste con Costa Rica, Nicaragua,
México, Honduras y Belice; con Cuba, Jamaica, República Dominicana y Puerto Rico en su
parte norte y con las Antillas Menores al sur. En el mar caribe se encuentra una gran cantidad de
El bloque fuerte sur, ubicado en el caribe colombiano justo frente a los departamentos de
Córdoba y Antioquia se encuentra dentro del golfo Darién. Dentro de este bloque se encuentran
El área del polígono formado por el bloque fuerte sur es de cuatrocientos setenta y tres mil
cuatrocientos noventa (473.490) hectáreas con cuatrocientos doce (412) metros cuadrados. El
Geográfico Agustín Codazzi, cuyas coordenadas planas GAUSS con origen Bogotá, Datum
Figura 10. Bloque fuerte sur. Fuente: Contrato de exploración y producción de hidrocarburos
ubicado en el bloque fuerte sur. Esta simulación servirá de ayuda para los campos cercanos como
Purple Angel y Gorgon. A continuación, se nombran las coordenadas de cada uno de ellos.
realizaron bajo los criterios de conservación ambiental de la zona en el mar caribe. Ver figura 11.
Se analizaron los mapas de zonas protegidas en el caribe colombiano (figura 12) y se realizó el
sondeo para la ruta de la línea y la profundidad de las zonas potenciales para la ubicación de la
plataforma. Las coordenadas de la plataforma propuesta son Lat 9° 22’58.68’’N – Long 76°
consideró la protección de los Corales del Rosario y de San Bernardo. En este se reporta una
Figura 12. Mapa zonas protegidas en el mar caribe Fuente: Invemar- Instituto de
Los datos del perfil batimétrico para el campo Kronos fueron obtenidos por medio del software
Google Earth PRO, con base en las coordenadas del campo. Dichas coordenadas corresponden a
la ubicación de los buques perforadores, las cuales fueron tomadas de la página de la Dirección
temperatura del agua del mar como función de la profundidad. La investigación realizada por
la temperatura del agua del mar en 3 grupos de latitudes, los cuales son identificados como 0S,
30S y 60S. El grupo 0S incluye latitudes de -15S hasta 15 N, el grupo 30S incluye latitudes de -
Ya que el pozo Kronos tiene una latitud de 9° 9’53.74’’ N, de acuerdo con la investigación de
Joseph (2010) pertenece al grupo de latitudes 0S. El autor plantea un conjunto de ecuaciones
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 80
para determinar la temperatura en el mar de acuerdo a la profundidad para este grupo, las cuales
se presentan a continuación:
-0.338+ 1.25
1.26
Donde:
S = Temperatura superficial del mar más 0.338, grados Celsius.
D = Profundidad, m.
se instalan sobre el lecho marino y a lo largo del tirante de agua con el fin de transportar los
hidrocarburos desde el pozo hasta la superficie; optimizando costos, de tal manera que confluya
se instala en la parte superior del pozo para tener seguridad y eficiencia en la producción de un
pozo de gas o de petróleo proporcionando una interface entre pozo y árbol submarino. En la
parte superior de dicho cabezal se instala el sistema de árbol submarino. (Fang y Menglan.
2014).
También se puede decir que un árbol submarino es la conexión o interface entre la cabeza del
Jumper o secciones de tubería: este equipo se extiende entre las estructuras submarinas para
servir como enlaces a través de los cuales se transmiten los fluidos. Estos segmentos de tubería
varían en longitud desde unos pocos metros hasta cientos de metros. (Varhaug, 2016)
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 82
varios pozos reduciendo así el número de líneas de producción que deberían ser necesarias.
(Varhaug, 2016)
Plet: estructura que conecta una línea de flujo con otra estructura submarina o se une a la
Esta es la infraestructura básica para los campos offshore según Verhaung, el campo Kronos y
los campos cercanos Purple Angel y Gorgon contarían con cabezales de pozo, arboles de
producción submarinos, jumpers, manifolds submarinos, plets y líneas de flujo los cuales están
conectados a la unidad de producción en superficie, en nuestro caso una plataforma fija a 294
metros.
Varhaug,2016.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 83
Para realizar la simulación de la línea de flujo es necesario contar con la composición del gas,
perfiles de presión y temperatura, entre otros. Datos que aún se desconocen ya que las cuencas
Sinú – San Jacinto y Sinú Offshore se encuentran en estado de exploración y evaluación. Por lo
tanto, la tasa de flujo y las composiciones son objeto de hipótesis hasta que comience la
producción de los pozos. Para obtener estos datos se revisó un conjunto amplio de trabajos
reportados en la literatura sobre diversos campos offshore a nivel mundial teniendo en cuenta
simulación como son la composición de gas, presión en cabeza y fondo de pozo, temperatura en
5.1.1 Composición del gas. Para la simulación se estudiaron tres diferentes clases de
composiciones de gas, las cuales son gas pobre, intermedio y rico debido a que no se conoce con
exactitud la composición del gas de campo Kronos. Las composiciones que se trabajaron en la
Tabla 11.
Composición de gases pobres reportados en la literatura.
Gases
CO2 N2 C1 C2 C3 iC4 nC4 iC5 nC5 C6 C7+ Autor
pobres
Baker
Mezcla 1 1.5 0.69 89.92 5.72 1.74 0.13 0.21 0.09 0.00 0.00 0.00
(2005)
Huang
Mezcla 2 2.33 0.54 91.37 2.67 1.74 0.43 0.30 0.14 0.07 0.41 0.00
(2010)
Casares
Mezcla 3 0.55 1.77 97.21 0.45 0.02 0.00 0.00 0.00 0.00 0.00 0.00
(1997)
Casares
Mezcla 4 1.68 1.00 96.47 0.75 0.09 0.01 0.00 0.00 0.00 0.00 0.00
(1997)
Casares
Mezcla 5 0.01 0.50 95.48 3.48 0.40 0.04 0.07 0.01 0.01 0.00 0.00
(1997)
Promedio 1.21 0.9 94.09 2.61 0.8 0.12 0.12 0.05 0.02 0.08 0.00
Tabla 12.
Composición de gases intermedios reportados en la literatura.
Gases
CO2 N2 C1 C2 C3 iC4 nC4 iC5 nC5 C6 C7+ Autor
interm.
Huang
Mezcla 1 3.10 1.51 84.91 4.44 2.28 0.88 0.62 0.29 0.20 1.77 0.00
(2010)
Casares
Mezcla 2 0.35 0.94 86.36 7.40 3.40 0.50 0.90 0.10 0.05 0.00 0.00
(1997)
Budd
Mezcla 3 0.00 0.40 87.2 7.60 3.10 0.50 0.80 0.20 0.20 0.00 0.00
(2004)
Nemat.
Mezcla 4 4.99 1.71 89.83 2.29 0.51 0.12 0.13 0.06 0.04 0.06 0.24
(1978)
Moshf.
Mezcla 5 1.36 0.03 83.52 7.80 3.81 0.63 0.85 0.44 0.27 0.20 1.09
(2012)
Promedio 1.96 0.92 86.36 5.91 2.62 0.53 0.66 0.22 0.15 0.40 0.27
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 85
Tabla 13.
Composición de gases ricos reportados en la literatura.
Gases
CO2 N2 C1 C2 C3 iC4 nC4 iC5 nC5 C6 C7+ Autor
ricos
Gaard
Mezcla 1 0.88 1.03 77.6 10.6 6.90 0.75 1.61 0.22 0.22 0.11 0.00
(2003)
Hank.
Mezcla 2 1.23 0.64 78.6 9.76 6.29 0.68 1.57 0.25 0.35 0.40 0.15
(1982)
Hank.
Mezcla 3 0.41 0.69 77.8 11.5 7.79 0.53 1.00 0.08 0.08 0.02 0.03
(1982)
Huang
Mezcla 4 0.57 0.72 77.9 3.31 0.94 0.73 0.55 0.27 0.91 14.0 0.00
(2010)
Kgogo
Mezcla 5 3.47 0.34 72.1 8.40 4.90 0.70 2.00 0.65 0.82 0.90 5.65
(2011)
Promedio 1.31 0.68 76.85 8.74 5.36 0.68 1.35 0.29 0.48 3.09 1.17
Con base en las composiciones promedio se realizó la envolvente de fase para cada gas, la cual
las características del yacimiento. Como es sabido las presiones de los campos ubicados en el
bloque fuerte sur aun no son conocidas, por lo anterior se realizó un estudio de los diferentes
campos offshore a nivel mundial que presentan características similares al campo Kronos en
cuanto a la profundidad tanto del yacimiento como de la lámina de agua para determinar los
valores de presión que se usaran en la simulación. Las figuras 16 y 17 presentan las presiones de
fondo y cabeza de pozo respectivamente para cada uno de los campos estudiados. Se determinó
que la presión de fondo varia en un rango de 6000 a 8000 psi, por lo que se manejará, una
Figura 18. Presión de cabeza de pozo de diferentes pozos a nivel mundial. Fuente: Jackson V
et al, (2001). Triolo D. et al. (2013). Lang P.P et al. (1998). Healy J et al. (2013).
5.1.3 Temperatura en cabeza y fondo de pozo. Al igual que la presión, la temperatura también
conocimiento del valor de la temperatura para el campo Kronos, por lo que el valor empleado en
igual a 170°F, por lo tanto, un valor aceptado para la temperatura de cabeza de pozo es 150°F
5.1.4 Caudal de gas y de agua. El caudal de gas es un parámetro importante a tener presente en
la simulación. Como no se conocen datos para el campo Kronos se tuvo en cuenta los caudales
expuestos en los trabajos revisados previamente. Se determinó que los siguientes caudales serian
evaluados en la simulación: 50, 100 y 150 MMSCFD. Además, es importante mencionar, que la
En cuanto al caudal de agua, Rivera et al. (2007) plantea que la producción de agua en
promedio para el campo Chuchupa es cercana a 0.2 bbl/MMSCF, este valor también se encontró
diámetros internos 8, 12, 16, 20 y 24 pulgadas para este estudio. La tubería empleada fue mild
HYSYS es un programa de simulación muy completo, diseñado por la firma de software AEA
temperatura.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 89
yacimientos de gas.
Para iniciar la simulación se debe introducir todas las corrientes y equipos a utilizar en el
proyecto con sus respectivas características como presión, temperatura, composiciones de los
gases y etc.
Los valores de temperatura ambiente fueron estimados con la ayuda de la figura 14, el
software determino el coeficiente global de trasferencia con la correlación PROFES y con los
Figura 21. Diseño por milla de temperatura ambiente y coeficiente de transferencia de calor
5.3.1 Gas pobre. La figura 22 presenta el perfil de la presión del fluido durante el trayecto de
49.2 km desde cabeza de pozo hasta la plataforma para el gas pobre a diferentes diámetros
internos de tubería. Se observa que la caída de presión disminuye a medida que aumenta el
La figura 23 presenta el hold up durante el trayecto de la tubería para el gas pobre a diferentes
diámetros internos. Se observa que el hold up aumenta a medida que el diámetro de la tubería es
mayor. Es importante mencionar que el aumento y la disminución del hold up se debe a los
La figura 24 presenta la velocidad del gas durante el trayecto de la tubería para el gas pobre a
diferentes diámetros internos. Se observa que la velocidad disminuye a medida que el diámetro
de la tubería es mayor.
con base en los resultados anteriores emplear un diámetro interno de operación de 16 pulgadas,
ya que genera una menor caída de presión y un menor hold up al compararse con las demás
tuberías.
La figura 25 presenta el perfil de la presión del fluido durante el trayecto de la tubería para el
gas pobre a diferentes caudales con un diámetro interno de tubería de 16 in. Se observa que la
caída de presión es mayor a medida que aumenta el caudal, sin embargo, la variación de la
presión no es mayor a 100 psi para los tres casos evaluados a condiciones de plataforma.
5.3.2 Gas Intermedio. La figura 26 presenta el perfil de la presión del fluido durante el trayecto
de 49.2 km desde cabeza de pozo hasta la plataforma para el gas intermedio evaluado a
diferentes diámetros internos de tubería. Se observa que la caída de presión disminuye a medida
diferentes diámetros internos. Se observa que el hold up aumenta a medida que el diámetro de la
tubería es mayor y de igual manera cambia a medida que varían las elevaciones.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 96
La figura 28 presenta la velocidad del gas durante el trayecto de la tubería para el gas intermedio
a diferentes diámetros internos. Se observa que la velocidad disminuye a medida que el diámetro
de la tubería es mayor.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 97
La figura 29 presenta el perfil de la presión del fluido durante el trayecto de la tubería para el
observa que la caída de presión es mayor a medida que aumenta el caudal, sin embargo, la
variación de la presión no es mayor a 100 psi para los tres casos evaluados a condiciones de
plataforma.
5.3.3 Gas rico. La figura 30 presenta el perfil de la presión del fluido durante el trayecto de 49.2
km desde cabeza de pozo hasta la plataforma para el gas rico evaluado a diferentes diámetros
internos de tubería. Se observa al igual que en los casos anteriores que la caída de presión
La figura 31 presenta el hold up durante el trayecto de la tubería para el gas rico a diferentes
Ya que existen inconsistencias en los datos obtenidos por parte del software cuando se simula
este diámetro. Al analizar los resultados se observó que a una distancia de 35 km se presentaba
un cambio de fase del fluido gaseoso de la tubería a liquido en un 100 %, lo cual es incoherente
al analizar las envolventes de fases y las presiones y Temperatura del fluido. Así mismo se
concluye para ninguno de los casos el hold up superior al 7 % de líquidos durante el trayecto de
la tubería.
La figura 32 presenta la velocidad del gas durante el trayecto de la tubería para el gas rico a
diferentes diámetros internos. Se observa al igual que en los casos anteriores que la velocidad
La figura 33 presenta el perfil de la presión del fluido durante el trayecto de la tubería para el
gas rico a diferentes caudales con un diámetro interno de tubería de 16 in. Se observa que la
5.3.4 Comparación de gases. Al comparar la caída de presión entre los gases pobre y rico para
los diferentes caudales simulados, ver figura 34, se observa que la caída de presión es menor para
el gas pobre, lo cual indica, que a mayor cantidad de componentes ricos la caída de presión es
mayor.
tubería es de 16 pulgadas.
Se analizó el hold up entre los gases pobre y rico para diferentes caudales simulados 50,100 y
150 MMSCFD, ver figura 35, se observa que el hold up es mayor para el gas rico, lo cual indica
que a mayor cantidad de componentes ricos se produce mayor formación de líquido. Así mismo
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 102
los cambios presentes en el hold up se debe a las variaciones en las elevaciones de la línea en la
trayectoria.
los caudales y diametros de la linea, para los escenarios estudiando, tal y como se presenta en la
figura 36.
figura 36, se puede observar que la temperatura cae abruptamente de 150 °F, la cual era la
temperatura del agua, que tiene valores muy bajos, cercanos a los 32 °F. Despues del primer
kilometro, la temperatura tiende a ser constante hasta los 30 km, punto en el cual empieza a
Se analizó la temperatura a la cual ocurre la formación de hidratos por medio de los datos
obtenidos por el simulador ver figura. En esta figura se observa que la temperatura a la cual se
forman los hidratos es cercana a los 80 °F a lo largo de toda la línea, al observar el perfil de la
temperatura del fluido, se observa que la temperatura del fluido se encuentra por debajo de la
temperatura de formación de hidratos a partir de 160 metros del pozo, lo cual indica que desde
este punto se formaran hidratos. Además, el simulador señala que se forman hidratos tipo II, sin
hidratos, las envolventes de cada gas y se presenta el perfil del fluido, es decir, la variación de la
presión y temperatura del fluido durante la línea, desde las condiciones iniciales (presión 5000 y
T 150) hasta las condiciones de plataforma. En esta grafica también se puede observar el cambio
Figura 38. Envolvente de las diferentes composiciones de gases. Perfil del fluido y de
formación de hidratos.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 105
A continuacion en la figura 39 se compara la velocidad del gas entre los gases pobres y rico para
los diferentes caudales simulados en el proyecto. Se puede analizar que se presenta mayor
velocidad del gas para mayores caudales, en este caso las composisiones del gas pobre y gas rico
Figura 39. Velocidad del gas vs longitud para diferentes composiciones a diferentes caudales.
Figura 40. Velocidad del líquido vs longitud para diferentes composiciones a diferentes
caudales.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 106
Con base en las velocidades tanto del gas como del liquido se puede determinar el patron de
flujo del fluido presente en la linea para cada uno de los escenarios propuestos haciendo uso del
mapa propuesto por Xiao et al. (1990) en su trabajo, Ver figura 41. Este mapa presenta los
siguientes flujos: estratigráfico liso (SS), flujo estratigráfico ondulado (SW), flujo intermitente
Al determinar el patrón de flujo con las velocidades de gas y liquido (figura 39 y 40) obtenidas
del simulador se determinó que el patrón de flujo presente en la línea de flujo del proyecto es
estratigráfico liso para todos los casos analizados, ya que la velocidad del gas obtuvo su máximo
valor en 1,2 m/s mientras que la velocidad del líquido tuvo como máximo valor 0.0009 m/s.
Figura 41. Mapa de patrón de flujo propuesto por Xiao et al. Fuente: Xiao et al .1990. A
comprehensive Mechanistic Model For Two-Phase Flow in Pipelines. Presentado en The 65th
Annual Technical Conference and Exhibition of the Society of Petroleum Engineers, New
6. Conclusiones
en la literatura concluyendo que el modelo más adecuado para la simulación fue Beggs &
Brill ya que puede ser aplicado en tuberías horizontales e inclinadas y puede determinar
Se estableció una ruta para la línea submarina desde el campo Kronos hasta la plataforma
Se estableció la ubicación de una plataforma fija a 49.2 km del pozo Kronos y a una
profundidad de 294 m. con el objetivo de que esta plataforma sirva para el desarrollo de
se conoce con exactitud la composición del campo Kronos. Se determinó que la caída de
medida que aumenta el diámetro de la tubería, sin embargo, para diámetros superiores a
genera menores costos de inversión, un menor hold up y una menor velocidad del gas al
Los resultados de la simulación indican que el patrón de flujo del fluido a lo largo de la
tubería fue un patrón estratificado liso donde el líquido fluye en el fondo de la tubería y el
gas en la parte superior, este resultado se presentó para todos los escenarios estudiados y
se analizó con base en el mapa de patrones de flujo propuesto por Xiao et al (1990).
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 109
7. Recomendaciones
Evaluar el impacto ambiental generado por la tubería sobre el lecho marino y la fauna y
Cartagena.
Realizar una simulación teniendo en cuenta los campos Purple Angel y Gorgon para
Referencias Bibliográficas
Kelleher, G. (2002). google books. Obtenido de Cómo evaluar una AMP: Manual de
indicadores naturales y sociales:
https://books.google.com.uy/books?id=ma3JFYljJwQC&pg=PA43&lpg=PA43&dq=Kellehe
r,+G.+2002+AREAS+MARINAS&source.
Kgogo, T.C.(2011). Well Test Analysis of Medium-Rich to Rich Gas Condensate Layered
Reservoirs. Presentado a SPE Europec/EAGE Annual Conference and Exhibition, Vienna ,
Austria, 23-26 Mayo. SPE 143621.
Lang P.P et al.(1998). Mensa Project: System Desing and Operation. Presentado en The
Offshore Technology Conference. Houston, Texas. 4-7 de Mayo.
Lares, M. C. (2017). Medidas para la protección del medio marino en la región del Gran
Caribe por daños ocasionados por la industria mar adentro. Revista de Derecho [online], n.
47, 207-233.
Lemay, M. H. (1998). Manejo de los recursos costeros y marinos en America Latina y el
Caribe. Washington, D.C.: Banco Interamericano de Desarrollo.
Ministerio de Minas y Energía. (11 de febrero de 2016). Ministerio Minas y Energía.
Obtenido de
https://www.minminas.gov.co/documents/10180/674559/PUBLICACI%C3%93N+REGLA
MENTACIO%CC%81N+OPERACIONES+COSTA+AFUERA+%40%2011-02-
2016.pdf/bbd144c1-878e-4235-89a0-2376c9053423
Moshfeghian, M. (2012). Transportation of Natural Gas in Dense Phase, PetroSkills, John
M. Campbell.
Natural Resources Conservation Authority. (2018). National Environment and Planning
Agency. Recuperado el 15 de 05 de 2018, de Guidelines for the construction, maintenance
and monitoring of underwater pipelines & cables in the coastal zone:
http://nepa.gov.jm/publications/guidelines/Underwater%20Cables%20&%20Pipelines/Under
sea%20Cables%20and%20Pipelines%20-%20part%201.html
Nematizadeh, F,. y Betpolice, A. (1978). Pars Offshore Gas Field Development. Presentado
en the European Offshore Petroleum Conference and Exhibition. London 24-27. Octubre,
SPE 8127.
MODELAMIENTO FLUJO MULTIFÁSICO OFFSHORE 114