Correlacion de Orkiszewski

Correlación de Orkiszewski
UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO

FACULTAD DE INGENIERÍA
FLUJO MULTIFÁSICO EN TUBERIAS.
GRUPO: 01
 Bernal Martínez Rodrigo

 Calva Romero Luis Mario
 Franco Jaramillo Uriel
 Fortunato Orozco Sayil Kristel
 Gonzáles Bandala Silvia
 Hernández Benítez Luis
Antonio
 Rivas Ruiz Blanca Estela
1
OBJETIVOS
● Conocer uno de los métodos para calcular caídas de presión en flujo multifásico
vertical.
● Conocer el desarrollo, importancia y uso de la correlación identificando sus
características y parámetros esenciales así como sus diferencias con otros métodos.
● Aplicar la correlación analíticamente a un caso práctico.
2
CONTENIDO
1. Introducción 6. Ejemplo
2. Flujo multifásico en tuberías verticales 7. Conclusiones
3. Correlación de Orkiszewski 8. Bibliografía
3.1 ¿Qué hizo Orkiszewski?
3.2 Observaciones
3.3 Selección y mejora
4. Patrones de flujo
4.1 Flujo Burbuja
4.2 Flujo Bache
4.3 Flujo Transitorio
4.4 Flujo Neblina
5. Procedimiento de cálculo
3
INTRODUCCIÓN
Durante el transporte de los fluidos desde el fondo del pozo y hasta los separadores, existen pérdidas
de energía a lo largo de la tubería e incluso en las líneas de flujo en superficie. Estas pérdidas se
deben principalmente a los efectos gravitacionales, de fricción y cambios en la energía cinética a lo
largo del recorrido que realizan los fluidos.
Mediante correlaciones de flujo multifásico en tuberías, estimamos la capacidad que tiene un pozo
para extraer y transportar fluidos del yacimiento.
El flujo multifásico se divide en:

Vertical
Horizontal
Inclinado
4
Entre las correlaciones para flujo multifásico se

encuentran:.
Para flujo horizontal:

Beegs & Brill, Eaton, Duckler, etc.
Para flujo vertical:

Hagedorn & Brown, Beggs & Brill, Ansari y
Orkiszewski.
5
FLUJO MULTIFÁSICO EN TUBERÍAS

VERTICALES
● El problema de predecir con precisión las caídas de presión en los pozos ha dado lugar a muchas
soluciones especializadas para condiciones limitadas. La razón de estas muchas soluciones es que
el flujo de dos fases es complejo y difícil de analizar incluso para condiciones limitadas de estudio.
● Al pasar los fluidos provenientes del yacimiento a través de la tubería de producción, se consume
la mayor parte de presión disponible para llevarlos del yacimiento a los separadores, por lo que es
de suma importancia realizar una evaluación precisa de la distribución de la presión a lo largo de
la tubería. Al hacerlo conjuntamente con un análisis integral del sistema de producción, es
posible...
6
a. Diseñar las tuberías de producción y

líneas de descarga
b. Proyectar aparejos de producción
artificial (neumático, mecánico, etc)
c. Obtener presión de fondo fluyendo, sin
intervenir los pozos
d. Calcular el efecto de los estranguladores
sobre el gasto
e. Determinar la vida fluyente de los pozos
f. Corroborar los datos obtenidos con las
correlaciones para su ajuste
7
CORRELACIÓN DE ORKISZEWSKI
● Orkiszewski presentó en 1967 su estudio “Predicting Two-Phase Pressure Drops in Vertical Pipe”
dónde proponía un método para predecir con una precisión de aproximadamente el 10% las
caídas de presión para el flujo de dos fases en tuberías verticales. El método es una extensión del
trabajo hecho por Griffith y Wallis y es considerado superior a otros cinco métodos publicados.
● La precisión del método se verificó al comparar los valores obtenidos de la correlación con 148
medidas de caídas de presión.
● La característica de este método con respecto a otros radica en que el colgamiento del líquido se
deriva del fenómeno físico observado y que el gradiente de presión está relacionado con la
distribución geométrica de las fases líquido y gas. Además de que provee una buena analogía de
lo que sucede dentro de la tubería.
● La correlación de Orkiszewski surge de la necesidad de contar con un método que pudiera ser
usado para condiciones de pozo más generales, ya que los métodos con los que se contaban
8
proporcionaban soluciones para condiciones particulares (limitadas).
¿QUÉ HIZO ORKISZEWSKI?

● Basado en la similitud de los conceptos teóricos y dividiendo los métodos publicados en 3
diferentes categorías, Orkiszewski comparó cinco métodos seleccionados y determinó la cantidad
de desviación entre las caídas de presión medidas y predichas. Orkiszewski observó que los
mejores resultados, bajo ciertas condiciones de flujo, se obtenían con los métodos de Griffith y
Wallis, y Duns y Ros, por lo que tomó estas correlaciones como base para desarrollar su método,
combinándolas para los diferentes patrones de flujo considerados.
.
● Las correlaciones de Griffith y Wallis, y Duns y Ros las clasificó dentro del Grupo III de acuerdo a
los criterios considerados para su desarrollo, los cuales son:
a) Se considera resbalamiento entre las fases

b) La densidad de la mezcla se calcula utilizando el efecto de colgamiento
c) El factor de fricción se correlaciona con las propiedades del fluido en la fase continua
d) Se distinguen diferentes patrones de flujo
9
OBSERVACIONES
● Orkiszewski observó que ninguno de los dos métodos demostró ser exacto en toda la gama de
condiciones utilizadas.
○ El método de Griffith y Wallis era confiable para velocidades de flujo bajas, sin embargo no
lo era así para el rango de velocidades altas.
○ El método de Duns y Ros exhibió el mismo comportamiento, excepto que también era
inexacto para los aceites de alta viscosidad en rangos de baja velocidad de flujo.
● Para Orkiszewski el método de Griffith y Wallis parecía proporcionar la mejor base para obtener
una solución general mejorada, aunque sus valores predichos tenían un mayor error que el método
de Duns y Ros.
10
OBSERVACIONES
● El núcleo del método de Griffith y Wallis es la predicción de la velocidad de resbalamiento
derivado de la observación física. Por otra parte, en el trabajo la caída de fricción fue insignificante
por lo que el método para predecir las pérdidas por fricción es una aproximación y, por lo tanto,
está abierto a mejoras.
●
El método de Duns y Ros es difícil relacionarlo con lo que ocurre dentro de la tubería debido a que
presenta un conjunto complejo de parámetros y ecuaciones interrelacionadas.
SELECCIÓN Y MEJORA
El trabajo de Griffith y Wallis se extendió para incluir el rango de alta velocidad de flujo. Al modificar el
método se desarrolló un parámetro para tener en cuenta:
1. La distribución del líquido entre el bache de líquido, la película de líquido y el líquido arrastrado
por la burbuja de gas.
2. El colgamiento de líquido a velocidades altas de flujo.
Este parámetro sirvió para aproximar mejor las pérdidas por fricción y la densidad del fluido.
El método de Griffith y Wallis modificado es suficientemente exacto (error promedio) y preciso

(desviación estándar) en el rango completo de condiciones.
12
PATRONES DE FLUJO
Las diferentes velocidades y geometrías presentes
en el flujo de dos fases influyen fuertemente en las
caídas de presión. Estos factores proveen las bases
para categorizar el flujo bifásico en: burbuja
(bubble), bache (slug), transición (annular-slug) y
neblina (annular-mist).
13
CARACTERÍSTICAS DEL FLUJO BURBUJA

(BUBBLE FLOW)
● La tubería está casi completamente llena de líquido y la fase gas (gas libre) es pequeña.
● El gas se presenta en pequeñas burbujas cuya distribución y diámetros varían aleatoriamente.
● Las burbujas de gas se mueven a diferentes velocidades dependiendo de sus respectivos
diámetros.
● El líquido sube por la tubería con una velocidad bastante uniforme.
● El gas tiene un pequeño efecto sobre el gradiente de presión.
14
CARACTERÍSTICAS DEL FLUJO BACHE

(SLUG FLOW)
● La fase gas es más pronunciada, sin embargo la fase líquida sigue siendo continua.
● Las burbujas de gas se juntan y se forman burbujas más estables de aproximadamente el mismo
tamaño y forma, además de que son casi del diámetro de la tubería.
● Las burbujas de gas se encuentran separadas por tramos o tapones de líquido, y rodeadas además
por una película de líquido.
● La velocidad de la burbuja es mayor que la velocidad del líquido.
● La velocidad del líquido no es constante.
● El bache de líquido siempre se desplaza hacia arriba (en la dirección del flujo), mientras que la
película de líquido puede moverse ya sea en la misma dirección del bache (posiblemente con una
velocidad menor) o incluso podría moverse hacia abajo.
15
CARACTERÍSTICAS DEL FLUJO BACHE

(SLUG FLOW)
● Las variaciones de velocidad del líquido podrían resultar no sólo en una variación de las pérdidas
por fricción en la pared sino también en el colgamiento del líquido, que está influenciado por la
densidad de la mezcla que fluye por la tubería.
● A altas velocidades de flujo el líquido puede ser arrastrado dentro de las burbujas de gas.
16
CARACTERÍSTICAS DEL FLUJO DE TRANSICIÓN

(TRANSITION FLOW)
● El cambio de la fase continua de líquido a la fase continua de gas ocurre en ésta región.
● Los tapones de líquido entre las burbujas de gas virtualmente desaparecen.
● Una cantidad significativa de líquido se arrastra en la fase gaseosa.
● Aunque los efectos de líquido son significativos, la fase gaseosa es más predominante.
17
CARACTERÍSTICAS DEL FLUJO NEBLINA

(ANNULAR-MIST FLOW)
● La fase gas es continua.
● La mayor parte del líquido es arrastrado y transportado en la fase gas.
● Una película de líquido se adhiere a la pared de la tubería pero sus efectos son secundarios
debido a que la fase gas es el factor predominante.
18
PROCEDIMIENTO DE CÁLCULO
1.- Se selecciona un punto donde empezar. Ya sea del fondo del pozo o en la cabeza del pozo.
cabeza de pozo
Fondo de pozo
19
2.- Determinar el gradiente de temperatura.
∆𝑇
𝑇𝑝𝑟𝑜𝑚 =
2
Donde:
𝑇𝑝𝑟𝑜𝑚 = Temperatura promedio

𝑇1 = Temperatura 1
𝑇2 = Temperatura 2
20
3. Encontrar la P promedio.
𝑃1 + 𝑃2
𝑃𝑃𝑟𝑜𝑚 = + 𝑃𝑎𝑡𝑚
2
Donde:
𝑃𝑃𝑟𝑜𝑚 = 𝑃𝑟𝑒𝑠𝑖ó𝑛 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜
𝑃1 = Presión 1
𝑃2 =Presión 2
𝑃𝑎𝑡𝑚 = 𝑃𝑟𝑒𝑠𝑖ó𝑛 𝑎𝑡𝑚𝑜𝑠𝑓é𝑟𝑖𝑐𝑎
21
4.- Suponer un Δh, para realizar los cálculos.
5.- Determinar las propiedades de los fluidos a P media y a la T.
• Densidad del liquido

1 1
𝜌𝐿 = 62.428 𝛾𝑜 + 𝛾𝑤
1+WOR 1+WOR
● Densidad del gas
520 Pprom 1
ρg = γg ∗ 0.0764
T+460 14.7 Z
22
● Velocidad superficial del liquido
5.6142 144 𝑞𝐿 1 𝑊𝑂𝑅

𝑉𝑠𝐿 = 𝐵𝑜 + 𝐵𝑊
86400∗𝐴𝑝 1+𝑊𝑂𝑅 1+𝑊𝑂𝑅
● Velocidad superficial del gas
1
144∗𝑞𝐿 𝑅−𝑅𝑠 1+𝑊𝑂𝑅 14.7 𝑇+460
𝑉𝑠𝑔 = 𝑍
86400𝐴𝑝 𝑃𝑝𝑟𝑜𝑚 520
● Velocidad de la mezcla
𝑉𝑚 = 𝑉𝑠𝐿 + 𝑉𝑠𝑔
23
● Viscosidad de la mezcla del líquido.
1 𝑊𝑂𝑅
𝜇𝐿 = 𝜇𝑜 + 𝜇𝑤
1+𝑊𝑂𝑅 1+𝑊𝑂𝑅
● Tensión superficial de la mezcla de líquido
1 𝑊𝑂𝑅
𝜎𝐿 = 𝜎𝑜 + 𝜎𝑤
1+𝑊𝑂𝑅 1+𝑊𝑂𝑅
24
● Número de velocidad del liquido 𝑁𝐿𝑣
𝜌𝐿 0.25
𝑁𝐿𝑣 = 1.938 ∗ 𝑉𝑠𝐿
𝜎𝐿
● Número de velocidad del gas 𝑁𝑔𝑣

𝜌𝐿 0.25
𝑁𝑔𝑣 = 1.938 ∗ 𝑉𝑠𝑔
𝜎𝐿
● Área de flujo.
𝜋(𝑑ℎ )2
𝐴𝑝 =
4
25
● Gasto gas
3.27∗10−7(𝑍)(𝑞𝑜 )(𝑅−𝑅𝑠 )(𝑇𝑝𝑟𝑜𝑚 +460)
𝑞𝑔 =
𝑃𝑝𝑟𝑜𝑚
● Gasto líquido
𝑞𝐿 = 6.49 ∗ 10−5 ( 𝑞𝑜 𝐵𝑜 + 𝑞𝑤 𝐵𝑤 )
● Gasto total.
𝑞𝑇 = 𝑞𝐿 + 𝑞𝑔
26
7. Determine 𝐿𝐵 , 𝐿𝑆 y 𝐿𝑀 para poder obtener el tipo de régimen de flujo de la siguiente tabla:
Limites de patrones de flujo
Limites Régimen de Flujo

𝑉𝑠𝑔 Burbuja
< 𝐿𝐵
𝑉𝑚
𝑉𝑠𝑔 Bache
 𝐿𝐵, 𝑁𝑔𝑣 < 𝐿𝑠 Donde 𝐿𝐵 , 𝐿𝑠 y 𝐿𝑀 son los limites
𝑉𝑚
adimensionales de:
Burbuja – bache
𝐿𝑀 𝑁𝑔𝑣 𝐿𝑠 Transición bache – niebla Bache – transición
𝑁𝑔𝑣 > 𝐿𝑀 Niebla Transición – niebla
27
Donde:
2
𝑉𝑚
𝐿𝐵 = 1.071 − 2.6616
𝑑ℎ
Con el limite 𝐿𝐵 0.13
𝑑ℎ = diámetro hidráulico de la tubería = área transversal/perímetro mojado
𝐿𝑠 = 50 + 39𝑁𝑙𝑣
0.75
𝐿𝑀 = 75 + 84 ∗ 𝑁𝐿𝑉
28
8. Basándose en el régimen de flujo obtenido en el paso 7 determinar el gradiente de elevación y
de fricción:
a) Régimen de burbuja:
∆𝑝 1
= (𝜌𝐿 ∗ 𝐻𝐿 ∗ 𝜌𝑔 (1 − 𝐻𝐿 ))
∆𝐿 𝑐 144
Donde 𝐻𝐿 se calcula con:
𝐶1 −𝐶2
𝐻𝐿 = 1 −
2
4 0.5
𝐶2 = 𝐶1 2 − 𝑉𝑠𝑔
0.8
𝑉𝑚
𝐶1 = 1 +
0.8
29
El gradiente de fricción se calcula con:
∆𝑝 1 12∗𝑓∗𝜌𝐿 ∗𝑉𝐿2
=
∆𝐿 𝑓 144 2∗𝑔∗𝑑ℎ
El factor de fricción se calcula utilizando un proceso iterativo
30
● b) Régimen de Bache
El gradiente de elevación se obtiene de acuerdo al procedimiento delineado por Griffith y
Wallis:
∆𝑝 1 𝐶3
= + 𝜌𝐿 ∗ 𝛿
∆𝐿 𝑒 144 𝑉𝑚 +𝑉𝑏
Donde:
𝐶3 = 𝜌𝐿 𝑉𝑠𝐿 + 𝑉𝑏 + 𝜌𝑔 ∗ 𝑉𝑠𝑔
31
● 𝛿 = Coeficiente de distribución del líquido. Se calcula como se indica en la tabla.
Tabla: Relación de la ecuación a aplicar, con la velocidad de la mezcla y la fase continua.
Las ecuaciones para obtener 𝛿 son:

𝑑ℎ −1.38 𝑑ℎ
𝛿 = −0.681 + 0.013 ∗ 𝑙𝑜𝑔𝜇𝐿 + 0.232 ∗ 𝑙𝑜𝑔𝑉𝑚 − 0.428 ∗ 𝑙𝑜𝑔 ……(3.88)
12 12
32
𝑑ℎ −0.799 𝑑ℎ
𝛿 = −0.709 + 0.0451 ∗ 𝑙𝑜𝑔𝜇𝐿 − 0.162 ∗ 𝑙𝑜𝑔𝑉𝑚 − 0.888 ∗ 𝑙𝑜𝑔 ……(3.89)
12 12
𝑑ℎ −1.415 𝑑ℎ
𝛿 = −0.284 + 0.0127 ∗ 𝑙𝑜𝑔(𝜇𝐿 +1) + 0.167 ∗ 𝑙𝑜𝑔𝑉𝑚 + 0.113 ∗ 𝑙𝑜𝑔 ……(3.90)
12 12
𝑑ℎ −1.371 𝑑ℎ 𝑑ℎ −1.571 𝑑ℎ
𝛿 = −0.161 + 0.0274 ∗ 𝑙𝑜𝑔(𝜇𝐿 +1) + 0.569 ∗ 𝑙𝑜𝑔 − 0.397 + 0.01 ∗ 𝑙𝑜𝑔(𝜇𝐿 +1) + 0.63 ∗ 𝑙𝑜𝑔 𝑙𝑜𝑔𝑉𝑚 … …(3.91)
12 12 12 12
33
El valor de 𝛿 debe estar dentro de los limites siguientes:
Para 𝑉𝑚 < 10
𝛿  −0.065∗ 𝑉𝑚
Para 𝑉𝑚 > 10
𝑉𝑏 𝜌𝑚
𝛿 − 1−
𝑉𝑚 +𝑉𝑏 𝜌𝐿
34
𝑉𝑏 Velocidad de bache, se determina con las siguientes ecuaciones:
124∗𝑑ℎ∗𝑉𝑚 ∗𝜌𝐿
𝑁𝑅𝐸𝐿 =
𝜇𝐿
124∗𝑑ℎ∗𝑉𝑏 ∗𝜌𝐿
𝑁𝑅𝐸𝑏 =
𝜇𝐿
El procedimiento para calcular la velocidad de bache se inicia suponiendo un valor 𝑉𝑏

= 1.75 pies/seg.
Después este valor se compara con los resultados de las siguientes ecuaciones:
35
Cuando 𝑁𝑅𝐸𝐿 > 6000 tenemos los siguientes tres casos para 𝑁𝑅𝐸𝑏 :
1. Si 𝑁𝑅𝐸𝑏  3000:
𝑔∗𝑑ℎ 0.5
𝑉𝑏 = (0.546 + 8.7410−6 ∗ 𝑁𝑅𝐸𝐿
12
2. Si 3000  𝑁𝑅𝐸𝑏  8000:

13.59∗𝜇𝐿 0.5
𝑉𝑏 = 0.5 ∗+ 2 + 𝑑 0.5
𝜌𝐿 ∗ 12ℎ
Donde:
𝑔∗𝑑ℎ 0.5
𝛼 = (0.251 + 8.74 × 10−6 ∗ 𝑁𝑅𝐸𝐿 )
12
36
3. Si 𝑁𝑅𝐸𝑏  8000:
𝑔∗𝑑ℎ 0.5
𝑉𝑏 = (0.350 + 8.7410−6 ∗ 𝑁𝑅𝐸𝐿
12
Cuando 𝑁𝑅𝐸𝐿  6000 y 𝑁𝑅𝐸𝑏  32.5 se utiliza la siguiente ecuación para determinar 𝑉𝑏 :
𝑔∗𝑑ℎ 0.5
𝑉𝑏 = 𝐶1 ∗ 𝐶2
12
Donde:
𝐶2 = 1.36+ 𝐶5 + 𝐶6 + 𝐶3 + 𝐶7 + 𝐶3 2
37
𝐶7 = −1.161 × 10−3 + 4.6 × 10−5 ∗ 𝐶4 + 2.954 × 10−3 ∗ 𝐶4 2 + 5.5 × 10−4 ∗ 𝐶4 3 − 6.67 × 10−4* 𝐶4 4
𝐶6 = 0.0413 − 0.0122 ∗ 𝐶4 +0.012 ∗ 𝐶4 2 −1.1 × 10−3 * 𝐶4 3 -1.118 × 10−3 𝐶4 4
𝐶5 = −0.220623 − 0.03408 ∗ 𝐶4 +9.549999 × 10−3 * 𝐶4 2 -8.283001 × 10−3 * 𝐶4 3 +2.645 × 10−3 * 𝐶4 4
𝑁𝑅𝐸𝑏 −5,500
𝐶4 =
1000
38
𝑁𝑅𝐸𝐿 −3000
𝐶3 =
1000
𝐶1 =0.013805+0.4246* 𝐶8 − 0.1753 ∗ 𝐶8 2 + 0.02363 ∗ 𝐶8 3
𝑁𝑅𝐸𝑏
𝐶8 =
10
Y si 𝑁𝑅𝐸𝐿 6000 y 𝑁𝑅𝐸𝑏 > 32.5; 𝐶1 = 0.351 y calculamos 𝐶2
39
El gradiente de fricción se obtiene de la ecuación:
∆𝑝 2 ∗𝜌
𝑓∗ 𝑉𝑚 𝑉𝑠𝐿 +𝑉𝑏
𝐿
= +𝛿
∆𝐿 𝑓 12∗ 2 ∗𝑔∗𝑑ℎ 𝑉𝑚 +𝑉𝑏
El factor de fricción se puede calcular con la ecuación
𝜀 2.514
𝑓 = −2𝑙𝑜𝑔 +
3.71∗𝑑 𝑓 ∗𝑁𝑅𝐸
Y con el numero de Reynolds

124∗𝑑ℎ∗𝑉𝑚 ∗𝜌𝐿
𝑁𝑅𝐸𝐿 =
𝜇𝐿
40
c) Régimen de transición bache – niebla
Orkiszewski adopto el método de interpolación propuesto por Duns y Ros que consiste en
calcular (∆𝑝 /∆𝐿)𝑒 y (∆𝑝 /𝐿∆)𝑓 en las fronteras para flujo bache y flujo niebla.
La zona de transición esta definida por:
𝐿𝑚 > 𝑁𝑔𝑣 > 𝐿𝑠
Donde:
𝐿𝑚 = 84 ∗ 𝑁𝐿𝑣 0,75 + 75
El valor del termino por elevación, esta dado por:

∆𝑝 𝐿𝑚 −𝑁𝑔𝑣 ∆𝑝 𝑁𝑔𝑣 −𝐿𝑠 ∆𝑝
𝑒 = 𝑒𝑏𝑎𝑐ℎ𝑒 +
∆𝐿 𝐿𝑚 −𝐿𝑠 ∆𝐿 𝐿𝑚 −𝐿𝑠 ∆𝐿 𝑒𝑛𝑖𝑒𝑏𝑙𝑎
41
Y el termino por fricción, por:
∆𝑝 𝐿𝑚 −𝑁𝑔𝑣 ∆𝑝 𝑁𝑔𝑣 −𝐿𝑠 ∆𝑝

𝑓 = 𝑓𝑏𝑎𝑐ℎ𝑒 +
∆𝐿 𝐿𝑚 −𝐿𝑠 ∆𝐿 𝐿𝑚 −𝐿𝑠 ∆𝐿 𝑓𝑛𝑖𝑒𝑏𝑙𝑎
Se obtiene un valor mas preciso del factor de fricción de la región de niebla, si el gasto de
has se obtiene con la siguiente ecuación:
𝜌𝐿 −0,25
𝑞𝑔 = 𝐴𝑝 ∗ 𝐿𝑚
𝑔∗𝜎𝐿
42
● d) Régimen de niebla
La región de niebla está definida para:
𝑁𝑔𝑣 > 𝐿𝑚
El gradiente por elevación, dado que el líquido va en suspensión dentro de la corriente

de gas y no existe diferencia de velocidad entre las fases, se calcula:
∆𝑃 1 𝜌𝐿 𝑉𝑠𝐿 +𝜌𝑔𝑉𝑠𝑔
=
∆𝐿 𝑒 144 𝑉𝑚
43
En el gradiente por fricción, se considera que la mayor parte de caídas de presión por fricción se
deben al flujo de gas por la tubería:
2 ∗𝜌
12∗𝑓∗𝑉𝑠𝑔
∆𝑃 𝑔
=
∆𝐿 𝑓 𝜇𝑔
En este caso la rugosidad relativa 𝜀/𝑑 se determina a través de una función del número de Weber
según los lineamientos de Duns y Ros, quienes señalan que sólo será significativo cuando su valor esté
comprendido entre 1 x 10-3 y 0.5. Entre estos límites, ε/d se calcula de la siguiente manera:
Si: 𝑁𝑤 < 0.005

Y si: 𝑁𝑤 ≥ 0.005
𝜀 0.8994∗𝜎𝐿 𝜀 4.6242∗𝜎𝐿 ∗𝑁𝑤 0.302 𝜌𝑔 𝑉𝑠𝑔 ∗𝜇𝐿 2

𝑑
= 2 ∗𝑑
𝜌𝑔 ∗𝑉𝑠𝑔 = 2 ∗𝑑
Donde: 𝑁𝑤 = 0.093
ℎ 𝑑 𝜌𝑔 ∗𝑉𝑠𝑔 ℎ 𝜌𝐿 𝜎𝐿
44
9.- Calcular ∆ℎ de la siguiente ecuación:
∆𝑃 ∆𝑃
+
∆𝑃 ∆ℎ 𝑒 ∆ℎ 𝑓
=
∆ℎ 𝑇 1−𝐸𝑘
Donde:
𝐸𝑘 es el término donde se incluyen las caídas de presión por aceleración
𝑊𝑚 ∗𝑉𝑠𝑔
𝐸𝑘 =
4637∗𝑃𝑝𝑟𝑜𝑚∗𝐴𝑝
45
10. Si la Δh calculada en el paso 9 no es igual a la supuesta en el paso 4, entonces repetir el
procedimiento a partir del paso 4 tomando como supuesta la Δh calculada en el paso 9. Este
procedimiento se repite hasta que Δhsup = Δhcal
11. En caso de que el punto 10 no se cumpla se debe realizar una nueva iteración utilizando
la ∆h obtenida en el paso 9. A Partir de este valor se calcula una nueva p y Z y se realizan
nuevamente todos los pasos.
12. Repetir el procedimiento del paso 3 hasta completar la profundidad total del pozo.
Nota: En caso de que el punto 10 no se cumpla se debe realizar una nueva iteración utilizando la ∆h
obtenida en el paso 9. A partir de este valor se calcula una nueva p y Z y se realizan nuevamente
todos los pasos. Esto hasta que se cumpla la condición del punto 10 46
EJEMPLO
Calcular la caída de presión a través de la TP con los siguientes datos:
47
48
49
50
51
52
53
54
55
56
57
58
59
▪ 11.-En caso de que el punto 10 no se cumpla se

debe realizar una nueva iteración utilizando la ∆h
obtenida en el paso 9. A Partir de este valor se
calcula una nueva p y Z y se realizan nuevamente
todos los pasos.
▪
12. Repetir el procedimiento del paso 3 hasta
completar la profundidad total del pozo.
CONCLUSIONES
Para el trabajo de ingeniería general el método de Griffith y Wallis modificado (Correlación de
Orkiszewski) predecirá caídas de presión con suficiente exactitud y precisión en una amplia gama de
condiciones de pozo. Por lo que se recomienda su uso. Sin embargo, el método debe usarse con
discreción para aquellas condiciones de pozo que no fueron suficientemente evaluadas en su desarrollo
(flujo en espacio anular y patrón de flujo neblina).
Para que un método de predicción sea general, debe expresarse en términos de los patrones de flujo y
de la distribución del líquido.
60
Bibliografía
● Orkiszewski, J. (1967). Predicting Two-Phase Pressure Drops in Vertical Pipe. Journal of Petroleum
Technology, 19(06), 829-838.
● Torres, J., & Trauwvitz E. (2008). Flujo multifásico en tuberías (tesis de pregrado). Universidad
Nacional Autónoma de México, Ciudad de México, México.

Correlacion de Orkiszewski

Cargado por

Copyright:

Formatos disponibles

Correlacion de Orkiszewski

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Correlacion de Orkiszewski

Cargado por

Copyright:

Formatos disponibles

Correlación de Orkiszewski

UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO

 Bernal Martínez Rodrigo

El flujo multifásico se divide en:

Entre las correlaciones para flujo multifásico se

Para flujo horizontal:

Para flujo vertical:

FLUJO MULTIFÁSICO EN TUBERÍAS

a. Diseñar las tuberías de producción y

¿QUÉ HIZO ORKISZEWSKI?

a) Se considera resbalamiento entre las fases

El método de Griffith y Wallis modificado es suficientemente exacto (error promedio) y preciso

CARACTERÍSTICAS DEL FLUJO BURBUJA

CARACTERÍSTICAS DEL FLUJO BACHE

CARACTERÍSTICAS DEL FLUJO BACHE

CARACTERÍSTICAS DEL FLUJO DE TRANSICIÓN

CARACTERÍSTICAS DEL FLUJO NEBLINA

𝑇𝑝𝑟𝑜𝑚 = Temperatura promedio

5.- Determinar las propiedades de los fluidos a P media y a la T.

• Densidad del liquido

● Densidad del gas

5.6142 144 𝑞𝐿 1 𝑊𝑂𝑅

● Velocidad superficial del gas

● Tensión superficial de la mezcla de líquido

● Número de velocidad del gas 𝑁𝑔𝑣

Limites de patrones de flujo

Limites Régimen de Flujo

Con el limite 𝐿𝐵 0.13

𝑑ℎ = diámetro hidráulico de la tubería = área transversal/perímetro mojado

Donde 𝐻𝐿 se calcula con:

El gradiente de fricción se calcula con:

El factor de fricción se calcula utilizando un proceso iterativo

Tabla: Relación de la ecuación a aplicar, con la velocidad de la mezcla y la fase continua.

Las ecuaciones para obtener 𝛿 son:

El procedimiento para calcular la velocidad de bache se inicia suponiendo un valor 𝑉𝑏

2. Si 3000  𝑁𝑅𝐸𝑏  8000:

𝐶6 = 0.0413 − 0.0122 ∗ 𝐶4 +0.012 ∗ 𝐶4 2 −1.1 × 10−3 * 𝐶4 3 -1.118 × 10−3 𝐶4 4

𝐶5 = −0.220623 − 0.03408 ∗ 𝐶4 +9.549999 × 10−3 * 𝐶4 2 -8.283001 × 10−3 * 𝐶4 3 +2.645 × 10−3 * 𝐶4 4

𝐶1 =0.013805+0.4246* 𝐶8 − 0.1753 ∗ 𝐶8 2 + 0.02363 ∗ 𝐶8 3

Y si 𝑁𝑅𝐸𝐿 6000 y 𝑁𝑅𝐸𝑏 > 32.5; 𝐶1 = 0.351 y calculamos 𝐶2

El factor de fricción se puede calcular con la ecuación

Y con el numero de Reynolds

𝐿𝑚 > 𝑁𝑔𝑣 > 𝐿𝑠

El valor del termino por elevación, esta dado por:

∆𝑝 𝐿𝑚 −𝑁𝑔𝑣 ∆𝑝 𝑁𝑔𝑣 −𝐿𝑠 ∆𝑝

La región de niebla está definida para:

El gradiente por elevación, dado que el líquido va en suspensión dentro de la corriente

Si: 𝑁𝑤 < 0.005

𝜀 0.8994∗𝜎𝐿 𝜀 4.6242∗𝜎𝐿 ∗𝑁𝑤 0.302 𝜌𝑔 𝑉𝑠𝑔 ∗𝜇𝐿 2

𝐸𝑘 es el término donde se incluyen las caídas de presión por aceleración

▪ 11.-En caso de que el punto 10 no se cumpla se

También podría gustarte