Produccio
Produccio
Produccio
1
controlada por las características de estos componentes y subcomponentes. Teniendo en
cuenta la interacción existente entre todos y cada uno de ellos, puesto que cualquier
cambio o alteración de las condiciones de flujo en alguno de ellos afectará en mayor o
menor grado las condiciones de flujo de los restantes, se puede inferir que todo sistema
de producción debe ser tratado y manejado de manera integral. De esta forma, su diseño
final estará basado tanto en los comportamientos esperados del yacimiento y del flujo
vertical y horizontal, como en los subcomponentes agregados tales como reductores,
válvulas, codos, etc. Así, los criterios de selección de las características, tamaño,
diámetro, etc., de los elementos del sistema deberán estar fundamentados en un análisis
físico riguroso, aunque generalmente aproximado, del sistema de flujo como un todo, de
manera integral, en correspondencia con las expectativas de potencial de producción y de
tasas de producción por pozo adecuadas para drenar el yacimiento de manera eficiente.
En la figura I.2 se presenta un diagrama ilustrativo del perfil de presión a través de toda la
trayectoria del flujo en un sistema de producción.
2
COMPRESORES
TANQUE DE ALMACENAMIENTO
SEPARADOR
MANIFOLD DE PRODUCCION
REGISTRADOR DE PRESIONES
3
1 - COMPORTAMIENTO DE INFLUJO
El componente más importante y de mayor relevancia económica en la industria de
los hidrocarburos es el yacimiento, fuente generadora del producto comercial. En torno a
su comportamiento giran los demás componentes de este complejo sistema industrial. De
sus características y condiciones dependen el éxito y la continuidad operativa de su
explotación comercial.
5
k dp
v (1.1)
dx
o, en término de tasa volumétrica de flujo
k A dp
qvA (1.2)
dx
Para flujo lineal, la ec. 1.2 puede ser integrada entre los puntos de entrada y salida
del flujo, a través de la longitud L, (fig. 1.1).
P2 k dP 1 L
P1
A 0
q dx (1.3)
q
Ck A
P1 P2 (1.5)
BL
En la ec. 1.4 se puede observar que un gráfico de P vs. L produciría una línea recta
de pendiente q B / k A . Esto es, la presión varía linealmente con la distancia.
Si el fluido es gas, entonces la tasa de flujo varía con la presión. Usando el hecho
de que la tasa de flujo de la masa, q, debe ser constante y expresando la densidad como
función de presión, temperatura y gravedad específica del gas, la ec. 1.4 se transforma
en:
8.93 Z T g L
P12 P22 q sc (1.6)
kg A
donde,
7
a.- Para petróleo
13
o Bo L 9.08 10 Bo2 o L
P1 P2 3
qo 2
qo2 (1.7)
1.127 10 ko A A
donde,
= Coeficiente de velocidad. (pies-1)
8.93 g Z T L 1247
. x10 16 Z T L
P1 P2 q sc
2 2 g
2
q sc 2 (18
. )
kg A A
donde,
2.329 x 1010 k 1.2 arenas consolidad as
(1.9)
1.47 x 10 7 k 0.55 arenas no consolidad as
Aunque en los yacimientos muy raramente existen condiciones de flujo lineal, las
ecuaciones precedentes podrán ser usadas en completaciones de pozos con empaques
de grava, con resultados satisfactorios.
8
2.- Flujo radial
La Ley de Darcy puede ser usada para calcular el flujo hacia un pozo, donde los
fluidos convergen radialmente. En este caso, el área abierta al flujo no es constante y, por
lo tanto, debe ser incluida en la integración de la ec. 1.2. El área abierta al flujo para
cualquier radio es A 2 k h ( ver fig. 1.2 ). Así, la ec. 1.2 se transforma en:
2 r h k p
q (1.10)
r
2 r h ko dP
qo (1.11)
o Bo dr
9
Pe ko re dr
2 h dP qo (1.12)
Pwf oBo rw r
0.00708 k o h P R Pwf
qo (1.14)
o Bo ln 0.472 re rw
La ec. 1.12 puede ser integrada entre cualquier radio arbitrario y el límite exterior del
área de drenaje, resultando la siguiente expresión en términos de presión en el radio
arbitrario:
141.2qo o Bo
P PR ln re 0.75 141.2qo o Bo ln(r ) (1.15)
ko h ko h
10
Un gráfico de P vs. Ln( r ) estará representado por una línea recta de pendiente:
141.2 qo o Bo
m (1.16)
ko h
5000
PR
4000
141.2 qo o Bo
Pwf
P wf1 m
ko h
3000 P wf2
P wf3
2000
0.1 1 10 100 1000 10000
r
Fig. 1.3 Distribución de Presiones en un Yacimiento
Esta pendiente variará cada vez que las condiciones de flujo sean alteradas. Es
decir, sí alguno de los factores del lado derecho de la ec. 1.15 es alterado; por ejemplo, la
tasa de producción, la pendiente tendrá un valor diferente. En la fig. 1.3 se muestra un
ejemplo ilustrativo de la distribución de presiones en un yacimiento para diferentes
condiciones de flujo.
11
q P Tsc 2 r h k g dp
qsc (1.19)
sc Psc T Z g dr
2 2 q sc g Z T Psc ln re rw
Pe Pwf (1.21)
k g hTsc
Modificando esta ecuación para flujo estabilizado y reemplazando las variables Psc y
Tsc (presión y temperatura a condiciones normales) por sus valores correspondientes de
14.7 lpca y 520 R, la ecuación resultante para la tasa de flujo de gas, expresada en
unidades de campo, será:
q sc
70310
6
2
k g h PR Pwf
2
(1.22)
g Z T ln 0.472re rw
La ec. 1.20 puede ser integrada entre cualquier radio arbitrario y el limite exterior del
área de drenaje, resultando la siguiente expresión:
q sc g Z T q sc g Z T
P2 PR
2
6
ln re 0.75 6
ln( r ) (1.23)
70310 kg h 70310 kg h
Un gráfico de P² vs. ln(r) estará representado por una línea recta, cuya pendiente es
dada por:
1422 q sc g Z T
m (1.24)
kg h
12
II.- RESTRICCION DEL FLUJO Y CONCEPTO DE DAÑO.
141.2qo o Bo
ln 0.472re rw
P R Pwf (1.25)
ko h
13
El resultado de todo esto es que la distribución de presiones para un pozo real
puede diferir notablemente de los cálculos para un pozo ideal. Esta diferencia es mayor
en las cercanías del pozo del pozo y ocurre en detrimento de su productividad.
PR
Pwf*
Pwf
rw ra re
r
Fig 1.4 Distribución de presiones
ó, en términos de Ps,
* 141.2 qo o Bo
PS Pwf Pwf S (1.27)
ko h
Combinando las ecuaciones 1.25 y 1.27 para expresar la caída de presión real en
términos del factor de daño, resulta,
14
141.2 qo o Bo
P R Pwf ln re rw 0.75 S (1.28)
ko h
Una ecuación similar para flujo de gas puede ser desarrollada fácilmente, y es la
siguiente:
q sc g Z T
2
P R Pwf
2
6
ln re rw 0.75 S (1.29)
70310 kg h
El factor de daño da una indicación del carácter del flujo en las vecindades del pozo,
con relación a un pozo ideal. Cuando es positivo indica daño o restricción al flujo. Si es
negativo indica estimulación o mejoramiento de las condiciones de flujo.
*
P R Pwf
Ef (1.31)
P R Pwf
ln re rw 0.75
Ef (1.32)
ln re rw 0.75 S
Generalmente, el término ln re rw está dentro del rango 6.5 – 8.5. Usando un
7
Ef (1.33)
7S
15
Otra expresión usada para cuantificar las condiciones de flujo es mediante la figura
de radio aparente del pozo, rwa, definido como:
rwa rw e S (1.34)
El cual, cuando se introduce en la ec. 1.28 sustituyendo a rw, ésta se transforma en :
141.2 qo o Bo
P R Pwf ln re rwa 0.75 (1.35)
ko h
Una ilustración del concepto de radio aparente del pozo se presenta en la figura 1.5.
Un pozo dañado se manifiesta por tener un radio aparente menor que su radio real.
Lo contrario es indicativo de que el pozo está estimulado. En este caso, el radio aparente
podría llegar a aproximarse al radio de drenaje.
16
Turbulencia y puede ser expresado analíticamente mediante la ecuación:
S S a D qo (1.36)
(2)
Craft y Hawkins propusieron la siguiente expresión analítica para calcular el daño
de la formación debido a la alteración de la permeabilidad en las vecindades del pozo,
excluyendo el efecto de turbulencia.
S a k k a 1 ln ra rw (1.37)
EFECTO CUANTIFICABLE
Estado S EF Ps rwa
Dañado >0 <1 >0 < rw
No alterado =0 =1 =0 = rw
Estimulado <0 >1 <0 > rw
Tabla 1.2. Características de pozos
Otro aspecto importante que conviene discutir es el referente a la restricción del flujo
causada por la penetración parcial o incompleta del pozo en la formación productora.
17
Desde los inicios de las operaciones de producción en campos petrolíferos, los
pozos han sido perforados o completados a través de solamente una fracción del espesor
neto petrolífero, a objeto de evitar o minimizar los efectos de conificación de agua y/o gas
y, en algunos casos, debido a dificultades operacionales para controlar la circulación del
fluido de perforación o completación en la zona objetivo. Obviamente, esta práctica,
conocida como penetración parcial, restringe el caudal de entrada de fluídos en el pozo.
Sin embargo, no debe verse como una operación en detrimento de la capacidad de
producción del pozo, puesto que más bien es un mecanismo de control para mejorar la
eficiencia de producción y prolongar la vida activa del pozo.
Varios trabajos relacionados a este tema han sido publicados a través de los años,
siendo Muskat (3) el pionero en presentarlo en forma analítica para cuantificar su efecto en
el comportamiento de producción de pozos. Sin embargo, aquí solamente se hará
referencia a 2 (dos) de ellos, ambos basados en el trabajo de Muskat, dada la simplicidad
analítica de sus ecuaciones y la aceptabilidad que han tenido en el campo de la
ingeniería de producción. Estos métodos se refieren al de Brons y Martins y el de Odeh.
En 1.961, Brons y Martins (4) sugirieron que el efecto de la penetración parcial puede
ser expresado como un factor de daño, adicional al hasta ahora discutido. Ellos
propusieron la siguiente ecuación:
18
(5)
En 1.980, Odeh publicó un artículo técnico sobre este tema. En él se formula una
relación empírica para calcular el efecto de daño debido a un intervalo parcial abierto al
flujo, hp, localizado arbitrariamente.
S c 1.351 b 1
0. 825
lnrw hD 7 1.95 0.49 0.1 lnrw hD lnrwc (1.41)
donde,
zm y hp / 2 (1.43)
y = Distancia desde el tope de la formación al tope del intervalo abierto
(pies)
qo
J (1.44)
P R Pwf
0.00708 k o h
J (1.45)
o Bo ln re rw 0.75 S
Arreglando la ec. 1.44 para Pwf P R qo J es notorio que un gráfico de Pwf vs. qo
19
Si las condiciones son tales que J permanece constante con cambios en la caída de
presión, una vez obtenido su valor de una prueba de producción representativa o
mediante la ec. 1.44, este valor podrá ser usado para predecir el comportamiento de
influjo bajo otras condiciones. Sin embargo, esto implica que la función de presión
f(p)= ko/oBo permanecerá constante, lo cual ocurre muy raramente. En la práctica, esta
opción solamente es usada en caso de flujo bajo condiciones de subsaturación. Así,
conviene expresar ocasionalmente el índice de productividad como:
0.00708 h PR ko
J
P R Pwf lnre rw 0.75 S wf o Bo dP
P
(1.46)
6
70310 kg h
Jg (1.48)
g Z T ln re rw 0.75 S
20
SISTEMA X SIST EMA X
re 0.966A1/ 2
1
rw rw
2
0.571A1/ 2 1.44A1/ 2
1
rw rw
2
0.604A1/ 2 1
1.925A1/ 2
rw rw
4
0.61A1/ 2 1
6.59A1/ 2
600
rw rw
4
0.678A1/ 2 1
9.36A1/ 2
1
3
rw rw
4
1
1
0.668A1/ 2 1.724A1/ 2
rw rw
2
1
1.368A1/ 2 1.794A1/ 2
1
4
rw rw
2
0.884A1/ 2 9.523A1/ 2
1
rw rw
2
1.485A1/ 2 10.135A1/ 2
rw rw
21
IV.- COMPORTAMIENTO DE IPR’S EN POZOS DE PETROLEO.
Un simple análisis de la ec. 1.28 permite inferir que varios parámetros dependientes
de presión afectan el flujo de petróleo hacia el pozo. Si todas estas variables pudieran ser
calculadas en cualquier momento de la vida productiva de un yacimiento, las ecuaciones
resultantes de la integración de la ecuación de Darcy podrían ser usadas para cuantificar
el IPR. Desafortunadamente, rara vez existe suficiente información que permita estos
cálculos con un grado de confiabilidad aceptable; por lo tanto, el Ingeniero de Producción
se apoya en métodos empíricos para llevar a cabo este proceso.
En esta sección se presentan los métodos empíricos mas usados para determinar el
comportamiento IPR en pozos productores. Algunos de ellos requieren de, al menos, una
prueba de producción en condiciones estables. Otros requieren varias pruebas con
medidas de Pwf y qo. También se presenta un procedimiento para estimar el IPR cuando
no se dispone de pruebas estabilizadas.
(7)
Vogel presentó un modelo empírico para calcular el comportamiento IPR de
pozos productores de petróleo en yacimientos saturados. El método fue basado en
resultados de modelos aplicados a yacimientos hipotéticos para diferentes tipos y
saturaciones de crudo y diferentes espaciados de pozos. Un total de 21 (veintiún)
condiciones de yacimientos fueron analizados.
observó que el comportamiento gráfico de estos dos parámetros era similar en todos los
casos estudiados, en base a lo cual desarrolló su ecuación empírica, dada a continuación,
mediante un modelo de regresión numérico.
22
2
qo Pwf Pwf
1 0.2 0.8 (1.49)
qo ( max) PR PR
En esta ecuación, las presiones son dadas en valores medidos, lpcm. En la fig. 1.7
se presenta un ejemplo ilustrativo del comportamiento IPR de Vogel para un yacimiento
bajo diferentes condiciones de flujo.
4000
3000
Pwf
2000
1000
0
0 1000 2000 3000 4000
qo
Fig. 1.7 Comportamiento IPR (según Vogel)
expresada como:
2
qo Pwf P
1.8 1 0.8 1 wf (1.51)
qo ( max) PR PR
23
Muchas veces será necesario calcular la presión de fondo fluyente que satisfaga
una tasa de producción deseada. En este caso, la ec. 1.50 se expresa en términos de ,
resultando:
0. 5
5.0 qo
1.125 0.5 5.0625
(1.52)
qo ( max)
y
Pwf P R 1 (1.53)
q L / q L( max) , donde q L qo q w .
que la presión fluyente sea mayor o igual que la presión de burbujeo, Pwf Pb , ó que
24
considerando que el Indice de Productividad es constante para Pwf Pb y suponiendo
que la ecuación de Vogel aplique para Pwf Pb . También se supone que la curva de IPR
es continua, esto es, las pendientes de los dos segmentos, por debajo y sobre el punto de
burbujeo, son idénticas a Pwf Pb . Ver gráfico ilustrativo en la fig. 1.8.
PR
Pb
Pwf
qb q o (m ax)
qo
Fig. 1.8 IPR para un yacimiento sub-saturado
qo q b
b 1.8 0.8 b (1.55)
qo ( max) qb
ó, en términos de qo,
Derivando esta expresión con respecto a Pwf se obtiene la pendiente recíproca, que
relaciona el cambio en la tasa de producción con respecto al cambio de Pwf.
25
Evaluando esta pendiente en el punto de burbujeo, a Pwf Pb , resulta:
J Pb
qo ( max) qb (1.60)
1.8
qo qb J Pb b 1 b (1.61)
2.25
Una vez calculado el valor de J a Pwf Pb , la ec. 1.61 puede usarse para generar
la curva IPR.
qb J P R Pb (1.62)
Muchas veces será necesario calcular la presión de fondo fluyente requerida para
obtener una tasa de producción deseada. En este caso, la ec. 1.61 se expresa en
términos de b, resultando:
9.0 qo qb
0.5
b 1.125 0.5 5.0625 (1.64)
J Pb
y,
Pwf Pb 1 b (1.65)
*
P P R Pwf q J* J
FE ideal * (1.66)
Preal P R Pwf q J J
*
P R Pwf PS P R Pwf
FE (1.67)
P R Pwf P R Pwf
27
Usando la ec. previa de eficiencia de flujo, y definiendo * 1 Pwf
*
P R , la
ecuación de Vogel puede ser expresada como:
q
q o
FE 1
* 1.8 0.8 * (1.68)
o ( max)
1
) = Tasa de producción (a Pwf 0 ) que se obtendría sí FE=1 ó S=0.
*
donde qoFE( max
Para relacionar los parámetros Pwf, Pwf* y FE se usa la ec. 1.67; de manera que,
Pwf P R FE P R Pwf
*
(1.69)
* FE (1.70)
Así, con las ecs. 1.68 y 1.70 es posible construir curvas de IPR adimensional para
diferentes valores de FE, como se muestra en la fig. 1.9. Para ello se calcula * de la
ec. 1.70 para cada valor de 1 Pwf P R y el valor estimado de FE. Luego se usa el
0.8
0.6
Pwf 2.0
1.8
PR 0.4
1.6
1.4
0.2 1.2
28
Combinando las ecuaciones 1.68 y 1.70 se obtiene la ecuación del comportamiento
IPR como función de la tasa máxima de producción para FE=1. En términos de tasa de
producción, esta ecuación es:
) FE 1.8 0.8 FE
1
qo qoFE( max (1.71)
La tasa máxima de producción para las condiciones reales del pozo ocurriría a =1;
es decir, para Pwf = 0. Sustituyendo esta condición en la ec. 1.71, resulta:
) FE 1.8 0.8 FE
1
qo( max) qoFE( max (1.72)
qo 1.8 0.8 FE
(1.73)
qo ( max) 1.8 0.8 FE
Una simple observación a esta ecuación permite inferir que su aplicabilidad está
condicionada a ciertos valores máximos de y FE; en consecuencia, la variable qo(max) no
representa necesariamente la tasa de producción máxima real para las condiciones de
flujo dadas, sino la tasa máxima de aplicabilidad de la ecuación. En términos generales,
estas condiciones limitantes son:
1.8 0.8 FE 0
o sea, FE 2.25
2 - La tasa de producción adimensional tiene que ser menor o igual que 1.0, puesto
que no pueden existir tasas de producción mayor que la máxima. Analíticamente,
esta ecuación es dada por:
29
1.8 0.8 FE
1.0
1.8 0.8 FE
2.25 2.25
2 1 0
FE FE
2.25
1 1 0
FE
o sea,
2.25
1 0 ; o sea,
FE
2.25
1
FE
2.25
Así, Pwf P R 2
FE
Las ecuaciones 1.71 1.73 aplican para yacimientos saturados, Pb P R , bajo las
30
qo q b 1.8 0.8 FE b
b (1.74)
qo ( max) qb 1.8 0.8 FE
qo ( max) qb
qo qb b 1.8 0.8 b FE (1.75)
1.8 0 .8 FE
Derivando respecto a Pwf, de igual manera como se hizo en la ec. 1.56, se obtiene:
d qo 1 d qo Pwf
J , con b 1
d Pwf Pb d b Pb
J Pb
qo ( max) qb 1.8 0.8 FE (1.76)
1.8
qo
J (1.78)
PR Pb Pb b 1 b FE 2.25
31
Para Pwf Pb , J debe ser calculado de las ecs. 1.44 ó 1.62
en:
qo qb J Pb 1 FE 2.25 (1.79)
qo q b 1 FE b / 2.25
b (1.80)
qo ( max) qb 1 FE / 2.25
(i) FE 2.25
2.25 2.25
(ii ) b 1 Pwf Pb 2
FE FE
32
2.25 1 qo 2 2 qo1
FE (1.81)
1 2 qo 2 2 2 qo1
(9)
Fetkovich propuso un método para calcular el comportamiento de influjo en
pozos petrolíferos usando el mismo tipo de ecuación que ha sido usada por mucho tiempo
para pozos de gas. Su validez fue comprobada con los resultados de numerosos análisis
de pruebas isocronales y de reflujo en yacimientos de variadas características: saturados
y subsaturados a condiciones iniciales, parcialmente agotados con saturación de gas por
encima de la crítica y permeabilidades en el rango 6 – 1000 md. En todos los casos se
encontró que el comportamiento de las curvas de presiones con tasas de producción
seguía la misma forma que la usada para expresar el comportamiento de influjo para
pozos de gas. Esto es:
2
qo C P R Pwf
2 n
(1.82)
donde,
C es el coeficiente de flujo.
n es el exponente dependiente de las características del pozo.
P R y Pwf vienen expresadas en lpca.
33
Como puede observarse, la ec. 1.82 tiene dos variables desconocidas, C y n; por lo
tanto, se requieren al menos dos pruebas de producción estabilizadas para poder aplicar
el método.
línea recta de pendiente igual a 1/n. Una vez conocido el valor de n es posible calcular el
valor de C mediante la aplicación de la ec. 1.82 en cualquier punto arbitrario de la recta.
En la fig. 1.10 se presenta un ejemplo ilustrativo del método.
qo
Fig. 1.10 Gráfico de Fetkovich
log qo1 log C n log P R Pwf
2 2
1
log qo3 log C n log P R Pwf
2 2
3
34
log qo1 log qo 2
n
log P R Pwf 2
2
log P
1
R
2
Pwf 2
2
(1.83)
qok
C
P R
2
Pwf 2
k
(1.84)
35
q4
q3
q2
qo
q1
q=0
PR
Pw f1
Pw f2
Pwf
Pw f3
Pw f4
tiempo
Fig. 1.11 Prueba de Reflujo Convencional
36
q3
qo q2
q1 q4 perm anente
q=0
PR
Pw f4
Pwf
Pw f1 Pw f2
Pw f3
tiempo
Fig. 1.12 Prueba de Reflujo isocronal norm al
graficados vs. qo en papel log-log, obteniendo así el valor de n, dado por la pendiente de
la línea recta resultante del gráfico. Para determinar el valor de C es conveniente efectuar
una prueba de producción en condiciones de flujo estabilizado. En la fig. 1.12 se muestra
un esquemático ilustrativo de una prueba isocronal.
En las pruebas isocronales pudiera darse el caso de que en los periodos de cierre
del pozo, entre cada período de producción constante, el tiempo de restauración de P R
sea excesivamente grande, lo que incrementaría considerablemente el costo del servicio.
Para evitar este problema, es frecuente usar pruebas isocronales modificadas para cerrar
el pozo, entre cada etapa de pruebas, por un tiempo igual al período de producción.
Obviamente, P R puede no ser alcanzado, pero es posible manejar el análisis mediante
un gráfico de log Pwsi Pwfi
2 2
vs. log q , el cual producirá una línea recta. Los P
o wsi, se
refieren al ultimo valor medido de presión estática en cada período de cierre, como se
ilustra en la fig. 1.13.
37
q3
qo
q2
q1 q4 perm anente
q=0
PR Pw S1
Pw S2 Pw S3
Pw f1 Pw f4
Pwf
Pw f2
Pw f3
tiempo
Fig. 1.13 Prueba de Reflujo isocronal m odificada
o sea,
qo ( max)
C 2n
(1.87)
PR
P 2 n
qo
1
wf
(1.88)
qo ( max) PR
Esta ecuación también presenta dos variables desconocidas, qo(max) y n, por lo que
se requieren dos pruebas de producción estabilizadas, por lo menos, para su aplicación.
38
El método de Fetkovich, usando cualquiera de las dos ecuaciones propuestas, 1.82
ó 1.88, para el análisis de comportamiento IPR en pozos productores de petróleo, ofrece
mayor confiabilidad que los métodos de Vogel y Standing por cuanto incorpora
intrínsecamente el efecto de altas velocidades de flujo (turbulencia, flujo no Darcy) a
través de la inclusión del exponente n.
Se puede apreciar fácilmente que las ecuaciones 1.49 (Vogel) y 1.88 (Fetkovich)
son parecidas en su forma para n=1. En este caso la ec. 1.88 se transforma en:
2
qo Pwf
1 (1.89)
qo ( max) PR
Como podrá observarse, esta ecuación, al igual que la ecuación de Vogel, contiene
solamente una variable desconocida, qo(max), y por lo tanto, solamente se requiere una
prueba de producción para su aplicación. Vale decir que los resultados obtenidos de la
aplicación de esta ecuación son generalmente mas conservadores que los de Vogel. Por
otra parte, es mas sencilla y más fácil de usar. En consecuencia, es recomendable su uso
antes que el método de Vogel cuando solamente se disponga de una prueba de
producción.
Para aplicar la ec. 1.89 en yacimientos subsaturados, será necesario proceder como
en el caso de Vogel, a aplicar la ecuación por debajo del punto de burbujeo mediante un
desplazamiento del eje de las ordenadas de manera que intersecte la curva IPR en el
punto de burbujeo. Esto es:
2
q o qb Pwf
1 (1.90)
qo ( max) qb Pb
ó, en términos de qo,
Pwf
2
qo qb qo ( max) qb 1 (1.91)
Pb
39
Evaluando esta pendiente en el punto de burbujeo, a Pwf = Pb, resulta:
2 qo ( max) qb
J (1.93)
Pb
ó, en términos de producción:
J Pb
qo ( max) qb (1.94)
2
J Pb Pwf
2
q o qb 1 (1.95)
2 Pb
ecuación siguiente:
qo
J (1.96)
P
2
Pwf
PR Pb b 1
2 Pb
2 qo
J (1.97)
Pwf
2
PR 1
Pb
Fetkovich propone agregar el exponente ―n‖ a las ecuaciones 1.96 y 1.97 cuando se
disponga de varias pruebas de producción. Esto es, respectivamente, en términos de qo:
40
n
J Pb Pwf
2
qo J P R Pb 1 (1.98)
2 Pb
n
J P R Pwf
2
qo 1 (1.99)
2 P R
141.2 qo o Bo
P R Pwf ln re rw 0.75 S a D qo (1.100)
ko h
1.6289 10 16 k o Bo o
D (1.101)
h rw o
Definiendo,
141.2 o Bo
A ln re rw 0.75 S a (1.102)
ko h
y,
141.2 o Bo
B D (1.103)
ko h
41
La ec. 1.100 puede ser expresada como:
2
P R Pwf A qo B qo (1.104)
P R Pwf
A B qo (1.105)
qo
42
Como puede observarse, la ecuación de Jones presenta 2 (dos) variables
desconocidas, A y B, por lo tanto, al igual que en el método de Fetkovich, se requieren de
por lo menos dos pruebas estabilizadas para su aplicabilidad. Resulta obvio que un
gráfico cartesiano de P R Pwf qo vs. qo generará una línea recta de pendiente igual a
La ec. 1.104 puede ser resuelta para expresarla en términos de qo, resultando:
qo
A A 2 4 B P R Pwf 0.5 (1.106)
2B
P R Pwf
qo B
qo
Fig. 1.14 Gráfico de Jones y asoc.
La tasa máxima de producción puede ser obtenida de la ec. 1.106, haciendo Pwf = 0.
Esto es,
qo ( max)
A A 2 4 B PR 0.5 (1.107)
2B
43
V.- COMPORTAMIENTO IPR EN POZOS DE GAS.
2
q sc C P R Pwf
2 n
(1.108)
Esta ecuación no es otra que la propuesta por Fetkovich, años mas tarde, para
pozos de petróleo.
aplicada. A altas velocidades, el término a tiene muy poca relevancia y puede ser omitido
de la ecuación. En estos casos, la caída de presión es proporcional al cuadrado de la
velocidad del fluido (análogo al flujo turbulento en tuberías). En flujo de transición, entre
baja y alta velocidad, será necesario tomar en cuenta ambos términos, ya que la omisión
de alguno de ellos podría conllevar a errores importantes en los cálculos del
(11)
comportamiento de influjo. En 1.952, Cornell y Katz propusieron expresar la ecuación
de Furchheimer en función de las propiedades del fluido fluyente y de la roca de
44
formación y proporcionaron la forma necesaria para aplicar la ecuación como una
herramienta predictiva. Los autores concluyeron que la constante a es definida por la ley
de Darcy (a = /k) y la constante b depende de la densidad del fluido y de una constante
empírica , usualmente conocida como coeficiente de velocidad. Así,
dP
2 (1.110)
dr k
d P C1 g Z T q sc C 2 g Z T q sc
2
(1.111)
dr kg r h P r 2 h2 P
Pe C1 g Z T q sc re dr C 2 g Z T q sc
2
re dr
Pwf
P dP
kg h rw r h2
rw r2
(1.112)
Integrando resulta:
C1 g Z T q sc r C 2 g Z T q sc
2
1 1
Pe2 Pwf2 ln e (1.113)
kg h rw h2 rw re
45
Para fluido pseudo-contínuo ó estabilizado, como generalmente se considera,
haciendo 1/rw – 1/re 1/rw, e introduciendo el factor de daño debido a la alteración de la
permeabilidad en el término correspondiente al flujo tipo Darcy, la ec. 1.113 se transforma
en:
2
P R Pwf2 A q sc B q sc2 (1.114)
Con,
1422 g Z T ln re rw 0.75 S a
A (1.115)
kg h
y
46
VI.- PREDICCION DEL COMPORTAMIENTO IPR.
A medida que la presión del yacimiento declina por efecto del agotamiento, el
caudal de flujo de fluidos en la formación también declinará. Así, resulta evidente que el
comportamiento de influjo de un pozo será diferente para cada cambio de condiciones del
yacimiento.
(8)
Standing publicó un procedimiento, basado en la ecuación de Vogel, que
puede ser usada para predecir los cambios de la tasa de producción máxima, qo(max), con
los cambios en las condiciones del yacimiento.
qo Pwf P
1 1.8 0.8 wf (1.117)
qo ( max) PR PR
resulta:
qo ( max) P
J 1 0.8 wf (1.118)
PR PR
47
Pwf
1 0.8
J P P RP qo ( max) P PR P
(1.119)
J F P RF q o ( max) F Pwf
1 0.8
PR F
P RF J F
qo ( max) F qo ( max) P (1.120)
P RP J P
Asumiendo que las condiciones de daño (factor S) y el radio de drenaje del pozo
permanecen constantes, la relación JF/JP puede ser expresada, de la ec. 1.45, como:
ko
JF
o o F
B
(1.121)
JP ko
o Bo P
ko
P RF
o o F
B
qo ( max) F qo ( max) P (1.122)
P RP ko
o Bo P
Una vez conocido el nuevo valor de qo(max) se puede aplicar la ec. De Vogel para
determinar el nuevo comportamiento IPR del pozo. Debe quedar claro que este
procedimiento requiere del conocimiento del comportamiento futuro del yacimiento.
(9)
Fetkovich propone el uso de su método para predecir el comportamiento IPR en
un pozo de petróleo mediante dos suposiciones básicas: 1.- el exponente n no cambia
con el tiempo y 2.- el coeficiente de flujo C es una función lineal de P R. Por lo tanto,
48
P RF
CF CP (1.123)
P RP
gas, los únicos parámetros que cambian con presión son la viscosidad (g) y el factor de
compresibilidad (Z).
Aplicando los mismos artificios y suposiciones que para el caso de petróleo, los
coeficientes A y B, de la ec. 1.114, y C, de la ec. 1.108, pueden ser calculadas para las
nuevas condiciones mediante las siguientes expresiones:
Z
AF AP
g
Z
F
(1.124)
g P
ZF
BF BP (1.125)
ZP
Z
CF CP
g
Z
F
(1.126)
g P
49
VII.- EJERCICIOS RESUELTOS.
Un pozo productor de un yacimiento con una presión promedio de 2.500 lpc fue
probado con una tasa de producción de 354 BPPD en condiciones estabilizadas. La
presión de fondo fluyente fue medida en 2.100 lpc. La presión de saturación es 2.600.
aplicando el método de Vogel, construir la curva IPR y calcular:
a) El potencial máximo de producción [qo(max)]
b) La tasa de producción esperada si la presión fluyente se logra reducir a 2000 lpc
c) La presión de fondo fluyente necesaria para obtener una tasa de producción de
700 BPD.
Solución.
a.- La ec. 1.50 de Vogel puede ser escrita como
qo
qo ( max)
1.8 0.8
De la prueba de producción, qo= 354 BPPD y Pwf = 2100 lpc, resulta:
Pw f
1 016
. qo (max) 1323 BFD
PR
3000
Pw f qL
2500
0 1 1323 2000
500 0.8 1228
Pwf
1500
1000 0.6 1048
1000
1500 0.4 783
2000 0.2 434 500
2500 0.0 0 0
0 250 500 750 1000 1250 1500
qL
50
b.- De la curva de IPR o aplicando la ecuación de Vogel se obtiene una tasa de
producción de 434 BPPD para una presión de fondo fluyente de 2000 lpc.
c.- Aplicando la ec. 1.52 para una tasa de producción de 700 BFD
5.0 700
0.5
1.125 0.5 5.0625 0.3477
1323
Pwf P R 1
qo 354
a.- q o ( max) 2
2
Pwf 2100
1 1
PR 2500
3000
Pwf qL
2500
0 1202 2000
500 1154
Pwf
1500
1000 1010
1000
1500 770
500
2000 433
2500 0 0
0 250 500 750 1000 1250 1500
qL
51
b.- De la curva IPR ó aplicando la ec. 1.89 se obtiene una tasa de producción de
433 BFD para una presión de fondo fluyente de 2000 lpc.
qo
Pw f PR 1
q o (max)
700
Pw f 2500 1
1202
Pw f 1616 lpc
Se desea :
1.- Calcular la tasa de producción correspondiente al punto de burbujeo.
2.- Construir la curva IPR.
3.- Determinar la tasa de producción correspondiente a una presión de fondo
fluyente de 1550 lpc
4.- Determinar la presión fluyente requerida para obtener una tasa de producción
de 780 BFD.
5.- Determinar la presión de fondo fluyente requerida para obtener una tasa de
producción de 2000 BPD.
52
qo 536
J 2.68 bpd / lpc
PR Pw f 2000 1800
qb J P R Pb
qb 804 BFD
Pwf b qL
2000 - 0
1800 - 536
1700 0 804
1500 0,1176 1312
1000 0,4118 2337
500 0,7059 3011
0 1,0 3335
Puede observarse que por encima de la presión de burbujeo la IPR es lineal, puesto
que J es constante. Esto es, en el intervalo P R Pwf Pb la curva IPR se basa en la ec.
1.44. Por debajo de la presión de burbujeo la curva IPR se construye mediante la ec. 1.61,
con:
Pw f
b 1 2500
Pb 2000
1500
Pwf
1000
500
0
0 1000 2000 3000 4000
qL
53
(c) – Para calcular qL a presiones fluyentes por debajo del punto de burbujeo (1550
lpc < Pb) se aplica la ec. 1.61, tal como se construyó la curva IPR anterior. Entonces:
Pwf
b 1 Pb
1 1550
1700
b 0.0882
q L qb J Pb b 1 b
2.25
Pwf P R q L
J
(e) – La tasa requerida de 2000 BFD es mayor que qb; por lo tanto, será necesario
aplicar las ecs. 1.64 y 1.65 para determinar la presión de fondo fluyente
correspondiente.
9.0 q L qb
0.5
b 1.125 0.5 5.0625
J Pb
b = 0.3034
Pwf Pb 1 b
54
Ejercicio Nº 1.4
Usando la siguiente información:
qL = 800 BFD
Pwf = 2100 lpc
P R = 2600 lpc
Pb = 2350 lpc
Construir la curva de comportamiento IPR
Pwf b qL
2600 - 0
Pb 2350 0 410 Calculado de la ec. 1.62
Pwf 2100 0.1064 800 Usando la ec. 1.61
1500 0.3617 1579
1000 0.5745 2058
500 0.7872 2381
0 1 2550
55
La tasa máxima ( 2550 BFD ) puede ser verificada mediante la ec. 1.60
J Pb
q L ( max) qb
1.8
1.639 2350
q L ( max) 410
1.8
3000
2500 Pb
2000 Prueba
Pwf
1500
1000
500
0
0 500 1000 1500 2000 2500 3000
qL
Ejercicio Nº 1.5
Resolver el problema anterior usando la ecuación normalizada de Fetkovich para
n=1
Solución
56
qo
J
P
2
Pwf
PR Pb b 1
2 Pb
800
J
2350 2100
2
2600 2350 1
2 2350
La IPR por encima del punto de burbujeo será calculada con J constante, mediante
la ec. 1.44. La tasa de producción en el punto de burbujeo es dada por la ec. 1.62 y
la curva IPR por debajo de Pb es obtenida mediante la ec. 1.95
3000
qL
57
q L qb
q L ( max) qb 2
Pwf
1
Pb
800 411
q L ( max) 411 2
2100
1
2350
qL(max) = 2342 BFD
Ejercicio Nº 1.6
Solución
Se requiere determinar el comportamiento IPR para las tres condiciones dadas:
dañado, no alterado y estimulado. La eficiencia de flujo para un factor de daño de
+3 es aproximadamente 0.7 (calculada mediante la ec. 1.33).
58
Pwf P R P R Pwf FE
Pwf P R FE P R Pwf
FE = 0.7 FE = 1 FE = 1.3
Pwf 1600 (dada) 1849 1983
Como las presiones de fondo fluyentes son menores que la presión del punto de
burbujeo para los tres casos, será necesario calcular el Indice de Productividad
constante ( por encima del punto de burbujeo ) mediante la ec. 1.78. Las tasas de
producción en el punto de burbujeo se obtienen de la ec. 1.62.
FE = 0.7 FE = 1 FE = 1.3
Pwf 1600 1849 1983
FE = 0.7 FE = 1 FE = 1.3
Pwf b qL qL qL
2430 - 0 0 0
2100 0 104 145 186
1983 0.0557 140 195 250
1849 0.1195 180 250 318
1600 0.2381 250 342 429
1300 0.3810 326 438 538
1000 0.5238 394 517 619
700 0.6667 454 579 672
400 0.8095 505 624 - (*)
0 1 560 659 - (*)
59
(*) No aplica por condiciones de límite b 2.25 FE 1
3000
2500
Pb = 2100 Lpc
Pwf ( Lpc )
2000
1600 Lpc
1500
EF
1000
=1
EF
.3
=1
.0
500
250
429
EF
342
=0
.7
0
0 100 200 300 400 500 600 700 800
qL ( BPD )
Ejercicio Nº 1.7
Solución
Aplicando la ec. 1.70 para ambas condiciones de flujo
FE1
FE1 1 FE 2 2 2 1
FE 2
60
FE = 0.7 FE = 1.3
0.15348 0.08264
IPR’s
Pwf qL qL
2085 0 0 0
1912.7 0.08264 111 202
1800 0.13669 181 323
1765 0.15348 202 359
1600 0.23261 298 517
1300 0.37650 459 756
1000 0.52038 603 934
700 0.66427 728 1051
561.3 0.73079 780 1084 LIMITE DE APLICACION
2500
2000
Pwf (Lpc)
1500
EF
1000 =1
.3
EF
=1
500
.0
0
0 200 400 600 800 1000 1200
qL (BPD)
61
Ejercicio Nº 1.8
Solución
La ec. 1.82, de Fetkovich, expresada en término de presiones es :
1 1
q n 1 n 1
P P L qL n
2 2
R wf
C C
Un gráfico de P R Pwf2 vs. qL en papel log-log resulta en una recta de pendiente igual a
2
log P R Pwf2
2
1
n
1
log qo log
n
1
C
n xy x y
m
n x 2 x
2
62
b
y m x
n
Resolviendo:
y= x=
log P R Pwf2
2
log (qL) xy x2
1 .E + 0 7
1 4 70.7263 10.5058 26.8903
m
n 4 27.6742 10.5058 10.5058 PR 2 - Pw f2
1
m 1.2342 n = 0.8102
n
1 1 26.8903 1.2342 10.5058
b log
n C 4
1 .E + 0 6
b = 3.4810 100 1000
QL
q L 0.001512 P R Pwf2
2
0.8102
63
Construcción de la IPR
Pwf qL
4000
3600 0
3500
3000 335 3000
(Lpc)
2500 514 2500 IP
2000
R
2000 649
Pwf
1500
1500 750 1000
Ejercicio Nº 1.9
Resolver el ejercicio anterior aplicando el método de Jones y asoc.
Solución
P R Pwf
La ecuación de Jones: A B q L puede ser ajustada por mínimos
qL
cuadrados, de la misma forma que la ecuación de Fetkovich.
y= x=
P R Pwf q L qL xy x2
qL
1.1813 1.3955 7.504 3600 Pwf 10 3 0. 5
3.752 10 3
Construcción de la IPR
Pwf qL 3
3600 0 2.5
(PR-Pwf)/qL
3000 332
2
2500 513
1.5
2000 661
1
1500 789
1000 904 0.5
500 1009 0
0 1106 0 200 400 600 800
qL
4000
3000
Pwf (Lpc)
2000
IP
R
1000
0
0 250 500 750 1000 1250
qL (BPD)
65
REFERENCIAS
1. - Muskat, M.: ― The Flow of Homogeneus Fluid Through Porous Media ―.
Publicado por McGraw Hill Book Company, 1937.
2. - Craft, B.C. y Hawking, M.F.: ― Applied Petroleum Reservoir Engineering ―.
Publicado por Prentice-Hall Book Co., 1959.
3. - Muskat, M.: ― Partially Penetrating Wells in Isotropic Formations; Potential
Distribution‖. Physics, 1932.
4. - Brons, F. y Martin, V.E.: ―The Effect of Restricted Fluid Entry on Well
Productivity ―. Actas del AIME, 1961.
5. - Odeh, A.S.: ― An Equation for Calculating Skin Factor Due to Restricted Entry ―.
JPT. Junio, 1980.
6. - Odeh, A.S.: ― Pseudo Steady-State Flow Equation and Productivity Index for a
Well with Non-circular Drainage Area ―. Mobil Res. and Dev. Corp.
7. - Vogel, J.V.: ― Inflow Performance Relationships for Solution Gas Drive Wells ―.
JPT. Enero, 1968.
8. - Standing, M.B.: ―Inflow Performance Relationships for Solution Gas Drive Wells
―. JPT. Noviembre, 1970.
9. - Fetkovich, M.J.: ― The Isochronal Testing of Oil Wells ―. SPE No. 4529, 1973.
10. - Jones, L.G., Blount, E.M. y Glaze, O.H.: ― Use of Short Term Multiple Rate
Flow Test to Predict Performance of Wells Having Turbulence ―. SPE No. 6133,
1976.
11. - Cornell, D. y Katz, D.L.: ― Flow of Gases through Consolidated Porous Media ―.
Ind. and Eng. Chem. Octubre, 1953.
BIBLIOGRAFIA
1 - Golan, M. y Whitson, C.H.: ― Well Performance ― Publicado por Prentice-Hall
Book Company. 1991.
2 - Beggs, M.D.: ― Production Optimization – Using Nodal Analysis ―. Publicado por
OGCI Publications. Tulsa, 1987.
3 – Nind,T.E.W.: ― Principles of Oil Well Production ―. Publicado por McGraw-Hill,
Inc. 1964.
66
2 - FLUJO EN TUBERIAS Y RESTRICCIONES
En este capitulo se discutirán las técnicas para calcular las pérdidas de presión
en las partes del sistema de producción correspondientes al flujo en el pozo (flujo
vertical) y al flujo en la superficie ( flujo horizontal), incluyendo las pérdidas en los
sub-componentes de esas partes del sistema, como: estranguladores de flujo,
válvulas de seguridad, etc.
67
Para diseñar y analizar un sistema de producción para flujo multifásico es
necesario entender claramente el fenómeno físico, así como conocer las bases
teóricas y las ecuaciones correspondientes a los diferentes métodos de cálculos
existentes. Todos estos métodos son empíricos y están basados en datos reales de
campo, experimentos de laboratorio o una combinación de ambos. Sus aplicaciones
requieren del conocimiento de ciertos parámetros físico-químicos y termodinámicos
involucrados en sus ecuaciones, tales como las propiedades PVT de los fluidos,
tensión superficial, masa fluyente, gradiente dinámico de temperatura, etc. En el
anexo A se presentan varias correlaciones empíricas para calcular algunos de estos
parámetros.
68
I.- ECUACION DE ENERGIA (FLUJO CONTINUO).
mv12 m g Z1 mv 22 m g Z2
U1 P1 V1 q w U 2 P2 V2 (2.1)
2 gc gc 2 gc gc
69
Donde q y w son la energía calórica añadida al fluido y el trabajo,
respectivamente.
dP g d dP
(2.2)
dL g c g c dL dL f
donde,
dP
= Gradiente de presión total,
dL
= Densidad del fluido responsable de la energía potencial,
d
= Energía cinética generada por el movimiento, y
g c dL
dP
= Gradiente de presión debido a las pérdidas por fricción o
dL f
fuerzas viscosas.
(1)
Beggs presenta una expresión analítica, desarrollada por Fanning, para las
pérdidas debidas al efecto de fricción. Esta expresión relaciona el esfuerzo cortante
de la pared de la tubería ( w ) y la energía cinética por unidad de volumen ( v /
2
70
El esfuerzo cortante de una tubería de diámetro t viene dado por:
t dP
w (2.4)
4 dL f
dP 2 f' 2
(2.5)
dL f g c t
dP f 2
(2.6)
dL f 2 g c t
t2 g c dP
(2.7)
32 dL f
Resultando en:
64
f (2.8)
t
t
N RE (2.9)
71
En unidades de campo,
t
N RE 120.6667 (2.10)
QSC g
N RE 2.009 *10 4 (2.10a)
t g
donde,
72
Varios estudios experimentales sobre perfiles de velocidad y gradientes de
presión, para flujo en tuberías bajo condiciones de turbulencia, han demostrado que
el efecto de la fricción depende de las condiciones de las paredes expuestas al flujo
con relación a su diámetro, además del Número de Reynolds. Varias ecuaciones han
sido desarrolladas para tuberías lisas, cada una de ellas válida para cierto rango de
Número de Reynolds, siendo la más usada la presentada por Drew, Koo y McAdams,
aplicable con muy buena aproximación para números de Reynolds en el rango de
3000 < NRE < 3x106. Esta ecuación es:
0. 32
f 0.0056 0.5 N RE (2.12)
1 2
1.74 2 log (2.13)
f t
73
Esta ecuación es:
1 2 18.7
1.74 2 log (2.14)
f t N RE f
2
1
f (2.15)
2 18.7
1.74 2 log
t N RE f
(2)
En 1.976 Jain presentó una correlación para el factor de fricción expresado
explícitamente, cuyos resultados comparan satisfactoriamente con los de Colebrook.
Jain demostró que para valores de rugosidad relativa en un rango entre 10-6 y 10-2 y
Números de Reynolds en un rango entre 5x103 y 108 el error obtenido es de 1.0 %
en relación a los resultados de la ec. de Colebrook. Su ecuación produce un error
máximo de 3 % para números de Reynolds del orden de 2000. Esta ecuación es:
1 21.25
1.14 2 log 0.9 (2.16)
f t N RE
74
Dada la excelente aplicabilidad de la ecuación de Jain, la carta universal de
factor de fricción presentada en la figura 2.1 está basada en dicha ecuación.
75
II.- VARIABLES DE FLUJO.
76
de un elemento volumétrico de tubería que es ocupado por líquido en cualquier
instante. Esto es,
Hg + HL = 1
Esto es,
qL
L (2.17)
qL qg
77
En términos de la fase gaseosa,
qg
g 1 L (2.18)
qL qg
qL
v SL (2.19)
A
qg
v Sg (2.20)
A
qL
vL (2.21)
A HL
qg
vg (2.22)
A Hg
qL qg
vm SL Sg (2.23)
A
78
Muchas veces es conveniente determinar el grado de deslizamiento y calcular
el factor de entrampamiento en función de la velocidad de deslizamiento, s, definida
como la diferencia entre la velocidad superficial del gas y la del líquido.
SL
L (2.24)
m
Por definición:
Sg SL
S g L (2.25)
Hg HL
S m S m 2 4 g SL
0. 5
HL (2.26)
2 S
qg 6.5 10 5 qo RGP Rs Bg
Sg (2.27)
A A
79
Con
0.00504 Z T 460
Bg (2.29)
P 14.7
Donde,
Sg = Velocidad superficial del gas, (pies/seg)
SL = Velocidad superficial del liquido, (pies/seg)
qO = Tasa de producción de petróleo, (BN/día)
RGP = Relación gas-petróleo de producción, (PCN / BN)
RS = Relación gas-petróleo en solución, (PCN / BN)
Bg = Factor volumétrico del gas, (Bls / PCN)
BO = Factor volumétrico del petróleo, (Bls / BN)
Bw = Factor volumétrico del agua, (Bls / BN)
RAP = Relación agua-petróleo de producción, (BN / BN)
A = Area seccional de la tubería, (pies2)
Z = Factor de compresibilidad del gas, (adim)
P = Presión, (lpc)
T = Temperatura, (oF)
80
n L L g g (2.30)
S LH g g
H
L
(2.31)
S L H L g H g (2.32)
L o f o w f w (2.33)
Las viscosidades del gas natural, del petróleo crudo y del agua pueden ser
calculadas mediante correlaciones empíricas (ver anexo A) si no se dispone de datos
de laboratorio.
L o fo w fw (2.34)
donde,
81
g.- densidad de los fluidos.
La densidad de los fluidos fluyentes es, tal vez, la variable de mas peso en
la ecuación general de pérdidas de presión en tuberías, principalmente en flujo
vertical, donde el gradiente de energía potencial corresponde al peso de la columna
de fluido. Las ecuaciones son las siguientes:
62.4 O 0.0136 g RS
O (2.35)
BO
g
g 0.0136 (2.36)
Bg
62.4 w
w (2.37)
Bw
Con,
Mg
g (2.38)
M aire
141.5
o (2.39)
131.5 API
Donde,
O = Densidad del petróleo y su gas en solución, lbs/pie3
g = Densidad del gas, Lbs/pie3
w = Densidad del agua, Lbs/pie3
O = Gravedad especifica del petróleo, adim.
g = Gravedad especifica del gas (aire = 1.0)
Mg = Peso molecular del gas, Lbs / Mol
Maire= Peso molecular del aire = 28.96 Lbs / Mol
API = Gravedad API del petróleo.
82
La densidad de la fase líquida se calcula en proporción al flujo fraccional de
petróleo y agua. Esto es,
L o fo w f w (2.40)
83
III.- COMPORTAMIENTO DE FLUJO EN TUBERIAS EDUCTORAS.
P2
144 V dP h W f 0 (2.42)
P1
Donde,
V = Volumen específico de la mezcla fluyente, pie3 / Lb.
P = Presión, lpca.
84
h = Diferencia de altura, pies.
Wf = Pérdida de energía debidas a irreversibilidad y otras causas.
4 f ' h
2
Wf (2.43)
2 g c t
1
2 g c t
Vm 1
V (2.45)
M
Donde,
M = masa total (petróleo, agua y gas) asociada a 1 BN de petróleo, lbs.
Vm = Volumen de la mezcla asociada a 1 BN de petróleo.
144
Vm dP
h M (2.46)
4 f '
2
1
2 g c t
V m dP V m P (2.47)
85
Sustituyendo esta expresión en la ec. 2.46 y arreglando en términos de
gradiente de presión, resulta:
2M f'
2
M
P g c t
144 (2.48)
h Vm
haciendo uso de la ec. 2.45 para introducir el término densidad, , la ecuación 2.48
puede ser expresada, en unidades de campo, como:
P 1 f ' Qo2 M 2
m 5
(2.50)
h 144 7.413 10 m t
10
Donde,
86
Poettman y Carpenter usaron la ec. 2.50 para calcular el factor de fricción de
correlación, f’, con la información de los 49 pozos mencionados.
144 P
m t5 m
f ' 7.413 1010 h (2.52)
2 2
Qo M
con
Q M
Log o
t 2.4158 (2.54)
2
87
POETTMAN & CARPENTER
FACTOR DE FRICCION
100
10
f'
0.1
0.01
0.001
4 5 6 7
10 10 10 10
Q M /
o t
88
b.- Baxendell y Thomas (5)
Con,
Qo M
0.14286 10 6 0.07143 (2.56)
t
R1 Rg K
K
5 RL
e 0. 011K (2.57)
89
BAXENDELL & THOMAS
FACTOR DE FRICCION
0.1
f'
0.01
0.001
0 1 2 3 4 5 6 7 8
Q M 10 /
o
6
t
90
Donde,
R1 = Nnúmero de Reynolds bifásico, adim..
Rg = Número de Reynolds de la fase gaseosa (gas libre), adim.
RL = Número de Reynolds de la fase líquida y gas en solución, adim.
K = Masa de gas libre/masa de líquido, basado en condiciones de
superficie.
0.0764 g RGP RS
K (2.58)
350.18 O RAP W
L L t
RL (2.58a)
L
g g t
Rg (2.58b)
g
En unidades de campo:
RL
Qo 350.18 O RAP W 0.0764 g RS (2.59)
45.6 L t
1.3 10 3 Qo g RGP RS
Rg (2.60)
g t
91
TEK
FACTOR DE FRICCION
100
10
f'
0.1
0.01
0.001
2 2.5 3 3.5 4 4.5 5 5.5 6
R2
92
d.- Hagedorn & Brown (7)
m2
m
P 1 f QL2 M L2 2 gc
(2.64)
h 144 2.965 10 m t h
m 11 5
P
m m2 2 g c (2.65)
h C 144 h
93
f 2 h
Wf (2.66)
2 g c t
L 62.4 O fW W O (2.68)
donde,
m L H L g 1 H L (2.69)
94
La ecuación 2.64 involucra dos incógnitas, f y HL. El factor de
entrampamiento de líquido es evaluado en función de los parámetros de flujo. Para
ello, Hagedorn y Brown emplearon ciertos grupos adimensionales desarrollados
previamente por Ros (8). Los 4 (cuatro) grupos usados son:
N L 1.938 SL L /
0.25
(2.70)
N g 1.938 Sg L /
0.25
(2.71)
N d 120.872 t L /
0.5
(2.72)
N L 0.15726 L 1 / L 3 0.25
(2.73)
donde,
N L P 0.1 C NL
= 0.575 (2.74)
N g 14.7 Nd
95
donde P es la presión promedio, expresada en lpca, y CNL es una función
del número de viscosidad del líquido, NL, y puede ser obtenido de la fig. 2.5 ó
mediante la siguiente ecuación de ajuste:
HL
0.0429 1.228 7.139 2 20.412 3 19.2475 4 5.696 5 (2.77)
Con,
Log 5.69897
(2.78)
3.69897
N g N L0.38
C2 (2.79)
N d2.14
Condicionado a C2 0.01
96
1 0.3423 15.3605 2 40.5783 3 41.082 4 14.695 5 (2.80)
Con,
12.5 C2 0.125 (2.81)
0.022 QL M L
N RE ( H .B ) (2.82)
t LH g1 H
L L
m1 L1 g1
m2 L2 g 2
m2 m2 2 m21 (2.83)
L O fW W O (2.84)
97
HAGEDORN & BROWN
Correlación para el COEFICIENTE DE VISCOSIDAD
0.1
0.01
CNL
0.001
0.001 0.01 0.1 1
NL
0.8
0.6
0.4
H L
0.2
0
-6 -5 -4 -3 -2
10 10 10 10 10
Fig. 2.6 Correlación de Entrampamiento de Líquido (Hagedorn & Brown)
98
HAGEDORN & BROWN
Factor de Corrección Secundario
2.0
1.8
1.6
1.4
1.2
1.0
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
C N N
2 g
0. 38
L
/N 2.14
d
99
e.- Orkiszewski (9)
Flujo transición: Este tipo de patrón de flujo fue definido para caracterizar la
condición transitoria en la cual la fase continua cambia de líquido a gas. El tapón de
100
líquido entre las burbujas desaparece y una cantidad significante de líquido entra en
la fase gaseosa. En este caso, aunque el efecto del líquido sea significativo, los
efectos de la fase gaseosa son predominantes.
WL QO 4.053 10 3 O RAP W 8.84 10 7 g RS (2.86)
101
Bg = Factor volumétrico del gas, BLS/PCN.
BO = Factor volumétrico del petróleo, BLS/BN
BW ·= Factor volumétrico del agua, BLS/BN
RAP = Relación agua-petróleo de producción, BN/BN
RGP = Relación gas-petróleo de producción, PCN/BN.
RS = Relación gas-petróleo en solución, PCN/BN.
qL = Tasa de flujo líquido, pies3/seg.
qg = Tasa de flujo de gas, pies3/seg.
QO = Tasa de producción de petróleo, BPD
Entonces,
WL
L (2.90)
qL
Wg
g (2.91)
qg
Orkiszewski limita los regímenes de flujo de acuerdo a las siguientes
condiciones:
Patrón de flujo Límite
Burbuja qg/qt < LB
donde,
qg
gD 1.938 L L 0.25 (2.92)
A
102
qL
LS 50 36 gD (2.94)
qg
0.75
q
LM 75 84 gD L (2.95)
q g
qt q L q g
t (2.96)
A A
A t2 (2.97)
4
2
4 qg
0.5
1 qt qt
E g 1 1 (2.99)
2 b A b A b A
f L L
2
f (2.100)
64.4 t
103
Donde L es la velocidad superficial de la fase líquida, en pies/seg, y es dada
por la ecuación:
qL
L
A 1 E g
(2.101)
1488 t L t
N Re (2.102)
L
L b Wt A
m L (2.103)
t b
y
f L t L b
2
f (2.104)
64.4 t t b
Con,
qL
L (2.105)
A
= Coeficiente de distribución del líquido, adim.
b C1 C2 gc t (2.106)
104
0.4
0.3
Factor C1
0.2
0.1
0
0 10 20 30 40 50
NRe(b)
La fig. 2.9 fue extrapolada por Nicklin y asoc.(11) para poder evaluar el
1.5
00
30
1.4
b) <
e(
NR
Factor C2
1.3
0<
00
40
=
e ( b) 0
1.2 Nr 500
6000
8000
1.1
1
0 2000 4000 6000 8000
NRe
Fig 2.9 Coeficiente C2 (Orkiszewski)
105
Estos valores extrapolados pueden ser calculados mediante el siguiente grupo
de ecuaciones:
Si N Re (b ) 3000
Si N Re (b ) 8000
0.5
1359 L
b (2.108)
2 L t0.5
Con,
1.4237 4.957 10 5 N Ret0.5 (2.109)
0.0127 log L 1
0.248 0.167 t 0.113 log t (2.111)
1.415
t
106
b – Si t < 10 pies/seg. y fw 0.5
0.013 log L
0.681 0.232 log t 0.428 log t (2.112)
t1.38
0.0274 log L 1
0.161 0.569 log t F (2.113)
t1.371
Con,
0.01 log L 1
F log t 0.397 0.63 log t (2.114)
t
1.571
d – Si t 10 pies/seg. y fw 0.5
0.045 log L
0.709 0.162 log t 0.88 log t (2.115)
t0.079
A q Wt L
Para t 10 b t (2.116b)
qt qt b A
1488 b L t
N Re (b ) (2.117)
L
107
El factor de fricción, f, de la ec. 2.104 puede ser calculado mediante las
ecuaciones de Jain o Colebrook, para el número de Reynolds del líquido calculado
mediante la ec. 2.102
qg
Eg (2.118)
qt
f g g
2
f (2.119)
64.4 t
1488 g g t
N Re ( g ) (2.120)
g
L L
m m (tapón) g m ( neblina)
M g S
(2.121)
L M LS L M LS
108
De igual manera se calcula el gradiente de fricción,
L g g LS
f M f (tapón) f ( neblina) (2.122)
LM L S L M LS
109
f.- Otras correlaciones.
Tapón NY 4 y LB > NX LS
NY > 4 y LB > NX 26.5
Transición NY 4 y LS > NX LM
Neblina NY 4 y NX > LM
NY > 4 y NX > 26.5
Con
g
1/ 3 0.25
L w
N X Sg (2.123)
0.0764 w L
0.25
N Y SL L w (2.124)
w L
LB 0.51100 NY
0.172
(2.125)
NY
LS 8.6 (2.126)
0.263
110
LM 70 100 NY
0.152
(2.127)
Donde,
Sg
Eg (2.128)
bf
(13)
Zuber y Findlay demostraron que la velocidad de ascenso de las burbujas
en una corriente de flujo turbulento puede ser calculada por:
bf 1.2 t bs (2.129)
g c t L g
0.5
bs C (2.130)
L
111
Con,
C 0.345 1 e 0.029N 1 e3.37 Eo / m (2.131)
8440.8 t
1.5
1 L 0.5
N
g
(2.132)
L
1.46 10 4 t L g
2
Eo (2.133)
L
m = 10 para N 250
m = 25 para N 18
2 f ' m t
2
f (2.134)
g c t
112
adicionales: flujo líquido, flujo espumoso y lento, flujo vapor y flujo transitorio entre
espumoso y vapor.
Por otra parte, Chierici y Ciucci concluyeron que las propiedades físicas de los
fluídos, las cuales deben ser calculadas en cada etapa de cálculo del gradiente de
presión en la corriente de flujo, tienen un efecto considerable en el resultado final de
la caída de presión calculada. En consecuencia, trabajaron siempre usando las
mismas correlaciones para el cálculo de estas propiedades.
113
Beggs y Brill definieron tres regímenes de flujo, a saber: Segregado,
intermitente y distribuido, con una zona de transición entre los flujos segregados e
intermitente. Para cada patrón de flujo correlacionaron el factor de entrampamiento
de líquido, calculando primero el entrampamiento que existiría si la tubería fuera
horizontal y, luego, corrigiendo de acuerdo al ángulo de inclinación de la tubería.
t 2 t
rh 4 (2.135)
t 4
114
4
i
2
o
2
i o
rh (2.136)
i o 4
Donde,
o = diámetro externo de la tubería eductora, pies.
i = diámetro interno del revestidor, pies.
rh = radio hidráulico, pies.
Igualando las ecuaciones 2.135 y 2.136 resulta una expresión para el diámetro
hidráulico.
h i o (2.137)
dP g d f 2
sen (2.138)
dL g c g c dL 2 g c t
115
i.- Flujo en pozos de gas.
Para pozos de gas hay varios métodos que permiten calcular las pérdidas de
presión, tanto en tuberías verticales e inclinadas como en tuberías horizontales.
Economides (29) presenta una derivación de la ecuación general 2.138 para el flujo de
un fluído newtoniano compresible en una sola fase, como es el caso del gas. Esta
ecuación es:
Pf e CS Pi 2685
2 2
i 5 Sen
f M ( Z TQSC ) 2 CS
e 1 (2.139)
Con
CS 0.0375 g L Sen Z T (2.140)
Donde,
Pi = Presión de entrada en la corriente de flujo, lpca.
Pf = Presión de salida en la corriente de flujo, lpca.
Qsc = Tasa de producción de gas, MMPC/día.
T = Temperatura promedio fluyente, ºR.
Z = Factor de compresibilidad promedio del gas.
t = Diámetro interno de la tubería, pulgadas
= Angulo de inclinación con la horizontal (+ hacia arriba), grados
L = Longitud de la tubería, pies.
fM = Factor de fricción.
La ec. 2.139, expresada en términos de Pwf ó Pwh debe ser resuelta por ensayo
y error, puesto que el factor Z es función de la presión. A pesar de que esta
ecuación fue desarrollada para flujo de una fase gaseosa, ella puede ser usada en
pozos produciendo pequeñas cantidades de líquido junto con el gas, con resultados
satisfactorios. Para ello se usa la gravedad de la mezcla, dada por:
116
g 4591 L RGP
m (2.141)
1 1123 RGP
Una suposición válida para muchos pozos de gas es que el flujo es turbulento.
De esta manera, el factor de fricción dependerá solamente de la rugosidad relativa
de la tubería, puesto que el número de Reynolds acusa valores muy grandes. Una
expresión para fM es:
f M 2 Log3.71 t
2
(2.142)
117
IV.- COMPORTAMIENTO DE FLUJO EN TUBERIAS HORIZONTALES.
SL
L (2.144)
t
L1 316 L
0.302
(2.145)
2.4684
L2 0.9252 10 3 L (2.146)
118
1.4516
L3 0.1 L (2.147)
6.738
L4 0.5 L (2.148)
Los límites de los regímenes de flujo horizontal para esta correlación son:
119
El factor de entrampamiento de líquido depende del régimen de flujo y viene
dado por la expresión siguiente:
H L H L ( 0) (2.151)
a bL
H L ( 0) c (2.152)
N FR
Patrón de flujo a b c
Segregado 0.98 0.4846 0.0868
Intermitente 0.845 0.5351 0.0173
Distribuido 1.065 0.5824 0.0609
Tabla 2.1 - Patrones de flujo (Begs & Brill)
HL(0) L
C 1 L ln d eL N LV
f g
N FR (2.154)
120
NLV es el número de velocidad del líquido, dado por la ec. 2.70. Las constantes
d, e, f y g para cada condición de flujo se dan en la tabla 2.2.
Una vez conocido el valor de HL, la densidad de la mezcla bifásica, m, puede
ser calculada mediante la ec. 2.69.
Patrón de flujo d e f g
P g
m Sen (2.155)
L elev. g c
P f n t2
(2.156)
L fricc. 62.4 t
Donde,
n L L g 1 L (2.157)
f
f f n (2.158)
fn
121
El factor de fricción sin resbalamiento, fn, es calculado del diagrama de Moody
ó mediante la ecuación de Jaín o Colebrook para un Número de Reynolds dado por:
1488 t n t
N Re (2.159)
n
Con,
n L L g 1 L (2.160)
f
eS (2.161)
fn
Con
ln X
S (2.162)
0.0523 3.182 ln X 0.8725 ln X 0.01853 ln X
2 3
y
L
X (2.163)
H L2
P m t Sg P
L (2.165)
L acel. g c P
122
Definiendo
m t Sg
EK (2.166)
32.2 P
f n t2
m Sen
P 1 64.4 t
(2.167)
L 144 1 K
E
144
P L m
0.5
L (2.168)
P L L
P L m
0.5
g (2.169)
P L g
P L L
0.5
(2.170)
P L g
123
Baker propuso la siguiente ecuación general para los diferentes patrones de
flujo, suponiendo que ambas fases fluyen en forma turbulenta
a b
m c (2.171)
L
Patrón de flujo a b c
Burbuja 14.2 0.75 0.1
Tapón 27.315 0.855 0.17
Estratificado 15400 1.0 0.8
Ondulado 15400 0.75 0.65
Disperso 4.8 – 0.3125 t 0.343 – 0.021 t 0
Tabla 2.3 Regímenes de flujo (Baker).
NY L3 y L1 > NY L2
Tapón
L2 > NY L4
Estratificado
L2 > NY > L4
Ondulado
Donde,
124
L1 10 5.87 0.667 Log N X (2.172)
Con,
WL
NX (2.176)
Wg
Wg
NY (2.177)
0.5
g L
(2.178)
0.0764 62.4
1
72
2 3
62.4
L (2.179)
L L
Donde
WL = Velocidad de la masa de líquido, lbs/hr.
Wg = Velocidad de la masa de gas, lbs/hr.
f B
WL 0.234 Qo Bo w w o f w w o (2.180)
1 fw
125
Las densidades de los fluidos a condiciones de flujo vienen dados por las
ecuaciones 2.35 – 2.37.
f ('L ) L QL'
2
P
L 3.38 10 11 5 (2.182)
L t
f ('g ) g Qg'
2
P
L 3.38 10 11 5 (2.183)
g t
donde,
t = Diámetro interno de la tubería, pies.
QL’, Qg’ =Tasas de producción de líquido y gas a condiciones de flujo, bls/día.
L, g = Densidades del líquido y del gas, lbm/pie3
f B
QL' Qo Bo w w (2.184)
1 fw
P P
L 2m (2.186)
f L g
Los factores de fricción de ambas fases, f’(L) y f’(g), son calculados mediante la
ecuación de Jaín o de Colebrook, a pesar de que los gradientes de presión son
expresados en la forma de Fanning.
126
Los números de Reynolds para cada fase en flujo simple son obtenidos
mediante:
L QL'
N RE ( L ) 0.1231 (1.187)
L t
g Qg'
N RE ( g ) 0.1231 (1.188)
g t
P m P
L 144 Sen L (2.189)
f
donde,
m = Densidad de la mezcla (ec. 2.69), Lbs/pie3.
= Angulo de inclinación (positivo hacia arriba), grados.
1
HL (2.190)
1 0.3264 Sg
1.006
127
1
0.9
0.8
0.7
0.6
HL
0.5
0.4
0.3
0.2
0.1
0
0 10 20 30 40 50
sg (pies/seg.)
Fig. 2.10 Correlación de entrampamiento de líquido (Flanigan)
128
c.- Eaton – Brown(19)
P 1 f Wt t2 WL L2 Wg g2
L 144 g c WL Wg (2.191)
WL W g
t L
L g L g
q
2
q
2
WL L W g g
P f Wt qt
1.75 10 4
HL 1 H L (2.192)
t qt L t
L 5 4
donde,
WL, Wg, Wt = Flujo de masa (líquido, gas y total), lbs/seg.
qL, qg, qt = Tasa de flujo volumétrico (líquido, gas y total), pies3/seg.
t = Diámetro interno de la tubería, pies
HL = Factor de entrampamiento de líquido, adim.
129
Las tasas de flujo de masa y volumétricas para líquido y gas pueden ser
calculadas mediante las ecuaciones 2.86 – 2.89.
N LV 0.575
P
0.05
NL
0.1
HL (2.193)
0.00226
0.0277
N gV N d 14.65
H L a0 a1 a2 2 a3 3 a4 4 a5 5 (2.194)
Con,
Log 2.7
(2.195)
3.7
y
a0 = 0.03159
a1 = 2.56935
a2 = - 18.72592
a3 = 50.02431
a4 = - 49.86117
a5 = 16.95090
130
1
0.8
0.6
HL
0.4
0.2
0
0.001 0.01 0.1 1 10
Fig. 2.11 Correlación de Entrampamiento de Líquido (Eaton & Brown)
Las velocidades promedio de las fases vienen dadas por las ecuaciones 2.21 y
2.22. La densidad de la mezcla puede ser calculada mediante la ec. 2.69.
Un grupo adimensional (LR) que relaciona la tasa másica de líquido con la tasa
másica total. Esto es,
ML
( LR) (2.196)
MT
Un grupo adimensional (GR) que relaciona la tasa másica de gas con la tasa
másica total:
MG
(GR) (2.197)
MT
131
g g t
N RE ( g ) (2.199)
g
10
1
f (LR) 0.1
0.1
2" H2O
4" H2O
17" Destilado
0.01
0.001
1.E+03 1.E+04 1.E+05 1.E+06 1.E+07
1. 25
0.08333 WT t
GR 0.5
t g
132
Mediante pruebas de ensayo y error determinaron que para 1 = 0.1 y 2 =
0.5 los valores correlacionaban satisfactoriamente para cada diámetro de tubería
analizado. En razón de ello introdujeron el efecto del diámetro de tubería en el eje de
las abscisas, resultando finalmente las ecuaciones presentada en la figura 2.12.
En las ordenadas:
0.1
W
F f L (2.200)
WT
En las abscisas:
0.5 1.25
W 0.08333 WT t
G (2.201)
WT t g
(21)
Lockhart y Martinelli fueron los primeros en presentar una correlación
práctica, desarrollada en la Universidad de California en 1.949. El estudio fue llevado
a efecto en tuberías pequeñas, con diámetros entre 0.0586‖ y 1‖.
133
Los autores propusieron cuatro correlaciones de pérdidas de presión para
cuatro mecanismos de flujo diferentes, a saber:
(19)
En 1.958, Orin Flanigan presentó un estudio sobre flujo en tuberías sobre
terrenos accidentados (flujo inclinado), concluyendo que toda la caída de presión
ocurre prácticamente en las secciones de tubería con flujo hacia arriba y que la caída
de presión total en la línea disminuye a medida que aumenta el flujo de gas.
P
L EH H (2.203)
144
Donde,
L = Densidad del líquido, lbs/pie3
H = Sumatoria de los cambios de nivel, pies.
134
e.- Flujo de gas en tuberías horizontales.
P12 P22 a3
1 a4
Qg a1 E C a2 ta5 (2.204)
ZT L g
Correlación a1 a2 a3 a4 a5
Weymouth 433.5 1.0 0.5 0.5 2.667
Panhandle A 435.87 1.0788 0.5394 0.4604 2.618
135
En los casos de flujo de gas con contenido de líquido, el factor de eficiencia
puede ser obtenido mediante la correlación presentada por Flanigan, como se
muestra en la fig. 2.13, donde:
Qg Z
g 3.1194 10 4 (2.205)
t2 P
QL BN
RLG (2.206)
Qg MMPCN
100
Factor de Eficiencia (E) %
10
0.1 1 10
g RLG 0.32
Fig. 2.13 Correlación de Eficiencia (Flanigan)
Log g RLG 1
0. 32
(2.208)
2
a0 = -0.97924
a1 = 1.98166
a2 = -5.45825
a3 = 15.93089
a4 = -19.84962
a5 = 8.29970
136
V.- CAIDA DE PRESION EN RESTRICCIONES.
137
A medida que P2 disminuye gradualmente, el flujo másico aumenta a un ritmo
diferencial decreciente, hasta que finalmente alcanza un nivel constante, que
representa el flujo máximo que puede pasar a través del orificio para una presión
corriente arriba dada, P1. La relación Rc = P2 / P1 es llamada relación de presión crítica
o condición de flujo crítico. Esta condición se presenta en todo tipo de flujo
compresible, incluyendo flujo de una mezcla gas-líquido, no así en flujo de líquido
incompresible. La relación de presión crítica en flujo de gas ha sido calculada
analíticamente y observada experimentalmente en rangos de valores muy cercanos a
0.5. En flujo bifásico, esta relación acusa valores entre 0.5 y 0.6, dependiendo del
tamaño del orificio y de las características de los fluidos fluyentes.
Tasa de Flujo Másico
FLUJO CRITICO FL
UJ
O
SU
B-
CR
I TI
CO
P2 = P1 /2 P2 = P1
Presión Corriente Abajo (P2)
Fig. 2.16 Concepto de Flujo Crítico
138
A.- Chokes de superficie.
Una ecuación general para flujo de gas a través de restricciones puede ser
derivada combinando la Ecuación de Bernoulli y la Ecuación de Estado de los Gases.
Las pérdidas irreversibles están implícitas en un coeficiente de descarga usado como
variable en la ecuación, el cual depende del tipo de restricción. Esta ecuación
general aplica tanto para flujo crítico como sub-crítico.
k 2/ k
q SC
C n P1 d 2
y y k 1 / k (2.209)
g T1 Z1 k 1
Con,
C s C d TSC
Cn (2.210)
PSC
Donde,
qSC = Tasa volumétrica de gas, MPCN/día.
Cn = Coeficiente de flujo, adim.
d = Diámetro del orificio, pulgadas.
k = Relación de calor específico = Cp/Cv, adim.
P1 = Presión corriente arriba, lpca.
P2 = Presión de descarga, lpca.
T1 = Temperatura de entrada, ºR.
Z1 = Factor de compresibilidad del gas a P1 y T1, adim.
Cs = Constante = 27.611
Cd = Coeficiente de descarga (empírico), adim.
139
El valor del coeficiente de descarga, Cd, depende de la longitud y forma de
apertura de la restricción, así como del número de Reynolds. Si no se dispone de
esta información, se recomienda usar Cd = 0.82.
k / k 1
2
Yc (2.211)
k 1
C n P1 d 2
q SC (2.212)
g T1
b QL RGL c
P1 (2.213)
da
Donde,
P1 = Presión corriente arriba, lpca.
QL = Tasa de flujo de líquido, BN/día.
RGL = Relación gas/líquido, PCN/BN.
d = Diámetro del choke, pulgs.
140
Las constantes a, b y c para las diferentes correlaciones se presentan en la
tabla 2.5.
a b c
Gilbert 1.89 3.86 x 10 -3 0.546
Baxendell 1.93 3.12 x 10-3 0.546
Ros 2.0 4.25 x 10-3 0.5
Achong 1.88 1.54 x 10-3 0.65
10 50 1"
1/8"
141
COMPORTAMIENTO DE ESTRANGULADORES (ROS)
(1) Conectar
2000 3/4"
2000 5/8"
500 1000
1000 1/2"
500 7/16"
200 3/8"
150 200 500
QL con RGL
5/16"
100 100
1/4"
200
50
50
100
20
10 50 1"
1/8"
Las figuras 2.17 y 2.18 representan cronogramas de la ecuación 2.213 para los
métodos de Gilbert y Ros, respectivamente. Ellos son de uso práctico y sencillo para
calcular cualquiera de las variables involucradas en dicha ecuación. Para calcular QL
ó RGL, conecte los ejes correspondientes a las variables PW y Choke. El punto de
intersección en el ejec Pivote conectarlo al eje de RGL para obtener QL ó viciversa.
Para calcular la presión del cabezal (PW) ó el diámetro de Choke se procede de
manera similar. Obtener un punto en el eje Pivote conectando los ejes QL y RGL.
Luego, conectar este punto al eje de Choke para obtener PW, ó a éste para obtener
el diámetro de Choke requerido.
142
Con,
g
N (2.215)
L
P1
NP (2.216)
L L / g c 0. 5
1
QD (2.217)
1 R
N D d L g c L
0. 5
(2.218)
d 2 P1
Qo
Bo f w 0. 5
T Z RGL Rs 151 P1 o 2.17 10 4 g Rs f w
T Z RGL Rs 111 P1 o 2.17 10 4 g RGL f w
(2.219)
La determinación del límite entre los flujos crítico y sub-crítico para casos de
flujo bifásico es mas complicado que para flujo de una fase compresible. La
velocidad sónica en una mezcla bifásica depende de las propiedades del líquido y del
gas. Una correlación para calcular la relación de presión crítica fue presentada por
Sachdeva (28) en 1.986. Esta condición de flujo crítico es calculada iterativamente.
Yc N D
k k 1
(2.220)
143
Donde,
k 1 X 1 g1 1 Y
N (2.221)
k 1 X 1 L1
n 1 X 1 g 2
2
k
D 1 (2.222)
k 1 2 X 1 L1
X 1 C P CV
n 1 (2.223)
X 1 CV 1 X 1 C L
La calidad del gas es la relación entre el flujo másico del gas y el flujo másico
total. En unidades de campo viene expresada como:
0.0764 g RGL 1 f w Rs
X (2.224)
0.0764 g RGL 1 f w Rs 5.6146 Bo o 1 f w f w Bw w
1 X 1 1 Y k X 1 1 Y k 1 k
0. 5
0.525 C D d 2
QL P1 m 2 (2.225)
C M2
L1 g1 k 1
144
Donde,
1
X1 1 X 1
m2 (2.226)
g1 Y
1k
L1
(29)
Díaz presentó una correlación para calcular el comportamiento de
producción restringida por choke de superficie. La correlación está basada en
condiciones de flujo crítico y relaciona la presión del cabezal, Pwh, o presión corriente
arriba, con las tasas de flujo másico de líquido y gas y el diámetro del choke. La
ecuación, en unidades prácticas, viene expresada como:
Pwh (2.228)
S 1. 873
donde,
Pwh = Presión del cabezal (corriente arriba), lpc
QL = Tasa de producción de líquido, BD.
RGL = Relación gas/líquido de producción, PCN/BN.
g = Gravedad del gas.
o = Gravedad del petróleo.
w = Gravedad del agua.
fw = Flujo fraccional de agua, fracción.
S = Diámetro del choke, (1/64‖).
145
Díaz también presentó una ecuación para calcular la caída de presión a través
de estranguladores, basado en los resultados de las pruebas de producción de 50
pozos en el oriente de Venezuela. Esta ecuación es:
0. 4562
4.2744 Pwh S 0. 6794
Pck Pwh (2.229)
C QL0.1446
con,
C g RGL o 1 f w fw w
0. 0891 0. 0555
(2.230)
La condición para que exista flujo crítico es que P2 / Pwh Yc. Suponiendo que
Yc 0.5, la condición de flujo crítico podrá ser expresada como Pck P2.
COMPORTAMIENTO DE ESTRANGULADORES
7
Log (6)
P2 d 2
6
P d
1 1
5
Pwh 2 / Pwh 1
0
0 1 2 3 4 5 6 7
d2 / d1
Fig. 2.19 Comportamiento de Estranguladores. Relación PW - Diámetros de Chokes
146
Este gráfico permite calcular con aceptable aproximación la presión de cabezal
esperada cuando se cambia el tamaño del reductor de flujo. Con esto es posible
predecir el nuevo comportamiento de producción usando cualquiera de las
ecuaciones presentadas para "comportamiento de estranguladores", si se asume que
las variaciones en las condiciones del influjo son despreciables; es decir, tanto la
relación gas líquido (RGL) como el corte de agua permanecen invariables en las
nuevas condiciones.
Pivote Q1 Q2 S1 S2
5000 5000 64 64
2000 2000 48 48
40 40
1000 1000
32 32
500 500 28 28
24 24
200 200
20 20
100 100
16 16
50 50
12 12
20 20
10 10 8 8
1.112
S
Q2 Q1 2
S1
Fig. 2.20 NOMOGRAMA. Cambio de Producción con cambio de Choke (ec. de Gilbert)
147
CAMBIO DE PRODUCCION CON CAMBIO DE CHOKE (Díaz)
Pivote Q1 Q2 S1 S2
5000 5000 64 64
2000 2000 48 48
40 40
1000 1000
32 32
500 500 28 28
24 24
200 200
20 20
100 100
16 16
50 50
12 12
20 20
10 10 8 8
1.67
S
Q2 Q1 2
S1
Fig. 2.21 NOMOGRAMA. Cambio de Producción con cambio de Choke (ec. de Díaz)
Para usar estos Nomogramas se procede como sigue: Conectar los ejes Q1 y
S1, correspondientes a tasa de producción de líquido y diámetro del choke ( 1/64"),
respectivamente. El punto de intersección en el eje Pivote se conecta con el nuevo
diámetro de choke en el eje S2. El valor leído en el eje Q2 representa la tasa de
producción esperada.
148
a.- Flujo de gas.
P1 P2
1.048 10 6 g Z1 T1 q sc2 1 4 (2.231)
P1 d 4 C D2 Y 2
con,
P P2
Y 1 0.41 0.35 4 1 (2.232)
k P1
donde,
P1 = Presión corriente arriba, lpca.
P2 = Presión corriente abajo, lpca.
Z1 = Factor de compresibilidad del gas a P1 y T1.
T1 = Temperatura corriente arriba, ºR.
qSC = Tasa de flujo de gas, MPCN / día.
= Relación d/t.
d = Diámetro de la restricción, pulgs.
t = Diámetro interno de la tubería, pulgs.
CD = Coeficiente de descarga 0.9
= Factor de expansión, adim.
k = Cp/Cv (relación de calores específicos del gas).
149
1.078 10 4 n vt2
P1 P2 (2.233)
CD
con,
C D 0.233 8.4 10 4 N v 6.672 11.661 2 (2.234)
Nv = qg/qL
= d/t
n = Densidad de la mezcla sin deslizamiento, lbs/pie3 (ec. 2.157)
vt = Velocidad de la mezcla a través del choke, pies/seg (ec. 2.23)
CD = Coeficiente de descarga.
150
VI.- EJERCICIOS RESUELTOS.
g = 0.817 = 0.00060
Obtener las curvas de gradiente de presión vertical por los métodos de:
Poettmann & Carpenter, Baxendell & Thomas y Tek.
SOLUCION
T f Ts
GT 0.00667 º F / pie
h
151
T 160 0.00667 400 162.67 º F
160 162.67
T 161.33 º F
2
Tr T Tc 1.533
Pr P Pc 0.319 Z = 1.013
QO M O
6.8 10 6
t
152
Q M
Log O O
t 2.4158 1.0
2
P 1 f QO2 M O2
m 5
h 144 7.413 10 m t
10
P
0.133 lpc / pie P = 53.2 lpc
h
P (calculada) = 220.2 lpc
POETTMAN - CARPENTER
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12)
MD TVD P* Pwf Pprom T Tprom Ppr Tpr Z Bg Bo
pies pies lpc lpc lpc °F °F (lpca) (°F) (BY/MPCN) (BY/BN)
0 0 167.00 167.00 167.0 160.0 160.0 0.281 1.530 1.018 17.502 1.05751
400 400 220.03 220.03 193.5 162.7 161.3 0.322 1.533 1.013 15.228 1.06041
800 800 269.13 269.13 244.6 165.3 164.0 0.401 1.539 1.002 12.157 1.06622
1500 1500 352.12 352.12 310.6 170.0 167.7 0.503 1.548 0.989 9.615 1.07413
3000 3000 536.83 536.83 444.5 180.0 175.0 0.709 1.567 0.962 6.704 1.09083
6000 6000 1006.17 1006.17 771.5 200.0 190.0 1.214 1.604 0.900 3.751 1.13302
(13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25)
Rs Uo Ug Uw TAUL ALFA f qL qg qt RHO prom Gfric. DELT P
(PCN/BN) (cps) (cps) (cps) dns/cm pc/seg pc/seg pc/seg Lbs/pie3 lpc
26 2.404 0.012 0.425 21.6 1.0 0.001817 0.339 3.9485 4.287 4.59 15.8197
31 2.316 0.012 0.420 21.1 1.0 0.001817 0.340 3.4113 3.751 5.25 13.8408 53.0
41 2.160 0.012 0.411 20.3 0.342 2.6851 3.027 6.51 11.1682 49.1
54 1.980 0.012 0.399 19.4 0.344 2.0834 2.428 8.11 8.9576 83.0
81 1.687 0.013 0.376 17.9 0.350 1.3936 1.743 11.30 6.4318 184.7
152 1.243 0.014 0.336 14.9 0.363 0.6945 1.058 18.63 3.9021 469.3
153
Método de Baxendell & Thomas
T f Ts
GT 0.00667 º F / pie
h
T 160 0.00667 400 162.67 º F
160 162.67
T 161.33 º F
2
Tr T Tc 1.533
Pr P Pc 0.319 Z = 1.013
QO M O
6.8 10 6
t
QO M O
0.14286 10 6 0.07143 0.901
t
P 1 f QO2 M O2
m 5
h 144 7.413 10 m t
10
P
0.2938 lpc / pie P = 117.5 lpc
h
155
BAXENDELL - THOMAS
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12)
MD TVD P* Pwf Pprom T Tprom Ppr Tpr Z Bg Bo
pies pies lpc lpc lpc °F °F (lpca) (°F) (BY/MPCN) (BY/BN)
0 0 167.00 167.00 167.0 160.0 160.0 0.281 1.530 1.018 17.502 1.05751
400 400 274.05 274.05 220.5 162.7 161.3 0.363 1.533 1.007 13.405 1.06272
800 800 361.59 361.59 317.8 165.3 164.0 0.514 1.539 0.987 9.331 1.07294
1500 1500 500.20 500.20 430.9 170.0 167.7 0.688 1.548 0.963 6.834 1.08586
3000 3000 793.38 793.38 646.8 180.0 175.0 1.022 1.567 0.918 4.443 1.11237
6000 6000 1499.12 1499.12 1146.3 200.0 190.0 1.793 1.604 0.827 2.335 1.17842
(13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25)
Método de Tek
T f Ts
GT 0.00667 º F / pie
h
T 160 0.00667 400 162.67 º F
160 162.67
T 161.33 º F
2
156
3) Calcular las propiedades de los fluidos
Tr T Tc 1.533
Pr P Pc 0.437 Z = 0.997
0.0764 g RGP RS
k 0.1426
350.38 O RAP W
RL
QO 350.38 O RAP W 0.0764 g RS 61272
45.6 L t
157
1.3 10 3 QO g RGP RS
RG 1.17 10 6
g t
k e 0.11 k
R1 RG 5k RL 76038
8) Calcular f’
R2 2.98
0.8942
2.11
P 1 f QO2 M O2
m 5
h 144 7.413 10 m t
10
P
0.5068 lpc / pie P = 202.73 lpc
h
158
TEK
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11)
MD TVD P* Pwf Pprom T Tprom Ppr Tpr Z Bg
pies pies lpc lpc lpc °F °F (lpca) (°F) (BY/MPCN)
(12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22)
Se recomienda resolver este ejercicio por los métodos de Hagedorn & Brown y
Orkiszewski como medio de desarrollar habilidades y conocimientos sobre estos dos
métodos, los cuales son los mas usados en la industria petrolera nacional.
159
HAGEDORN - BROWN
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12)
MD TVD P* Pwf Pprom T Tprom Ppr Tpr Z Bg Bo
pies pies lpc lpc lpc °F °F (lpca) (°F) (BY/MPCN) (BY/BN)
0 0 167.00 167.00 167.0 160.0 160.0 0.281 1.530 1.018 17.502 1.05751
400 400 258.50 258.50 212.8 162.7 161.3 0.351 1.533 1.009 13.886 1.06205
800 800 349.94 349.94 304.2 165.3 164.0 0.493 1.539 0.990 9.758 1.07166
1500 1500 507.43 507.43 428.7 170.0 167.7 0.685 1.548 0.963 6.872 1.08564
3000 3000 857.09 857.09 682.3 180.0 175.0 1.077 1.567 0.911 4.181 1.11636
6000 6000 1636.49 1636.49 1246.8 200.0 190.0 1.949 1.604 0.809 2.101 1.19146
(13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24)
Rs Uo Ug Uw TAUL NL coef1 CNL VSL NLV VSG NGV
(PCN/BN) (cps) (cps) (cps) dns/cm pies/seg pies/seg
26 2.404 0.012 0.425 21.6 0.0140 0.352 0.00247 6.903 16.77 80.44 195.46
35 2.268 0.012 0.420 20.8 0.0135 0.346 0.00245 6.933 17.00 63.03 154.54
53 2.029 0.012 0.411 19.5 0.0127 0.335 0.00241 6.996 17.43 43.13 107.47
79 1.766 0.013 0.399 18.0 0.0117 0.320 0.00236 7.087 18.01 29.19 74.18
136 1.389 0.013 0.376 15.6 0.0103 0.297 0.00229 7.287 19.20 16.21 42.69
271 0.946 0.015 0.336 11.6 0.0088 0.267 0.00223 7.778 22.07 6.29 17.85
(25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36)
Nd PHI a1 HL/PHI C2 a2 PSI HL (NRe)2 FASES f DENSgas DENSMEZCLA
Lbs/pie3 Lbs/pie3
47.5 0.00005 0.387 0.250 0.0100 -0.0006 1.000 0.250 3,312,585 0.02465 0.635 13.791
48.4 0.00006 0.404 0.268 0.0075 -0.0315 1.000 0.268 3,036,942 0.02466 0.800 14.886
50.0 0.00008 0.429 0.299 0.0047 -0.0657 1.000 0.299 2,646,988 0.02467 1.139 16.742
52.0 0.00010 0.455 0.333 0.0029 -0.0885 1.000 0.333 2,297,811 0.02468 1.617 18.839
55.8 0.00013 0.494 0.389 0.0014 -0.1079 1.000 0.389 1,858,472 0.02470 2.657 22.360
64.8 0.00023 0.555 0.488 0.0004 -0.1201 1.000 0.488 1,315,658 0.02475 5.289 28.758
160
ORKISZEWSKI
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13)
MD TVD P* Pwf Pprom T Tprom Ppr Tpr Z Bg Bo Rs
pies pies lpc lpc lpc °F °F (lpca) (°F) (BY/MPCN) (BY/BN) (PCN/BN)
0 0 167.00 167.00 167.0 160.0 160.0 0.281 1.530 1.018 17.502 1.05751 26
400 400 288.49 288.49 227.7 162.7 161.3 0.375 1.533 1.005 12.987 1.06335 38
800 800 385.21 385.21 336.9 165.3 164.0 0.543 1.539 0.982 8.789 1.07475 60
1500 1500 537.34 537.34 461.3 170.0 167.7 0.735 1.548 0.956 6.353 1.08897 87
3000 3000 858.60 858.60 698.0 180.0 175.0 1.101 1.567 0.907 4.074 1.11815 140
6000 6000 1626.81 1626.81 1242.7 200.0 190.0 1.942 1.604 0.810 2.110 1.19093 270
(14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26)
Uo Ug Uw TAUL qL qg qt WL Wg Wt DENSL DENSg Vt
3 3 3
(cps) (cps) (cps) dns/cm pie /seg pie /seg pie /seg Lbs/seg Lbs/seg Lbs/seg Lbs/pie3 Lbs/pie3 pies/seg
2.404 0.012 0.425 21.6 0.339 3.9485 4.287 17.182 2.510 19.692 50.705 0.636 87.34
2.231 0.012 0.420 20.6 0.341 2.8813 3.222 17.223 2.469 19.692 50.548 0.857 65.64
1.962 0.012 0.411 19.1 0.344 1.8872 2.232 17.303 2.389 19.692 50.242 1.266 45.46
1.714 0.013 0.399 17.7 0.349 1.3098 1.659 17.398 2.294 19.692 49.859 1.751 33.79
1.373 0.013 0.376 15.5 0.358 0.7701 1.128 17.588 2.103 19.692 49.089 2.731 22.99
0.948 0.015 0.336 11.6 0.382 0.3108 0.692 18.052 1.639 19.691 47.304 5.275 14.11
B U R B U J A
(27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39)
VgD (L)B (L)B (L)S (L)M FLUJO NRe(L) f Eg VL RHOprom Gfricción DELT P
adim corr. TIPO pies/seg Lbs/pie3 lpc
193.0 -6766.9 0.13 646.25 764.64 TAPON 685362 0.0249 0.92 86.61 4.63 588.11
142.4 -3821.3 0.13 656.19 773.24 TAPON 553327 0.0250 0.89 64.92 6.17 330.39 956.2
94.9 -1832.6 0.13 673.32 787.99 TAPON 433001 0.0251 0.84 44.79 8.94 156.87 465.2
67.0 -1012.0 0.13 692.56 804.44 TAPON 365711 0.0251 0.79 33.16 12.06 85.62 477.3
40.6 -467.8 0.13 729.44 835.62 TAPON 305656 0.0252 0.67 22.45 17.80 38.77 590.5
17.4 -175.5 0.13 820.28 910.68 TAPON 261942 0.0253 0.43 13.76 29.02 14.09 898.6
T A P O N N I E B L A
(40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52)
Vb NReB DELTA DELTA RHOprom Gfric. DELT P qg (corr) Eg Vg (corr) NRe gas f RHOprom
corr. Lbs/pie3 lpc Lbs/pie3
18.0 141154 -0.210 -0.155 4.59 48.457 72.2225 0.9953 80.82 6257200 0.0246 0.87
14.7 124040 -0.202 -0.161 6.11 36.650 121.5 72.2350 0.9953 58.97 4455510 0.0246 1.09
11.7 111726 -0.193 -0.169 8.82 25.650 96.7 72.3567 0.9953 38.63 2779731 0.0247 1.50
10.1 108892 -0.188 -0.175 11.87 19.267 152.1 72.6030 0.9952 26.81 1813209 0.0247 1.98
8.6 113984 -0.182 -0.175 17.45 13.330 321.3 73.2346 0.9951 15.77 933047 0.0248 2.96
7.5 139064 -0.178 -0.138 28.44 8.419 768.2 74.9827 0.9949 6.36 257352 0.0253 5.49
TRANSICION
(53) (54) (55) (56) (57) Vt < 10 Vt > 10
Gfric. DELT P RHOprom Gfric. DELT P ALFA Vb(1) NReB Vb(2) NReB DELTA DELTA F
lpc Lbs/pie3 lpc >= >=
6.348 22.38 209.67 17.71 20.96 164471 17.99 141154 -5.677 -0.155 0.126
4.559 35.6 31.73 177.51 594.5 14.43 18.81 158556 14.71 124040 -4.266 -0.161 0.114
2.894 19.6 49.37 140.44 532.5 11.45 16.57 157815 11.73 111726 -2.955 -0.169 0.098
1.932 26.3 70.72 116.20 913.2 9.78 15.05 162872 10.06 108892 -2.196 -0.175 0.086
1.046 51.2 114.87 93.03 2169.8 8.29 13.33 177230 8.57 113984 -1.494 -0.175 0.069
0.336 135.8 235.09 80.20 6571.5 7.21 11.46 212738 7.49 139064 -0.917 -0.138 0.050
161
En la gráfica siguiente se presentan las curvas de gradientes correspondientes
a los cinco métodos considerados. Como puede observarse, para las condiciones de
flujo dadas en este problema, los gradientes de Hagedorn & Brown y Orkiszewski
son bastante aproximados. El gradiente de Baxendell & thomas se muestra
ligeramente inferior a los dos anteriores. Sin embargo, los gradientes de Tek y
Poettman & Carpenter muestran importantes desviaciones en comparación a los dos
primeros. Tek es demasiado pesimista en el sentido de que sus pérdidas de presión
son demasiado elevadas, mientras que para Poettman & carpenter las pérdidas de
presión son muy pequeñas.
Qo = 4931 BPD
RGP = 730 PCN/BN
Tw f = 200 ºF
Tw h = 160 ºF
% A&S = 0
API = 34
Pw h = 167 lpc
Grv. Gas = 0.817
Diám . = 3''
PRESION FLUYENTE, lpc
1500
1000
500
0
0
500
1000
1500
2000
2500
3000
3500
4000
4500
5000
5500
6000
162
Ejercicio No. 2.2
= 0.00060 Twh = 90 ºF
Construir las curvas de gradiente de presión horizontal por los métodos de Ovid
Baker y Beggs&Brill.
1) Estimar P y calcular P
Tr T Tc 1.457
Pr P Pc 0.64 Z = 0.958
163
Bg = 6.252 x 10-3 Bls/PCN
RS = 96 PCN/BN (Standing)
BO = 1.0469 Bls/BN (Standing)
O = 1.661 cps
g = 0.012 cps
L = 17.6 dinas/cm.
f B
QL' Qo Bo w w 7475 bls / dia
1 fw
f B
WL 0.234 Qo Bo w w o f w w o 90053 lbs / hr
1 fw
0.5
g L
4.002
0.0764 62.4
1
72
2 3
62.4
L 5.622
L
L
WL
NX 33.87
Wg
164
Wg
NY 14948
5) Calcular los factores de fricción y los gradientes de fricción para cada fase,
considerando que ellas fluyen solas en la tubería.
L QL'
N RE ( L ) 0.1231 27731
L t
g Qg'
N RE ( g ) 0.1231 2724516
g t
f '( g ) 0.0044
165
f ('L ) L QL'
2
P
L 0.000053
L 3.38 10 11 t5
f ('g ) g Qg'
2
P
L 3.38 10 11 5 0.000558
g t
P L L
0.5
0.307
P L g
4.4875 0 . 322
m 3.069
L0
P P
L 2m 0.0053 lpc / pie
f L g
5.6146 Qg
Vsg 13.869 pies / seg
86400 A
1
HL 0.179
1 0.3264 Sg
1.006
166
c – Calcular la densidad de la mezcla
P
L m Sen 0.0040 lpc / pie
densidd 144
P P P
L L L 0.0093 lpc / pie
f
1) Estimar P y calcular P
167
Tr T Tc 1.457
Pr P Pc 0.589 Z = 0.966
f B
QL' Qo Bo w w 7447 bls / dia
1 fw
f B
WL 0.234 Qo Bo w w o f w w o 89722 lbs / hr
1 fw
0.5
g L
3.825
0.0764 62.4
1
72
2 3
62.4
L 5.567
L L
168
WL
NX 31.85
Wg
Wg
NY 15680
5) Calcular los factores de fricción y los gradientes de fricción para cada fase,
considerando que ellas fluyen solas en la tubería.
L QL'
N RE ( L ) 0.1231 26671
L t
g Qg'
N RE ( g ) 0.1231 2745106
g t
f '( g ) 0.0044
169
f ('L ) L QL'
2
P
L 0.000053
L 3.38 10 11 t5
f ('g ) g Qg'
2
P
L 3.38 10 11 5 0.000615
g t
P L L
0.5
0.293
P L g
4.4875 0 . 322
m 3.022
L0
P P
L 2m 0.0056 lpc / pie
f L g
5.6146 Qg
Vsg 15.252 pies / seg
86400 A
1
HL 0.165
1 0.3264 Sg
1.006
170
c – Calcular la densidad de la mezcla
P
L m Sen 0.0066 lpc / pie
densidd 144
P P P
L L L 0.00094 lpc / pie
f
172
Cálculos para el primer tramo de la línea de flujo ( 5280 pies )
1) Estimar P y calcular P
Tr T Tc 1.457
Pr P Pc 0.64 Z = 0.958
f B
q L 0.65 10 4 Qo Bo w w 0.486 pies 3 / seg
1 fw
173
qL
SL 0.618 pies / seg
A
qg
Sg 13.847 pies / seg
A
SL
L 0.0428
t
t2
N FR 6.5
32.2 t
L
0. 25
N LV 1.938 SL 1.557
L1 121.97
L2 2.216
L3 9.711
L4 8.38 10 8
174
C(seg) = 5.662 (ec. 2.154)
(seg) = 1.576 (ec. 2.153)
0.98 0L.4846
H L ( seg ) 0.0868
( seg ) 0.2850
N FR
0.845 0L.5351
H L (int) 0.0173
(int) 0.1530
N FR
L3 N FR
A 0.4286
L3 L2
B 1 A 0.5714
62.4 O 0.0136 g RS
O 50.060 Lbs / pie 3
BO
2.7024 P g
g 1.525 Lbs / pie 3
ZT
L
X 0.97
H L2
S = 0.195 (ec. 2.162)
f f ( n) e S 0.0218
P
0.00638 lpc / pie
L
P 33.7 lpc
P2 P1 P 425 33.7
P2 391.3 lpc
177
Un pozo productor está fluyendo con 3.0 MMPCN de gas natural con una
gravedad específica de 0.72 (aire = 1) a través de una tubería vertical de 12000 pies
de longitud 2-7/8" O.D. (i = 2.259"). Las temperaturas de cabezal y de fondo son
95°F y 300°F, respectivamente. La rugosidad de la tubería es 0.00065 (este valor se
usa comunmente en tuberías nuevas). Se desea obtener los gradientes de presión,
con valores cada 1000 pies, para las siguientes condiciones:
1 - Presión en el cabezal de 1600 lpca
2 - Presión de fondo fluyente de 3000 lpca
Solución caso 1:
T f Ts
GT 0.01708 º F / pie
h
T 95 0.01708 1200 115.5 º F
95 115.5
T 105.3 º F
2
3) Calcular el factor de desviación del gas, Z.
178
Tr (T 460) Tc 1.478
Pr ( P 14.7) Pc 2.482 Z = 0.687
g = 0.012 cps
g = 8.217 Lbs/pie3
3.0 * 0.72
N RE 2.009 *10 4 1.16 *10 6
2.259 * 0.01654
fM = 0.0181
5) Calcular el nuevo nivel de presión a 1200'. Para ello, la ecuación 2.139 debe
ser arreglada en función de Pi.
Pi e CS Pf 2685
2 2 f M ( Z TQSC ) 2
i 5 Sen
1 e C S
Pi 1697.8 lpcm
6) Hacer una nueva iteración estimando una presión de 1697.8 lpcm. De esta
manera se obtendrá como resultado una presión igual a la estimada de
1697.8 lpcm, lo que significa que se obtuvo convergencia en la segunda
iteración.
179
7) Estimar una nueva profundidad ( 2400' para este caso) y repetir el
procedimiento según los pasos 1 - 6. De esta manera se procederá hasta
cubrir la profundidad total de 12000', donde la presión calculada será de
2577.7 lpcm.
La siguiente tabla contiene una hoja de cálculo EXCEL, cuyo uso facilita el
proceso de cálculo de esta metodología.
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14)
LONG. Tx Tprom P* Pprom. P Ppr Tpr A B C D Z Cs
pies °F °F lpcm lpcm (lpca) (°F)
0 95.0 - 1600.0 - 1600.0
1200 115.5 105.3 1697.8 1648.9 1697.8 2.518 1.478 0.405 1.152 0.078 0.154 0.683 -0.084
2400 136.0 125.8 1795.6 1746.7 1795.6 2.666 1.531 0.435 1.152 0.073 0.136 0.696 -0.079
3600 156.5 146.3 1893.5 1844.6 1893.5 2.814 1.585 0.462 1.148 0.068 0.119 0.710 -0.075
4800 177.0 166.8 1991.4 1942.5 1991.4 2.963 1.639 0.487 1.140 0.063 0.104 0.722 -0.072
6000 197.5 187.3 2089.3 2040.4 2089.3 3.111 1.692 0.511 1.127 0.059 0.091 0.735 -0.068
7200 218.0 207.8 2187.2 2138.3 2187.2 3.259 1.746 0.534 1.110 0.055 0.079 0.747 -0.065
8400 238.5 228.3 2285.0 2236.1 2285.0 3.407 1.799 0.555 1.088 0.050 0.069 0.759 -0.062
9600 259.0 248.8 2382.7 2333.9 2382.7 3.555 1.853 0.575 1.061 0.046 0.060 0.772 -0.059
10800 279.5 269.3 2480.3 2431.5 2480.3 3.703 1.907 0.593 1.029 0.042 0.052 0.784 -0.057
12000 300.0 289.8 2577.7 2529.0 2577.7 3.851 1.960 0.611 0.992 0.038 0.045 0.796 -0.054
180
Solución caso 2:
Tr (T 460) Tc 1.96
Pr ( P 14.7) Pc 4.526 Z = 0.765
g = 0.022 cps
g = 10.14 Lbs/pie3
3.0 * 0.72
N RE 2.009 *10 4 8.67 *10 5
2.259 * 0.02215
fM = 0.0182
181
5) Calcular el nuevo nivel de presión a 10800', usando la ecuación 2.139
Pf e CS Pi 2685
2 2 f M ( Z TQSC ) 2
i 5 Sen
e CS
1
Pi 2892.9 lpcm
6) Hacer una nueva iteración estimando una presión de 2892.9 lpcm. De esta
manera se obtendrá como resultado una presión de 2893 lpcm, que será
la presión convergente en la tercera iteración.
182
REFERENCIAS
1. – Beggs, H.D.: ― Production Optimization – using Nodal Analysis ―. Publicado por
OGCI Publications. Tulsa, 1987.
2. – Jain, A.K.: ― Accurate Explicit Equation for Friction Factor ―. J. Hyd. Div. ASCE.
Mayo, 1976.
3. – Moody,L.F.: ― Friction Factors for Pipe Flow ‖. Actas del ASME, vol. 66, 1944.
4. – Poettman, F.H. y Carpenter, P.G.: ― The Multiphase Flow of Gas, Oil and Water
Through Vertical Flow Strings with Application to the Design of Gas-Lift
Installations ―. Drilling and Production Practice, 1952.
5. – Baxendell, P.B. y Thomas, R.: ― The Calculation of Pressure Gradients in High-
Rate Flowing Wells ‖. JPT. Octubre, 1961.
6. – Tek, M.R.: ― Multiphase Flow of Water, Oil and Natural Gas Through Vertical
Flow String ―. JPT. Octubre, 1961.
7. – Hagedorn, A.R. y Brown, K.E.: ― Experimental Study of Pressure Gradients
Ocurring during Continuos Two-Phase Flow in Small Diameter Vertical Conduits ―.
Actas del AIME. Abril, 1965.
8. – Ros, N.C.J.: ― Simultaneous Flow of Gas and Liquid as Encountered in Well
Tubing ―. JPT. Octubre, 1961.
9. – Orkiszewski, J.: ― Predicting Two-Phase Pressure Drops in Vertical Pipe ―. JPT.
Junio, 1967.
10. – Griffith, P. y Wallis, G.B.: ― Two-Phase Slug Flow ―. JHT. Agosto, 1961.
11. – Nicklin, D.J., Wilkes, T.O. y Davidson, J.F.: ― Two-Phase Flow in Vertical Tubes
―. Actas del Inst. Chem. Engrs. 1962.
12. – Aziz, K., Govier, G. y Fogarasi, M.: ― Pressure Drop in Wells Producing Oil and
Gas ―. Journal of Pet. Canadian. Septiembre, 1972.
13. – Zuber, N. y Findlay, J.A.: ― Average Volumetric Concentration in Two-Phase
Systems ―. Actas del ASME. Noviembre, 1965.
14. – Chierici, G.L. y Ciucci, G.M.: ― Two-Phase Flow in Oil Wells, Prediction of
Pressure Drop ―. SPE No. 4316. Abril, 1973.
15. – Beggs, H.D. y Brill, J.P.: ― A Study of Two-Phase Flow in Inclined Pipes ―. JPT.
Mayo, 1973.
16. – Katz, D. y otros.: ― Handbook of Natural Gas Engineering ―. Publicado por
McGraw-Hill Book Company. 1959.
183
17. – Baker, Ovid: ― Multiphase Flow in Pipe Lines ―. Oil and Gas Journal. Noviembre,
1958.
18. – Lockhart, R.W. y Martinelli, R.C.: ― Proposed Correlation of Data for Isothermal
Two-Phase Two Component Flow in Pipe Lines ―. Chem. Eng. Progress. 1949.
19. – Flanigan, O.: ― Effect of Uphill Flow on Pressure Drop in Design of Two-Phase
Gathering Systems ―. OGJ. Marzo, 1958.
20. – Eaton, B.A. y Brown, K.E.: ― The Prediction of Flow Pattern, Liquid Holdup and
Pressure Losses Ocurring during Continuos Two-Phase Flow in Horizontal
Pipelines ‖. Reporte Técnico de la Universidad de Texas. Octubre, 1965.
21. – Gilbert, W.E.: ― Flowing and Gas-Lift Well Performance ―. Drilling and Production
Practice. API. 1954.
22. – Baxendell, P.B.: ― Bean Performance – Lake Wells ―. Internal Report. Octubre,
1957.
23. – Ros, N.C.J.: ― An Analysis of Critical Simultaneous Gas/Liquid Flow through a
Restriction and its Application to Flow Metering ―. Applied Science Research.
Febrero, 1961.
24. – Achong, I.: ― Revised Bean Performance Formula for Lake Maracaibo Wells ―.
Internal Report. Octubre, 1961.
25. – Omana, R.A., Houssiere, C.R., Brill, J.P. y Brown, K.E.: ― Multiphase Flow
Through Chokes‖. SPE No. 2682. Septiembre, 1969.
26. - Ashford, F.E.: ― An Evaluation of Critical Multiphase Flow Performance Through
Wellhead Chokes ―. JPT. Agosto, 1974.
27. - Sachdeva, R., Schmidt, Z., Brill, J.P. y Blais, R.M.: ― Two-Phase Flow Through
Chokes ―. SPE No. 15657. 1986.
28. – Díaz, C.: ― Comportamiento de Flujo Bifásico a través de Estranguladores en
Condiciones de Flujo Crítico ―. Reporte Interno. 1996.
29. - Economides, M.J., Daniel Hill, A. Y Ehlig-Economides, C.: "Petroleum
Production System". Publicado por Prentice Hall PTR. 1993.
184
BIBLIOGRAFIA
1. - Golan, M. y Whitson, C.H.: ― Well Performance ― Publicado por Prentice-Hall
Book Company. 1991.
2. - Beggs, M.D.: ― Production Optimization – Using Nodal Analysis ―. Publicado por
OGCI Publications. Tulsa, 1987.
3. - Brown, K.E.: ― Gas Lift Theory and Practice – Including a Review of Petroleum
Engineering Fundamentals ―. Publicado por Prentice - Hall, INC. 1937.
4. – Nind, T.E.W.: ― Principles of Oil Well Production ―. Publicado por McGraw-
Hill,Inc. 1964.
185
3 - METODOS DE PRODUCCION
Pwh1
R
1
PR
TP
Pwf1 R2 Pwh2
TP
(lpc)
Pwf2
PRESION
IP
R
Pwh1 WP
R
Pwh2
Q1 Q2 Qm ax
186
Otros dos términos que relacionan tasas de producción y presión se refieren al
comportamiento del flujo a través de estranguladores (chokes) y en la línea de flujo.
El CPR (Choke Performance Relation) establece la relación entre la presión del
cabezal del pozo y la tasa de producción para un tamaño de estrangulador o reductor
de flujo dado y el FPR (Flowline Performance Relation) establece la misma relación
para una determinada presión en la línea de flujo. Esta presión de la línea se refiere
generalmente a la presión de salida o corriente abajo del estrangulador. Obviamente,
la presión de la línea será igual a la presión del cabezal cuando no haya choke en el
cabezal del pozo. Una de estas dos relaciones debe ser usada conjuntamente con el
WPR para analizar el comportamiento del pozo en cuanto a las condiciones en el
cabezal. En la fig. 3.2 se presentan ejemplos ilustrativos de estas relaciones.
Las relaciones CPR y FPR pueden combinarse para obtener la relación del
comportamiento del flujo en superficie, SPR, que relaciona la presión del cabezal con
la tasa de producción para determinadas condiciones de choke y nivel de presión de
separación.
WP CK1
WP
R
PRESION DE CABEZAL
PRESION DE CABEZAL
R PL 1
Pwh 1
CK2>CK1 Pwh 1
FPR1
Pwh 2 PL 2< PL 1
R1
Pwh 2
CP
FPR2
2
PR
C
Q1 Q2 Q1 Q2
TASA DE PRODUCCION TASA DE PRODUCCION
187
INFLUJO o FLUJO ENTRANTE . Se refiere a las condiciones de un punto
seleccionado arbritrariamente en el sistema de producción, conocido como punto de
referencia o punto de balance (NODO), calculadas en la dirección del flujo. Por
ejemplo, si el nodo es la profundidad del punto medio de las perforaciones, la presión
del nodo será Pwf (presión de fondo fluyente) y el INFLUJO coincidirá con el IPR del
pozo. Si el nodo es el cabezal del pozo, el INFLUJO será dado por el
comportamiento del flujo en el yacimiento y el comportamiento del flujo en la tubería
vertical, calculado en la dirección del flujo (hacia arriba en pozos productores) y la
presión del nodo será Pwh. En general, un nodo localizado en cualquier punto del
sistema de producción tendrá el siguiente balance de INFLUJO
Debe quedar claro que este balance de INFLUJO no tiene nada que ver con el
COMPORTAMIENTO de INFLUJO de un pozo. No debe haber confusión al
respecto.
188
I - PRODUCCION POR FLUJO NATURAL
PUNTO DE FLUJO
NATURAL INESTABLE
PUNTO DE FLUJO
NATURAL ESTABLE
IPR
Pwf
TPR
qL
Fig. 3.3 Condición de Flujo Natural
189
mientras que el de la izquierda es el punto de flujo inestable. El punto de flujo
natural estable es expresado analíticamente como el punto de intersección donde las
pendientes de ambas curvas, IPR y TPR, son de signos opuestos. Como corolario
de esto se establece que si ambas pendientes son de signos iguales en un punto de
intersección, el mismo es el punto de flujo natural inestable y cualquier pequeño
cambio o fluctuación del flujo provocará un cambio en el estado de equilibrio del
sistema; es decir, o el pozo muere o se torna en condiciones de flujo estable..
En esta gráfica, la curva TPR construida para una presión de cabezal requerida
y una relación gas/líquido esperada no intersecta la curva del IPR. En el caso (a) se
logra una disminución de las pérdidas de presión mediante la inyección de un
volumen de gas adicional en el eductor, lo que se traduce en una reducción de la
presión de fondo fluyente. En el caso (b) se suple una energía adicional al sistema;
es decir, se crea una ―ganancia‖ de presión, proporcionada por una bomba de
subsuelo, que permite restaurar las condiciones de flujo en el punto de equilibrio P.
P
IP
R
IPR
(a) (b)
qL qL
190
El parámetro que gobierna el comportamiento del flujo en la superficie en
función del comportamiento de la tubería eductora, para una determinada presión de
separación, es la presión del cabezal del pozo, Pwh. El punto de equilibrio de flujo
será, en este caso, el punto de intersección de las curvas WPR y SPR, como se
muestra en la fig. 3.5, y representa las condiciones requeridas para transportar el
fluido desde el cabezal del pozo hasta el separador
de producción o hasta el múltiple de WP
R
producción general en la estación
Pwh
recolectora, según sea el caso. SPR
Si estas curvas no se intersectan
significa que la presión requerida es qL
menor que la disponible en el cabezal Fig. 3.5 Equilibrio del Flujo en Superficie
191
En la fig. 3.6 se ilustran los efectos de la variación de los parámetros de flujo
en las diferentes curvas de comportamiento, las cuales se explican por sí solas.
IPR TPR
AC
TU
AL d1
d2 > d1
Pwf
Pwf
d3 > d2
ACTUAL
DAÑADO
d4 > d3
IP
FUTURO
R
DEPLETADO
qL qL
TPR TPR
RGL1
Pwh1 RGL2 > RGL1
Pwh2<Pwh1
RGL3 > RGL2
Pwf
Pwf
Pwh3<Pwh2
IP
IP
R
R
qL qL
SPR SPR
CK1
Ps1
CK2 > CK1 Ps2 < Ps1
d1
CK3 > CK2 Ps2
d1
Pwh
Pwh
d2 > d1
W
W
PR
PR
qL qL
192
II - PRODUCCION POR LEVANTAMIENTO ARTIFICIAL CON GAS
(GAS LIFT)
193
Dos puntos son de relevante importancia para el ingeniero de producción en
torno a este tema: (1) diseñar la completación de un pozo para retornarlo a
producción o mejorar su productividad y (2) analizar el comportamiento de un pozo
productor para mejorar su eficiencia.
194
2500
0
50
=
P
TH
0
2000
40
0
30
Presión, LPC
1500 TPR
00
=5
L
1000 RG 00
10 IPR
00
WPR
8"
00
15
3/
20
2"
ck
ck
1/
500
no
SPR
0
0 200 400 600 800 1000 1200
qL, BPD
195
Pwf Pwh Gav H v Gbv H f H v (3.1)
Donde,
En la ec. 3.1 hay dos parámetros que pueden ser cambiados o manipulados
por el analista: la profundidad del punto de inyección y el gradiente de presión por
encima de ese punto. De esta manipulación dependerá la eficiencia del diseño en
cuanto a la respuesta en términos de tasa de producción y presión de fondo fluyente.
La presión del anular viene dada por el gradiente de la columna de gas. Puede
ser calculada para diámetros superiores a 3 - 1/2", despreciando el efecto de fricción,
por la ecuación:
0.0375 g H / Z T
PH Pso 14.7e 14.7 (3.3)
196
Fig. 3.9 Diagrama de Presión de un Pozo de LAG
Donde,
PH = Presión del anular a la profundidad H (vertical), lpc
T = Temperatura promedio, ºR
P = Presión promedio, lpca
197
Usualmente, la profundidad del punto de inyección es suministrada como
información básica para el diseño y debe ser la máxima posible a fin de generar la
menor presión de fondo fluyente y maximizar la capacidad de producción del pozo.
Sin embargo, esta profundidad podría estar limitada por las presiones disponibles en
el sistema de distribución de gas en superficie, en cuyo caso debe ser calculada por
el analista. Para ello, se determina el punto de balance de presión en el sistema
eductor-anular, donde ambas presiones son iguales, dado por la intersección de las
curvas de gradiente correspondientes. En términos prácticos, éste sería el punto de
inyección si ésta se llevara a cabo sin restricciones; es decir, si no hubiese caída de
presión en la válvula de inyección. Normalmente, la presión diferencial de la válvula
está en el rango 50 – 150 lpca.
Pv
Hb Hv (3.6)
Gbv Gc
198
La presión de balance, Pb , puede ser calculada mediante la ec. 3.3 y el
gradiente de flujo promedio en la tubería vertical por encima de la válvula de
operación, Gav , se obtiene despejando de la ec. 3.1
Pwf Pb
Gbv (3.7)
Hf Hb
f Hb
Hb Hb
f 'Hb
Pso 14.7 C Hb
f 'Hb 1 C e
Gbv
(3.9)
199
Con,
C 0.01875 g / Z T (3.10)
Luego, calcular Hb iterativamente hasta que ambos valores, Hb y Hb , coincidan
en un grado de tolerancia establecido. Las demás variables desconocidas pueden
ser calculadas como se explicó anteriormente, re-ordenando algunas ecuaciones
cuando se requiera.
200
Hv( k 1) Gs Gu Pso Pdh
Hv( k ) (3.11)
Gs Gc
Con
Pko Gs Ho Pdhi
Hv(1) (3.12)
Gs Gc
Ptv Pdh
Gu (3.13)
Hv
Las presiones del eductor y del anular en cada válvula son dadas por:
Pspring
Hv( k ) HV( k 1) (3.16)
Gs
Con
Pko
Hv(1) Ho (3.17)
Gs
201
P RES I O N
Pdhi Pw h Pdh Pt1 Pt2 Pt3 Pso Pko
N IV E L D E F L U ID O E S T A T IC O
GRADIENTE
Ho
GR A D I EN T E
GR
ES T A T I C O D E
D ES C A R GA
AD
P
IEN
DE OPERACIÓ
R
TE
O
DE
F
LF
H1 Pc1
U
LU
IDO
N
N
D
DE
I H2 Pc2
D DE
A SC
H3 A
RG Pc3
D GR A D I EN T E A F
F L U YEN T E D E
LU
P R OD U C C I ON
YE
NT
E
202
La veracidad de esta información es de extrema importancia por cuanto la
profundidad de la primera válvula depende de ella y el espaciado de las válvulas
restantes depende de la ubicación de la primera. Como factor de seguridad se
recomienda que la primera válvula se diseñe con una presión de arranque en el
cabezal igual a la presión del separador de prueba.
Pwh es la presión del cabezal del pozo en estado de producción. Esta presión
no tiene ninguna inherencia en el diseño de las válvulas de descarga, a no ser que
se considere que los fluidos de descarga y producción sean similares y que la
descarga se realice bajo las mismas condiciones de flujo establecidas para la
producción del pozo, en cuyo caso el gradiente fluyente de descarga coincidirá con el
gradiente fluyente de producción.
203
El gradiente fluyente de descarga es el usado generalmente para diseñar el
espaciado de las válvulas de descarga por debajo de la primera. Usualmente se
denomina gradiente de diseño.
204
Pc1, Pc2, Pc3, etc. corresponden a las presiones en el espacio anular a la
profundidad de la válvula correspondiente. Esta presiones son calculadas en función
del gradiente de la columna de gas en el anular, la primera de ellas, Pc1, usando la
presión de ―kick-off‖. Las siguientes, usando la presión de operación.
f c Pd Ab (3.22)
sobre el vástago del pistón asentado sobre la superficie del orificio que lo comunica
con el eductor, A p . Esto es,
205
f a Pcv Ab Ap Ptv Ap (3.23)
Pd Ptv R
Pcv (3.24)
1 R
Con
R Ap / Ab , adim.
La ec. 3.24 aplica tanto para inyección continua como para inyección
intermitente. La única diferencia es que, en válvulas para inyección intermitente, el
elemento difusor (―spread‖), definido como el diferencial entre la presión del anular
necesaria para abrir la válvula y la presión del anular a la cual la válvula cierra, debe
ser controlado para evitar el uso de altos volúmenes de gas inyectado requerido para
levantar la columa de líquido en la tubería.
206
El ―spread‖ es la característica
mas importante de las válvulas
operadas por presión y es
causado, básicamente, por la
existencia de áreas no balan -
ceadas en la válvula.
F Pc1 Ab A p Pt A p Fig. 3.11 “Spread” de válvula operada por presión
Cuando la válvula abre se tiene que F Pc1 Ab . Para que la válvula cierre se tiene que
“SPREAD”. Pc1 es la presión del anular necesaria para abrir la válvula y P c2 es la presión del
anular a la cual la válvula cierra.
207
expuesta al flujo mas la presión ejercida por el gas del espacio anular sobre el orificio
de asiento del vástago. O sea,
Pd Ptv 1 R
Pcv (3.27)
R
por:
208
INFLUJO
Ptv P R P( formación) P( debajode la válvula) (3.28)
EXFLUJO
Ptv Pwh P( sobrela válvula) (3.29)
RGL excesiva
L
RG 1
Ptv
EXFLUJO
INF
GL 1 LU
>R
RGL2 JO
GL 2
>R
L3
RG
qL
209
de inyección de gas podría estar limitado por la disponibilidad en el sistema de gas
lift del campo. Entonces, es conveniente determinar la RGL de inyección mas
efectiva en base a un análisis del comportamiento esperado de tasa de producción
con tasa de inyección de gas, como se ilustra en la fig. 3.13.
qL máxima teórica
qL permisible
qL
Qgi
Fig 3.13 Comportamiento qL vs. Qgi en Pozos de LAG
Este tipo de análisis también es usado a menudo para determinar cuál válvula
está en funcionamiento. Para la condición de flujo del pozo constante, se construyen
curvas de INFLUJO y EXFLUJO tomando como nodos la profundidad de la válvula
del punto de inyección y de dos o tres válvulas por encima de ella, como se muestra
en la fig. 3.14.
Ptv1 Hv1
Ptv3
Hv3
Ptv
Hv2
IN
FL
Hv3
UJ
O
EXFLUJO
q3 q2 q1
qL
Fig. 3.14 Efecto de la Profundidad de Inyección en Pozos de LAG
210
Cada punto de intersección de ambas curvas para cada profundidad de
inyección representa la tasa de producción esperada para cada presión del eductor a
la profundidad de la válvula. Conociendo la tasa de producción real del pozo se
podrá determinar cuál válvula está operando,
Hv3 q3
o cuales válvulas funcionan en forma
PROFUNDIDAD
alternada. Un gráfico de profundidad
Hv2 q2
versus presión para cada tasa de
producción esperada, como el ilustrado Hv1 q1
Hv 3 < Hv 2
Hv 2 < Hv 1
EX
Hv 1
Pwf
FL
UJ
IN ( IPR
O
FL )
UJ
O
qL
3.16 Efecto de la Profundidad de Inyección en Pozos de LAG
211
III - PRODUCCION POR BOMBEO DE SUBSUELO
212
abre para admitir la entrada de fluido al cilindro, mientras que la válvula de trabajo se
cierra. Este ciclo de succión termina cuando el pistón se posiciona en el tope del
cilindro o cámara de compresión, el cual en ese momento está cargado con fluido de
producción. Durante el movimiento descendente del pistón, o ciclo de compresión, la
acción de las válvulas se invierte. La válvula fija se cierra para retener el fluido de
producción en el cilindro, mientras que la válvula viajera se abre debido a la fuerza
de compresión ejercida por el pistón en su movimiento hacia abajo, permitiendo que
el fluido contenido en el cilindro se impulse hacia arriba por la fuerza de compresión
generada, con una presión mucho mayor que la de succión, provocando el
movimiento ascendente de la columna de fluido en la tubería eductora. La longitud
del recorrido del pistón, usualmente denominado ―carrera del pistón‖ (S) es
controlada en superficie mediante ajustes del recorrido de una cabilla especial, ―barra
pulida‖, colocada en el tope de la sarta de cabillas. Este parámetro, conjuntamente
con su frecuencia o velocidad, generalmente dada en carreras por minuto (SPM),
determinan la capacidad de producción de un pozo para ciertas características de la
Tubería
eductora
Cabillas
Ancla superior
Tubo de la bomba
Pistón
Válvula viajera
Válvula fija
Tubo de extensión
Ancla de gas
(A) (B)
Fig. 3.17 Bombas Operadas con Cabillas de Succión
A - Bajada con la sarta de tubería
B - Bajada con la sarta de cabillas
213
bomba. Esta velocidad también es controlada en superficie ajustando la potencia del
motor que acciona a la unidad de bombeo. En la fig. 3.17 se muestra un diagrama
esquemático de una bomba de tubería y una de cabillas.
A B
214
(2) Bombas Reciprocantes Hidráulicas. En este tipo de bombas de
desplazamiento positivo, la potencia es transmitida a un motor hidráulico acoplado a
la bomba. Estas pueden ser del tipo insertada, en la tubería o en el revestidor de
producción, o del tipo libre. Estas últimas son de dos clases: bomba libre con sartas
paralelas de tubería y bomba libre con bombeo a través del espacio anular. En la fig.
3.18 se muestra un diagrama esquemático de los cuatro tipos de bombas.
La unidad de bombeo hidráulico está formada por dos secciones principales. Una es
la bomba, que usa un pistón reciprocante de vaivén, en cierto modo similar a la
bomba de cabillas. La otra sección corresponde a la parte motriz, que también tiene
un pistón reciprocante conectado directamente al pistón de la bomba. La parte motriz
es accionada por la potencia generada por un fluido inyectado desde la superficie,
fluido motriz, que generalmente es petróleo extraído del pozo y filtrado
cuidadosamente para eliminar la presencia de sólidos que pudieran taponar el
sistema o dañar el motor de la bomba.
215
B. Bombas de Desplazamiento Dinámico. Aquí se hace referencia a los dos tipos
de bombas de desplazamiento dinámico mas usados en pozos petrolíferos.
Motor
Cabezal
Prensaestopas
Barra pulida
Revestidor de producción
Tubería eductora
Cabillas
Centralizadores
Bomba
Niple de asiento
Empacadura
216
Fluido
Producción
motriz
(mezcla)
Cable eléctrico
Fluido motriz
Tubería de la bomba
Orificio
Garganta Succión
Difusor
Bomba
Gargant a
Revestidor
Rejilla de succión
Difusor
Producción
(mezcla)
Perforaciones
Motor eléctrico
(2) Bombas de Chorro. Son bombas operadas por una corriente de fluido
motriz a alta presión, que converge en forma de chorro en la succión de la
bomba después de atravesar el orificio eyector. El chorro entra a la succión
a alta velocidad y baja presión, mezclándose con el fluido de formación.
Esta mezcla es expandida en un elemento difusor provocando una caída
abrupta de la velocidad de flujo. A medida que esta velocidad baja se
genera una ganancia de energía potencial que permite levantar la columna
de fluido del pozo (fluido motriz y de formación) y transportarla a través del
sistema de producción. En la fig. 3.21 se presenta un diagrama
esquemático de este tipo de bomba.
217
La diferencia fundamental entre el comportamiento de bombas de
desplazamiento positivo y bombas de desplazamiento mecánico estriba en la
relación tasa/presión de bombeo para ambos tipos de bombas. En las primeras, la
tasa de bombeo es independiente de la presión de bombeo. La tasa de bombeo
puede ser controlada desde la superficie mediante simples ajustes del sistema de
desplazamiento de la bomba; bien sea, variando la carrera del pistón o su velocidad.
En bombas de desplazamiento dinámico la tasa de bombeo depende de la presión
de bombeo. La tasa es baja a alta presión y viceversa. La relación tasa/presión de
bombeo se conoce como ―característica de la bomba‖, y es especificada por el
fabricante. Esta condición no puede ser controlada por el operador.
218
TIPO DE BOMBA PROFUNDIDAD CONDICIONES MINIMAS DE SUCCION
MACANICAS
DE TUBERIA HASTA 8000’ Nivel de fluido 10’ sobre la succión (qL<20 bpd)
Nivel de fluido 100’ sobre la succión (qL>50 bpd)
DE CABILLAS HASTA 6000’ Idem
219
Tomando como punto de balance (nodo) la profundidad de la bomba, la presión
de succión será dada por el INFLUJO y la presión de descarga es determinada por la
curva de EXFLUJO. La presión de fondo fluyente, Pwf, es dada por el
comportamiento IPR del pozo, de acuerdo a la tasa de producción deseada. Este
comportamiento IPR puede ser obtenido de una prueba de producción reciente o a
partir de un índice de productividad estimado. En ambos casos debe conocerse la
presión promedio del yacimiento; o en su defecto, la presión estática en el punto
medio del intervalo perforado, después de cierto tiempo de cierre del pozo; o sea,
con un nivel de fluido estable.
Donde,
PFLS = Presión estática de la columna de gas sobre el nivel de fluido, lpc
Las ecuaciones 3.30 y 2.40 deben ser resueltas de manera iterativa, por cuanto
Z y L son funciones de presión y temperatura.
220
Suponiendo un gradiente lineal de flujo por debajo de la succión de la bomba,
éste vendrá dado por:
O sea,
PS Pwf G fb H F H B (3.33)
Donde,
H B = Profundidad de la bomba, pies
PS = Presión de succión de la bomba, lpc
221
INFLUJO
EXFLUJO
qL Pw h Gradiente
P fluyente por
R encim a de la
O bom ba
HFL F
U
PFLS
N
D Gradiente
HB I PB fluyente por
D PD debajo de la
PS
A bom ba
D
HF Pw f PR
PRESION
qL
qL IPR
222
PD EXFLUJO
PRESION (nodo)
Presión de la
bomba, PB
PS
INF
LU
JO
qL
Fig. 3.23 Requerimiento de Bombeo de un Pozo
223
Lo discutido previamente aplica para determinar los requerimientos de bombeo
para cualquier tipo de bomba que se desee usar. Sin embargo, para el análisis del
comportamiento de pozos productores con equipo de bombeo, la aplicación
solamente es válida para bombas centrífugas, ya que en bombas reciprocantes la
tasa de producción y la caída de presión en la tubería de producción son
independientes de la presión de fondo fluyente.
q L 0.11655 S P N P2 (3.38)
Donde,
q L = Tasa de producción líquida máxima, BFD
S P = Carrera del pistón, pulgadas
N = Frecuencia o velocidad del pistón, carreras/minuto
P = Diámetro del pistón, pulgadas
S N
Wmax W fe Wcs 1 P H B
(3.39)
5400
S N2
Wmin W fe Wcs 1 P 0.1275 Wcs H B (3.40)
70500
224
Con,
W fe 0.7854 P2 GL (3.41)
GL 0.433 o w o f w (3.42)
CARGA DE TRABAJO
WT f T E c2 (3.44)
Con,
fT c
0.79360 5/8‖
0.78560 3/4‖
0.78498 7/8‖
0.78540 1‖
3 – De la barra pulida
225
4 – Por el motor
E = Esfuerzo, lpc
c = Diámetro de las cabillas, pulgadas
226
V - USO DEL DINAMOMETRO.
A1
CARRERA DEL PISTON (pulg)
Fig. 3.24 Dinagrama de un pozo funcionando normalmente
227
energía retorna al sistema. La magnitud de esta energía es suficiente para reducir las
cargas en el movimiento ascendente, resultando en una disminución de la carga de
la sarta de cabillas. A partir de E comienza a aumentar la velocidad del movimiento
ascendente y, en consecuencia, a incrementar la carga hasta un valor máximo
representado por el punto F. A partir de este punto y hasta A, la carga disminuye
como consecuencia de una desaceleración de la barra pulida.
228
VI.- EJERCICIOS RESUELTOS.
Se requiere diseñar una completación de Gas Lift para una tasa de producción de
600 BFPD, conociendo que un pozo vecino tiene una prueba reciente de 383 BFPD,
con una presión de fondo fluyente (Pwf) de 1850 lpc. Se conoce, además, que la
presión del yacimiento es del orden de 2000 lpc.
Solución
229
383
QL ( max ) 2935 BFPD
0.075 1.8 0.8 0.075
0 .5
5.0 QL
1.125 0.5 5.0625
QL ( max)
= 0.12
Pwf P R 1
H b 2439.9 e
0 .0121875 Hb / Z T 460 2605.6 0
PROCEDIMIENTO
(a) Estimar Hb
(b) Calcular PB ec. 3.7
(c) Calcular P , T , Z
(d) Calcular Hb
PB Pwf Gbv H f H b
TB 0.01093 H f 100
230
Tsc 159.6 331.1 g 374.8
Psc 639.2 35.55 g 616.1
TB 100
T
2
PB 870
P
2
231
(c) Calcular la presión del casing en el punto de inyección ( ecuación 3.5 )
Pcv Pso Gc H b
Pcv 870 0.0232 5081
Pcv 988 lpc
Hv = 5081 pies
Pcv = 988 lpc
Ptv = 888 lpc
Tv = 156 ºF
232
SOLUCION ALTERNA.
RGL1 RGL2
Z1 T1 Z 2 T 2
T 1 140 ºF
T 2 (156 + 100)/2 = 128 ºF
Z1 Z 2
233
Qgi 45.0 MPCN / DIA
(METODO PROPUESTO)
Considerando que las válvulas serán operadas por presión o por fluido.
234
4) Calcular el espaciado y la temperatura de las válvulas restantes. Tabular.
(ecuación 3.11)
1 1720 119
2 2645 129
3 3373 137
4 3946 143
5 4397 148
6 4752 152
7 5031 155
235
REFERENCIAS
1. – Brown, K.E.: ― Gas Lift Theory and Practice – Including a Review of Petroleum
Engineering Fundamentals ―. Publicado por Prentice - Hall, INC. 1937.
2. - Beggs, M.D.: ― Production Optimization – Using Nodal Analysis ―. Publicado por
OGCI Publications. Tulsa, 1987.
3. - Golan, M. y Whitson, C.H.: ― Well Performance ― Publicado por Prentice-Hall
Book Company. 1991.
BIBLIOGRAFIA
1. – Zaba, J. y Doherty, W.T.: ― Practical Petroleum Engineers Handbook ―.
Publicado por Gulf Publishing Company. 1956.
2. – Fagg, L.W.: ― Dinamometer Charts and Well Weighing ―. Petroleum
Transactions, AIME. 1950.
3. – Mene Grande Oil Co.: ― Pozos de Bombeo ―. Curso de Producción de Petróleo y
Gas. 1960.
4. – Kudu Industries.: ― Bombas de Cavidad Progresiva para Petróleo ―. Publicación
Interna.
5. – Winkler, H.W. y Smith, S.S.: ― Gas Lift Manual ―. Publicado por CAMCO,
Incorporated. Texas, 1962.
236
ANEXO A:
CORRELACIONES EMPIRICAS
(4)
En 1985, Sutton derivó sus ecuaciones de propiedades pseudo-críticas del
gas natural ( fig. A.2 ) basadas en medidas del factor de compresibilidad, Z, de 634
muestras de gas. El usó la ecuación de Dranchuk Abou-Kassem para el factor Z y
los factores de ajustes propuestos por Aziz para presencia de componentes no
hidrocarburos. Las ecuaciones de Sutton son:
237
700
650
Psc ( lpca )
600
550
500
450
Tsc ( ºR )
400
350
300
0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2
g
Fig. A.1 Propiedades Pseudo-Críticas del Gas Natural ( según Gatlin )
700 700
Temperatura pseudo-crítica (ºR)
600 650
500 600
400 550
300 500
0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7
g
Fig. A.2 Propiedades Pseudo-Críticas del Gas Natural (según Sutton)
238
Fig. A.3 Gráfico Comparativo de las Propiedades Pseudo-Reducidas
según Gatlin, Beggs y Sutton
239
II - FACTOR DE COMPRESIBILIDAD DEL GAS NATURAL
Z A 1 A e B C PprD (A2.)
con,
De
0.715 1.128 Tpr 0.42 Tpr2 (A.2d)
1.1 1.8
TPR = 2.5
1 05 1.7
2 1.
=
FACTOR DE COMPRESIBILIDAD, Z
0.3 1.05 1
0.2 0.9
0.1 0.8
0 1 2 3 4 5 6 7 7 8 9 10 11 12 13 14 15
PRESION PSEUDO-REDUCIDA, Psr PRESI ON PSEUDO-REDUCI DA, Psr
Fig. A.4 Factor de Compresibilidad del Gas ( Z ) como Función de las Propiedades
Pseudo-Reducidas
240
III - FACTOR VOLUMETRICO DEL PETROLEO
1 – Standing (5)
g
0. 5
con, F RS
1.25 T (A.3a)
o
2 – Frick (6)
con,
API 30 API 30
C1 = 4.677 x 10-4 4.670 x 10-4
C2 = 1.751 x 10-5 1.100 x 10-5
C3 = -1.811 x 10-8 1.337 x 10-9
gc g 1.0 5.912 10 5 API T Log P 114.7 ( A.6)
BO 1.0 10 F (A.7)
con,
g
0. 526
BO RS
0.968 T ( A.7c)
o
241
5 – Al - Marhoun (9)
BO
2.49 RS g 0.1046 ( A.9)
P 0. 0526 10 0. 48 O
con,
1.69
O ( A.9a)
P 0. 0429 10 0. 000796 T
1 – Standing (5)
1. 20482
0.054945 P 1.4
RS g 0. 00091 T 0. 0125 API ( A.11)
10
242
2 – Frick (6) (Standing ajustado)
1. 20482
P
RS g 0. 00091 T 0. 0125 API
( A.12)
1810
con,
API 30 API 30
C1 = 0.0362 0.0178
C2 = 1.0937 1.1870
C3 = 25.7240 23.9310
1. 22549
P API 0.989
RS g
( A.14)
T 0.172
con,
5 – Al - Marhoun (9)
RS 185.8432P 1g.87764 O 3.1437 T 460
1. 32657 1. 398441
( A.15)
243
6 – Mannucci y Rosales (10)
1.88679
P
RS g 0. 000922 T 0. 0072 API
( A.16)
84.88 10
RS C1
c2
P
C
3
10
C4 API T 460 ( A.17)
g
con,
API 30 API 30
C1 = 0.05958 0.0315
C2 = 0.7972 0.7567
C3 = 1.0014 1.0937
C4 = 13.1405 11.2895
244
VI - VISCOSIDAD DEL PETROLEO MUERTO
1. 163
OD 10 X T 1.0 ( A.19)
con,
O A OD
B
(A.21)
con,
245
VIII - VISCOSIDAD DEL PETROLEO SUB-SATURADO
O Ob P Pb B (A.22)
con,
11. 513 8.98 10 5 P
B 2.6 P 1.187 e
( A.22a)
g A 10 4 e
B g 62.4
C
( A.23)
con,
986.4
B 3.3448 0.01009 ( A.23b)
T 460
M g 28.96 g ( A.23d )
g 0.0136 g Bg ( A.23e)
246
X - VISCOSIDAD DEL AGUA
T 68 68 100
O 1.0 0.024 P 0. 45 68 ( A.26)
32
con,
68 39 0.2571 API ( A.26a)
247
XII - TENSION INTERFACIAL GAS/AGUA
CW f C ( A B T C T 2 ) ( A.28)
con,
248
Mg = Peso molecular del gas, Lbs/mol
P = Presión, Lpca
PC = Presión crítica, Lpca
Ppr = Presión pseudo-reducida, fracción
RS = Solubilidad del gas en petróleo, PCN/BN
S = Salinidad del agua, ppm cl-
T = Temperatura, ºF
TC = Temperatura crítica, ºF
Tpr = temperatura pseudo-reducida, fracción
Z = Factor de compresibilidad del gas, adim.
249
REFERENCIAS
1. – Gatlin, C.: ― Petroleum Engineering, Drilling and Well Completion ―. Publicado
por Prentice - Hall, INC. 1960.
2. – Brown, K.E.: ― Gas Lift Theory and Practice – Including a Review of Petroleum
Engineering Fundamentals ―. Publicado por Prentice – Hall, Inc. 1967.
3. – Beggs, H.D.: ― Production Optimization – Using Nodal Analysis ―. Publicado por
OGCI Publications. 1987.
4. – Sutton, R.P.: ― Compressibility Factors for High – Molecular Weigth Reservoir
Gases ―. SPE Nº 14265. Septiembre, 1985.
5. – Standing, M.B.: ― Volumetric and Phase Behavior of Oil Field Hydrocarbon
Systems ―. Publicado por Reinhold Publishing Corp. 1952.
6. – Frick, T.C.: ― Petroleum Production Handbook ―. Publicado por McGraw – Hill
Book Company. 1962.
7. – Vasquez, M. y Beggs, H.D.: ― Correlations for Fluid Physical Property
Prediction ―. J.P.T. Mayo, 1980.
8. – Glaso, Oistein: ― Generalized Pressure – Volume – Temperature Correlations ―.
J.P.T. Mayo,1980.
9. – Al – Marhoun, M.A.: ― PVT Correlations for Middle East Crude Oils ―. J.P.T.
Mayo, 1988.
10. – Mannucci, J.E. y Rosales, E.E.: ― Correlaciones de Presión de Burbujeo y
Factor Volumétrico del Petróleo para Crudos del Oriente de Venezuela ―.
Documento Interno CVP. Septiembre, 1968.
11. – Kartoatmodjo, F. y Schmidt, Z.: ― Large Data Bank Improves Crude Oil
Physical Property Correlations ―. O.G.J. Julio, 1994.
12. – Beggs, H.D. y Robinson, J.R.: ― Estimating the Viscosity of Crude Oil Systems
―. J.P.T. Septiembre, 1975.
13. – Ng, J.T.H. y Egbogah, E.O.: ― An Improved Temperature – Viscosity
Correlation for Crude Oil Systems ―. Pet. Society of CIM. Mayo, 1983.
14. – Lee, A.L., Gonzales, M.H. y Eakin, B.E.: ― The Viscosity of Natural Gases ―.
J.P.T. Agosto, 1966.
15. - Dodson y Standing, M.B.: " A Correlation for Water Compressibility ". Petroleum
Engineer International. Noviembre, 1980
250
ANEXO B:
CONSTANTES UNIVERSALES
FACTORES DE CONVERSION
TABLAS
I - CONSTANTES UNIVERSALES
251
II - FACTORES DE CONVERSION
LONGITUD
1 Pie = 0.3048 Metros
= 12 Pulgs
= 0.3333333 Yardas
= 1.8939 x 10-4 Millas (US)
AREA
1 Pie2 = 0.092903 Metros2
= 144 2
Pulgs
= 2.2956841x10-5 Acres
= 9.2903 x 10-6 Hectáreas
VOLUMEN
1 Pie3 = 0.1781 Bls (US)
= 7.4805 Gals (US)
= 28.317 Litros
-5
= 2.2957 x 10 Acres – pies
MASA
1 Lbm = 453.59 Gramos
= 7000 Granos
VELOCIDAD
1 Pie/seg. = 0.3048 Metros/seg.
= 1.09728 Kms/hr.
= 0.68182 Millas/hr.
= 12 Pulgs./seg.
ACELERACION
1 Pie/seg2 = 30.48 Cms/seg2
252
DENSIDAD
1 Lb/pie3 = 5.787 x 10-4 Lbs/pulg3
= 0.13368 Lbs/gals
= 0.016018 Grs/cm3
= 5.6146 Lbs/bbl
= 6.944 x 10-3 Lpc/pie
FUERZA
1 Lb.F = 4.4482 x 105 Dinas
= 453.59 Grms.F
= 4.4482 Newtons
POTENCIA
1 HP = 550 Lbs–pies/seg
= 745.7 Watts
= 7.457 x 109 Ergs/seg
PRESION
2
1 Lb/pulg = 0.068046 Atm.
= 2.036 Pulgs H2O a 32ºF
= 6.8948 x 104 Dinas/cm2
= 0.0703 Kgs.F/cm
2
PERMEABILIDAD
1 md = 0.001 Darcies
= 1.0623 x 10-14 Pies2
TASA de FLUJO
1 Bbl/dia. = 0.159 Mts3/dia
= 6.4984 x 10-5 Pies3/seg.
= 42 Gals/dia
253
VISCOSIDAD
1 Cps = 0.01 .Poises
= 9.9 x 10-9 Atm.-seg.
TENSION
SUPERFICIAL
1 Dina/cm. = 2.2 x 10-3 Lbm/seg.
TEMPERATURA
1 ºF = 1.8 x ºC + 32
= ºR - 460
= 1.8 x ºK - 460
254
CONSTANTES FISICAS DE LOS GASES
M g Pc Tc
Peso Gravedad Presión Temperatura
Molecular del Gas Crítica Crítica
Grms/mol Aire = 1 lpca lpca
Componente Símbolo
Metano CH4 16.042 0.554 673.1 343.5
255
VALVULAS OPERADAS POR PRESION - CAMCO
ESPECIFICACIONES
Ab Av Av/Ab 1 - Av/Ab T.E.F. = St *
Area efect. Diám. Del Area del Av/Ab/(1-Av/Ab) Factor de efecto
Tipo de del domo Orificio Orificio Factor de efecto del RESORTE
2 2
Válvula ( Pulg ) ( Pulgs ) ( Pulg ) del TUBING ( Lpc )
BK, BP-1, 0.3109 3/16 0.0291 0.094 0.906 0.103 SIN RESORTE
PK-1, J-40 1/4 0.0511 0.164 0.836 0.197
5/16 0.0792 0.255 0.745 0.342
CP-1, RP-1, 0.7096 1/4 0.0511 0.072 0.928 0.078 SIN RESORTE
CK, CKF 5/16 0.0792 0.112 0.888 0.126
3/8 0.1134 0.160 0.840 0.190
J-20 ó R-20 0.7650 1/4 0.0511 0.067 0.933 0.072 SIN RESORTE
5/16 0.0792 0.103 0.897 0.115
3/8 0.1134 0.148 0.852 0.174
7/16 0.1538 0.201 0.799 0.252
1/2 0.2003 0.262 0.738 0.355
o
* Para válvulas con presión de apertura (Pvo) a 60 F < 1000 Lpc
** Válvulas parcialmente balanceadas para reducir la relación A v/Ab
BAKER
0.7650 3/16 0.0291 0.038 0.962 0.040 SIN RESORTE
1/4 0.0511 0.067 0.933 0.072
5/16 0.0792 0.103 0.897 0.115
3/8 0.1134 0.148 0.852 0.174
7/16 0.1538 0.201 0.799 0.252
1/2 0.2003 0.262 0.738 0.355
256
ESPECIFICACIONES API DE TUBERIAS
TUBING - UPSET
Diámetro Peso DIAMETROS Coupling CAPACIDADES
3
Nominal Lbs/pie Externo Interno Drift O.D. pies/bl Bls/100' pies /100'
TUBING - NO UPSET
3/4" 1.14 1.050 0.824 0.730 1.313 1510.00 0.066 0.37
257
ESPECIFICACIONES API DE CASING
258
ANEXO C
CURVAS DE
GRADIENTE
VERTICAL
259
CURVAS DE GRADIENTE VERTICAL
P R E S I O N (Lpc)
PARAMETROS DE FLUJO
QL = 50 BPD
A&S = 0 %
-2000 API = 20
Ggas = 0.7
Dtbg = 1.995"
-4000
P R O F U N D I D A D (Pies)
-6000
-8000
RGL
0
-10000
-12000
100
-14000
50
25 0
15 0
10 0
8000
60
40
30
20
0
0
0
0
0
260
CURVAS DE GRADIENTE VERTICAL
P R E S I O N (Lpc)
PARAMETROS DE FLUJO
QL = 50 BPD
A&S = 25 %
API = 20
-2000
Ggas = 0.7
Dtbg = 1.995"
-4000
P R O F U N D I D A D (Pies)
-6000
-8000
RGL
0
-10000
-12000
100
-14000
25 00
60
30
20
75
50 0
15 0
10 0
80
40
0
0
0
0
00
0
261
CURVAS DE GRADIENTE VERTICAL
P R E S I O N (Lpc)
-4000
P R O F U N D I D A D (Pies)
-6000
-8000
RGL
-10000 0
-12000
100
-14000
15 0
60
20
50
2500
10
80
40
30
00
0
0
00
0
262
CURVAS DE GRADIENTE VERTICAL
P R E S I O N (Lpc)
PARAMETROS DE FLUJO
QL = 100 BPD
A&S = 25 %
API = 20
-2000
Ggas = 0.7
Dtbg = 1.995"
-4000
P R O F U N D I D A D (Pies)
-6000
-8000
RGL
0
-10000
-12000 100
200
-14000
50
25 0
15 0
10 0
80
60
40
30
0
0
0
00
0
0
263
CURVAS DE GRADIENTE VERTICAL
P R E S I O N (Lpc)
PARAMETROS DE FLUJO
QL = 200 BPD
A&S = 0 %
API = 20
-2000
Ggas = 0.7
Dtbg = 1.995"
-4000
P R O F U N D I D A D (Pies)
-6000
-8000
RGL
-10000 0
-12000 100
200
-14000
60
25
15
10
80
40
30
0
00
00
00
0
264
CURVAS DE GRADIENTE VERTICAL
P R E S I O N (Lpc)
PARAMETROS DE FLUJO
QL = 300 BPD
A&S = 0 %
API = 20
-2000
Ggas = 0.7
Dtbg = 1.995"
-4000
P R O F U N D I D A D (Pies)
-6000
-8000
RGL
0
-10000
100
-12000
200
-14000
10
40
30
15
8000
60
0
00
0
0
265
CURVAS DE GRADIENTE VERTICAL
P R E S I O N (Lpc)
PARAMETROS DE FLUJO
QL = 300 BPD
A&S = 25 %
API = 20
-2000
Ggas = 0.7
Dtbg = 1.995"
-4000
P R O F U N D I D A D (Pies)
-6000
-8000
RGL
0
-10000
100
-12000
200
-14000 300
80
40
15
10
60
0
0
00
00
266
CURVAS DE GRADIENTE VERTICAL
P R E S I O N (Lpc)
CONDICIONES DE FLUJO
QL = 400 BPD
A&S = 0 %
API = 20
-2000
Ggas = 0.7
Dtbg = 1.995"
-4000
P R O F U N D I D A D (Pies)
-6000
-8000
RGL
-10000 0
100
-12000
200
-14000
10 00
40
15
8000
60
30
0
0
0
267
CURVAS DE GRADIENTE VERTICAL
P R E S I O N (Lpc)
CONDICIONES DE FLUJO
QL = 400 BPD
A&S = 25 %
API = 20
-2000
Ggas = 0.7
Dtbg = 1.995"
-4000
P R O F U N D I D A D (Pies)
-6000
-8000 RGL
-10000
100
-12000 200
300
-14000
15
10 0
80 0
60
40
0
0
0
0
268