Formulario Producción I-Flujo Multifásico-2024-2

FORMULARIO DE PRODUCCIÓN I
EFECTOS DE LA COMPLETACIÓN
POR CAÑONEO
∆𝑷𝒄 = 𝑷𝒘𝒇𝒔 − 𝑷𝒘𝒇 = 𝒂𝒒𝟐 + 𝒃𝒒
𝟏 𝟏 𝒓𝒄
𝟐.𝟑𝟎𝒙𝟏𝟎−𝟏𝟒 𝜷𝑩𝒐 𝟐 𝝆𝒐 (𝒓𝒑−𝒓𝒄) 𝝁𝒐 𝑩𝒐 𝑳𝒏 (𝒓𝒑)
𝒂= 𝒃 = 𝟎.𝟎𝟎𝟕𝟎𝟖𝑳
𝑳𝒑 𝟐 𝒑 𝑲𝒑
𝟐.𝟑𝟑𝒙𝟏𝟎𝟏𝟎 𝑸
𝜷= 𝟏.𝟐𝟎𝟏 𝒒=𝑯
𝑲𝒑 𝒑 ∗𝒅𝒆𝒏𝒔𝒊𝒅𝒂𝒅 𝒅𝒆𝒍 𝒕𝒊𝒓𝒐
POR GRAVA
∆𝑷𝒄 = 𝑷𝒘𝒇𝒔 − 𝑷𝒘𝒇 = 𝒂𝒒𝟐 + 𝒃𝒒
𝟗.𝟎𝟖𝒙𝟏𝟎−𝟏𝟑 𝜷𝑩𝒐 𝝆𝒐 𝑳 𝒐 𝒐 𝝁 𝑩 𝑳 𝟏.𝟒𝟕𝒙𝟏𝟎𝟕

𝒂= 𝒃 = 𝟏.𝟏𝟐𝟕𝒙𝟏𝟎 𝜷=
𝑨𝟐 −𝟑 𝑲
𝒈𝑨 𝑲𝒈 𝟎.𝟓𝟓
𝟏𝟒𝟏.𝟓
𝜸𝒐 = 𝟏𝟑𝟏.𝟓+𝑨𝑷𝑰 𝑳 = ∅ 𝑯𝒐𝒚𝒐 − ∅𝑻𝒖𝒃𝒊𝒏𝒈 𝑨 = 𝑨∅ ∗ 𝑻𝒑𝒑 ∗ 𝑯𝒑
FLUJO MULTIFÁSICO DE TUBERÍAS

Consideraciones teóricas del flujo monofásico y multifásico en tuberías
𝟏𝟒𝟖𝟖𝒅𝑽𝝆
𝑵𝑹𝒆 = 𝝁
Flujo laminar
𝟔𝟒𝝁 𝟔𝟒
𝒇 = 𝒇𝒎 = 𝒅𝑽𝝆 = 𝑵 ;
𝑹𝒆
Flujo turbulento lisas
𝒇 = 𝟎. 𝟑𝟏𝟔𝑵𝑹𝒆 −𝟎.𝟐𝟓 ⟶ 𝑵𝑹𝒆 < 𝟏𝟎𝟓
𝒇 = 𝟎. 𝟎𝟎𝟓𝟔 + 𝟎. 𝟓𝑵𝑹𝒆 −𝟎.𝟑𝟐 → 𝟑𝟎𝟎𝟎 < 𝑵𝑹𝒆 < 𝟑𝒙𝟏𝟎𝟔
Flujo turbulento rugosas
−𝟐
𝟐𝝃 𝟏𝟖. 𝟕
𝒇𝒄 = [𝟏. 𝟕𝟒 − 𝟐𝒍𝒐𝒈 ( + )]
𝒅 𝑵𝑹𝒆 √𝒇𝒔
𝝃 𝟐𝟏.𝟐𝟓 −𝟐
𝒇 = [𝟏. 𝟏𝟒 − 𝟐𝒍𝒐𝒈 (𝒅 + 𝟎.𝟗 )]
𝑵𝑹𝒆
𝝃 𝝃
−𝟐
𝟓.𝟎𝟐
𝒅 𝒅 𝟏𝟑
𝒇= [−𝟐𝒍𝒐𝒈 [ 𝟑.𝟕 − 𝑵 𝑳𝒐𝒈 ( 𝟑.𝟕 + 𝑵 )]]
𝑹𝒆 𝑹𝒆
Hagedorn & Brown

∆𝑷 𝟏 𝒇 𝒎 𝝆 𝒇 𝑽𝒎 𝟐
( ) = (𝝆 + )
∆𝑯 𝑻 𝟏𝟒𝟒 𝒎 𝟐𝒈𝒄 𝒅
𝟐
𝝆𝒏 𝟐 [𝝆𝑳 𝝀𝑳 +𝝆𝒈 (𝟏−𝝀𝑳 )]
𝝆𝒎 = 𝝆𝑳 𝑯𝑳 + 𝝆𝒈 (𝟏 − 𝑯𝑳 ) ; 𝝆𝒏 = 𝝆𝑳 𝝀𝑳 + 𝝆𝒈 (𝟏 − 𝝀𝑳 ) ; 𝝆𝒇 = 𝝆𝒎
= 𝝆𝑳 𝑯𝑳 +𝝆𝒈 (𝟏−𝑯𝑳 )
𝑽𝒔𝑳 𝑽𝒔𝒈
𝝀𝑳 = ; 𝝀𝒈 = 𝟏 − 𝝀𝑳 ; 𝝀𝒈 =
𝑽𝒎 𝑽𝒎
𝟎.𝟐𝟐𝟏𝟖𝑽𝒎 𝟐
𝑳𝑩 = 𝟏. 𝟎𝟕𝟏 − 𝒅
𝑽𝒎 = 𝑽𝒔𝒈 + 𝑽𝒔𝑳
𝑽𝒎 𝑽𝒎 𝟐 𝑽𝒔𝒈 𝒇𝒕
𝑯𝑳 = 𝟏 − 𝟎. 𝟓 [𝟏 + − √(𝟏 − ) − 𝟒 ( 𝑽 )] ; 𝑽𝒔 = 𝟎. 𝟖
𝑽𝒔 𝑽𝒔 𝒔 𝒔
𝝆𝑳 = 𝝆𝒐 𝑭𝒐 + 𝝆𝒘 𝑭𝒘 ; 𝝈𝑳 = 𝝈𝒐 𝑭𝒐 + 𝝈𝒘 𝑭𝒘
𝟏
𝝆𝑳 𝟒
𝑵𝑽𝑳 = 𝟏. 𝟗𝟑𝟖𝑽𝒔𝑳 ( )
𝝈𝑳
𝟏
𝝆 𝟒
𝑵𝑮𝑽 = 𝟏. 𝟗𝟑𝟖𝑽𝒔𝒈 (𝝈𝑳 )
𝑳
𝟏
𝝆𝑳 𝟐
𝑵𝑫 = 𝟏𝟐𝟎. 𝟖𝟕𝟐𝒅 ( )
𝝈𝑳
𝟏
𝟏 𝟒
𝑵𝑳 = 𝟎. 𝟏𝟓𝟕𝟐𝟔𝝁𝑳 (𝝈 𝟑𝝆 )
𝑳 𝑳
𝑵𝑳𝑽 𝑷 𝟎.𝟏 𝑪𝑵𝑳 𝑵𝑮𝑽∗𝑵𝑳𝟎.𝟑𝟖𝟎
(𝑵𝑮𝑽𝟎.𝟓𝟕𝟓) (𝟏𝟒.𝟕) ( 𝑵𝑫 ) ;
𝑵𝑫𝟐.𝟏𝟒
𝑯𝑳
𝑯𝑳 = ( ) 𝝍
𝝍
𝟏𝟒𝟖𝟖𝒅𝝆𝒏 𝑽𝒎
𝑵𝑹𝒆 = ; 𝝁𝒎 = 𝝁𝑳 𝑯𝑳 ∗ 𝝁𝒈 (𝟏−𝑯𝑳)
𝝁𝒎
Duns & Ros

∆𝑷 ∆𝑷
(∆𝑯) + (∆𝑯)
∆𝑷 𝑬 𝒇
( ) =
∆𝑯 𝑻 𝟏 − 𝑬𝑲
𝑽𝒔𝒈 𝝆𝒏𝒔 𝑽𝒎
𝑬𝑲 = 𝟏𝟒𝟒𝒈𝒄 𝑷
; 𝝆𝒏𝒔 = 𝝆𝑳 𝝀 + 𝝆𝒈 (𝟏 − 𝝀)
𝟏
𝝆 𝟐
𝑵𝑫 = 𝟏𝟐𝟎. 𝟖𝟕𝟐𝒅 (𝝈𝑳 ) ; 𝑳𝑺 = 𝟓𝟎 + 𝟑𝟔(𝑵𝑳𝑽) ; 𝑳𝑴 = 𝟕𝟓 + 𝟖𝟒(𝑵𝑳𝑽)𝟎.𝟕𝟓
𝑳
𝟏 𝟏
𝝆 𝟒 𝝆 𝟒
Con 𝑵𝑳𝑽 = 𝟏. 𝟗𝟑𝟖𝑽𝒔𝑳 (𝝈𝑳 ) y 𝑵𝑮𝑽 = 𝟏. 𝟗𝟑𝟖𝑽𝒔𝒈 (𝝈𝑳 )
𝑳 𝑳
Patrón burbuja
∆𝑷 ∆𝑷 ∆𝑷
( ) =( ) +( )
∆𝑯 𝑻 ∆𝑯 𝑬 ∆𝑯 𝒇
∆𝑷 𝒎 𝝆
(∆𝑯) = 𝟏𝟒𝟒 ; 𝝆𝒎 = 𝝆𝑳 𝑯𝑳 + 𝝆𝒈 (𝟏 − 𝑯𝑳 )
𝑬
−(𝑽𝒎 − 𝑽𝒔 ) + √(𝑽𝒎 − 𝑽𝒔 )𝟐 + 𝟒𝑽𝒔 𝑽𝒔𝑳

𝑯𝑳 =
𝟐𝑽𝒔
𝑺 𝑵𝑮𝑽
𝑽𝒔 = 𝟏 ; 𝑺 = 𝑭𝟏 + (𝑭𝟐 𝑵𝑳𝑽) + [𝑭′𝟑 ((𝟏+𝑵𝑳𝑽)𝟐 )]
𝝆𝑳 𝟒
𝟏.𝟗𝟑𝟖( )
𝝈𝑳
𝟏 𝟏
𝑭𝟒 𝝆 𝟐 𝟏 𝟒
𝑭′𝟑 = 𝑭𝟑 − 𝑵𝑫 ; 𝑵𝑫 = 𝟏𝟐𝟎. 𝟖𝟕𝟐𝒅 (𝝈𝑳 ) ; 𝑵𝑳 = 𝟎. 𝟏𝟓𝟕𝟐𝟔𝝁𝑳 (𝝈 𝟑𝝆 )
𝑳 𝑳 𝑳
∆𝑷 𝒇𝒎 𝝆𝑳 𝑽𝒎 𝑽𝒔𝑳 𝒇 𝟏𝟒𝟖𝟖𝒅𝝆𝑳 𝑽𝒔𝑳
(∆𝑯) = ; 𝒇𝒎 = 𝒇𝟏 (𝒇𝟐 ) ; 𝑵𝑹𝒆 𝒃 =
𝒇 𝟏𝟒𝟒∗𝟐𝒈𝒄 𝒅 𝟑 𝝁𝑳
𝒇𝟏 =factor de fricción de Moody
𝟐
𝒇𝟏 𝑽𝒔𝒈 𝑵𝑫𝟑 𝑽𝒔𝒈
𝒇𝟐 obtiene de la grafica ; 𝒇𝟑 = 𝟏 + 𝒇𝟏 √𝟓𝟎𝑽
𝑽𝒔𝑳 𝒔𝑳
Patrón Tapón
∆𝑷 ∆𝑷 ∆𝑷
( ) =( ) +( )
∆𝑷 𝒎 𝝆
(∆𝑯) = 𝟏𝟒𝟒 con 𝝆𝒎 = 𝝆𝑳 𝑯𝑳 + 𝝆𝒈 (𝟏 − 𝑯𝑳 )
𝑬
−(𝑽𝒎 − 𝑽𝒔 ) + √(𝑽𝒎 − 𝑽𝒔 )𝟐 + 𝟒𝑽𝒔 𝑽𝒔𝑳

𝑯𝑳 =
𝟐𝑽𝒔
𝑺 (𝟏+𝑭𝟓 )(𝑵𝑮𝑽𝟎.𝟗𝟖𝟐 +𝑭′𝟔 )

𝑽𝒔 = 𝟏 ; 𝑺= ; 𝑭′𝟔 = 𝟎. 𝟎𝟎𝟐𝟗𝑵𝑫 + 𝑭𝟔
(𝟏+𝑭𝟕 𝑵𝑳𝑽)𝟐
𝝆𝑳 𝟒
𝟏.𝟗𝟑𝟖(𝝈 )
𝑳
𝟏 𝟏
𝝆𝑳 𝟐 𝟏 𝟒
𝑵𝑫 = 𝟏𝟐𝟎. 𝟖𝟕𝟐𝒅 ( ) ; 𝑵𝑳 = 𝟎. 𝟏𝟓𝟕𝟐𝟔𝝁𝑳 ( 𝟑 )
𝝈𝑳 𝝈𝑳 𝝆𝑳
∆𝑷 𝒇𝒎 𝝆𝑳 𝑽𝒎 𝑽𝒔𝑳 𝒇 𝟏𝟒𝟖𝟖𝒅𝝆𝑳 𝑽𝒔𝑳
(∆𝑯) = ; 𝒇𝒎 = 𝒇𝟏 (𝒇𝟐 ) ; 𝑵𝑹𝒆 𝒃 =
𝒇 𝟏𝟒𝟒∗𝟐𝒈𝒄 𝒅 𝟑 𝝁𝑳
𝒇𝟏 =factor de fricción de Moody
𝟐
𝒇𝟏 𝑽𝒔𝒈 𝑵𝑫𝟑 𝑽𝒔𝒈
𝒇𝟐 obtiene de la grafica ; 𝒇𝟑 = 𝟏 + 𝒇𝟏 √𝟓𝟎𝑽
𝑽𝒔𝑳 𝒔𝑳
Patrón Neblina
∆𝑷 ∆𝑷
(∆𝑯) + (∆𝑯)
∆𝑷 𝑬 𝒇
( ) =
∆𝑯 𝑻 𝟏𝟒𝟒(𝟏 − 𝑬𝑲)
∆𝑷
( ) = 𝝆𝒎 = 𝝆𝑳 𝝀 + 𝝆𝒈 (𝟏 − 𝝀)
∆𝑯 𝑬
∆𝑷 𝒇𝝆𝒈 𝑽𝒔𝒈 𝟐
( ) =
∆𝑯 𝒇 𝟐𝒈𝒄 𝒅
𝟏𝟒𝟖𝟖𝒅𝝆𝒈 𝑽𝒔𝒈
𝑵𝑹𝒆 =
𝝁𝒈
𝝆𝒈 𝑽𝒔𝒈 𝟐
𝑵𝒘𝒆 = 𝟒𝟓𝟒 ∗ 𝝃 ∗ ( )
𝝈𝑳
𝝃 𝟎.𝟎𝟕𝟒𝟗𝝈𝑳
Si 𝑵𝒘𝒆 ∗ 𝑵𝑳 < 𝟎. 𝟎𝟎𝟓, entonces: =
𝒅 𝒅𝝆𝒈 𝑽𝒔𝒈 𝟐
𝝃 𝟎.𝟑𝟖𝟔𝝈𝑳 (𝑵𝒘𝒆 ∗𝑵𝑳)𝟐
Si 𝑵𝒘𝒆 ∗ 𝑵𝑳 ≥ 𝟎. 𝟎𝟎𝟓, entonces:𝒅 =
𝒅𝝆𝒈 𝑽𝒔𝒈 𝟐
𝝃
3. Si 𝟏𝟎−𝟑 < 𝒅 < 𝟎. 𝟎𝟓, 𝒇 se obtiene de Moody.
𝝃 𝟏 𝝃 𝟏.𝟕𝟑
Si 𝒅 ≥ 𝟎. 𝟎𝟓, 𝒇 = 𝝃 𝟐
+ 𝟎. 𝟐𝟔𝟕 (𝒅)
𝟒𝒍𝒐𝒈(𝟎.𝟐𝟕 )
𝒅
(𝝆𝑳 𝑽𝒔𝑳 +𝝆𝒈 𝑽𝒔𝒈 )
𝑬𝑲 = 𝑽𝒔𝒈
𝟏𝟒𝟒 ∗ 𝟐𝒈𝒄 𝑷
Patrón Transición
∆𝑷 ∆𝑷 ∆𝑷
( ) =( ) +( )
∆𝑷 ∆𝑷 ∆𝑷
( ) = 𝑨( ) + 𝑩( )
∆𝑯 𝑬 ∆𝑯 𝑬𝑻𝒂𝒑ó𝒏 ∆𝑯 𝑬𝑵𝒆𝒃𝒍𝒊𝒏𝒂
𝑳𝑺 −𝑵𝑮𝑽 𝑵𝑮𝑽−𝑳𝑺
𝑨= ; 𝑩=
𝑳𝒎 −𝑳𝑺 𝑳𝒎 −𝑳𝑺
∆𝑷 ∆𝑷 ∆𝑷
( ) = 𝑨( ) + 𝑩( )
∆𝑯 𝒇 ∆𝑯 𝒇𝑻𝒂𝒑ó𝒏 ∆𝑯 𝒇𝑵𝒆𝒃𝒍𝒊𝒏𝒂
Orkiszewski
Patrón de burbuja
∆𝑷 ∆𝑷
(∆𝑯) + (∆𝑯)
∆𝑷 𝑬 𝒇
( ) =
∆𝑯 𝑻 𝟏𝟒𝟒
𝑽𝒔𝒈
< 𝑳𝑩
𝑽𝒎
𝟎.𝟐𝟐𝟏𝟖𝑽𝒎 𝟐
𝑳𝑩 = 𝟏. 𝟎𝟕𝟏 − ; 𝑽𝒎 = 𝑽𝒔𝒈 + 𝑽𝒔𝑳
𝒅
𝟓.𝟔𝟏𝟓(𝒒𝒐 𝒇𝒐 𝑩𝒐 +𝒒𝒘 𝒇𝒘 𝑩𝒘 ) (𝒒𝑳 𝑹𝑮𝑳−𝒒𝒐 𝑹𝒔 )𝑩𝒈

𝑽𝒔𝑳 = ; 𝑽𝒔𝒈 =
𝟖𝟔𝟒𝟎𝟎𝑨𝒕 𝟖𝟔𝟒𝟎𝟎𝑨𝒕
∆𝑷
(∆𝑯) = 𝝆𝒎 ; 𝝆𝒎 = 𝝆𝑳 𝑯𝑳 + 𝝆𝒈 (𝟏 − 𝑯𝑳 )
𝑬
𝑽𝒎 𝑽𝒎 𝟐 𝑽𝒔𝒈
𝑯𝑳 = 𝟏 − 𝟎. 𝟓 [𝟏 + − √(𝟏 + ) − 𝟒 ( 𝑽 )] ; 𝑽𝒔 = 𝟎. 𝟖𝒇𝒕/𝒔
𝑽𝒔 𝑽𝒔 𝒔
𝑽 𝟐
∆𝑷 𝒇𝝆𝑳 ( 𝑯𝒔𝑳 ) 𝟏𝟒𝟖𝟖𝒅𝝆𝑳 𝑽𝒔𝑳
𝑳
( ) = ; 𝑵𝑹𝒆 =
∆𝑯 𝒇 𝟐𝒈𝒄 𝒅 𝝁𝑳 𝑯 𝑳
𝒇 es el factor de fricción de Moody
Patrón de tapón
∆𝑷 ∆𝑷 ∆𝑷
( ) =( ) +( )
𝑽𝒔𝒈
> 𝑳𝑩 y 𝑵𝑮𝑽 < 𝑳𝑺
𝑽𝒎
𝟎.𝟐𝟐𝟏𝟖𝑽𝒎 𝟐
𝑳𝑩 = 𝟏. 𝟎𝟕𝟏 − ; 𝑳𝑺 = 𝟓𝟎 + 𝟑𝟔𝑵𝑳𝑽
𝒅
∆𝑷 𝝆𝒎
( ) =
∆𝑯 𝑬 𝟏𝟒𝟒
𝝆𝑳 (𝑽𝒔𝑳 + 𝑽𝒃 ) + 𝝆𝒈 𝑽𝒔𝒈
𝝆𝒎 = + 𝜹𝝆𝑳
𝑽𝒎 + 𝑽𝒃
𝑽𝒃 = 𝑪𝟏 𝑪𝟐 √𝒈𝒅 ;
𝑽𝒃𝒂 = 𝟎. 𝟓√𝒈𝒅
𝟏𝟒𝟖𝟖𝑽𝒃𝒂 𝒅𝝆𝑳 𝟏𝟒𝟖𝟖𝑽𝒎 𝒅𝝆𝑳
𝑵𝑹𝒆 𝒃 = ; 𝑵𝑹𝒆 𝑳 =
𝝁𝑳 𝝁𝑳
Si 𝑵𝑹𝒆 𝒃 ≤ 𝟑𝟎𝟎𝟎, entonces:

𝑽𝒃𝒄 = (𝟎. 𝟓𝟒𝟔 + 𝟖. 𝟕𝟒𝒙𝟏𝟎−𝟔 𝑵𝑹𝒆 𝑳 )√𝒈𝒅
Si 𝑵𝑹𝒆 𝒃 ≥ 𝟖𝟎𝟎𝟎, entonces:
𝑽𝒃𝒄 = (𝟎. 𝟑𝟓 + 𝟖. 𝟕𝟒𝒙𝟏𝟎−𝟔 𝑵𝑹𝒆 𝑳 )√𝒈𝒅
Si 𝟑𝟎𝟎𝟎 < 𝑵𝑹𝒆 𝒃 < 𝟖𝟎𝟎𝟎, entonces:
𝟏 𝟏𝟑.𝟓𝟗𝝁𝑳
𝑽𝒃𝒄 = 𝟐 [𝜽 + √𝜽𝟐 + ] ; 𝜽 = (𝟎. 𝟐𝟓𝟏𝟖. 𝟕𝟒𝒙𝟏𝟎−𝟔 𝑵𝑹𝒆 𝑳 )√𝒈𝒅
𝝆𝑳 √𝒅
El proceso termina cuando |𝑽𝒃𝒄 − 𝑽𝒃𝒂 | < 𝟎. 𝟎𝟐, si |𝑽𝒃𝒄 − 𝑽𝒃𝒂 | > 𝟎. 𝟎𝟐 se reemplaza 𝑽𝒃𝒄 por
𝑽𝒃𝒂 para obtener un nuevo 𝑵𝑹𝒆 𝒃 , y si sigue siendo 𝑵𝑹𝒆 𝒃 ≤ 𝟑𝟎𝟎𝟎, entonces se toma el valor 𝑽𝒃𝒄
como 𝑽𝒃 .
𝟎. 𝟎𝟏𝟑𝒍𝒐𝒈(𝝁𝑳 )
𝜹= − 𝟎. 𝟔𝟖𝟏 + 𝟎. 𝟐𝟑𝟐𝒍𝒐𝒈(𝑽𝒎 ) − 𝟎. 𝟒𝟐𝟖𝒍𝒐𝒈(𝒅)
𝒅𝟏.𝟑𝟖
𝟎. 𝟎𝟒𝟓𝒍𝒐𝒈(𝝁𝑳 )
𝜹= − 𝟎. 𝟕𝟎𝟗 − 𝟎. 𝟏𝟔𝟐𝒍𝒐𝒈(𝑽𝒎 ) − 𝟎. 𝟖𝟖𝟖𝒍𝒐𝒈(𝒅)
𝒅𝟎.𝟕𝟗𝟗
𝟎.𝟎𝟏𝟐𝟕𝒍𝒐𝒈(𝝁𝑳 +𝟏)
𝜹= − 𝟎. 𝟐𝟖𝟒 + 𝟎. 𝟏𝟔𝟕𝒍𝒐𝒈(𝑽𝒎 ) + 𝟎. 𝟏𝟏𝟑𝒍𝒐𝒈(𝒅)
𝒅𝟏.𝟒𝟏𝟓
𝟎. 𝟎𝟐𝟕𝟒𝒍𝒐𝒈(𝝁𝑳 + 𝟏)
𝜹= − 𝟎. 𝟏𝟔𝟏 + 𝟎. 𝟓𝟔𝟗𝒍𝒐𝒈(𝒅) + 𝑿
𝒅𝟏.𝟑𝟕𝟏
𝟎. 𝟎𝟏𝒍𝒐𝒈(𝝁𝑳 + 𝟏)
𝑿 = −𝒍𝒐𝒈(𝑽𝒎 ) ( ) + 𝟎. 𝟑𝟗𝟕 + 𝟎. 𝟔𝟑𝒍𝒐𝒈(𝒅)
𝒅𝟏.𝟓𝟕𝟏
Si 𝑽𝒎 < 𝟏𝟎, y 𝜹 ≥ −𝟎. 𝟎𝟔𝟓𝑽𝒎 , entonces 𝜹 es el correcto, sino 𝜹 = −𝟎. 𝟎𝟔𝟓𝑽𝒎

𝑽𝒃 𝝆 𝑽𝒃 𝝆
Si 𝑽𝒎 ≥ 𝟏𝟎, y 𝜹 ≥ − 𝑽 +𝑽 (𝟏 − 𝝆𝒎 ), entonces 𝜹 es el correcto, sino 𝜹 = − 𝑽 +𝑽 (𝟏 − 𝝆𝒎)
𝒎 𝒃 𝑳 𝒎 𝒃 𝑳
∆𝑷 𝒇𝝆𝑳 𝑽𝒎 𝟐 𝑽𝒔𝑳 + 𝑽𝒃
( ) = ( + 𝜹)
∆𝑯 𝒇 𝟏𝟒𝟒 ∗ 𝟐𝒈𝒄 𝒅 𝑽𝒎 + 𝑽𝒃
𝒇 es el factor de fricción de Moody
𝟏𝟒𝟖𝟖𝒅𝝆𝑳 𝑽𝒎
𝑵𝑹𝒆 𝒃 =
𝝁𝑳
Patrón de Neblina y Transición: los gradientes de presión se calculan de la misma forma que
en la correlación de Duns & Ros.
Beggs & Brill
∆𝑷 ∆𝑷
∆𝑷 (∆𝑯) + (∆𝑯)
𝑬 𝑭
( ) =
∆𝑯 𝑻 (𝟏 − 𝑬𝑲)
∆𝑷 𝝆𝒕𝒑 𝒔𝒆𝒏(𝜽)
( ) =
∆𝑯 𝑬 𝟏𝟒𝟒
𝝆𝒏𝒔 = 𝝆𝑳 𝝀𝑳 + 𝝆𝒈 (𝟏 − 𝝀𝑳 )
𝑽𝒔𝑳
𝝀𝑳 = 𝑽𝒎
𝝆𝒕𝒑 = 𝝆𝑳 𝑯𝑳 (𝝓) + 𝝆𝒈 (𝟏 − 𝑯𝑳 (𝝓))
𝑳𝟏 = 𝟑𝟏𝟔𝝀𝑳 𝟎.𝟑𝟎𝟐 ; 𝑳𝟐 = 𝟎. 𝟎𝟎𝟎𝟗𝟐𝟓𝟐𝝀𝑳 −𝟐.𝟒𝟔𝟖𝟒
𝑳𝟑 = 𝟎. 𝟏𝟎𝝀𝑳 −𝟏.𝟒𝟓𝟏𝟔 ; 𝑳𝟒 = 𝟎. 𝟓𝝀𝑳 −𝟔.𝟕𝟑𝟖

𝑽𝒎 𝟐
𝑵𝑭𝑹 =
𝒈𝒅
Segregado:
𝝀𝑳 < 𝟎. 𝟎𝟏 y 𝑵𝑭𝑹 < 𝑳𝟏 ó 𝝀𝑳 ≥ 𝟎. 𝟎𝟏 y 𝑵𝑭𝑹 < 𝑳𝟐
𝟎. 𝟗𝟖𝝀𝑳 𝟎.𝟒𝟖𝟒𝟔
𝑯𝑳 (𝟎) =
𝑵𝑭𝑹𝟎.𝟎𝟖𝟔𝟖
Intermitente:
𝟎. 𝟎𝟏 ≤ 𝝀𝑳 < 𝟎. 𝟒 y 𝑳𝟑 < 𝑵𝑭𝑹 ≤ 𝑳𝟏 ó 𝝀𝑳 ≥ 𝟎. 𝟒 y 𝑳𝟑 < 𝑵𝑭𝑹 < 𝑳𝟒
𝟎. 𝟖𝟓𝟒𝝀𝑳 𝟎.𝟓𝟑𝟓𝟏
𝑯𝑳 (𝟎) =
𝑵𝑭𝑹𝟎.𝟎𝟏𝟕𝟑
Distribuido:
𝝀𝑳 < 𝟎. 𝟒 y 𝑵𝑭𝑹 ≥ 𝑳𝟏 ó 𝝀𝑳 ≥ 𝟎. 𝟒 y 𝑵𝑭𝑹 > 𝑳𝟒
𝟏. 𝟎𝟔𝟓𝝀𝑳 𝟎.𝟓𝟖𝟐𝟒
𝑯𝑳 (𝟎) =
𝑵𝑭𝑹𝟎.𝟎𝟔𝟎𝟗
Transición:
𝝀𝑳 ≥ 𝟎. 𝟎𝟏 y 𝑳𝟐 ≤ 𝑵𝑭𝑹 ≤ 𝑳𝟑
𝑯𝑳 (𝟎)𝑻𝒓𝒂𝒏𝒔𝒊𝒄𝒊ó𝒏 = 𝑨𝑯𝑳 (𝟎)𝑺𝒆𝒈𝒓𝒆𝒈𝒂𝒅𝒐 + 𝑩𝑯𝑳 (𝟎)𝑰𝒏𝒕𝒆𝒓𝒎𝒊𝒕𝒆𝒏𝒕𝒆
𝑳𝟑 −𝑵𝑭𝑹
𝑨= ; 𝑩=𝟏−𝑨
𝑳𝟑 −𝑳𝟐
𝑯𝑳 (𝝓) = 𝑯𝑳 (𝟎)𝝍
𝝍 = 𝟏 + 𝑪[𝒔𝒆𝒏(𝟏. 𝟖𝝓) − 𝟎. 𝟑𝟑𝟑𝒔𝒆𝒏𝟑 (𝟏. 𝟖𝝓)]
𝝍 = 𝟏 + 𝟎. 𝟑𝑪 ; Para tubería vertical hacia arriba 𝝓 = 𝟗𝟎°
𝑪 = (𝟏 − 𝝀𝑳 )𝑳𝒏(𝑫𝝀𝑳 𝑬 𝑵𝑳𝑽𝑭 𝑵𝑭𝑹𝑮 )
∆𝑷 𝒇𝒕𝒑 𝝆𝒏𝒔 𝑽𝒎 𝟐
( ) =
∆𝑯 𝑭 𝟏𝟒𝟒 ∗ 𝟐𝒈𝒄 𝒅
𝒇𝒕𝒑 = 𝒇𝒏𝒔 ∗ 𝒆𝑺
−𝟐
𝑵𝑹𝒆
𝒇𝒏𝒔 = [𝟐𝑳𝒐𝒈 ( )]
𝟒. 𝟓𝟐𝟐𝟑𝑳𝒐𝒈(𝑵𝑹𝒆 ) − 𝟑. 𝟖𝟐𝟏𝟓
𝟏𝟒𝟖𝟖𝒅𝝆𝒏𝒔 𝑽𝒎
𝑵𝑹𝒆 =
𝝁𝒏𝒔
𝝆𝒏𝒔 = 𝝆𝑳 𝝀𝑳 + 𝝆𝒈 (𝟏 − 𝝀𝑳 ) ; 𝝁𝒏𝒔 = 𝝁𝑳 𝝀𝑳 + 𝝁𝒈 (𝟏 − 𝝀𝑳 )
𝑿
𝑺=
−𝟎. 𝟎𝟓𝟐𝟑 + 𝟑. 𝟏𝟖𝟐𝑿 − 𝟎. 𝟖𝟕𝟐𝟓𝑿𝟐 + 𝟎. 𝟎𝟏𝟖𝟓𝟑𝑿𝟒
𝝀𝑳
𝑿 = 𝑳𝒏(𝒀) ; 𝒀 = [𝑯 𝟐 ; 𝑺 = 𝑳𝒏(𝟐. 𝟐𝒀 − 𝟏. 𝟐)
𝑳 (𝝓)]
𝝆𝒕𝒑 𝑽𝒔𝒈 𝑽𝒎
𝑬𝑲 =
𝟏𝟒𝟒𝒈𝒄 𝑷

Formulario Producción I-Flujo Multifásico-2024-2

Cargado por

Copyright:

Formatos disponibles

Formulario Producción I-Flujo Multifásico-2024-2

Cargado por

Información del documento

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

Formulario Producción I-Flujo Multifásico-2024-2

Cargado por

Copyright:

Formatos disponibles

FORMULARIO DE PRODUCCIÓN I

∆𝑷𝒄 = 𝑷𝒘𝒇𝒔 − 𝑷𝒘𝒇 = 𝒂𝒒𝟐 + 𝒃𝒒

𝟗.𝟎𝟖𝒙𝟏𝟎−𝟏𝟑 𝜷𝑩𝒐 𝝆𝒐 𝑳 𝒐 𝒐 𝝁 𝑩 𝑳 𝟏.𝟒𝟕𝒙𝟏𝟎𝟕

FLUJO MULTIFÁSICO DE TUBERÍAS

Hagedorn & Brown

Duns & Ros

−(𝑽𝒎 − 𝑽𝒔 ) + √(𝑽𝒎 − 𝑽𝒔 )𝟐 + 𝟒𝑽𝒔 𝑽𝒔𝑳

−(𝑽𝒎 − 𝑽𝒔 ) + √(𝑽𝒎 − 𝑽𝒔 )𝟐 + 𝟒𝑽𝒔 𝑽𝒔𝑳

𝑺 (𝟏+𝑭𝟓 )(𝑵𝑮𝑽𝟎.𝟗𝟖𝟐 +𝑭′𝟔 )

𝟓.𝟔𝟏𝟓(𝒒𝒐 𝒇𝒐 𝑩𝒐 +𝒒𝒘 𝒇𝒘 𝑩𝒘 ) (𝒒𝑳 𝑹𝑮𝑳−𝒒𝒐 𝑹𝒔 )𝑩𝒈

𝒇 es el factor de fricción de Moody

Si 𝑵𝑹𝒆 𝒃 ≤ 𝟑𝟎𝟎𝟎, entonces:

Si 𝑽𝒎 < 𝟏𝟎, y 𝜹 ≥ −𝟎. 𝟎𝟔𝟓𝑽𝒎 , entonces 𝜹 es el correcto, sino 𝜹 = −𝟎. 𝟎𝟔𝟓𝑽𝒎

𝑳𝟏 = 𝟑𝟏𝟔𝝀𝑳 𝟎.𝟑𝟎𝟐 ; 𝑳𝟐 = 𝟎. 𝟎𝟎𝟎𝟗𝟐𝟓𝟐𝝀𝑳 −𝟐.𝟒𝟔𝟖𝟒

𝑳𝟑 = 𝟎. 𝟏𝟎𝝀𝑳 −𝟏.𝟒𝟓𝟏𝟔 ; 𝑳𝟒 = 𝟎. 𝟓𝝀𝑳 −𝟔.𝟕𝟑𝟖

𝝀𝑳 < 𝟎. 𝟎𝟏 y 𝑵𝑭𝑹 < 𝑳𝟏 ó 𝝀𝑳 ≥ 𝟎. 𝟎𝟏 y 𝑵𝑭𝑹 < 𝑳𝟐

𝟎. 𝟎𝟏 ≤ 𝝀𝑳 < 𝟎. 𝟒 y 𝑳𝟑 < 𝑵𝑭𝑹 ≤ 𝑳𝟏 ó 𝝀𝑳 ≥ 𝟎. 𝟒 y 𝑳𝟑 < 𝑵𝑭𝑹 < 𝑳𝟒

𝝀𝑳 < 𝟎. 𝟒 y 𝑵𝑭𝑹 ≥ 𝑳𝟏 ó 𝝀𝑳 ≥ 𝟎. 𝟒 y 𝑵𝑭𝑹 > 𝑳𝟒

𝑯𝑳 (𝟎)𝑻𝒓𝒂𝒏𝒔𝒊𝒄𝒊ó𝒏 = 𝑨𝑯𝑳 (𝟎)𝑺𝒆𝒈𝒓𝒆𝒈𝒂𝒅𝒐 + 𝑩𝑯𝑳 (𝟎)𝑰𝒏𝒕𝒆𝒓𝒎𝒊𝒕𝒆𝒏𝒕𝒆

𝝍 = 𝟏 + 𝑪[𝒔𝒆𝒏(𝟏. 𝟖𝝓) − 𝟎. 𝟑𝟑𝟑𝒔𝒆𝒏𝟑 (𝟏. 𝟖𝝓)]

𝝍 = 𝟏 + 𝟎. 𝟑𝑪 ; Para tubería vertical hacia arriba 𝝓 = 𝟗𝟎°

𝑪 = (𝟏 − 𝝀𝑳 )𝑳𝒏(𝑫𝝀𝑳 𝑬 𝑵𝑳𝑽𝑭 𝑵𝑭𝑹𝑮 )

También podría gustarte