Analisis Nodal
Analisis Nodal
Analisis Nodal
5.1 Introducción
Todas las perdidas de carga son función del caudal de producción y de las
propiedades o característica de los componentes del sistema. En el caso de un flujo en fase
simple, sea este liquido o gas, la perdida de presión puede ser calculada fácilmente, sin
embargo durante la producción de un pozo se tiene un fluido multifasico lo que complica el
calculo de dichas perdidas de presión. De los tres términos que contribuyen a la perdida de
carga total en un sistema total Fig. 5.2 son Componentes de Fricción, Aceleración y
Elevación
La base teórica para la mayoría de las ecuaciones de flujo de fluidos es la ecuación general
de energía, la cual es una expresión del balance o conservación de energía entre dos puntos
en un sistema. Considerando un sistema de estado estable, el balance de energía se puede
escribir como:
mu12 mgZ1 mu22 mgZ 2
U1 p1V1 Q w U 2 p2V2 eq. (5.1)
2 gc gc 2 gc gc
donde:
U1 = Energía interna
PV = Energía de expansión o compresión.
2
mu
= Energía cinética.
2 gc
mgZ
= Energía potencial.
gc
Q = Transferencia de calor.
w = Trabajo desenvuelto por el fluido.
Dividiendo la ecuación 5.1 por m para obtener un balance de energía por unidad de
masa y escribiendo la ecuación resultante en forma diferencial.
p u du g
dU d dZ dQ dw 0 eq. (5.2)
gc gc
La ecuación 5.2 se puede convertir en un balance de energía mecánica usando la
siguiente relación termodinámica:
dp
dh Tds eq. (5.3)
y
p
dU dh d eq. (5.4)
o
dp p
dU Tds d eq. (5.5)
Donde:
h = Entalpía.
s = Entropía.
T = Temperatura.
dQ
ds eq. (5.7)
T
Tds=-dQ + d(lw) eq. (5.8)
Donde
lw = Perdida de trabajo, debido a la irreversibilidad.
dZ d lw dw 0
dp u du g
eq. (5.9)
gc gc
dZ d lw 0
dp u du g
eq. (5.10)
gc gc
Considere una tubería inclinada por donde escurra un fluido en régimen permanente,
conforme se muestra en la figura 5.2, ya que dZ 0 dL sen, la ecuación de energía se
vuelve:
dL sen d lw 0
dp u du g
eq. (5.11)
gc gc
dp u du g d l
sen w 0 eq. (5.12)
dL g c dL g c dL
Considerando la caída de presión positiva en la dirección del flujo, la ecuación 5.12 puede
ser escrita como:
dp g dp u du
sen eq. (5.13)
dL gc dL f gc dL
Donde el gradiente de presión debido al esfuerzo viscoso o las perdidas por atracción
son expresados como:
dp d l
w eq. (5.14)
dL f dL
w 2 w g c
f¡ eq. (5.15)
u / 2 gc
2
u2
Donde
w = Esfuerzo cortante.
u 2
= Energía cinética por unidad de volumen.
2 gc
dp 2 f 'u2
eq. (5.18)
dL f gc d
dp f u 2
eq. (5.19)
dL f 2 g c d
Donde
u = Velocidad de flujo, ft / seg.
gc = Factor de conversión gravitacional 32.17 lbmft / lbf seg2
d = diámetro de la tubería, ft
f = Factor de fricción Moody.
El gradiente de presión total puede ser considerado como la composición de los tres
dp
componentes del sistema. El componente debido al cambio de elevación , debido
dL elevf
dp
a la pérdida por fricción , y el último componente debido a los cambios en la
dL friccion
dp
energía cinética o aceleración , Rescribiendo la ecuación en función al gradiente
dL acele
de presión total se tiene:
dp dp dp dp
Eq. (5.20)
dL dL dL f dL
dp
, es la perdida de presión debido al peso de la columna de fluido
dL elevf
hidrostático esta es cero si se trata de un flujo horizontal. Es la perdida de
presión debido al cambio de la energía potencial entre las posiciones 1 y 2 de
la figura 5.2 la expresión viene dada:
dp g . * sen
Eq. (5.21)
dL elev gc
dp
, es la perdida de presión por fricción, que se origina por el
dL friccion
movimiento del fluido contra las paredes de la tubería de producción y están
en función del diámetro de la tubería y a las propiedades del fluido, como se
muestra en la ecuación siguiente:
dp f u 2
Eq. (5.22)
dL f 2 g c d
dp
, es la caída de presión que resulta de un cambio en la velocidad del
dL acele
fluido entre las posiciones 1 y 2 de la figura 5.2 la ecuación para la caída de
presión por aceleración es la siguiente:
dp vdv
Eq. (5.23)
dL accl gc dL
dp m um dum g fmmvm
m sen Eq. (5.24)
dL g c dL gc 2 gcD
Corregir ec (3.40) de la tesis
du
N Re Eq. (5.25)
Donde
d = Diámetro de la tubería.
= Densidad del fluido, lbm / ft3.
u = Velocidad del fluido, ft / seg.
μ = Viscosidad del fluido, lbm / ft seg.
du
N Re 1488 Eq. (5.26)
En términos de caudal de gas o número de Reynolds para condiciones bases de
temperatura y presión se puede escribir como:
20 q g
N Re Eq. (5.27)
d
Donde
e = Rugosidad absoluta, ft o in.
d = Diámetro interno, ft o in.
Los valores de flujo normalmente se encuentran entre dos extremos: flujo laminar y
flujo turbulento, dentro de estos valores se distinguen cuatro zonas: laminar, crítico,
transición y totalmente turbulento, para cada zona las ecuaciones varían en función del
número de Reynolds y la rugosidad (ver figura 5.4).
dp 32 u
Eq. (5.29)
dL f gc d 2
fu 2 32 u
Eq. (5.30)
2 gc d gc d 2
y
64
f Eq. (5.31)
N Re
Las ecuaciones que se presentan son válidas para valores específicos de número de
Reynolds, como la propuesta por Drew, Koo y McAdams en 1930, utilizado para intervalos
de 3x103 < NRe < 3x106.
1 2e 18.7
1.74 2 log Eq. (5.35)
f d N f e
Re
1 e 21.25
1.44 2 log 0.9 Eq. (5.36)
f d N Re
dp g f u 2 u du
sen Eq. (5.37)
dL g c 2 gc d g c dL
Donde el factor de fricción, f, es una función del número de Reynolds y la rugosidad
de la tubería. Esta relación se muestra en el diagrama de Moody (figura 5.4). La gradiente
de la presión total puede ser considerada para ser compuesta de tres componentes distintos:
dp dp dp dp
Eq. (5.38)
dL dL el dL f dL ace
Donde:
La ecuación 5.37 se aplica para cualquier fluido en estado estable, un flujo dimensional
para el cual f , , y u pueden ser definidos.
1
Production Obtimization Usin Nodal Análisis, H Dale Beggs, 1991, pag 62-64
2
Gas Production Operations, H. Dale Beggs, 1984, pag 100
transiente en las tuberías verticales o inclinadas. Para flujo descendente el seno del ángulo
es negativo, y la presión hidrostática aumenta en dirección del flujo.
Introduciendo una segunda fase dentro de una corriente de flujo complica el análisis
de la ecuación de gradiente de presión. El gradiente de presión es incrementado para la
misma masa de flujo, y el flujo podría desarrollar pulsaciones naturales. Los fluidos pueden
separarse debido a la diferencia de densidades y pueden fluir en velocidades diferentes en la
tubería. Un interfase de separación puede existir entre la fase de líquido y gas.
Analizaremos algunas de las propiedades mas importantes, las cuales deben ser
estudiadas antes de adaptar las ecuaciones de gradiente de presión para las condiciones de
dos fases.
3
Production Optimization Using Nodal Análisis, H. Dale Beggs, 1991, pag 64-68
Volumen de liquido en un elemento de la tubería
HL
Volumen de elemento de tubería
El valor de escurrimiento del líquido varía desde cero para flujo de gas de una sola
fase, a uno para flujo de líquido de una sola fase. El escurrimiento de líquido puede ser
medido experimentalmente por varios métodos, tales como pruebas de resistividad o
capacidad, densitómetros nucleares, o por entrampar un segmento de corriente de flujo
entre válvulas de echado rápido y midiendo el volumen de líquido apañado.
Hg = 1 - HL
qL
L Eq. (5.39)
qL q g
donde
qL = Suma de caudales de petróleo y agua in – situ.
qg = Caudal de flujo de gas in – situ.
qg
g 1 L
qL q g
5.8.1.3 Densidad
Toda ecuación de flujo de fluido requiere que un valor de densidad de fluido este
disponible. La densidad esta envuelta en la evaluación de los cambios de energía debido a
la energía potencial y los cambios de energía cinética. Para calcular los cambios de
densidad con los cambios de presión y temperatura, se necesita tener una ecuación de
estado para el líquido sobre consideraciones. Las ecuaciones de estado son fácilmente
disponibles para fluidos de fase simple.
L O fO w f w Eq. (5.40)
donde
qO
fO Eq. (5.41)
qO qW
y
fw = 1 – fO Eq. (5.42)
S L H L g H g Eq. (5.43)
L 2L g 2g
K Eq. (5.45)
HL Hg
5.8.1.4 Velocidad
El área real por la cual los flujos de gas están reducidos por la presencia de líquido al
A Hg.
qg
vg Eq. (5.47)
AH g
Donde A es el área de tubería.
qL
vL Eq. (5.49)
AH L
Desde que Hg y HL son menores que uno, las velocidades reales son mayores que las
velocidades superficiales.
La velocidad de dos fases o mezcla es calculada basada en los caudales de flujo instantáneo
total de la ecuación:
qL q g
vm vsL vsg Eq. (5.50)
A
vsg vsL
vs vg vL Eq. (5.51)
Hg HL
HL s m
v v vm vs 4vs vsL
2
1/ 2
Eq. (5.53)
2vs
5.8.1.5 Viscosidad
s LH g
Hg
L
Eq. (5.55)
n L H L g H g Eq. (5.56)
n L fO w f w Eq. (5.57)
L O fO w f w Eq. (5.58)
Donde
L = Tensión superficial de petróleo.
w = Tensión superficial de agua.
La ecuación de gradiente de presión, se aplica para cualquier flujo de fluido en una línea
inclinada dado un ángulo de horizontal, dado previamente como:
dp dp dp dp
Eq. (5.59)
dL dL el dL f dL ace
5.8.3.-Modelo simplificado para predecir velocidad mínima del gas para remover líquido
del fondo y la velocidad erosional
Las velocidades críticas encontradas por Turner para el modelo de caída resulta una
ecuación simple de la velocidad hacia arriba ascendente para evitar que las gotas más
grandes caigan, lo cual es considerada como un esfuerzo de corte. Las ecuaciones de Turner
fueron desarrolladas de la siguiente manera:
g CdgApVt
( l g )Vol Eq. (5.60)
gc 2 gc
Donde
Resolviendo
Vl 6.55 * ( L g ) D Eq.(5.61)
Cdg
Vt ^ 2 gD
Nwe Eq.(5.62)
.gc
Donde:
Nwe = numero de weber
σ = tensión superficial lb/pie
Usando la relación anterior, Tarner recomendó que el número de Weber sea igual a 30 para
caídas grandes, modificando la ecuación 5.61.
1.59 ^1 / 4( l g )^1 / 4
Vl Eq. (5.63)
g ^1 / 2
3.06 * Vg * A * P
Qg.. min Eq .(5.66)
TZ
Qg min = caudal mínimo del gas para remover el liquido del fondo MMPCD
Vg =velocidad del gas pie/seg
A =Área conducto o tubería pc
P = Presión psi
Z = Factor se compresibilidad
5.8.4.- Cambio del Componente de elevación
dp g
s sen Eq. (5.67)
dL el g c
dp
f u2
f
Eq. (5.68)
dL f 2 gc d
dp f u2
L L sL Eq. (5.69)
dL f 2 gc d
dp f g g usg2
Eq. (5.70)
dL f 2 gc d
dp f tp f um2
Eq. (5.71)
dL f 2 gc d
dp u du K
Eq. (5.72)
dL ace gc d L
Cuando dos fluidos con diferentes propiedades físicas simultáneamente están en una
tubería, se tiene una gama amplia de posibles modelos de flujo. Que hacen referencia a la
distribución de cada fase en la tubería. Muchos autores ponen su atención en el pronóstico
de modelo de flujo que podría existir para varias condiciones fijadas, a los cuales se tiene
diferentes nombres a cada uno de los modelos.
Govier presento una serie de descripción de los modelos de flujo horizontal aire –
agua, y para flujo vertical aire – agua. Estos son mostrados en las figuras 5.6 y 5.7, las
cuales dependen de la extensión de algunas magnitudes relativas u sL y u sg., cuando el flujo
ocurre en una tubería inclinada con un cierto ángulo el modelo de flujo toma otra forma.
Para flujos ascendentes inclinado en pozos, los modelos son casi siempre tapón o niebla, el
efecto de gravedad en líquido evita la estratificación. Para flujo descendente inclinado el
modelo es normalmente estratificado, niebla o anular.
udu g
dh dZ dq 0 Eq. (5.73)
gc gc
udu g
Jdh dZ Jdq 0 eq. (5.74)
gc gc
dq
El calor adicionado al sistema por unidad de longitud es negativo por tanto que se
dL
pierde calor cuando la temperatura del fluido es mayor que la temperatura que se encuentra
circundante. El gradiente de pérdida de presión de calor puede ser expresado como:
dq U d
T Tg eq. (5.76)
dL wT
Donde:
T = Temperatura media del fluido.
Tg = Temperatura media circundante.
U = Coeficiente de transferencia de calor.
wT = Caudal de flujo másico.
5.10.- Flujo en pozos de Gas
Varios métodos están disponibles para calcular la caída de presión estática y fluyente
en pozos de gas. El método mas utilizado es el de Cullender y Smith, y presión y Temp.
Media para gas seco (sistema monofasico) y Grey para Gas condensado (sistema
multifasico). Todos los métodos inician con la ecuación 5.13, con modificaciones para la
geometría de flujo. La mayoría de los casos, el gradiente de aceleración es ignorada.
Frecuentemente es necesario calcular la presión estática de fondo en un pozo de gas para
cualquier cálculo de pronóstico y capacidad de entrega.
La cual
gM
Pws Pwh EXP Eq. (5.80)
gc R Z T
Esta ecuación sostiene cualquier combinación consistente de unidades. Para las
unidades convencionales de campo,
Pws Pwh EXP 0.01875 g H / T Z Eq. (5.81)
Donde
Pws = Presión estática de fondo pozo, psia.
Pwh = Presión en cabeza de pozo, psia.
g = Gravedad del gas (aire = 1)
H = Profundidad del pozo, ft.
T = Temperatura media en tubería, ºR..
Z = Factor de compresibilidad, evaluado a las
Condiciones medias de presión y temperatura
P ( Pws Pwh ) / 2 .
Pws Pwh 1 2.5 105 H Eq. (5.82)
Este método considera las variaciones del factor de compresibilidad con los cambios de
presión y temperatura, así mismo como las alteraciones de temperatura con la profundidad.
De la ecuación 5.79
H
TZ M
p dp R 0 0.01875 g dH Eq. (5.83)
ZT
I
p
Integrando el lado derecho de la ecuación se tiene:
Pws
Idp 0.01875
Pwh
g H Eq. (5.84)
Pws
2 Idp Pm Pwh I ms I wh Pws Pms I ws I ms eq. (5.85)
Pwh
Donde:
Pwh = Presión a la profundidad cero, H = 0, psia.
Pms = Presión a la profundidad media, H / 2, psia.
Pws = Presión a profundidad total, H, psia.
Iwh = Integral evaluada a Pwh y Tsup.
Ims = Integral evaluada a Pms y Tmed.
Iws = Profundidad a la cual se quiere evaluar, Th, ft.
0.01875 g H
Pws Pwh Eq. (5.86)
I ms I ts
0.01875 g H
Pws Pms Eq. (5.87)
I ms I ws
Conociendo la temperatura de superficie y de fondo de pozo se puede determinar la
temperatura a cualquier profundidad utilizando la relación:
T f Ts
Th Ts h Eq. (5.88)
H
Donde
Th = Temperatura a cualquier profundidad, ºR.
Tf = Temperatura final, ºR.
Ts = Temperatura en superficie, ºR.
H = Profundidad final, pie.
h = Profundidad a la cual quiere evaluarse Th, pies.
Paso B
1. Determinar el valor inicial de presión estática, Pms, aplicando la ecuación 5.82.
Previamente sustituir Pwh por Pms y H por H / 2 en lugar de H se debe que el
análisis del paso B parte de la profundidad final del paso A, que dice, dar la
profundidad media del pozo.
2. Obtener Z a Pms y Tf.
3. Encontrar Iws.
4. Recalcular Pws como los resultados de los puntos 2 y 3, si Pws Pms el
procedimiento finaliza obteniendo de esta manera la presión estática de fondo
procurada, si no repetir hasta encontrar dos Pws próximas.
El pozo es puesto en producción, con una presión de fondo fluyente que impulsa el
fluido hasta superficie. Si la energía es suficiente como para vencer la resistencia
encontrada durante la trayectoria de flujo dice se que el pozo fluye por surgencia natural; en
caso de ser insuficiente el pozo es sometido a métodos de recuperación o levantamiento
artificial y obtener la presión adecuada para elevar el fluido.
dp g f u2
cos Eq. (5.89)
dL g c 2 gc d
dp pM f u2
(cos Eq. (5.90)
dL ZRT 2 gc d
La integración de la ecuación 5.90 asumiendo una temperatura media en la sarta de
flujo y evaluando Z a las condiciones medias de presión y temperatura,
Para una mejor comprensión en la aplicación del método, se mostrará una aplicación
práctica.
La derivación del método de Cullender y Smith para pozos fluyentes comienza con la
ecuación 5.90. Las siguientes substituciones son hechas para la velocidad:
q
u
A
Psc T Z
q qsc
Tsc P Zsc
La cual da:
dp pM cos MTZ Psc2 f qsc2
dL ZRT R p Tsc2 2 gc d A2
o
p dp M p
2
cos C
ZT dh R ZT
Donde
8 psc2 qsc2 f
C
Tsc2 g c 2 d 5
ZT MD
donde
0.667 f qsc2
F2 Eq. (5.94)
d5
y
TVD
cos Eq. (5.95)
MD
Para abreviar la escrita ecuación 5.93 y dividiendo el comportamiento total del pozo
en dos secciones, H / 2.
Para encima de H / 2 la ecuación se simplifica a la siguiente expresión:
0,10796 q
F d 4,277 in. Eq. (5.99)
d 2,612
0,10337 q
F d 4,277 in. Eq. (5.100)
d 2,582
Aplicando la regla de Simpson para obtener un resultado de presión más exacto, esta es:
18,75 g ( MD) 2
Pwf Ptf Itf 4 Im f Iwf Eq. (5.101)
3
1. Determinar el valor del lado izquierdo de la ecuación 6.82 para encima de H /2.
2. Determinar F2 de la ecuación 5.99 o 5.100.
3. determinar Itf de la ecuación 5.98 y las condiciones de cabeza del pozo.
4. Asumir Imf = Itf para condiciones medias de profundidad del pozo o el punto medio del
tuno de producción.
5. Determinar Pmf de la ecuación 5.96.
6. Utilizando el valor de Pwf determinado en el paso 5 y la temperatura media aritmética
Tmf determinar el valor de Imf de la ecuación 5.98.
7. Recalcular Pmf de la ecuación 5.96. Si este valor recalculado no es menos de 1 psi de
Pmf calculado en el paso 5, repita los pasos 6 y 7 hasta que el criterio anterior este
satisfecho.
8. Asuma Iwf = Imf para las condiciones de fondo del tubo de producción.
9. Repita los pasos 5 hasta el 7, utilizando la ecuación 5.97 para debajo de H /2 del tubo
de producción, y obtener un valor de la presión de fondo de pozo, Pwf.
10. Aplicar la regla de Simpson como expresa la ecuación 5.101, obtenemos un valor de
presión de fondo de pozo mas exacto.
Para un mejor entendimiento mostraremos una aplicación práctica con los mismos
datos del ejemplo 2 para ver la diferencia entre métodos.
5.10.2.3.- Método de Grey
De la ecuación 5.89 tenemos la ecuación básica de energía para una línea de producción
vertical
dp g f u2
cos Eq. (5.89)
dL g c 2 gc d
Pasos a seguir para la determinación de las pérdidas de presión por fricción y hidrostática
qg qsc Bg
Velocidad superficial del gas (pie/seg.) Vsg u sg Eq.(5.102)
A A
5.6144 qo Bo
Velocidad superficial del Cond.(pie/seg.) Vsc uso Eq.(5.103)
A
Fracción de líquido L qc qw
L eq. (5.106)
(qc qw) qg
qc
Fracción de Condensado c c eq. (5.108)
(qc qw)
Fracción de Agua w w 1 c eq.(5.109)
ns ^ 2Vm ^ 4
N1 eq. (5.113)
g ( l g )
gD ^ 2( l g )
N2 eq. (5.114)
1 730 Rv
N3 0.08141 0.0554 ln eq. (5.115)
Rv 1
Vsl
Relación de Velocidades Rv RV eq.(5.116)
Vsg
H l 2.314N11 ^ N 3
Fracción de Holup Hl 205 eq. (5.117)
N2
En base a las ecuaciones presentadas las perdidas de presión por fricción y hidrostática son:
mgdL
Phh eq. (5.123)
144 gc
2 ftpGm ^ 2dL
Pfric eq. (5.124)
144 gcDns
La línea que conecta la boca del pozo con el separador adquiere mayor importancia
cuando esta tiene una longitud considerable donde la caída de presión es relevante para
determinar la capacidad productiva del pozo.
Los factores básicos envueltos en un flujo horizontal son los mismos que se aplicarán
en flujo vertical. En ambos sistemas la caída de presión es la suma de las perdidas por
fricción y aceleración. La principal diferencia radica en las consideraciones del balance de
energía debido a la posición de la tubería.
dp f u
2
p M f u2
Eq. (5.125)
dL 2 gC d Z R T 2 gC d
0.5
CT P12 P22
q b d
2.5
eq. (5.127)
Pb g f T Z L
el valor de C depende de la unidad utilizada en la ecuación. En la tabla 5.3 están los valores
de C para varias combinaciones de unidades.
Un factor llamado eficiencia es usado algunas veces en la ecuación 5.127 para justificar de
echo que las tuberías frecuentemente entregan menos gas que el calculado. El factor de
eficiencia normalmente varía entre 0.7 y 0.92 y es normalmente obtenido de la experiencia.
TABLA N° 5.3
Valores de C para Diversas unidades
P T D L q C
psia ºR in. Mi scfd 77,54
psia ºR in. Ft scfd 5634
psia ºR in ft MMscfd 5.634x10-3
psia ºR in M M3/d 1,149x106
Utilizando estas relaciones para el factor de fricción en la ecuación 5.127, la forma general
de la ecuación de flujo de las líneas de surgencia es:
a4
Tb P12 P22 1 a5
2a a 3