Teoria - Ejemplo de Aplicación (Pozo Guaricho 12) PDF

UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO
FACULTAD DE INGENIERÍA
“Aplicación del Análisis Nodal para

Incrementar la Productividad de un Pozo”
TESIS PROFESIONAL
PARA OBTENER EL TÍTULO DE:
I NG E N I E R O PETROLERO
P R E S E N T A :
ALVAREZ RAYA ROGELIO
DIRECTOR DE TESIS: M. en I. Tomás Eduardo Pérez García.
MÉXICO, D.F., CIUDAD UNIVERSITARIA, “2008”

CONTENIDO
RESUMEN --------------------------------------------------------------------------------------------------- 1
INTRODUCCIÓN ----------------------------------------------------------------------------------------- 2
I COMPORTAMIENTO DE FLUJO EN EL YACIMIENTO
I.1 Flujo de aceite

I.1.1 Ley de Darcy ------------------------------------------------------------------------------- 5
I.1.2 Ecuaciones de flujo ---------------------------------------------------------------------- 5
a) Estado pseudo-estacionario
b) Estado estacionario
I.2 Índice de productividad y el IPR
I.2.1 Definición del índice de productividad (J) ------------------------------------------ 6
I.2.2 Definición de IPR (Inflow Performance Relationship) --------------------------- 8
I.2.3 Índice de productividad para yacimientos bajo saturados -------------------- 8
I.2.4 Índice de productividad para yacimientos saturados -------------------------- 9
I.3 Flujo de dos fases: Aceite y Gas

I.3.1 Ecuación de Vogel ---------------------------------------------------------------------- 11
I.3.2 Procedimiento de Standing para EF ≠ 1 ------------------------------------------ 14
I.3.3 Ecuación de Fetkovich ----------------------------------------------------------------- 18
I.3.4 Ecuación de Jones, Blount y Glaze ------------------------------------------------- 20
I.3.5 Ecuación de Klins y Majcher --------------------------------------------------------- 20
I.3.6 Ecuación de Sukarno ------------------------------------------------------------------ 20
I.4 Flujo combinado: Una fase (Pwf > Pb) y dos fases (Pwf < Pb)
I.4.1 Ecuación General ----------------------------------------------------------------------- 21
I.4.2 Ecuación de Hasan --------------------------------------------------------------------- 23
I.4.3 Ecuación de Fetkovich ----------------------------------------------------------------- 23
I.5 Curvas de IPR para pozos de gas

I.5.1 Ecuaciones de flujo -------------------------------------------------------------------- 24
a) Estado estacionario
b) Estado pseudo-estacionario
I.5.2 Pruebas de potencial (Pruebas multigastos) ------------------------------------ 25
I.5.3 IPR para pozos horizontales de gas ---------------------------------------------- 28
I.6 IPR de dos fases para pozos direccionales ----------------------------------------- 30

I.7 IPR de dos fases para pozos horizontales ------------------------------------------- 32
I.8 IPR para flujo de tres fases (aceite, gas y agua) --------------------------------- 35
I.9 Curvas de IPR futuras

I.9.1 Método de Standing -------------------------------------------------------------------- 39
I.9.2 Método de Fetkovich ------------------------------------------------------------------- 40
I.9.3 Método de Eckmier --------------------------------------------------------------------- 40
I.9.4 Método de Couto ------------------------------------------------------------------------ 40
I.9.5 IPR futura para pozos de gas ------------------------------------------------------ 41
II FLUJO EN TUBERÍAS
II.1 Fundamentos de flujo a través de tuberías

II.1.1 Ecuación de balance de energía --------------------------------------------------- 43
II.1.2 Pérdidas de presión por fricción --------------------------------------------------- 45
a) Ecuación de Darcy
b) Factor de fricción
c) Número de Reynolds
II.2 Flujo de una sola fase

II.2.1 Flujo de líquido -------------------------------------------------------------------------- 48
II.2.2 Flujo de gas ------------------------------------------------------------------------------ 49
II.3 Flujo multifásico en tuberías

II.3.1 Propiedades de los fluidos ----------------------------------------------------------- 53
a) Definiciones --------------------------------------------------------------------- 53
b) Factores de volumen Bo y Bg ----------------------------------------------- 54
c) Densidad del gas disuelto -------------------------------------------------- 55
d) Densidad del aceite ----------------------------------------------------------- 56
e) Densidad del gas -------------------------------------------------------------- 56
f) Densidad del agua saturada ----------------------------------------------- 56
g) Viscosidad del aceite --------------------------------------------------------- 57
h) Relación de solubilidad Rs -------------------------------------------------- 57
II.3.2 Patrones de flujo ----------------------------------------------------------------------- 58
II.3.3 Colgamiento ------------------------------------------------------------------------------ 59
II.3.4 Velocidades superficiales ------------------------------------------------------------- 60
II.3.5 Velocidad real --------------------------------------------------------------------------- 61
II.3.6 Densidad de la mezcla ----------------------------------------------------------------- 61
II.3.7 Gasto másico ----------------------------------------------------------------------------- 63
II.3.8 Viscosidad de la mezcla --------------------------------------------------------------- 63
II.4 Flujo multifásico en tuberías verticales
II.4.1 Correlación de Poettmann y Carpenter ------------------------------------------- 64
• Procedimiento de cálculo
II.4.2 Correlación de Beggs y Brill ---------------------------------------------------------- 66
II.4.3 Curvas de gradiente de presión ---------------------------------------------------- 70
• Procedimiento para utilizar las curvas de gradiente
II.5 Flujo multifásico en tuberías horizontales

II.5.1 Correlación de Bertuzzi, Tek y Poettmann --------------------------------------- 72
II.5.2 Correlación de Beggs y Brill ---------------------------------------------------------- 75
II.5.3 Curvas de gradiente de presión ---------------------------------------------------- 77
II.6 Flujo por espacio anular -------------------------------------------------------------------- 78
II.7 Flujo a través de restricciones

II.7.1 Flujo multifásico a través de estranguladores ---------------------------------- 79
III BASES PARA EL ANÁLISIS NODAL
III.1 ¿Qué es el análisis nodal? ------------------------------------------------------------------ 81
III.2 Elección del nodo solución

III.2.1 Fondo del pozo como nodo solución --------------------------------------------- 83
III.2.2 Cabeza del pozo como nodo solución ------------------------------------------- 86
III.2.3 Separador como nodo solución ---------------------------------------------------- 88
III.3 Análisis nodal para pozos con restricciones

III.3.1 Estranguladores ------------------------------------------------------------------------- 90
III.3.2 Válvulas de tormenta ------------------------------------------------------------------ 92
III.4 Análisis nodal para tuberías telescopiadas ----------------------------------------- 94

• Procedimiento de cálculo.
III.5 Análisis nodal para sistemas multipozos --------------------------------------------- 96

IV EJEMPLO DE APLICACIÓN (POZO GUARICHO 12) ----------------------------- 101
CONCLUSIONES ---------------------------------------------------------------------------------------------- 113
NOMENCLATURA ------------------------------------------------------------------------------------------- 115
BIBLIOGRAFÍA ------------------------------------------------------------------------------------------------ 118

Resumen
RESUMEN
La planeación de la explotación de un campo petrolero debe tomar en cuenta todos los factores
técnicos, operativos, ecológicos y económicos que de alguna manera puedan afectar las
condiciones existentes. Desde el punto de vista técnico, es de vital importancia disponer de
información de buena calidad y efectuar el análisis del sistema integral de producción (SIP), que
permita identificar con precisión los problemas que ocasionan la baja productividad de los pozos y
diseñar las intervenciones o tratamientos más recomendables para optimizar la explotación.
En esta tesis se describen los conceptos que definen el comportamiento de flujo de los fluidos
desde el yacimiento hasta el separador, pasando a través de posibles restricciones que se puedan
presentar tanto en la tubería de producción (T.P) como en la línea de descarga (L.D).
También se describen los procedimientos que se tienen que realizar a un pozo o sistemas
multipozos para efectuar el análisis nodal, con esta técnica se busca optimizar las instalaciones de
producción con el objetivo de producir el gasto óptimo con la menor inversión económica posible.
Se aplicó la técnica del análisis nodal al Pozo Guaricho 12. Para este pozo se evaluaron dos
posibilidades viables económicamente hablando para aumentar la producción, las cuales son:
• Cambio del diámetro del estrangulador
• Producir por espacio anular
Del estudio del análisis nodal se concluyó lo siguiente:
• Las caídas de presión debido a la fricción, son mayores cuando se produce por la tubería
de producción (T.P) que cuando se produce por espacio anular.
• Para mantener el flujo crítico cuando se tiene flujo por tubería de producción (T.P) el
máximo estrangulador es de 32 / 64 pg y la producción sería de 1,415 bpd.
• Para mantener el flujo crítico cuando se tiene flujo por espacio anular el máximo
estrangulador es de 40 / 64 pg y la producción sería de 2,238 bpd.
• Por lo tanto la mejor opción es producir el pozo Guaricho 12 por espacio anular.
1
Resumen
INTRODUCCIÓN
La planeación de la explotación de un pozo petrolero debe tomar en cuenta los factores técnicos,
ecológicos y económicos, que de alguna manera puedan modificar las condiciones existentes.
La razón fundamental de someter un sistema integral de producción (SIP) a la técnica de análisis

nodal es porque este involucra en sus cálculos todos los elementos del sistema, permite
determinar el efecto de su variación de la capacidad de transporte y tener una imagen en conjunto
del comportamiento del pozo. Desde la perspectiva de evaluación esto es posible; sin embargo en
condiciones de diseño (sin restricciones económicas) es factible dimensionar el SIP en su conjunto
para obtener la capacidad de transporte requerida, o bien, la capacidad de transporte óptima,
teniendo como límite tan sólo la capacidad de afluencia del yacimiento.
Para saber si un pozo produce en forma apropiada, es necesario conocer su potencial, el cual se
define como el gasto máximo que aportaría el pozo si se le impusieran las mejores condiciones,
por lo que si un pozo no produce de la forma esperada, las causas de su baja productividad deben
ser determinadas para establecer el método correctivo adecuado. La producción puede ser
restringida por el comportamiento de un sólo componente del sistema, el efecto de cada
componente sobre el sistema puede ser analizado en forma aislada y así optimizar
económicamente el sistema.
Los objetivos de someter un pozo a la técnica del análisis nodal son los siguientes:
1) Determinar el gasto al cual un pozo productor de aceite o gas, producirá considerando la

geometría del pozo y las limitaciones de la terminación.
2) Determinar bajo qué condiciones de flujo, el pozo dejará de producir naturalmente (pozo
muerto), lo cual está relacionado con el tiempo de explotación.
3) Seleccionar el momento más adecuado para la instalación de sistemas artificiales de
producción.
4) Optimizar el sistema para producir un gasto determinado.
5) Evaluar cada componente del sistema para determinar si alguno de ellos está restringiendo
la producción innecesariamente.
Una recomendación importante para lograr la producción óptima de un pozo, es que él ingeniero
de producción, necesita entender y comprender en forma clara y completa los principios que rigen
el movimiento de aceite, gas y agua desde el yacimiento, hasta el separador. Únicamente si se
comprende lo anterior, él ingeniero puede aplicar correctamente las técnicas de ingeniería
disponibles y decidir las especificaciones precisas de las instalaciones de producción que se
deben usar para cualquier pozo.
2
Resumen
En la actualidad, se puede aplicar la técnica del análisis nodal utilizando simuladores del
comportamiento de flujo en todo el sistema integral de producción (SIP), el principal beneficio del
uso de estos simuladores es que se pueden generar distintos escenarios con los que se puede
observar y determinar cuáles son las mejores condiciones de producción de un pozo, un sistema
multipozos o hasta un campo petrolero. Con los simuladores se puede generar un análisis de
sensibilidad para diferentes condiciones del sistema integral de producción (SIP), desde la
variación de propiedades de los fluidos producidos, variación de los diámetros de la tubería de
producción (T.P), línea de descarga (L.D) y estranguladores, variación de algunas condiciones del
yacimiento, etc.
Otro beneficio del uso de los simuladores es la disminución del tiempo de realización de un
análisis nodal, ya que como se verá en esta tesis, la elaboración de las curvas de inflow y outflow
manualmente requiere de mucho tiempo, como en el caso del cálculo de las caídas de presión en
tuberías para flujo multifásico debido a que este es un proceso iterativo.
En general, el uso de los simuladores es de mucha utilidad para los ingenieros de producción,
pero éstos deben de comprender que los resultados arrojados por el simulador estarán en función
de la calidad de los datos de entrada.
3
Capítulo I
CAPITULO 1
COMPORTAMIENTO DE FLUJO EN EL YACIMIENTO
El yacimiento es el espacio físico en donde se encuentran almacenados los fluidos petroleros

(aceite, gas y agua), una vez que se perfora el primer pozo los fluidos que se encontraban en un
estado estático cambian a un estado dinámico y comienzan a fluir del yacimiento hacia el pozo.
Es por esta razón que un primer paso para la optimización de la producción es analizar el
comportamiento de flujo del yacimiento hacia al pozo.
El comportamiento de flujo desde el yacimiento hasta el pozo, depende de factores como el tipo de
mecanismo de producción del yacimiento y de variables como; la presión media del yacimiento, la
permeabilidad, la saturación de fluidos etc. Además de otros factores que no son parte de las
características petrofísicas del yacimiento tales como; daño a la formación y tipo de terminación
del pozo.
En este capítulo se verán algunas herramientas analíticas útiles para el estudio del
comportamiento de flujo en el yacimiento, como son: Índice de productividad que es expresado
mediante la letra (J) y la curva de IPR (Inflow performance relationship). El IPR y el índice de
productividad (J) no son equivalentes, el IPR es la relación que existe entre el gasto producido (q)
y la presión de fondo fluyendo (Pwf), mientras que el índice de productividad (J) es la relación entre
el gasto producido (q) y la caída de presión ∆P (∆P =Pr − Pwf).
Donde:
Pr =Presión media del yacimiento
Pwf =Presión de fondo fluyendo
La elaboración de la curvas de IPR para pozos de aceite y gas son extremadamente importantes
en el análisis del sistema integral de producción, ya que si no se cuenta con alguna idea de la
capacidad productiva de un pozo, el diseño y la optimización del aparejo de producción llega a ser
mucho más complicado. Debemos tener siempre en cuenta que mientras más información
disponible se tenga de las propiedades del sistema roca-fluidos, los resultados de las curvas de
IPR se aproximarán más a la realidad.
En la actualidad se cuenta con buenos simuladores de yacimientos que pueden pronosticar de

manera más exacta la capacidad productiva de un pozo, si se tiene acceso a simuladores, éstos
deben de utilizarse en lugar de los procedimientos analíticos que se verán en este capítulo, pero
como no siempre se tiene acceso a los simuladores, las soluciones analíticas son herramientas
útiles y económicas que pueden dar resultados muy aproximados a lo real.
4
Capítulo I
El comportamiento de flujo en el yacimiento empezó con el desarrollo de una ecuación que nos
permitiera relacionar las dos variables que son de principal importancia en la optimización de la
producción, las cuales son: Presión y gasto. La ecuación que nos relaciona estas dos variables
es la ecuación de Darcy, esta ecuación es la base fundamental para entender el comportamiento
de flujo en el yacimiento.
I.1 Flujo de aceite
I.1.1 Ley de Darcy
La ley de Darcy es una ecuación general que indica el comportamiento de flujo en un medio
poroso y permeable, es utilizada para la predicción de flujo del yacimiento al pozo. Partiendo de la
ecuación general de Darcy para flujo lineal, sin considerar las fuerzas gravitacionales ni capilares,
tenemos la siguiente ecuación.
k dP
q=A × ---------------1.1
μ dL
Donde:
q =Gasto del fluido (cm3 / s)

A =Área expuesta al flujo (cm2)
k =Permeabilidad de la roca (darcys)
μ =Viscosidad del fluido (cp)
dP =Diferencial de Presión (atm)

dL =Diferencial de Longitud (cm)
I.1.2 Ecuaciones de flujo
Haciendo el desarrollo de la ecuación de Darcy (ecuación 1.1) para flujo radial y en unidades de
campo obtenemos lo siguiente:
7.08 × 10 −3 (k oh)(Pr − Pwf )

qo = (Estado pseudo-estacionario) ---------------1.2
⎡ ⎛ 0.472 re ⎞ ⎤
μo Bo ⎢In ⎜⎜ ⎟⎟ + S⎥
⎣⎢ ⎝ rw ⎠ ⎦⎥
7.08 × 10 −3 (k oh)(Pr − Pwf )

qo = (Estado estacionario) ---------------1.3
⎡ ⎛ re ⎞ ⎤
μo Bo ⎢In ⎜⎜ ⎟⎟ + S⎥
⎢⎣ ⎝ rw ⎠ ⎥⎦
5
Capítulo I
Donde:
qo =Gasto de aceite (bpd @ c.s)
k o =Permeabilidad efectiva al aceite (md)
h =Espesor de la zona productora (pie)

re =Radio de drene (pie)
rw =Radio del pozo (pie)
S =Daño total
Pr =Presión media del yacimiento (lb / pg2)
Pwf =Presión de fondo fluyendo a la altura media de los disparos (lb / pg2)
μo =Viscosidad del aceite (cp) medido a una presión media (Pr + Pwf ) / 2
Bo =Factor de volumen del aceite medido a una presión media
Con las ecuaciones anteriores podemos obtener el gasto de aceite que está aportando el
yacimiento, pero necesitamos conocer los valores de las propiedades del sistema roca-fluidos, el
radio de drene y pozo, así como la caída de presión en el yacimiento.
I.2 ÍNDICE DE PRODUCTIVIDAD Y EL IPR
I.2.1 Definición del índice de productividad (J)
Un parámetro del yacimiento comúnmente usado es el índice de productividad (J), el cual

representa una relación gasto-∆P, donde J es definido como sigue:
qo
J= ---------------1.4
(Pr − Pwf )
Esta ecuación relaciona el gasto aportado por un pozo dada una caída de presión, desde la
presión media del yacimiento hasta la presión de fondo fluyendo. Observemos la figura 1.1, la
cual nos indica el comportamiento del índice de productividad (J) con respecto a la producción
acumulada para diferentes mecanismos de producción.
6
Capítulo I
Figura 1.1 J vs. Np para diferentes mecanismos de producción
Tipos de mecanismos naturales de producción del yacimiento.
Expansión del sistema roca-fluidos.

Este mecanismo de producción se debe a que tanto la roca como el fluido contenido dentro de ella
se expanden cuando se genera una caída de presión en el yacimiento.
Expansión por gas disuelto.

Este mecanismo de producción se debe a que el gas que está disuelto en el aceite se libera y
debido a su alto valor de compresibilidad se expande ocupando el volumen del aceite, es
importante decir que el gas que se está liberando es producido junto con el aceite, es decir, el gas
liberado no forma un casquete de gas.
Empuje de acuífero activo.

Este mecanismo de producción se debe al empuje hidráulico de un acuífero activo.
Expansión del casquete de gas.

Este mecanismo de producción se debe al que el gas libre en el yacimiento (casquete de gas) se
expande desplazando del espacio poroso al aceite.
Segregación Gravitacional.
Este mecanismo de producción se debe, como su nombre lo dice, a la segregación de las dos
fases: Aceite y gas debido a la densidad de cada uno. Este mecanismo de producción se da
cuando el gas que se libera del aceite no es producido, si no que se incorpora al casquete de gas.
En el diseño de sistemas artificiales de producción, el índice de productividad es expresado en bpd

de líquido total, es decir (qo + qw) y es definido como sigue:
7
Capítulo I
qo + qw
J= ---------------1.5
Pr − Pwf
Donde:
qo =Gasto de aceite (bpd @ c.s)

qw =Gasto de agua (bpd @ c.s)
Otro término a veces usado es el índice de productividad específico, este es sólo el valor de J
dividido por los pies de espesor neto de la zona productora y es comúnmente usado para
comparar diferentes pozos de un mismo campo.
J ⎡ bpd ⎤
Js = =⎢ ⎥ ---------------1.6
h ⎣ (lb / pg ) ⋅ pie ⎦
2
I.2.2 Definición de IPR
El IPR (Inflow Performance Relationship), es una relación del comportamiento del gasto de líquido
o gas respecto a la presión de fondo fluyendo (q vs. Pwf). Cuando hablamos de la curva de IPR
debemos tomar en cuenta el tipo de mecanismo de producción del yacimiento, la forma de la
curva, así como también recordar que el índice de productividad (J) de un pozo es muy probable
que cambie con el tiempo y la producción acumulada.
I.2.3 Índice de productividad para yacimientos bajo saturados
El índice de productividad (J) puede ser constante o variable, esto dependerá si la presión de
fondo fluyendo se encuentra por arriba o por debajo de la presión de burbuja. Un yacimiento de
aceite bajo saturado se tiene cuando la Pwf > Pb y la curva de IPR que representa el
comportamiento del flujo en el yacimiento (relación entre q y Pwf), se comporta como en la
siguiente gráfica (figura 1.2).
Figura 1.2 Curva de IPR para yacimiento bajo saturado
8
Capítulo I
OB
Donde Tan θ = =J
OA
Otra forma de obtener el índice de productividad (J) cuando Pwf > Pb, teniendo tanto gasto de
aceite como de agua es a partir de la ecuación de Darcy, obteniendo lo siguiente:
7.08 × 10 −3 h ⎡ k o kw ⎤
J= ⎢ + ⎥ ---------------1.7
In (re / rw ) ⎣ Bo μ o B w μ w ⎦
Ejemplo:
Tenemos un yacimiento de aceite bajo saturado; de una prueba de producción realizada a un pozo
se obtuvieron los siguientes datos:
Pr = 2,800 (lb / pg2 )

qo = 200 bpd
Pwf = 2,400 (lb / pg2 )
Calcular el índice de productividad (J) de dicho pozo:
qo 200 ⎡ bpd ⎤
J= = = 0 .5 ⎢ 2 ⎥
(Pr − Pwf ) (2,800 − 2,400 ) ⎣ (lb / pg ) ⎦
El valor de J =0.5 se mantendrá constante hasta que el valor de Pwf alcance el valor de Pb, por lo
tanto con la ecuación anterior podemos calcular el gasto para cualquier Pwf mientras se cumpla la
condición mencionada, como se hace a continuación:
Calcular el gasto de aceite para una Pwf de 2,000 (lb / pg2).
qo = J (Pr − Pwf ) = 0.5 (2,800 − 2,000 ) = 400 bpd
I.2.4 Índice de productividad para yacimientos saturados
En el caso de que Pwf < Pb, la curva de IPR se comporta como en la figura 1.3.
− dq
Donde J = , el signo negativo indica un decremento de J para un incremento de gasto.
dPwf
9
Capítulo I
Figura 1.3 Curva de IPR para yacimientos saturados
Para un yacimiento de aceite saturado se tienen dos fases moviéndose dentro del yacimiento
(aceite y gas) en el caso dado de que la saturación de agua (Sw) sea menor o igual a la irreducible.
Como se puede observar en la figura 1.3, el índice de productividad no es constante a lo largo del
tiempo de explotación de este tipo de yacimientos, por lo que no podemos recurrir a la
ecuación 1.4 para pronosticar gastos futuros de aceite (qo) para diferentes Pwf. Muchos
investigadores se dieron a la tarea de obtener ecuaciones que nos permitieran resolver este
problema hasta que en 1968 Vogel propuso su ecuación, la cual es una de las más utilizadas en la
industria petrolera debido a su simplicidad y eficacia.
I.3 Flujo de dos fases: Aceite y gas
Si nos regresamos a la ecuación de Darcy para flujo de líquido con frontera cerrada, asumiendo
que conocemos el valor de la presión media del yacimiento (Pr) y que el valor de S =0, tenemos la
siguiente ecuación:
7.08 × 10 −3 kh
Pr
k ro
qo = ∫
⎡ ⎛ 0.472 re ⎞⎤ Pwf μo Bo
dp ---------------1.8
⎢In ⎜⎜ ⎟⎟⎥
⎢⎣ ⎝ rw ⎠⎥⎦
k ro
El término es una función de la presión y a su vez el valor de la permeabilidad relativa al
μ o Bo
aceite ( k ro ) es una función de la saturación de aceite ( So ). Una gráfica típica de esto se muestra a
continuación (figura1.4).
10
Capítulo I
k ro
Figura 1.4 vs. Presión
μ o Bo
La integral puede ser en dividida en dos partes, como se muestra a continuación:
Pr Pr Pb
k ro 1 k ro
∫
Pwf μ o Bo
dp = ∫
Pb μ o Bo
dp + ∫
Pwf μ o Bo
dp
Aunque estrictamente de la presión media del yacimiento (Pr) a la presión de burbuja (Pb) los
valores de μo y Bo no son constantes, su variación es mínima con respecto a la variación de
presión, por lo que la gráfica de (q vs. Pwf) es una línea recta para esta etapa (yacimiento de
aceite bajo saturado). Para la segunda etapa (yacimiento de aceite saturado) vemos que está
involucrado el valor de permeabilidad relativa al aceite (kro). Al caer la Pwf por debajo de la de Pb
se comienza a liberar gas con lo que la permeabilidad relativa a éste aumenta y la del aceite
comienza a disminuir, además que los valores de μo y Bo son más significativos porque al liberarse
los componentes ligeros del aceite, el valor de μo empieza a aumentar y el del Bo a disminuir.
Estas son las principales causas por las cuales las curvas de IPR para yacimientos de aceite bajo
saturado y saturado se comportan de diferente manera.
I.3.1 Ecuación de Vogel
Una solución simplificada para el flujo de dos fases en el yacimiento fue propuesta por Vogel,
quien en 1968 dio la siguiente ecuación general para determinar la curva de IPR para yacimientos
con empuje de gas disuelto y flujo por debajo de la Pb.
2
qo ⎛P ⎞ ⎛P ⎞
= 1 − 0.2 ⎜⎜ wf ⎟⎟ − 0.8 ⎜⎜ wf ⎟⎟ ---------------1.9
qo max ⎝ Pr ⎠ ⎝ Pr ⎠
11
Capítulo I
La curva general del trabajo de Vogel se muestra a continuación:
Figura 1.5 Curva de IPR para yacimientos con empuje de gas disuelto (Vogel)
De esta gráfica con datos dados de una prueba de producción qo y Pwf y conociendo la Pr
podemos obtener el gasto máximo que aportaría la formación, este es un gasto teórico porque se
obtiene cuando Pwf =0 y en la práctica esta condición no se presenta sin que antes el pozo deje
de fluir. Con el valor del gasto máximo podemos calcular los diferentes gastos esperados para
diferentes Pwf recordando que esta gráfica es para una eficiencia de flujo igual a uno (EF =1).
Para aplicar esta ecuación se deben tomar en cuenta las siguientes observaciones:
1) Este procedimiento no debe ser considerado correcto para otro tipo de mecanismo de
producción, sólo para expansión por gas disuelto. Aunque se ha visto que para otros
tipos de mecanismos puede dar resultados aceptables.
2) Este procedimiento es para dos fases: Aceite y gas únicamente, no debe ser considerado
valido cuando tres fases están fluyendo: Aceite, gas y agua.
12
Capítulo I
3) La gráfica de Vogel sólo puede ser utilizada para pozos sin daño, es decir cuando la
eficiencia de flujo EF =1, por lo que si el pozo está dañado o estimulado debemos usar otra
gráfica.
Si se requiere, la ecuación de Vogel puede ser modificada para encontrar el valor de Pwf dado un
gasto determinado:
[
Pwf = 0.125 Pr − 1 + 81 − 80 (qo / qo max ) ] ---------------1.10
La ecuación de Vogel, es la más utilizada en la industria, ya que sólo es necesario un par de datos
de una prueba de producción (qo y Pwf) para poder aplicar esta correlación con muy buenos
resultados dada las restricciones anteriores.
Ejemplo:
Se tienen los siguientes datos:
Pr = 2,400 (lb / pg2 )

qo = 100 bpd
Pwf = 1,800 (lb / pg2 )
Determinar el qo para una Pwf de 800 (lb / pg2).
Solución:
1) Determinar el gasto máximo aportado por dicho pozo utilizando la ecuación de Vogel:
qo 100
qo max = 2
= 2
= 250 bpd
⎛P ⎞ ⎛P ⎞ ⎛ 1,800 ⎞ ⎛ 1,800 ⎞
1 − 0.2 ⎜⎜ wf ⎟⎟ − 0.8 ⎜⎜ wf ⎟⎟ 1 − 0. 2 ⎜ ⎟ − 0 .8 ⎜ ⎟
⎝ Pr ⎠ ⎝ Pr ⎠ ⎝ 2,400 ⎠ ⎝ 2,400 ⎠
2) Obteniendo el valor de qo para la Pwf de 800 (lb / pg2).
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤ ⎡ ⎛ 800 ⎞ ⎛ 800 ⎞ ⎤
2
qo = qo max ⎢1 − 0.2 ⎜⎜ wf ⎟⎟ − 0.8 ⎜⎜ wf ⎟⎟ ⎥ = 250 ⎢1 − 0.2 ⎜ ⎟ − 0. 8 ⎜ ⎟ ⎥ = 211 bpd

⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦ ⎣⎢ ⎝ 2,400 ⎠ ⎝ 2,400 ⎠ ⎦⎥
Este procedimiento también se puede realizar gráficamente debido a que conocemos el valor de la
relación (Pwf / Pr) y el valor de qo. Para nuestro ejemplo (Pwf / Pr) =0.75, si partimos de este valor
en el eje de las ordenadas hasta intersectar la curva y leemos el valor en el eje de las abscisas
observamos que (qo / qomax) =0.40.
13
Capítulo I
Por lo tanto:
qo 100
qo max = = = 250 bpd
0.40 0.40
Para obtener gráficamente el qo para una Pwf de 800 (lb / pg2), primero se determina el valor de la
relación de presiones (Pwf / Pr)
Pwf 800
= = 0.333
Pr 2,400
De la gráfica se observa que para el valor de (Pwf / Pr) =0.333, se tiene un valor de
(qo / qomax) =0.84, por lo que:
qo = 0.84qo max = 0.84 × 250 = 210 bpd
Este valor difiere por un barril de la solución analítica debido a la imprecisión inherente a la lectura
de valores en la gráfica. Los puntos citados se muestran en la siguiente gráfica:
I.3.2 Procedimiento de Standing para EF ≠ 1
Standing extendió el trabajo de Vogel para aquellos casos donde la EF ≠ 1. La eficiencia de flujo
(EF) puede ser definida de la siguiente forma.
14
Capítulo I
Pr − Pwf*
EF = ---------------1.11
Pr − Pwf
Donde:
Pwf* =Presión de fondo fluyendo si el pozo no tuviera daño
Pwf =Presión de fondo fluyendo del pozo
Despejando el valor de Pwf* tenemos:
Pwf = Pr − (Pr − Pwf ) EF

*
---------------1.12
Standing presentó la figura 1.6 para usar eficiencias de flujo (EF) entre valores de 0.5 hasta 1.5.
Un paso necesario para el uso de esta gráfica es determinar el qomax para una EF =1, después
cualquier gasto puede ser determinado para cualquier EF. A continuación se muestra un ejemplo
para mostrar el procedimiento de cómo se utiliza la gráfica de Standing.
Figura 1.6 Correlación de Standing para pozos con valores de EF ≠ 1
15
Capítulo I
Ejemplo:
Tenemos los siguientes datos:
Pr = 2,600 (lb / pg2 ) < Pb

qo = 500 bpd para una Pwf = 1,800 (lb / pg2 )
EF = 0.6
Determinar:
1) qomax para una EF =0.6

2) qo para una Pwf =1,300 (lb / pg2) y una EF =0.6
Solución:
1) Determinar el qomax para una EF =1. De los datos dados obtener el valor de (Pwf / Pr)
Pwf 1,800
= = 0.692
Pr 2,600
De la figura 1.6 leemos en la curva para EF =0.6 que [qo / (qomax EF=1)] =0.3, de este valor
despejamos el valor de (qomax EF=1) y obtenemos lo siguiente:
qo 500
(qo max EF = 1) = = = 1,667 bpd
0 . 3 0 .3
2) Determinar qomax para una EF =0.6. De la figura 1.6 leemos el valor de [qo / qomax EF=1] para la
curva de EF =0.6 cuando (Pwf / Pr) =0. Donde qo =(qomax EF=0.6)
⎡ qo max EF = 0.6 ⎤
⎢ ⎥ = 0.77
⎣ qo max EF = 1 ⎦
(qo max EF = 0.6) = 0.77 × (qo max EF = 1) = 0.77 × 1,667
(qo max EF = 0.6) = 1,283 bpd
3) Determinar qo para una Pwf =1,300 (lb / pg2) y una EF =0.6. De los datos dados obtener el valor
de (Pwf / Pr)
16
Capítulo I
Pwf 1,300
= = 0. 5
Pr 2,600
De la figura 1.6 leemos en la curva de EF =0.6 que [qo / qomax EF=1] =0.462, de este valor
despejamos el valor de qo y obtenemos lo siguiente:
qo = 0.462 × 1,667 = 770 bpd
Existen ciertas condiciones de flujo en la que no podemos utilizar la ecuación de Vogel o la gráfica
de Standing. Estas condiciones se dan cuando tenemos altos valores de eficiencia de flujo
(EF >1.5) y bajos valores de Pwf, un ejemplo de esto es el siguiente:
Ejemplo:
Pr = 2,000 (lb / pg2 )

Pwf = 500 (lb / pg2 )
EF = 2
Solución:
Pwf = Pr − (Pr − Pwf ) EF

*
Pwf = 2,000 − (2,000 − 500 ) 2 = −1,000 (lb / pg2 )

*
Estos valores negativos causan conflictos en la ecuación de Vogel y mostrarán una reducción del
gasto. Un resultado típico de esto es mostrado en la figura 1.7
Figura1.7 Error en la extrapolación del trabajo de Standing
17
Capítulo I
I.3.3 Ecuación de Fetkovich
Una ecuación más exacta para el comportamiento de un yacimiento saturado, no importando altos
valores de EF y bajos valores de Pwf fue propuesta por Fetkovich(15) y es la siguiente:
(
qo = J'o Pr2 − Pwf2 ) n
---------------1.13
Donde:
J'o =Índice de productividad (coeficiente de la curva de contrapresión)
n =Exponente de la curva de contrapresión
Fetkovich presentó un método de análisis de pruebas multigastos (isocronales y pruebas de flujo

continuo) aplicado a pozos de aceite. Fetkovich encontró que los pozos produciendo por debajo de
la Pb deben realmente comportarse más como pozos de gas, esto significa que (Pr2 − Pwf2) vs. qo
dará una línea recta en una gráfica Log-Log como se muestra en la figura 1.8.
Figura 1.8 Prueba multigastos aplicada a un pozo de aceite (yacimiento de aceite saturado)
Esta curva también es llamada “curva de contra presión”. En donde J’o es la intersección en el
eje de las abscisas (qo) cuando (Pr2 − Pwf2) =1 y n =(1 / m), donde m es la pendiente de la recta
que puede ser encontrada con dos puntos. Al igual que la ecuación de Vogel, la ecuación de
Fetkovich es únicamente para yacimientos de aceite saturado, es decir que la Pr sea menor o igual
que la Pb (Pr ≤ Pb).
18
Capítulo I
Con la ecuación de Fetkovich, Harrison realizó una gráfica general para ser usada cuando los
valores de EF están fuera de rango de la gráfica presentada por Standing, la gráfica es mostrada a
continuación (figura 1.9).
Figura 1.9 Amplificación del trabajo de Standing incluyendo otros valores de EF (Harrison)
Hasta este momento sólo hemos visto dos ecuaciones generales que describen el comportamiento
de flujo para dos fases. Sin embrago, existen otras ecuaciones para flujo de dos fases para
yacimientos cuyo mecanismo de producción es por expansión de gas disuelto y diferentes
condiciones de flujo. A continuación se hace una recopilación de algunas de ellas.
Vogel
2
qo ⎛P ⎞ ⎛P ⎞
= 1 − 0.2 ⎜⎜ wf ⎟⎟ − 0.8 ⎜⎜ wf ⎟⎟
19
Capítulo I
Fetkovich
n
⎡ ⎛P ⎞
2
⎤
(
o
qo
qo max
= ⎢1 − ⎜⎜ wf ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎥⎦
I.3.4 Ecuación de Jones, Blount y Glaze
Jones, Blount y Glaze(13) propusieron la siguiente relación entre presión y gasto.
Pr − Pwf
= C + Dqo
qo
Pr − Pwf
En un plano cartesiano graficamos vs. qo lo cual nos dará una línea recta con ordenada al
qo
origen C y pendiente D, el termino C representa el coeficiente de flujo laminar y D es el coeficiente
de turbulencia; mínimo es necesario conocer dos puntos de la recta para poder determinar los
coeficientes mencionados. Una vez que se obtienen los valores de los coeficientes, el gasto puede
ser determinado a cualquier Pwf siguiendo la siguiente relación.
−C+ C2 + 4D (Pr − Pwf )

qo =
2D
I.3.5 Ecuación de Klins y Majcher
Klins y Majcher(13) desarrollaron un IPR incorporando el valor de Pb.
d
qo ⎛P ⎞ ⎛P ⎞
= 1 − 0.295 ⎜⎜ wf ⎟⎟ − 0.705 ⎜⎜ wf ⎟⎟
Donde:
⎛ P ⎞
d = ⎜⎜ 0.28 + 0.72 r ⎟⎟ (1.235 + 0.001Pb )
⎝ Pb ⎠
I.3.6 Ecuación de Sukarno
Sukarno(13) desarrolló una ecuación de IPR basado en los resultados de una simulación que
pretende tomar en cuenta la eficiencia de flujo (EF), debido a que el gasto depende tanto del daño
como del cambio de Pwf. Sukarno desarrolló la siguiente ecuación:
20
Capítulo I
qo ⎡ ⎛P ⎞ ⎛P ⎞
2
⎛P ⎞
3
⎤
= EF ⎢1 − 0.1489 ⎜⎜ wf ⎟⎟ − 0.4418 ⎜⎜ wf ⎟⎟ − 0.4093 ⎜⎜ wf ⎟⎟ ⎥
qomax@S = 0 ⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
Donde:
2 3
⎛P ⎞ ⎛P ⎞ ⎛P ⎞
EF = a0 + a1 ⎜⎜ wf ⎟⎟ + a2 ⎜⎜ wf ⎟⎟ + a3 ⎜⎜ wf ⎟⎟ y an = b0 + b1S + b2S2 + b3S3
⎝ Pr ⎠ ⎝ Pr ⎠ ⎝ Pr ⎠
S es el factor de daño y los valores de b necesarios para calcular las an son dados en la siguiente
tabla:
Constantes para IPR de Sukarno
b0 b1 b2 b3
a0 1.0394 0.12657 0.0135 -0.00062
a1 0.01668 -0.00385 0.00217 -0.0001
a2 -0.858 0.00201 -0.00456 0.0002
a3 0.00952 -0.00391 0.0019 -0.00001
I.4 Flujo combinado: Una fase (pwf > pb) y dos fases (pwf < pb)
Aunque estrictamente por arriba de la Pb los valores de μo y Bo no son constantes, su variación es

mínima con respecto a la variación de presión, por lo que cuando Pwf > Pb la curva de IPR es un
línea recta, cuando se llega al punto en que la presión del yacimiento Pr ≤ Pb tenemos que utilizar
la ecuación de Vogel u otras ecuaciones para flujo de dos fases. La ecuación de Vogel parte de
que la Pr ≤ Pb y por lo tanto si tenemos un yacimiento cuya Pr > Pb pero Pwf < Pb , tendremos que
hacer el análisis tanto para flujo en una sola fase, como para flujo en dos fases para obtener el
gasto que aporta dicho pozo dada una Pwf.
I.4.1 Ecuación General
Eickmeier y Neeley encontraron la siguiente correlación:
Si partimos de la ecuación general de Vogel.
2
qo ⎛P ⎞ ⎛P ⎞
= 1 − 0.2 ⎜⎜ wf ⎟⎟ − 0.8 ⎜⎜ wf ⎟⎟ ---------------1.14
Y tomando la derivada respecto a Pwf tenemos:
dq ⎡ 0.2 1.6 Pwf ⎤

− = qo max ⎢ + 2 ⎥ ---------------1.15
dp wf P
⎣ r Pr ⎦
21
Capítulo I
Recordando que para utilizar la ecuación de Vogel partimos de que Pr ≤ Pb, si hacemos el análisis
justo en el punto donde Pr = Pb y tomando la pendiente o el valor de J cuando Pwf →Pb se obtiene
lo siguiente:
1.8 qo max
J= ---------------1.16
Pb
Por lo tanto el gasto máximo aportado por la formación será la suma de los gastos máximos de las
2 etapas (bajo saturado y saturado) como se muestra a continuación (figura 1.10).
bajo saturado saturado

JP
qb = J (Pr − Pb ) ---------------1.17 q omax = b ---------------1.18
VOGEL 1.8
Figura 1.10 Combinación de los casos J constante y comportamiento de Vogel
Sumando las dos ecuaciones obtenemos el gasto máximo para dicho yacimiento.
J Pb
qo max = qb + ---------------1.19
1 .8
Por lo tanto, si queremos obtener un gasto a una Pwf < Pb, será la suma de los gastos de las dos
etapas qo= qb +qvogel y tenemos la siguiente ecuación:
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤
qo = qb + [qo max − qb ] ⎢1 − 0.2 ⎜⎜ wf ⎟⎟ − 0.8 ⎜⎜ wf ⎟⎟ ⎥ ---------------1.20
⎢⎣ ⎝ Pb ⎠ ⎝ Pb ⎠ ⎥⎦
22
Capítulo I
Si conocemos el valor de J en el punto donde Pwf = Pb podemos determinar el gasto máximo

(qo max) que aportaría el pozo. Si la prueba es tomada por debajo de la presión de burbuja,
encontramos el valor de J de la siguiente ecuación:
qo
J= ---------------1.21
Pb ⎡ ⎛ Pwf ⎞ ⎤
2
⎛ Pwf ⎞
Pr − Pb + ⎢1 − 0.2 ⎜⎜ ⎟ − 0.8 ⎜⎜ ⎟ ⎥
1. 8 ⎢ ⎝ Pb ⎟⎠ ⎝ Pb ⎟⎠ ⎥
⎣ ⎦
Si la prueba es tomada cuando Pwf > Pb, utilizamos la ecuación 1.17 o podemos utilizar la ecuación
de Darcy si contamos con los datos necesarios.
I.4.2 Ecuación de Hasan
Hasan extendió la ecuación de Vogel para la combinación de flujo en una y dos fases, la ecuación
a la que llegó es la siguiente.
J Pb ⎡ ⎛P ⎞ ⎛ Pwf ⎞ ⎛P ⎞
2
⎤
= ⎢1 . 8 ⎜ r
⎟ − 0 . 8 − 0 . 2 ⎜⎜ ⎟⎟ − 0.8 ⎜⎜ wf ⎟⎟ ⎥ ---------------1.22
J* 1.8 (Pr − Pwf ) ⎢ ⎜⎝ Pb ⎟⎠ ⎝ Pb ⎠ ⎝ Pb ⎠ ⎥⎦
⎣
Donde:
Pr =Presión media del yacimiento (Por arriba de la Pb)
J =Índice de Productividad a condiciones de flujo por debajo de la Pb

J * =Índice de Productividad cuando Pwf →Pb
El valor de J* se puede obtener con la ecuación de Darcy o de la ecuación 1.18. Con la ecuación
1.22, podemos conocer el índice de productividad J cuando se tiene flujo de dos fases y con este
qo
valor recurriendo a la ecuación general de índice de productividad J = , podemos obtener
(Pr − Pwf )
el gasto futuro (qo) para una determinada presión de fondo fluyendo (Pwf).
I.4.3 Ecuación de Fetkovich
Fetkovich propuso otra ecuación para flujo combinado, la cual es la siguiente:
(
+ J (Pr − Pb ) ---------------1.23
Con las ecuaciones anteriores podemos calcular los gastos para flujo combinado (una y dos fases
en el yacimiento), es decir cuando Pwf < Pb, pero Pr > Pb.
23
Capítulo I
I.5 Curvas de IPR para pozos de gas
I.5.1 Ecuaciones de flujo
De la ley de Darcy para flujo radial tenemos la siguiente relación:
k (2πrh) ⎛ dp ⎞
q= ⎜ ⎟ ---------------1.24
μ ⎝ dr ⎠
La ecuación anterior es una ecuación diferencial y debe ser integrada para su aplicación. Antes de
la integración, la ecuación de flujo debe ser combinada con una ecuación de estado y una
ecuación de continuidad para ser aplicable al flujo de gas. La ecuación de continuidad es la
siguiente:
ρ1q1 = ρ2q2 = Constante ---------------1.25
Donde:
ρ =Densidad del gas

q =Gasto del gas
Y la ecuación de estado es la ecuación de los gases reales y es la siguiente:
pM
ρ= ---------------1.26
ZRT
Combinando las ecuaciones de estado y continuidad con la ley de Darcy obtenemos el

comportamiento de flujo radial para pozos de gas. Para unidades convencionales de campo la
ecuación queda de la siguiente manera:
qg =
(
703 × 10 −6 kh Pr2 − Pwf2 ) (Estado estacionario) ---------------1.27
T ⋅ μg ⋅ Z ⋅ In (re / rw )
qg =
(
703 × 10 −6 kh Pr2 − Pwf2 ) (Estado pseudo-estacionario) ---------------1.28
T ⋅ μg ⋅ Z ⋅ In (0.472 re / rw )
Donde:
qg =Gasto de Gas (Mpcd)
k =Permeabilidad en (md)
24
Capítulo I
h =Espesor de la formación en (pie)

T =Temperatura del yacimiento en (°R)
Z =Factor de compresibilidad del gas medido a una presión media
re =Radio de drene (pie)
rw =Radio del pozo (pie)
Pr =Presión en el radio de drene en (lb / pg2)
Pwf =Presión de fondo fluyendo en (lb / pg2)
μg =Viscosidad del gas en (cp) medido a una presión media
Las ecuaciones de flujo para estado estacionario y pseudo-estacionario incluyendo el factor de

daño y el efecto de turbulencia son:
qg =
(
703 × 10 −6 kh Pr2 − Pwf2 )
(Estado estacionario) ---------------1.29
T ⋅ μ ⋅ Z [In (re / rw ) + S + Dqg ]
qg =
(
703 × 10 −6 kh Pr2 − Pwf2 ) (Estado pseudo-estacionario) ---------------1.30
T ⋅ μ ⋅ Z [In (0.472 re /rw ) + S + Dqg ]
Además.
2.22 × 10 −15 ⋅ γ g ⋅ k 2.33 × 1010

D= β y β=
μ ⋅ h ⋅ rw k1.2
Donde:
S =Factor de daño (adimensional)

D =Coeficiente de turbulencia
γ g =Densidad relativa del gas
I.5.2 Pruebas de potencial (pruebas multigastos)
El potencial o capacidad de flujo de un pozo de gas, es requerido para la planeación de las

operaciones de cualquier campo de gas. La capacidad de flujo debe ser determinada para
diferentes presiones de fondo fluyendo en cualquier momento de la vida del yacimiento, también
debe ser considerado el cambio en la capacidad de flujo con respecto al cambio en la presión
media del yacimiento.
25
Capítulo I
El método más común para determinar la potencialidad de un pozo de gas son las llamadas
pruebas multigastos, en las cuales un pozo es producido a varios gastos (usualmente 4) y se mide
ese gasto (qg) y la presión de fondo fluyendo correspondiente (Pwf). Con estas pruebas una
ecuación de IPR de gas puede ser desarrollada. Hay básicamente dos tipos de pruebas que
pueden ser realizadas. Una prueba es llamada de flujo continuo y la otra es llamada prueba
isocronal. La teoría detrás de las pruebas está basada en ecuaciones de flujo en estado pseudo-
estacionario.
La ecuación más utilizada en la industria petrolera para encontrar la curva de IPR para pozos de
gas es la siguiente, esta puede ser encontrada a partir de las pruebas antes mencionadas.
(
qg = C Pr2 − Pwf2 ) n
---------------1.31
Donde:
qg =Gasto de gas en unidades consistentes respecto a C.
C =Una constante que depende de las características del pozo y de las propiedades del gas, el
valor de C también puede ser obtenido analíticamente con la siguiente ecuación:
703 × 10 −6 kh
C=
⎡ ⎛ 0.472 re ⎞ ⎤
T ⋅ μg ⋅ Z ⎢In ⎜⎜ ⎟⎟ + S⎥
⎢⎣ ⎝ rw ⎠ ⎥⎦
En la ecuación 1.31, n es el coeficiente de turbulencia, el cual tiene un valor que varia entre
(0.5 y 1). Para valores cercanos a 0.5 indica flujo turbulento, y para valores cercanos a 1 indica
flujo laminar.
Los valores de C y n son determinados de las pruebas de pozo. Al menos dos pruebas son
requeridas debido a que hay dos variables desconocidas en la ecuación, cuatro pruebas son
recomendadas para minimizar los efectos de error por medición.
Si más de dos pruebas están disponibles, los datos pueden ser graficados en una gráfica
Log-Log y los datos serán ajustados a una línea recta para determinar los valores de C y n. Sobre
la gráfica Log-Log graficamos (Pr2 − Pwf2) vs. qg, 1 / n es la pendiente de la línea graficada y el
valor de C =qg cuando (Pr2 − Pwf2) =1. Una vez que el valor de n ha sido determinado de la gráfica,
el valor de C también puede ser calculado usando los datos de una de las pruebas que este sobre
la línea recta graficada, esto es:
qg
C=
(Pr
2
− Pwf2 ) n
26
Capítulo I
Ejemplo:
De una prueba de potencial aplicada a un pozo de gas se obtuvieron los siguientes datos:
Prueba qg [Mpcd] Pwf (lb / pg2) (Pr^2 − Pwf^2) x 10^ 3 [(lb / pg2)^2]
0 201 (Pr) 40.4 (Pr^2)
1 2,730 196 1.985
2 3,970 195 2.376
3 4,440 193 3.152
4 5,550 190 4.301
Si graficamos los datos en una gráfica Log-Log se obtiene la siguiente gráfica.
Como en la gráfica se observa que las pruebas 1 y 4 caen dentro de la línea recta, podemos usar
estos puntos para obtener el valor de n, el procedimiento es el siguiente:
log q1 − log q4
n=
log Δ(P2 )1 − log Δ(P2 )4
log (2,730 ) − log (5,550 )

n= = 0.92
log (1.985 × 10 3 ) − log ( 4.301× 10 3 )
Utilizando la prueba 4 para calcular el valor de C:
qg 5,550 ⎛ Mpcd ⎞
C= = = 2.52 ⎜⎜ ⎟
(Pr
2
− Pwf2 ) n
(4.301× 10 )3 0.92 2 1.84 ⎟
⎝ (lb / pg ) ⎠
27
Capítulo I
Por lo tanto la ecuación que representa el comportamiento de este pozo es la siguiente:
(
qg = 2.52 Pr2 − Pwf2 ) 0.92
Donde:
qg =Gasto de gas (Mpcd)
Pr =Presión media del yacimiento (lb / pg2)
Pwf =Presión de fondo fluyendo (lb / pg2)
Ejemplo:
Calcular el qg para una Pwf de 193 (lb / pg2).
Utilizando la ecuación que representa el comportamiento de este pozo se tiene lo siguiente:
(
qg = 2.52 Pr2 − Pwf2 ) 0.92
(
= 2.52 2012 − 193 2 ) 0.92
= 4,170 Mpcd
I.5.3 IPR para pozos horizontales de gas
En un estudio presentado por N.Akhimiona y M.L Wiggins(10) se dieron las siguientes ecuaciones
para curvas de IPR, los autores utilizaron un simulador comercial para esta investigación. Las
ecuaciones de IPR fueron generadas en términos de presión y presión al cuadrado.
IPR en Términos de Presión.
2
qg ⎛P ⎞ ⎛P ⎞
= 1 − 0.811⎜⎜ wf ⎟⎟ − 0.189 ⎜⎜ wf ⎟⎟ ---------------1.32
qg max ⎝ Pr ⎠ ⎝ Pr ⎠
IPR en Términos de Presión al Cuadrado.
2
qg ⎛ Pwf2 ⎞ ⎛ P2 ⎞
= 1 − 1.867 ⎜⎜ 2 ⎟⎟ + 0.867 ⎜⎜ wf2 ⎟⎟ ---------------1.33
IPR Basadas en el Grado de Depresionamiento.
Se ha notado que el comportamiento presión-producción para un pozo horizontal, puede estar en

función del depresionamiento. Esto permite que IPR individuales sean desarrolladas para cada
etapa de explotación:
28
Capítulo I
2
qg ⎛ P2 ⎞ ⎛ P2 ⎞
= 1 − (1 + d) ⎜⎜ wf2 ⎟⎟ + d ⎜⎜ wf2 ⎟⎟ ---------------1.34
Donde d es un valor que depende de la etapa de explotación y se puede encontrar por el uso de
tablas que a continuación se presentan.
Valores del coeficiente d
Etapa de
Promdio Absoluto de Coeficiente
depresionamiento Coeficiente d
Diferencia de Ajuste R2
%
10 1.0369 8.72 0.9991
20 0.8404 2.44 0.9999
30 0.6531 14.16 0.918
40 0.6413 5 0.9894
50 0.6928 18.7 0.9063
60 0.7231 15.39 0.9172
70 0.7415 11.22 0.9282
Ejemplo:
Si tomamos los datos del ejemplo anterior para pozos verticales de gas y despejando de la
ecuación 1.32 el término de qg max se obtiene:
qg
qg max = 2
⎛P ⎞ ⎛P ⎞
1 − 0.811 ⎜⎜ wf ⎟⎟ − 0.189 ⎜⎜ wf ⎟⎟
⎝ Pr ⎠ ⎝ Pr ⎠
Tomando los datos de la prueba 4 para el cálculo qg max obtenemos lo siguiente:
5,550
qg max = 2
= 86,042 Mpcd
⎛ 190 ⎞ ⎛ 190 ⎞
1 − 0.811⎜ ⎟ − 0.189 ⎜ ⎟
⎝ 201 ⎠ ⎝ 201 ⎠
Realizando el mismo procedimiento pero para la ecuación 1.33 se obtiene:
qg
qg max = 2
⎛ Pwf2 ⎞ ⎛ P2 ⎞
1 − 1.867 ⎜⎜ 2 ⎟⎟ + 0.867 ⎜⎜ wf2 ⎟⎟
⎝ Pr ⎠ ⎝ Pr ⎠
Tomando los datos de la prueba 4 para el cálculo qg max obtenemos lo siguiente:
29
Capítulo I
5,550
qg max = 2
= 231,396 Mpcd
⎛ 190 2 ⎞ ⎛ 190 2 ⎞
1 − 1.867 ⎜⎜ 2 ⎟
⎟ + 0.867 ⎜⎜ ⎟
2 ⎟
⎝ 201 ⎠ ⎝ 201 ⎠
Por lo que las ecuaciones que representan el comportamiento de flujo en términos de presión y
presión al cuadrado respectivamente son:
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤
qg = 86,042 ⎢1 − 0.811 ⎜⎜ wf ⎟⎟ − 0.189 ⎜⎜ wf ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
⎡ ⎛ Pwf2 ⎞ ⎛ Pwf2 ⎞
2
⎤
qg = 231,396 ⎢1 − 1.867 ⎜⎜ 2 ⎟⎟ + 0.867 ⎜⎜ 2 ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
Calcular el gasto de gas para una Pwf de 193 (lb / pg2).
1) En términos de presión se obtiene lo siguiente:
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤ ⎡ ⎛ 193 ⎞ ⎛ 193 ⎞ ⎤
2
qg = 86,042 ⎢1 − 0.811⎜⎜ wf ⎟⎟ − 0.189 ⎜⎜ wf ⎟⎟ ⎥ = 86,042 ⎢1 − 0.811⎜ ⎟ − 0 . 189 ⎜ ⎟ ⎥

⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦ ⎣⎢ ⎝ 201 ⎠ ⎝ 201 ⎠ ⎦⎥
qg = 4,046 Mpcd
2) En términos de presión al cuadrado se obtiene lo siguiente:
⎡ ⎛ Pwf2 ⎞ ⎛ P2 ⎞
2
⎤ ⎡ ⎛ 193 2 ⎞ ⎛ 193 2 ⎞ ⎤
2
qg = 231,396 1 − 1.867 ⎜⎜ 2
⎢ ⎟⎟ + 0.867 ⎜⎜ wf2 ⎟⎟ ⎥ = 231,396 ⎢1 − 1.867 ⎜⎜ 2 ⎟
⎟ + 0.867 ⎜⎜ 2 ⎟
⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦ ⎢⎣ ⎝ 201 ⎠ ⎝ 201 ⎠ ⎥⎦
qg = 3,622 Mpcd
I.6 IPR de dos fases para pozos direccionales
La curva de Vogel fue desarrollada para pozos verticales y no es valida para pozos direccionales u
horizontales. El desarrollo de curvas de IPR para pozos direccionales y horizontales fue obtenido
por medio de un simulador, varios puntos importantes fueron observados:
Las curvas de IPR para pozos direccionales u horizontales tienen el mismo comportamiento que la
gráfica de Vogel.
30
Capítulo I
Los datos de IPR generados para pozos direccionales y horizontales pueden diferir hasta un 22%
y 27% respectivamente de los datos obtenidos usando la ecuación de Vogel para pozos verticales.
Ecuaciones de IPR para pozos direccionales y horizontales.
Angulo Ecuación.
q′ = 1 − 0.2 (P′) − 0.8 (P′)

2
0° (Vertical)
q′ = 0.9998 − 0.2210 (P′) − 0.7783 (P′)
2
15°
q′ = 0.9969 − 0.1254 (P′) − 0.8682 (P′)
2
30°
q′ = 0.9946 − 0.0221 (P′) − 0.9663 (P′)
2
45°
q′ = 0.9926 + 0.0549 (P′) − 1.0395 (P′)
2
60°
q′ = 0.9915 + 0.1002 (P′) − 1.0829 (P′)
2
75°
q′ = 0.9915 + 0.1120 (P′) − 1.0942 (P′)
2
85° (Horizontales)
q′ = 0.9885 + 0.2055 (P′) − 1.1818 (P′)
2
90°
Donde:
Pwf
P' =
Pr
qo
q' =
qo max
Ejemplo:
De una prueba de producción se obtienen los siguientes datos:
Pr = 2,400 (lb / pg2 )

qo = 100 bpd
Pwf = 1,800 (lb / pg2 )
Inclinación del pozo θ = 60°
Determinar el gasto de aceite para una Pwf de 800 (lb / pg2)
Solución:
1) Primero obtenemos el qomax para una Pwf de 1,800 (lb / pg2) y una Pr de 2,400 (lb / pg2).
Utilizaremos la ecuación para una inclinación de 60°.
31
Capítulo I
qo 100
qo max = 2
= 2
⎛P ⎞ ⎛P ⎞ ⎛ 1,800 ⎞ ⎛ 1,800 ⎞
0.9926 + 0.0549 ⎜⎜ wf ⎟⎟ − 1.0395 ⎜⎜ wf ⎟⎟ 0.9926 + 0.0549 ⎜ ⎟ − 1.0395 ⎜ ⎟
⎝ Pr ⎠ ⎝ Pr ⎠ ⎝ 2,400 ⎠ ⎝ 2,400 ⎠
100
qo max = = 223 bpd
0.449
2) Obteniendo el gasto de aceite para una Pwf de 800 (lb / pg2).
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤
qo = qo max ⎢0.9926 + 0.0549 ⎜⎜ wf ⎟⎟ − 1.0395 ⎜⎜ wf ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
⎡ ⎛ 800 ⎞ ⎛ 800 ⎞ ⎤
2
qo = 223 ⎢0.9926 + 0.0549 ⎜ ⎟ − 1.0395 ⎜ ⎟ ⎥ = 200 bpd

⎢⎣ ⎝ 2,400 ⎠ ⎝ 2,400 ⎠ ⎥⎦
Si hacemos el mismo desarrollo utilizando la ecuación de Vogel obtenemos un valor de

qo =211 bpd. Por lo que estaríamos subestimando la producción por 11 barriles diarios.
I.7 IPR de dos fases para pozos horizontales
Los pozos horizontales de aceite y gas han tenido un gran auge en la industria desde mediados de
los años 80. La mayor ventaja de estos pozos es que incrementan el área de contacto con el
yacimiento y por lo tanto mejoran la producción o inyectividad de un pozo.
Los pozos horizontales tal vez sean una viable alternativa económica para yacimientos
naturalmente fracturados, yacimientos con baja permeabilidad y en yacimientos con problemas de
conificación.
Retnanto y Economides investigaron el desarrollo de las curvas de IPR para pozos horizontales
produciendo bajo el mecanismo de expansión de gas disuelto, ellos utilizaron un simulador y
determinaron la siguiente ecuación:
n
qo ⎛P ⎞ ⎛P ⎞
= 1 − 0.25 ⎜⎜ wf ⎟⎟ − 0.75 ⎜⎜ wf ⎟⎟ ---------------1.35
Donde n está en función de la presión de burbuja.
⎛ ⎛ Pr ⎞ ⎛ Pr ⎞ ⎞⎟
2
⎜
n = − 0.27 + 1.46 ⎜⎜ ⎟⎟ − 0.96 ⎜⎜ ⎟⎟ (4 + 0.00166 Pb ) ---------------1.36
⎜ ⎝ Pb ⎠ ⎝ Pb ⎠ ⎟⎠
⎝
32
Capítulo I
Wang y Wiggins(14) propusieron 3 ecuaciones empíricas de curvas IPR para pozos horizontales,
la primera es una ecuación general para IPR, la segunda es una ecuación de IPR que está en
función del grado de depresionamiento o de recuperación Np y la última es una ecuación que nos
permite calcular IPR futuras. La ecuación general de IPR para pozos horizontales que encontraron
estos autores fue:
2
qo ⎛P ⎞ ⎛P ⎞
= 1 − 0.4533 ⎜⎜ wf ⎟⎟ − 0.5467 ⎜⎜ wf ⎟⎟ ---------------1.37
La segunda IPR fue desarrollada basada en diferentes etapas de recuperación. La ecuación a la

cual los autores llegaron es:
2
qo ⎛P ⎞ ⎛P ⎞
= 1.0 − d ⎜⎜ wf ⎟⎟ − (1 − d) ⎜⎜ wf ⎟⎟ ---------------1.38
Donde d está en función del factor de recuperación. La siguiente tabla muestra los valores de d
para las diferentes etapas.
Promedio
Coeficiente Coeficiente
Recuperación de Error
d de ajuste R2
Absoluto
0.01 0.675 0.0781 0.985
0.02 0.475 0.0845 0.985
0.04 0.25 0.0959 0.984
0.06 0.125 0.1031 0.981
0.08 0.2 0.1173 0.975
0.1 0.225 0.1577 0.957
0.12 0.25 0.1388 0.952
0.14 0.275 0.1931 0.907
Por ultimo desarrollaron una ecuación que permite el cálculo de IPR futuras:
2 3 4
qo max,F ⎛ Pr ⎞ ⎛ Pr ⎞ ⎛ Pr ⎞ ⎛ Pr ⎞
= 1.0 − 7.36 ⎜ ,F ⎟ + 20.32 ⎜ ,F ⎟ − 23.92 ⎜ ,F ⎟ + 10.96 ⎜ ,F ⎟ ---------------1.39
qo max,P ⎜P ⎟ ⎜P ⎟ ⎜P ⎟ ⎜P ⎟
⎝ r ,P ⎠ ⎝ r ,P ⎠ ⎝ r ,P ⎠ ⎝ r ,P ⎠
Para utilizar esta ecuación, primero se debe de encontrar el gasto máximo presente con cualquiera
de las dos ecuaciones (1.37 o 1.38) dadas anteriormente y con este dato, obtener el gasto máximo
futuro de la ecuación 1.39. Con el gasto máximo futuro regresar a las ecuaciones (1.37 o 1.38)
para poder graficar la curva de IPR futura.
33
Capítulo I
Ejemplo:
Pr = 2,400 (lb / pg2 )

qo = 300 bpd
Pwf = 1,800 (lb / pg2 )
a) Determinar el gasto de aceite para una Pwf de 800 (lb / pg2)

b) Determinar el gasto de aceite para una Pwf de 800 (lb / pg2) y una Pr de 2,000 (lb / pg2)
Solución:
a.1) Determinar el qo max para una Pwf de 1,800 (lb / pg2) y una Pr de 2,400 (lb / pg2)
qo 300
qo max = 2
= 2
⎛P ⎞ ⎛P ⎞ ⎛ 1,800 ⎞ ⎛ 1,800 ⎞
1 − 0.4533 ⎜⎜ wf ⎟⎟ − 0.5467 ⎜⎜ wf ⎟⎟ 1 − 0.4533 ⎜ ⎟ − 0.5467 ⎜ ⎟
⎝ Pr ⎠ ⎝ Pr ⎠ ⎝ 2,400 ⎠ ⎝ 2,400 ⎠
300
qo max = = 851 bpd
0.3525
a.2) Obteniendo el qo para una Pwf de 800 (lb / pg2)
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤
qo = qo max ⎢1 − 0.4533 ⎜⎜ wf ⎟⎟ − 0.5467 ⎜⎜ wf ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
⎡ ⎛ 800 ⎞ ⎛ 800 ⎞ ⎤
2
qo = 851 ⎢1 − 0.4533 ⎜ ⎟ − 0.5467 ⎜ ⎟ ⎥ = 671 bpd

⎢⎣ ⎝ 2,400 ⎠ ⎝ 2,400 ⎠ ⎥⎦
Si hacemos el mismo desarrollo utilizando la ecuación de Vogel obtenemos un valor de

qo =633 bpd. Por lo que estaríamos subestimando la producción por 38 barriles diarios.
b.1) Determinar el gasto máximo de aceite futuro cuando Pr cae hasta 2,000 (lb / pg2)
⎡ ⎛ Pr ,F ⎞ ⎛ Pr ,F ⎞
2
⎛ Pr ,F ⎞
3
⎛ Pr ,F ⎞ ⎤
4
qo max,F = qo max,P ⎢1.0 − 7.36 ⎜ ⎟ + 20.32 ⎜ ⎟ − 23.92 ⎜ ⎟ + 10.96 ⎜ ⎟ ⎥

⎢ ⎜P ⎟ ⎜P ⎟ ⎜P ⎟ ⎜P ⎟ ⎥
⎣ ⎝ r ,P ⎠ ⎝ r ,P ⎠ ⎝ r ,P ⎠ ⎝ r ,P ⎠ ⎦
34
Capítulo I
⎡ ⎛ 2,000 ⎞ ⎛ 2,000 ⎞
2
⎛ 2,000 ⎞
3
⎛ 2,000 ⎞ ⎤
4
qo max,F = 851 ⎢1.0 − 7.36 ⎜ ⎟ + 20.32 ⎜ ⎟ − 23.92 ⎜ ⎟ + 10.96 ⎜ ⎟ ⎥

⎢⎣ ⎝ 2,400 ⎠ ⎝ 2,400 ⎠ ⎝ 2,400 ⎠ ⎝ 2,400 ⎠ ⎥⎦
qo max,F = 358 bpd
b.2) Obteniendo el qo para una Pwf de 800 (lb / pg2) cuando la Pr cae hasta 2,000 (lb / pg2)
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤
qo = qo max ⎢1 − 0.4533 ⎜⎜ wf ⎟⎟ − 0.5467 ⎜⎜ wf ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
⎡ ⎛ 800 ⎞ ⎛ 800 ⎞ ⎤
2
qo = 358 ⎢1 − 0.4533 ⎜ ⎟ − 0.5467 ⎜ ⎟ ⎥ = 262 bpd

⎢⎣ ⎝ 2,000 ⎠ ⎝ 2,000 ⎠ ⎥⎦
I.8 IPR para flujo de tres fases (aceite, gas y agua)
M.L Wiggins(12), con ayuda de un simulador presentó ecuaciones para calcular curvas IPR de tres
fases. Para desarrollar la curva de IPR general para tres fases, realizó un procedimiento parecido
al de Vogel y obtuvo las siguientes ecuaciones.
2
qo ⎛P ⎞ ⎛P ⎞
= 1 − 0.519167 ⎜⎜ wf ⎟⎟ − 0.481092 ⎜⎜ wf ⎟⎟ ---------------1.40
2
qw ⎛P ⎞ ⎛P ⎞
= 1 − 0.722235 ⎜⎜ wf ⎟⎟ − 0.284777 ⎜⎜ wf ⎟⎟ ---------------1.41
qw max ⎝ Pr ⎠ ⎝ Pr ⎠
El método propuesto asume que podemos tratar cada fase por separado. Como se observa en las
ecuaciones anteriores necesitamos una simple prueba de producción para encontrar los gastos
máximos de aceite y agua, con lo cual se podrá predecir el comportamiento del gasto (q) con
respecto a una presión de fondo fluyendo (Pwf). M.L Wiggins también desarrolló una ecuación para
predecir el gasto máximo de aceite y agua futuros, obteniendo las siguientes ecuaciones:
2
qo max F ⎛ Pr ⎞ ⎛ Pr ⎞
= 0.15376309 ⎜ ,F ⎟ + 0.83516299 ⎜ ,F ⎟ ---------------1.42
qo max P ⎜P ⎟ ⎜P ⎟
⎝ r ,P ⎠ ⎝ r ,P ⎠
2
qw max F ⎛ Pr ⎞ ⎛ Pr ⎞
= 0.59245433 ⎜ ,F ⎟ + 0.36479178 ⎜ ,F ⎟ ---------------1.43
qw max P ⎜P ⎟ ⎜P ⎟
⎝ r ,P ⎠ ⎝ r ,P ⎠
35
Capítulo I
Para utilizar las ecuaciones anteriores, primero se debe determinar los gastos máximos de aceite y
agua a las condiciones presentes (ecuaciones 1.40 y 1.41) y después calcular el gasto máximo
futuro de cada fase para la presión del yacimiento futura (ecuaciones 1.42 y 1.43).
Estas ecuaciones arrojan un mayor error cuando el grado de depresionamiento es mayor, Para
que estas ecuaciones no tengan un error muy significativo se recomienda que los valores de
gastos futuros se realicen cada 6 meses.
Existen otros métodos como el presentado por PETROBRAS, el cual es un procedimiento analítico
para las curvas de IPR de tres fases.
Ejemplo:
qt = 500 bpd
Pr = 2,550 (lb / pg2 )
Pwf = 2,300 (lb / pg2 )
Fw = 0.25
Determinar:
a) El gasto de aceite y agua para una Pwf de 2,000 (lb / pg2)

b) El gasto de aceite y agua para una Pwf de 1,500 (lb / pg2) y una Pr de 2,100 (lb / pg2)
Solución:
a.1) Determinar el gasto de aceite y agua de la prueba de producción
qo = 500 × 0.75 = 375 bpd

qw = 500 × 0.25 = 125 bpd
a.2) Determinar los gastos máximos de aceite y agua
qo 375
qo max = 2
= 2
⎛P ⎞ ⎛P ⎞ ⎛ 2,300 ⎞ ⎛ 2,300 ⎞
1 − 0.519167 ⎜⎜ wf ⎟⎟ − 0.481092 ⎜⎜ wf ⎟⎟ 1 − 0.519167 ⎜ ⎟ − 0.481092 ⎜ ⎟
⎝ Pr ⎠ ⎝ Pr ⎠ ⎝ 2,550 ⎠ ⎝ 2,550 ⎠
375
qo max = = 2,672 bpd
0.14
36
Capítulo I
qw 125
qw max = 2
= 2
⎛P ⎞ ⎛P ⎞ ⎛ 2,300 ⎞ ⎛ 2,300 ⎞
1 − 0.722235 ⎜⎜ wf ⎟⎟ − 0.284777 ⎜⎜ wf ⎟⎟ 1 − 0.722235 ⎜ ⎟ − 0.284777 ⎜ ⎟
⎝ Pr ⎠ ⎝ Pr ⎠ ⎝ 2,550 ⎠ ⎝ 2,550 ⎠
125
qw max = = 1,069 bpd
0.116
a.3) Obteniendo el gasto de aceite y agua para una Pwf de 2,000 (lb / pg2)
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤
qo = qo max ⎢1 − 0.519167 ⎜⎜ wf ⎟⎟ − 0.481092 ⎜⎜ wf ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
⎡ ⎛ 2,000 ⎞ ⎛ 2,000 ⎞ ⎤
2
qo = 2,672 ⎢1 − 0.519167 ⎜ ⎟ − 0.481092 ⎜ ⎟ ⎥ = 793 bpd

⎢⎣ ⎝ 2,550 ⎠ ⎝ 2,550 ⎠ ⎥⎦
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤
qw = qwmax ⎢1 − 0.722235 ⎜⎜ wf ⎟⎟ − 0.284777 ⎜⎜ wf ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
⎡ ⎛ 2,000 ⎞ ⎛ 2,000 ⎞ ⎤
2
qw = 1,069 ⎢1 − 0.722235 ⎜ ⎟ − 0.284777 ⎜ ⎟ ⎥ = 276 bpd

⎣⎢ ⎝ 2,550 ⎠ ⎝ 2,550 ⎠ ⎦⎥
b.1) Determinar los gastos máximos de aceite y agua futuros para una Pr de 2,100 (lb / pg2)
2
qo max F ⎛ Pr ,F ⎞ ⎛ Pr ,F ⎞ ⎛ 2,100 ⎞ ⎛ 2,100 ⎞
2
= 0.15376309 ⎜ ⎟ + 0.83516299 ⎜ ⎟ = 0.15376309 ⎜ ⎟ + 0.83516299 ⎜ ⎟

qo max P ⎜P ⎟ ⎜P ⎟ ⎝ 2,550 ⎠ ⎝ 2,550 ⎠
⎝ r ,P ⎠ ⎝ r ,P ⎠
qo max F
= 0.693 ∴ qo max F = 0.693 × 2,672 = 1,852 bpd
qo max P
2
qw max F ⎛ Pr ⎞ ⎛ Pr ⎞ ⎛ 2,100 ⎞ ⎛ 2,100 ⎞
= 0.59245433 ⎜ ,F ⎟ + 0.36479178 ⎜ ,F ⎟ = 0.59245433 ⎜ ⎟ + 0.36479178 ⎜ ⎟
qw max P ⎜ ⎟ ⎜ ⎟
⎝ Pr ,P ⎠ ⎝ Pr ,P ⎠ ⎝ 2,550 ⎠ ⎝ 2,550 ⎠
37
Capítulo I
qw max F
= 0.735 ∴ qw max F = 0.735 × 1,069 = 786 bpd
qw max P
b.2) Obteniendo el qo y qw para una Pwf de 1,500 (lb / pg2) y una Pr de 2,100 (lb / pg2)
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤
qo = qo max ⎢1 − 0.519167 ⎜⎜ wf ⎟⎟ − 0.481092 ⎜⎜ wf ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
⎡ ⎛ 1,500 ⎞ ⎛ 1,500 ⎞ ⎤
2
qo = 1,852 ⎢1 − 0.519167 ⎜ ⎟ − 0.481092 ⎜ ⎟ ⎥ = 711 bpd

⎢⎣ ⎝ 2,100 ⎠ ⎝ 2,100 ⎠ ⎥⎦
⎡ ⎛P ⎞ ⎛P ⎞
2
⎤
qw = qw max ⎢1 − 0.722235 ⎜⎜ wf ⎟⎟ − 0.284777 ⎜⎜ wf ⎟⎟ ⎥
⎢⎣ ⎝ Pr ⎠ ⎝ Pr ⎠ ⎥⎦
⎡ ⎛ 1,500 ⎞ ⎛ 1,500 ⎞ ⎤
2
qw = 786 ⎢1 − 0.722235 ⎜ ⎟ − 0.284777 ⎜ ⎟ ⎥ = 266 bpd

⎣⎢ ⎝ 2,100 ⎠ ⎝ 2,100 ⎠ ⎦⎥
I.9 Curvas de IPR futuras
La elaboración de curvas de IPR futuras es muy importante en la determinación de gastos futuros

de producción para pozos fluyentes o con sistema artificial de producción. Estas curvas también
representan una manera de determinar cuando un pozo debe ser equipado con un sistema
artificial de producción. La siguiente gráfica (figura 1.11) representa un conjunto típico de este tipo
de curvas.
Figura 1.11 Curvas de IPR para presiones de yacimiento presentes y futuras
38
Capítulo I
I.9.1 Método de Standing
Muskat notó que el índice de productividad de un pozo en el tiempo 1 puede ser relacionado con
el índice de productividad del mismo pozo al tiempo 2 mediante la siguiente relación:
⎛ k ro ⎞
⎜ ⎟
J1 ⎜⎝ μo Bo ⎟⎠1
= ---------------1.44
J2 ⎛ k ro ⎞
⎜⎜ ⎟⎟
⎝ μ o Bo ⎠ 2
Donde los subíndices 1 y 2 indican los tiempos respectivamente.
Standing utilizó la ecuación de Vogel en combinación con la relación de Muskat para obtener un
método que determine la curva de IPR futura. Él obtuvo la siguiente ecuación para determinar la
curva de IPR futura de un pozo:
J * Pr,F ⎡ ⎛ Pwf ⎞ ⎤
2
⎛ Pwf ⎞
qo = ⎢1 − 0.2 ⎜ ⎟ ⎜ ⎟ ⎥
1 .8 ⎢ ⎜ P ⎟ − 0 .8 ⎜ P ⎟ ⎥ ---------------1.45
⎣ ⎝ r ,F ⎠ ⎝ r ,F ⎠ ⎦
En la ecuación J* es el valor del índice de productividad a la Pr para la cual se quiere conocer la

nueva curva de IPR, es decir el valor futuro de la presión media del yacimiento.
Donde:
1.8 qo max,F
J* =
Pr ,F
Si J*P es el valor al día presente y J*F el valor futuro, la relación de índice de productividad es:
⎛ K ro ⎞
⎜⎜ ⎟
⎝ μo Bo ⎟⎠F
J * F = J *P
⎛ K ro ⎞
⎜⎜ ⎟⎟
⎝ μ o B o ⎠P
Procedimiento:
1) Calcular el valor del índice de productividad presente J*P.

2) Calcular el valor del índice de productividad futuro J*F de la relación de Muskat.
3) Sustituir el valor de J*F en la ecuación 1.45 y construir la gráfica de IPR futura para una Pr
dada.
39
Capítulo I
I.9.2 Método de Fetkovich
Fetkovich propuso la siguiente ecuación para el comportamiento de la curva de IPR futura.
⎛P ⎞ 2
(
qo = J'o,P ⎜⎜ r,F ⎟⎟ Pr ,F − Pwf2 ) n
---------------1.46
⎝ Pr ,P ⎠
Donde recordemos que el valor de J‘o,P y n son los valores de ordenada al origen y el inverso de la
pendiente respectivamente para la gráfica Log-Log de producción al tiempo 1, es decir a la Pr
(presente). qo y Pr,F son el gasto y la presión media del yacimiento en el tiempo 2 (futuro).
I.9.3 Método de Eckmier
Eckmier notó que si tomamos las ecuaciones de Fetkovich para presión media del yacimiento al
tiempo 1 (presente) y para la presión media del yacimiento al tiempo 2 (futuro), llegaremos a una
ecuación para determinar el qomax,F, después éste lo podemos usar directamente en la ecuación de
Vogel y así preparar la curva de IPR futura.
3
qo max,P ⎛ Pr ,P ⎞
=⎜ ⎟ ---------------1.47
qo max,F ⎜⎝ Pr ,F ⎟⎠
I.9.4 Método de Couto
El procedimiento de Couto puede ser usado tanto para curvas de IPR presentes como futuras,
esto sólo si se conocen los valores de permeabilidad relativa al aceite (kro), la viscosidad del aceite
(μo) y el factor de volumen (Bo) a ciertas condiciones de presión y temperatura para la cual se
quiere saber la curva de IPR futura. La ecuación dada por Couto es:
kh ⎛ k ⎞
qo = 3.49(A ) Pr ⎜⎜ ro ⎟⎟ (EF )(1 − P' ) [1.8 − 0.8 (EF )(1 − P')] ---------------1.48
In (0.472 re / rw ) ⎝ μo Bo ⎠
Donde:
A =Factor de conversión =0.001127, para (pie, lb / pg2, cp, md y bpd)

Pwf
P' =
Pr
EF =Eficiencia de flujo
La ecuación que se utilice para determinar las curvas de IPR futuras depende de los datos con los
que se cuente, si queremos utilizar el método de Standing o el de Couto necesitamos conocer las
propiedades del sistema roca-fluidos a la presión media del yacimiento para la cual se quiere
evaluar la curva IPR futura.
40
Capítulo I
Si tenemos 4 pruebas de producción podemos ocupar el método de Fetkovich. También podemos

utilizar el método de Eckmier si sólo contamos con el dato de gasto máximo al tiempo presente.
Por lo tanto la ecuación a utilizar para conocer la curva de IPR futura dependerá de la cantidad de
información con la que contemos.
I.9.5 IPR futura para pozos de gas
Para predecir el comportamiento de las curvas de IPR para pozos de gas en un tiempo futuro se
puede utilizar la siguiente aproximación.
(
qg = CF Pr2,F − Pwf2 ) n
---------------1.49
Comparando esta ecuación con la ecuación 1.30 puede ser visto que los efectos por turbulencia
Dqg son incluidos en el valor del exponente n, y que el coeficiente C contiene varios parámetros
sujetos a cambio, como se muestra a continuación.
703 × 10 −6 kh
Debido a que C =
⎡ ⎛ 0.472 re ⎞ ⎤
T ⋅ μg ⋅ Z ⎢In ⎜⎜ ⎟⎟ + S⎥
⎣⎢ ⎝ rw ⎠ ⎦⎥
En la ecuación anterior las propiedades que cambian fundamentalmente con el cambio en la

presión del yacimiento son la viscosidad y el factor de compresibilidad, por lo que se tiene la
siguiente relación:
C1 (μg Z )2
= , por lo tanto para nuestro caso tenemos que:
C2 (μg Z )1
(μ Z )
g p
CF = Cp ---------------1.50
(μ Z )
g F
Donde:
qg =Gasto futuro de gas
CP =Coeficiente al tiempo presente consistente con las unidades de gasto y presión
CF =Coeficiente al tiempo futuro consistente con las unidades de gasto y presión
Pr =Presión media del yacimiento a un tiempo dado (lb / pg2)
Pwf =Presión de fondo fluyendo (lb / pg2)
El exponente n se asume como constante y los subíndices P Y F indican presente y futuro

respectivamente.
41
Capítulo II
CAPITULO II
FLUJO EN TUBERÍAS
Para poder explotar un campo petrolero se necesita perforar varios pozos, un pozo es un agujero
que se hace a través de la roca hasta el yacimiento, en el cual se instalan tuberías y otros
elementos con el fin de transportar controladamente los fluidos desde el yacimiento hasta la
superficie.
Los fluidos al entrar al pozo son transportados hasta la superficie por la tubería de producción
(T.P), de la superficie hasta el separador son transportados por una tubería llamada línea de
descarga (L.D). Es por esta razón, que un segundo paso para la optimización de la producción es
analizar el comportamiento de flujo en tuberías, este flujo puede ser de líquido, gas o multifásico.
Las principales caídas de presión en el sistema integral de producción se dan en la tubería de

producción (T.P), y en la línea de descarga (L.D) cuando esta es demasiado extensa, por lo tanto
es de vital importancia para los ingenieros de producción contar con las herramientas adecuadas
para predecir las caídas de presión y de esta manera optimizar el diseño de los pozos, con el
objetivo de maximizar la vida productiva de cada uno de ellos. Para cumplir con este objetivo
muchos investigadores se han dado a la tarea de desarrollar correlaciones de flujo multifásico con
base en la experimentación, en este capítulo sólo se verán dos correlaciones para flujo multifásico
horizontal y dos para flujo vertical.
Los tres principales componentes que gobiernan las caídas de presión en tuberías para cualquier
tipo de fluido y cualquier dirección de flujo son:
• Pérdidas por elevación

• Pérdidas por fricción
• Pérdidas por aceleración
Para flujo vertical o direccional las caídas de presión por elevación son las más significativas y
para flujo horizontal las caídas de presión por fricción son las que predominan.
El desarrollo de todas las correlaciones de flujo es a partir de la ecuación general de balance de

energía, por lo tanto empezaremos este capítulo explicando brevemente su desarrollo.
42
Capítulo II
II.1 Fundamentos de flujo a través de tuberías
II.1.1 Ecuación de balance de energía
La ecuación general que gobierna el flujo de fluidos a través de una tubería, se obtiene a partir de
un balance microscópico de la energía asociada a la unidad de masa de un fluido que pasa a
través de un elemento aislado de sistema, como se muestra en la figura 2.1.
Figura 2.1 Diagrama de flujo para un conducto aislado
De acuerdo con la ley de la conservación de la energía:
E1 + ∆Wf + ∆Ws =E2 -----------------2.1
∆Wf =Pérdidas de energía por fricción, estas pérdidas corresponden a la fricción interna del fluido
(viscosidad) y a la fricción del fluido con las paredes rugosas de la tubería.
∆Ws =Pérdidas de energía por trabajo externo.
⎛ lb − pie ⎞
E1 =Energía por unidad de masa en la posición uno ⎜⎜ f ⎟⎟
⎝ lbm ⎠
⎛ lb − pie ⎞
E2 =Energía por unidad de masa en la posición dos ⎜⎜ f ⎟⎟
⎝ lbm ⎠
Energía de expansión
La energía de expansión está dada por la siguiente expresión:
43
Capítulo II
⎛ lb f − pie ⎞ ⎛ lb f ⎞ ⎛ pie3 ⎞
Ee ⎜⎜ ⎟⎟ = P ⎜⎜ ⎟ V ⎜⎜
2 ⎟
⎟⎟ = PV -----------------2.2
⎝ lbm ⎠ ⎝ pie ⎠ ⎝ lbm ⎠
Donde:
V =Volumen especifico (1 / ρ).
Energía potencial
La energía potencial está dada por la siguiente expresión:
⎛ lb − pie ⎞ ⎛ pie ⎞ 1 ⎛ lb f − s2 ⎞ g
Ep ⎜⎜ f ⎟⎟ = g ⎜ 2 ⎟ ⎜⎜ ⎟⎟ h(pie ) = h -----------------2.3
⎝ lbm ⎠ ⎝ s ⎠ gc ⎝ lbm − pie ⎠ gc
Energía cinética
La energía cinética está dada por la siguiente expresión:
⎛ lb − pie ⎞ v 2 ⎛ pie2 ⎞ 1 ⎛ lb f − s2 ⎞ v 2
Ec ⎜⎜ f ⎟⎟ = ⎜⎜ 2 ⎟⎟ ⎜⎜ ⎟⎟ = -----------------2.4
⎝ lbm ⎠ 2 ⎝ s ⎠ gc ⎝ lbm − pie ⎠ 2gc
Al sustituir las energías correspondientes a las posiciones 1 y 2 en la ecuación 2.1 se obtiene:
g v2 g v2
P1V1 + h1 + 1 + ΔWf + ΔWS = P2 V2 + h2 + 2 -----------------2.5
gc 2gc gc 2gc
Donde:
V =Volumen especifico (1 / ρ)
Multiplicando la ecuación 2.5 por (ρ / ∆L) y considerando despreciables las pérdidas de energía
por trabajo externo se tiene:
ΔP gΔh Δv 2 ΔWf
=ρ +ρ +ρ
ΔL gc ΔL 2gc ΔL ΔL
A esta ecuación se le acostumbra escribir de la siguiente forma:
⎛ ΔP ⎞ ⎛ ΔP ⎞ ⎛ ΔP ⎞ ⎛ ΔP ⎞
⎜ ⎟ =⎜ ⎟ +⎜ ⎟ +⎜ ⎟ -----------------2.6
⎝ ΔL ⎠T ⎝ ΔL ⎠E ⎝ ΔL ⎠ AC ⎝ ΔL ⎠F
44
Capítulo II
Donde:
⎛ ΔP ⎞
⎜ ⎟ =Gradiente de presión total
⎝ ΔL ⎠ T
⎛ ΔP ⎞
⎜ ⎟ =Gradiente de presión debido a la elevación
⎝ ΔL ⎠E
⎛ ΔP ⎞
⎜ ⎟ =Gradiente de presión debido a la aceleración
⎝ ΔL ⎠ AC
⎛ ΔP ⎞
⎜ ⎟ =Gradiente de presión debido a la fricción
⎝ ΔL ⎠F
II.1.2 Pérdidas de presión por fricción
Las pérdidas de presión por fricción en conductos circulares de diámetro constante han sido
determinadas experimentalmente por varios investigadores. Los resultados de esta
experimentación utilizando tuberías de diversos materiales, constituyen la base de las fórmulas
que actualmente se usan.
Ecuación de Darcy
Darcy, Weisbach y otros, dedujeron experimentalmente la siguiente ecuación expresada en

unidades consistentes para obtener el gradiente de presión debido a la fricción:
⎛ ΔP ⎞ fρv 2
⎜ ⎟ = -----------------2.7
⎝ ΔL ⎠F 2gc d
Factor de fricción
El valor del factor de fricción (f) es función de la rugosidad de la tubería (ε) y del número de
Reynolds (Nre), por lo tanto:
f = F(ε, Nre )
El Número de Reynolds (adimensional) se define como:
dvρ
Nre = -----------------2.8
μ
Donde:
45
Capítulo II
d =Diámetro de la tubería (pie)

v =Velocidad del líquido (pie / seg)
ρ =Densidad del líquido (lbm / pie3)
μ =Viscosidad del líquido (lbm / pie-seg)
La rugosidad (ε) de una tubería es una característica de su superficie, la cual está constituida por
pliegues o crestas unidas. Los valores más comúnmente empleados en la industria son:
Tipo de Tubería ε (pg)

Tubería estriada 0.00006
Tubería de producción o
0.0006
perforación
Tubería de escurrimiento 0.0007
Tubería galvanizada 0.006
Para calcular el valor de (f) es necesario determinar el régimen de flujo. En conductos, los fluidos
se mueven de acuerdo a cualquiera de los siguientes regímenes de flujo: Laminar, transitorio o
turbulento.
El flujo laminar ocurre cuando las partículas del fluido se mueven en líneas rectas paralelas al eje
del conducto. A velocidades mayores las partículas se mueven de una manera caótica, en este
caso el flujo es turbulento. Osborne Reynolds estableció experimentalmente un parámetro para
determinar el régimen de flujo en tuberías, a este parámetro le llamo número de Reynolds. El
flujo laminar se presenta cuando Nre < 2300, y el flujo turbulento cuando Nre > 3100.
Para flujo laminar de una sola fase (Nre < 2300), (f) depende exclusivamente del valor del (Nre)
y está dado por:
64
f= -----------------2.9
Nre
Para flujo turbulento (Nre > 3100), (f) está dado por la ecuación de Colebrook y White:
−2
⎡ ⎛ ε 2.514 ⎞⎤
f = ⎢− 2 log ⎜⎜ + ⎟⎥ ---------------2.10
⎣⎢ ⎝ 3.715 d f Nre ⎟⎠⎦⎥
Se observa que para calcular (f) se requiere un proceso iterativo, debido a que (f) se encuentra en
ambos lados de la igualdad. Basándose en la ecuación 2.10, Moody preparó un diagrama para
determinar el factor de fricción en tuberías de rugosidad comercial (figura 2.2).
El diagrama de Moody es el siguiente:
46
Capítulo II
Figura 2.2 Factor de fricción para cualquier tipo de tubería comercial
a) Se observa que para valores de (Nre < 2300) “flujo laminar”, (f) depende exclusivamente del
número de Reynolds.
b) A partir de (Nre > 3100) se inicia la zona de flujo turbulento. Dentro de esta, (f) depende tanto
del (Nre) como del valor de rugosidad relativa (ε / d), el cual es un valor adimensional.
c) La zona francamente turbulenta se inicia a diferentes valores de (Nre), dependiendo del valor de
la rugosidad relativa (ε / d). En esta zona (f) es independiente del (Nre) y varía únicamente con el
valor de (ε / d), el valor de (f) en esta zona es calculado por la siguiente ecuación:
f = (− 2 log ( ε / 3.715 d))

−2
---------------2.11
d) Cuando el flujo es crítico (2300< Nre <3100), (f) se puede aproximar con la siguiente expresión:
⎡ ⎤
⎢ ⎥
Nre − 2300 ⎢ 1.3521 ⎥
f= ×⎢ 2⎥
+ 0.032 ---------------2.12
2300 ⎢ ⎛⎜ ⎛ ε 2.514 ⎞⎞ ⎥
⎢ ⎜ 2.3026 log ⎜⎜ 3.715 d + 3100 f ⎟⎟ ⎟ ⎥
⎟
⎣⎝ ⎝ ⎠⎠ ⎦
47
Capítulo II
De esta manera se pueden calcular los valores del factor de fricción (f) dependiendo el tipo de flujo
que tengamos. Posteriormente el valor de (f) será utilizado para calcular las caídas de presión por
fricción.
II.2 Flujo de una sola fase
II.2.1 Flujo de líquido
De la ecuación 2.6, el gradiente de presión total despreciando el efecto por aceleración se puede
escribir de la siguiente manera:
ΔPT = ΔPE + ΔPF
La caída de presión por elevación se calcula de la siguiente manera:
ΔPE = 0.433 λL Δh ---------------2.13
Donde:
ΔPE =Caída de presión por elevación (lb / pg2)
λL =Densidad relativa del líquido
Δh =Diferencia de alturas (pie)
La pérdida de presión por fricción en unidades prácticas, se obtiene a partir de la ecuación de

Darcy, quedando de la siguiente forma:
f ⋅ λL ⋅ q2 ⋅ L
ΔPF = 0.06056 ---------------2.14
d5
Donde:
q =Gasto de líquido (bpd)
L =Longitud de la tubería (millas)

d =Diámetro de la tubería (pg)
ΔPF =Caída de presión por fricción (lb / pg2)
Para obtener la caída de presión total a lo largo de una tubería, sumamos las caídas de presión
debido a la elevación y a la fricción. La ecuación final queda de la siguiente manera:
48
Capítulo II
f ⋅ λL ⋅ q2 ⋅ L
ΔPT = 0.433 λL Δh + 0.06056 ---------------2.15
d5
Como ya vimos, necesitamos obtener el valor del número de Reynolds para encontrar el valor del
factor de fricción (f), por lo que también es importante encontrar el factor de conversión que nos
permita calcular el (Nre) en unidades prácticas, el cual queda de la siguiente forma:
q λL
Nre = 92.2 ---------------2.16
dμ
Donde:
q =Gasto de líquido (bpd)

μ =Viscosidad del líquido (cp)
λL =Densidad relativa del líquido
II.2.2 Flujo de gas
Como ya vimos anteriormente, la caída de presión total es la suma de la caída de presión por
elevación y por fricción. Para flujo de gas, la forma de calcular la caída de presión debido a la
fricción es partiendo de la ecuación de Darcy, haciendo las sustituciones convenientes la ecuación
final queda de la siguiente forma:
⎛ f ⎞ ⎛ P0 ⎞ q2 ⋅ λ g ⋅ Z ⋅ (T + 460 ) ⋅ L
ΔPF = ⎜ ⎟ ⎜⎜ ⎟⎟
⎝ 519.98 ⎠ ⎝ T0 + 460 ⎠ d5 ⋅ P
Observamos en la ecuación anterior que la caída de presión por fricción está en función del valor
de la presión media P y es definida como P =(P1 + P2) / 2, y además ∆PF puede ser definida como
∆PF =(P1 − P2). Por lo tanto, si multiplicamos ambos miembros de la ecuación por P y buscamos
el factor de conversión para unidades prácticas, obtenemos lo siguiente:
⎛ f ⎞ ⎛ P0 ⎞ q2 ⋅ λ g ⋅ Z ⋅ (T + 460 ) ⋅ L
P12 − P22 = ⎜ 2 ⎟⎜
⎜ ⎟⎟ ---------------2.17
⎝ 461.346 ⎠ ⎝ T0 + 460 ⎠ d5
Donde:
q =Gasto de gas (pcd)
L =Longitud de la tubería (millas)

49
Capítulo II
f =Factor de Fricción
P =Presión media (P1 + P2) / 2 (lb / pg2)
T =Temperatura media (T1 + T2) / 2 (°F)
λ g =Densidad relativa del gas
Z =Factor de compresibilidad del gas medido a P y T

P1 y P2 =Presiones al inicio y al final P1 > P2 (lb / pg2)
P0 y T0 =Presión y Temperatura @ C.S, normalmente 14.7 (lb / pg2) y 60°F
La caída de presión debido al cambio de elevación se expresa de la siguiente forma:
ΔPE = ρg h
Por lo que las pérdidas por elevación en unidades prácticas queda de la siguiente forma:
⎛ 2.7044 ⎞ P ⋅ λ g ⋅ h
ΔPE = ⎜ ⎟
⎝ 144 ⎠ Z ⋅ (T + 460 )
Nuevamente podemos sustituir el valor de ∆PE por (P1 – P2) y multiplicar ambos miembros de la
ecuación por (P1 + P2) / 2, obteniendo lo siguiente:
P 2 ⋅ λg ⋅ h
P12 − P´22 = 0.03756 ---------------2.18
Z ⋅ (T + 460 )
Donde:
h = Elevación con respecto a un nivel base (pie)
Es importante señalar que el valor de h se toma respecto a un nivel base que normalmente es la
superficie, el valor de h será positivo hacia arriba y negativo en caso contrario. Otra forma de ver el
signo de h es el siguiente: Si el flujo va en contra de la fuerza de gravedad, el valor de la diferencia
de alturas (h) será positivo, por otro lado si el flujo va en el mismo sentido de la fuerza de
gravedad, el valor de la diferencia de alturas (h) será negativo.
La ecuación 2.18 es la que nos permite evaluar exclusivamente las pérdidas de presión por
elevación. Como se observa, a diferencia de la ecuación de pérdidas por fricción, es necesario
suponer un valor de P2 para obtener una P y así calcular las pérdidas de presión por elevación.
Por lo tanto se necesita un proceso iterativo hasta alcanzar que la P2 supuesta sea
aproximadamente igual a la P2 calculada.
50
Capítulo II
Sumando las ecuaciones 2.17 y 2.18 podemos calcular P1 o P2 según sea el caso, tomando en
cuenta las siguientes suposiciones:
-Flujo en régimen permanente.

-Flujo de una sola fase (gas).
-Pérdidas por aceleración despreciables.
⎛ f ⎞ ⎛ P0 ⎞ q2 ⋅ λ g ⋅ Z ⋅ (T + 460 ) L P 2 ⋅ λg ⋅ h
P −P = ⎜
2 2
2 ⎟⎜
⎜ ⎟⎟ + 0.03756 ---------------2.19
Z ⋅ (T + 460 )
1 2
⎝ 461.346 ⎠ ⎝ T0 + 460 ⎠ d5
Como se puede observar en la ecuación anterior, se debe suponer un valor de P2 para poder
obtener el valor de P y con éste calcular el valor de Z .
La ecuación 2.19 puede ser despejada para obtener un gasto determinado dadas unas
condiciones de presión de entrada y salida, así como el diámetro de la tubería.
0. 5
C ⎡ 2 ⎛ 0.03756 P 2 ⋅ λ g ⋅ h ⎞⎤
qg = 0.5 ⎢P1 − P2 − ⎜⎜ ⎟⎥
2
d2.5 ---------------2.20
L ⎢⎣ ⎝ Z ⋅ (T + 460 ) ⎟⎥
⎠⎦
Donde:
0. 5
⎡⎛ 461.3462 ⎞ ⎛ T + 460 ⎞ ⎛ 1 ⎞⎤
C = ⎢⎜⎜ ⎟⎟ ⎜⎜ 0 ⎟⎟ ⎜ ⎟⎥ ---------------2.21
⎢⎣⎝ f ⎠⎝ P0
⎜
⎠⎝ gλ ⋅ Z ⋅ (T + 460 ) ⎟⎥
⎠⎦
O bien el diámetro dadas unas condiciones de presión y gasto:
0. 2
⎡⎛ L ⎞ ⎛ 1 ⎞ 2⎤
d = ⎢⎜ ⎟ ⎜ 2 ⎟ qg ⎥ ---------------2.22
2
(
⎣⎢⎝ C ⎠ ⎝ P1 − P2 − 0.03756 P ⋅ λ g ⋅ h / Z ⋅ (T + 460 )
⎜ 2
) ⎟ ⎥
⎠ ⎦
Al igual que el flujo de líquido por tuberías, es conveniente obtener una ecuación del número de
Reynolds para flujo de gas en las que sus factores estén en unidades prácticas. Haciendo las
respectivas sustituciones el (Nre) en unidades de campo queda de la siguiente forma:
qg λ g
Nre = 0.0201056 ---------------2.23
d μg
Donde:
qg =Gasto de gas (pcd)
51
Capítulo II
λ g =Densidad relativa del gas

μg =Viscosidad del gas (cp)
Como la distribución de presiones en un gasoducto no es lineal, se recomienda el uso de la

siguiente ecuación para obtener la presión media.
2 ⎡ P13 − P23 ⎤
P= ⎢ 2 ⎥ ---------------2.24
3 ⎣ P1 − P22 ⎦
Al igual que en el flujo de líquidos por tuberías, la eficiencia es un factor de ajuste para compensar
los efectos de corrosión, erosión, rugosidad e incrustaciones que no se consideran en la
deducción de las ecuaciones de flujo, de tal forma que los resultados que se obtengan deben ser
corregidos para obtener el gasto real. Los valores más comunes de EF (Eficiencia de Flujo) para
flujo de gas, según Ikoku son:
Contenido
de liquido
Linea EF
gal /
MMpie3
GAS SECO 0.1 0.92
GAS HUMEDO 7.2 0.77
GAS Y CONDENSADO 800 0.6
II.3 Flujo multifásico en tuberías
El flujo multifásico en tuberías es definido como el movimiento concurrente de gas libre y líquido
en las tuberías. El gas y el líquido pueden existir como una mezcla homogénea, o el líquido puede
estar en baches con el gas empujando detrás de el, u otras combinaciones de patrones de flujo
pueden estar presentes. El gas puede estar fluyendo con dos líquidos (aceite y agua) y existe la
posibilidad de que los dos líquidos puedan estar emulsificados. Los problemas de flujo multifásico
pueden ser divididos en 4 categorías:
1) Flujo multifásico vertical

2) Flujo multifásico horizontal
3) Flujo multifásico inclinado
4) Flujo multifásico direccional.
Los fluidos después de entrar al pozo viajan a través de la tubería de producción o por el espacio
anular hacia la superficie, esto implica problemas de flujo multifásico vertical o direccional. Al llegar
a superficie los fluidos pueden o no pasar a través de un estrangulador, y de ahí fluir a través de
52
Capítulo II
tuberías llamadas líneas de descarga, este flujo multifásico a través de las líneas de descarga
implica problemas de flujo multifásico horizontal o inclinado.
La aplicación de correlaciones de flujo multifásico para predecir las caídas de presión en tuberías
es extremadamente importante para la industria petrolera. Estos son algunos de sus usos:
1) Flujo natural: Cuando los fluidos son producidos desde el yacimiento y a través de todo aparejo
de producción la energía es disipada. Para prolongar la vida fluyente de los pozos tanto como sea
posible hay una necesidad de minimizar las pérdidas de energía (presión), esto se puede lograr
teniendo un óptimo diseño del aparejo de producción.
2) Instalaciones de sistemas artificiales: Las pérdidas de presión en tuberías verticales son

necesarias para el propio diseño de muchos sistemas artificiales de producción. Por ejemplo en
bombeo neumático, ya que al inyectar gas se tiene flujo multifásico en la tubería de producción.
3) Diseños de tuberías que transportan flujo multifásico: Determinar los diámetros y longitudes
óptimas de tuberías que transportan flujo multifásico. Un ejemplo de estos son tuberías de
producción, líneas de descarga, líneas de transporte de condensados, gasoductos cuando algún
químico está siendo inyectado para prevenir el congelamiento etc.
Como se ha estado mencionando el conocer las caídas de presión para flujo multifásico en
tuberías sirve para optimizar las producciones de los pozos petroleros y así obtener una mayor
rentabilidad de los proyectos.
II.3.1 Propiedades de los fluidos
Otro aspecto que es importante tratar cuando se habla de flujo multifásico en tuberías es el de las
propiedades de los fluidos (aceite, gas y agua), ya que en todo momento del transporte de fluidos,
la presión y la temperatura están variando por lo que las propiedades de los fluidos también
varían.
Algunas propiedades de los fluidos son las siguientes:
Factores de volumen: Bo, Bg, Bw, Bt.

Relaciones gas-aceite: R, Rs.
Densidades: ρo, ρg, ρw.
Viscosidades: μo, μg, μw.
Compresibilidades: Co, Cg, Cw.
Algunas definiciones de las propiedades de los fluidos son las siguientes:
Factor de volumen del aceite Bo: Es la relación de volumen de aceite más gas disuelto medido a
condiciones de yacimiento, entre el volumen de la misma cantidad de aceite medido a condiciones
estándar.
53
Capítulo II
Factor de volumen del gas Bg: Es la relación del volumen de una masa de gas medido a
condiciones de yacimiento o escurrimiento entre el volumen de la misma masa de gas medido a
condiciones estándar.
Condiciones estándar: Las condiciones entandar dependen de los reglamentos de los países.
Normalmente para los cálculos se consideran P =14.7 (lb / pg2) y T =60 °F.
Fase: Es la parte de un sistema que difiere en sus propiedades intensivas, de la otra parte del
sistema. Los sistemas de hidrocarburos generalmente se presentan en dos fases: Gaseosa y
líquida.
Gas disuelto: Es el conjunto de hidrocarburos que a condiciones atmosféricas constituyen un gas,

pero que forman parte de la fase líquida a condiciones de yacimiento ó de flujo.
Relación gas-aceite R: Son los pies cúbicos de gas producido por cada barril de aceite producido,
medidos ambos a condiciones estándar.
Relación de solubilidad Rs: Son los pies cúbicos de gas disuelto en el aceite a ciertas condiciones
de presión y temperatura, por cada barril de aceite en el tanque, medidos ambos volúmenes a
condiciones estándar.
Factor de compresibilidad Z: Se denomina también factor de desviación. Es un factor que se

introduce a la ley de los gases ideales para tomar en cuenta la desviación que experimenta un gas
real con respecto a un gas ideal.
Encogimiento: Es la disminución de volumen que experimenta una fase líquida por efecto de la
liberación del gas disuelto y por su contracción térmica. El factor de encogimiento es el reciproco
del factor de volumen ó de formación.
Viscosidad: Es la resistencia que presenta un fluido a fluir; a presión constante si la temperatura

aumenta la viscosidad del líquido disminuye, mientras que la viscosidad del gas aumenta.
A continuación se presentan algunas correlaciones que nos permiten calcular las propiedades de
los fluidos que más influyen en el cálculo de las caídas de presión para flujo multifásico.
Factor de volumen del aceite y gas (Bo y Bg)
Standing.
Bo = 0.972 + 0.000147 F 1.175( )

Donde:
⎛ γ gd ⎞
F = Rs ⎜⎜ ⎟⎟ + 1.25 T
γ
⎝ o ⎠
54
Capítulo II
Vazquez.
Bo = 1 + C1 Rs + (T − 60 ) (API / γ gs ) (C2 + C3 Rs )
Donde γ gs es el valor de la densidad relativa del gas normalizado a una presión de separación de
100 (lb / pg2) y se obtiene de la siguiente ecuación:
( )
γ gs = γ gp 1 + 5.912 × 10 − 5 (API)(Ts ) log (Ps / 114.7 )
Donde los coeficientes dependen de la densidad API del aceite
Coeficiente API < 30° API > 30°

C1 (4.677 x 10^ -4) (4.67x 10^ -4)
C2 (1.751 x 10^ -5) (1.1 x 10^ -5)
C3 (-1.811x 10^ -8) (1.337x 10^ -9)
De la ecuación de los gases reales se obtienen el factor de volumen del gas:
0.02825 Z ⋅ (T + 460 )
Bg =
P
Para estos casos se tiene la siguiente nomenclatura:
γ gd = Densidad relativa del gas disuelto
γ o = Densidad relativa del aceite
Bo = Factor de volumen del aceite
Bg = Factor de volumen del gas
Rs = Relación de solubilidad
Z = Factor de compresibilidad del gas

P = Presión
T = Temperatura
API = Densidad del aceite en grados API
Densidad del gas disuelto
γ gd = 0.25 + 0.02 API + R s × 10 −6 (0.6874 − 3.5864 API)
55
Capítulo II
Densidad del aceite
Saturado.
62.4 γ o + 0.01362 Rs ⋅ γ gd
ρo =
Bo
Bajo saturado.
⎛ ⎞
Co ⎜⎜ P −P ⎟⎟
b⎠
ρo = ρob ⋅ e ⎝
Donde:
Co =Compresibilidad del aceite
Pb =Presión de burbuja
ρo =Densidad del aceite
ρob =Densidad del aceite a la presión de burbuja
Densidad del gas
2.7044 P ⋅ γ gf
ρg =
Z ⋅ (T + 460 )
Donde:
γ gf =Densidad relativa del gas libre
ρg =Densidad del gas
Densidad del agua saturada
62.43
ρw =
Bw
Donde:
ρ w =Densidad del agua
B w =Factor del volumen del agua
56
Capítulo II
Viscosidad del aceite
Saturado
μo = a μbom
a = 10.715 (Rs + 100 )

− 0.515
b = 5.44 (Rs + 150 )

− 0.388
μom = 10D − 1
(
D = E T −1.163 )
E = 10F
F = 3.0324 − 0.02023 API
Bajo saturado
m
⎛P⎞
μo = μob ⎜⎜ ⎟⎟
⎝ Pb ⎠
⎡
⎢ −11.513 + ⎛⎜ − 8.98 × 10 − 5 P ⎞⎟⎤⎥
⎢
m = 2.6 P 1.187 ⋅ e ⎣ ⎝ ⎠⎥⎦
Donde:
μo =Viscosidad del aceite
Relación de solubilidad
Standing
1
⎡⎛ P ⎞ (0.0125 API−0.00091 T ) ⎤ 0.83
Rs = γ gd ⎢⎜ ⎟ 10 ⎥
⎣⎝ 18 ⎠ ⎦
Vazquez
Rs = C1 γ gs ⋅ PC 2 ⋅ exp [C3 API / (T + 460 )]
Donde:
[
γ gs = γ g 1 + 5.912 × 10 − 5 (API) (Ts ) log (Ps / 114.7 ) ]
57
Capítulo II
Los valores de los coeficientes se muestran en la siguiente tabla.
Coeficiente API < 30° API > 30°

C1 0.0362 0.0178
C2 1.0937 1.187
C3 25.724 23.931
II.3.2 Patrones de flujo
Al fluir dos fases simultáneamente en una tubería, las dos fases pueden fluir de diferentes formas,
cada una de estas formas presenta una distribución relativa de una fase con respecto a la otra,
constituyendo un patrón o tipo de flujo. Las siguientes figuras muestran los diferentes patrones de
flujo observados en tuberías horizontales y verticales por Beggs y Orkiszewsky figuras 2.3 y 2.4
respectivamente.
Figura 2.3 Patrones en flujo horizontal, observados por Beggs
58
Capítulo II
Figura 2.4 Patrones en flujo vertical, observados por Orkiszewsky
II.3.3 Colgamiento
El colgamiento se define como la relación entre el volumen de líquido existente en una sección de
la tubería a las condiciones de flujo y el volumen de la sección aludida. Esta relación de
volúmenes depende de la cantidad de líquido y gas que fluye simultáneamente en la tubería
(figura 2.5). Generalmente, la velocidad con que fluye el gas es diferente de la velocidad con que
fluye el líquido, propiciando un “resbalamiento” entre las fases.
Figura 2.5 Fenómenos del colgamiento y resbalamiento en una tubería
El término “resbalamiento” se usa para describir el fenómeno natural del flujo a mayor velocidad de
una de las dos fases. Las causas de este fenómeno son diversas; La resistencia al flujo por
fricción es mucho menor en la fase gaseosa que en la fase liquida. La diferencia de
compresibilidad entre el gas y el líquido, hace que el gas en expansión viaje a mayor velocidad
que el líquido. Cuando el flujo es ascendente o descendente, actúa la segregación gravitacional
ocasionando que el líquido viaje a menor velocidad que el gas en el primer caso y a mayor en el
segundo.
Para calcular las pérdidas de presión por elevación (carga hidrostática), es necesario predecir con
precisión el colgamiento considerando el resbalamiento entre las fases.
Las expresiones establecidas por Mukherjee y Brill son:
59
Capítulo II
⎡ NgvC 5 ⎤
(
YL = exp ⎢ C1 + C2 senθ + C3 sen2θ + C4 NLμ
2
⎥
NLvC 6 ⎥⎦
) ---------------2.25
⎢⎣
Donde:
0.25
⎛ 1 ⎞
N = 0.15726 μL ⎜⎜
2
Lμ
⎟
3 ⎟
---------------2.26
⎝ ρL σ ⎠
0.25
⎛ρ ⎞
NLv = 1.938 v sL ⎜ L ⎟ ---------------2.27
⎝σ⎠
0.25
⎛ρ ⎞
Ngv = 1.938 v sg ⎜ L ⎟ ---------------2.28
⎝σ⎠
σ = Tensión superficial del líquido
En el siguiente cuadro aparecen lo coeficientes de la ecuación 2.25 para los diferentes patrones
de flujo.
Dirección del
Tipo de Flujo C1 C2 C3 C4 C5 C6
Flujo
Horizontal o
Todos -0.38011 0.12988 -0.11979 2.34323 0.47569 0.28866
ascendente
Descendente Estratificado -1.33028 4.80814 4.17158 56.26227 0.07995 0.50489
Otros -0.51664 0.78981 0.55163 15.51921 0.37177 0.39395
Otro concepto que se usa con frecuencia en los cálculos de gradientes para flujo multifásico, es el
colgamiento sin resbalamiento (λ). Se define en la misma forma que (YL) pero se calcula a partir
de las condiciones de presión y temperatura de flujo existentes, considerando las producciones
obtenidas en superficie (qo y R):
qL 1
λ= = ---------------2.29
qL + qg qo ⋅ (R − R s ) ⋅ Bg
1+
5.615 (qo Bo + qw B w )
II.3.4 Velocidades superficiales
Es la velocidad que tendría cualquiera de las fases si se ocupara toda la tubería. Se define por las
expresiones siguientes:
qL 0.01191 (qo Bo + qw B w )
v sL = = ---------------2.30
Ap d2
60
Capítulo II
qg 0.002122 qo ⋅ (R − Rs ) ⋅ Bg
v sg = = ---------------2.31
Ap d2
qL + qg
vm = = v sL + v sg ---------------2.32
Ap
Donde:
Ap es el área de la sección transversal de la tubería. De estas ecuaciones se observa que:
v sL
λ= ---------------2.33
vm
II.3.5 Velocidad real
Aplicando el concepto de colgamiento, se puede obtener la velocidad real correspondiente a cada

fase:
qL qL v
vL = = = sL ---------------2.34
A L A p YL YL
qg qg v sg
vg = = = ---------------2.35
Ag A p (1 − YL ) (1 − YL )
Donde:
qL y qg =Gastos de líquido y gas
A L y A g =Áreas por donde pasa el líquido y el gas respectivamente
A p =Área total de la sección transversal
III.3.6 Densidad de la mezcla
La densidad real de la mezcla de fluidos se obtiene a partir del colgamiento con:
ρm = ρL YL + ρg (1 − YL ) ---------------2.36
Algunos autores calculan la densidad de la mezcla sin considerar el resbalamiento entre las fases,
esto es:
61
Capítulo II
ρns = ρL λ + ρg (1 − λ ) ---------------2.37
También la densidad puede obtenerse a partir de la siguiente expresión:
M
ρns = ---------------2.38
Vm
Donde:
M =Masa de la mezcla @ c.s por barril de aceite producido @ c.s

Vm =Volumen de la mezcla @ c.s por barril de aceite producido @ c.s
Los valores de M y Vm se obtienen de las ecuaciones siguientes:
M = Mo + Mg + Mw ---------------2.39
⎛ lb ⎞
Mo = 350.5 λ o ⎜⎜ o ⎟⎟ ---------------2.40
⎝ blso ⎠
⎛ lbg ⎞
Mg = 0.0764 λ g ⋅ R ⎜⎜ ⎟⎟ ---------------2.41
⎝ blso ⎠
⎛ lb ⎞
Mw = 350.5 λ w ⋅ Fw ⎜⎜ w ⎟⎟ ---------------2.42
⎝ blso ⎠
Sustituyendo las ecuaciones 2.40, 2.41 y 2.42 en la ecuación 2.39 se obtiene:
⎛ lb ⎞
M = 350.5 (λ o + λ w Fw ) + 0.0764 R ⋅ λ g ⎜⎜ m ⎟⎟ ---------------2.43
⎝ blso ⎠
Cálculo de Vm (pie3o+g+w @ c.s / blso @ c.s)
Vmo = 5.615 Bo ---------------2.44
Vmg = (R − R s ) Bg ---------------2.45
Vmw = Fw × 5.615 B w ---------------2.46

Vm = 5.615 (Bo + B w Fw ) + (R − Rs ) Bg ---------------2.47
Sustituyendo las ecuaciones 2.43 y 2.47 en la ecuación 2.38 se obtiene:
62
Capítulo II
350.5 (λ o + λ w Fw ) + 0.0764 R ⋅ λ g ⎛ lbm ⎞

ρns = ⎜ ⎟ ---------------2.48
5.615 (Bo + B w Fw ) + (R − Rs ) Bg ⎜⎝ pie3 ⎟⎠
III.3.7 Gasto másico
Se define por la siguiente expresión: wm = (lbm de líquido / seg). Y puede obtenerse con cualquiera
de las siguientes ecuaciones:
qo M
wm = ---------------2.49
86400
wm = wo + wg + w w ---------------2.50
Donde:
w o = ρo qo Bo / 15388
w w = ρ w qw B w / 15388
w g = ρg qo (R − Rs ) Bg / 86400
III.3.8 Viscosidad de la mezcla
Dependiendo del método que se aplica, se usan las siguientes ecuaciones para obtener la
viscosidad de la mezcla de fluidos:
λ (1− λ )
μns = μL μ g ---------------2.51
(1− YL )
μm = μLYL μ g ---------------2.52
Donde:
μns =Viscosidad de la mezcla sin resbalamiento
μm =Viscosidad de la mezcla con resbalamiento
Además:
μL = μo Fo + μ w Fw ---------------2.53
qo Bo Bo
Fo = = ---------------2.54
qw B w + qo Bo Bo + WOR B w
Fw = 1 − Fo ---------------2.55
63
Capítulo II
III.4 Flujo multifásico en tuberías verticales
El flujo multifásico vertical es encontrado en prácticamente toda el aparejo de producción de los

pozos de aceite, por lo tanto es necesario que él ingeniero de producción sea capaz de predecir
las caídas de presión de flujo multifásico vertical para seleccionar correctamente el aparejo de
producción por el cual producirá el pozo, predecir el gasto aportado por dicho pozo y optimizar el
diseño de sistemas artificiales de producción.
Mucho de los progresos encaminados hacia la solución del problema han sido hechos desde la
publicación del artículo de Poettmann y Carpenter en 1952. Muchas de las correlaciones
establecidas usan alguna forma de la ecuación general de energía.
El gradiente de presión para flujo vertical es la suma de los tres factores de contribución:
Gradiente de presión hidrostática, gradiente de presión por fricción y gradiente de presión por
aceleración. Los efectos de las reacciones químicas entre las fases son insignificantes; sin
embargo factores como viscosidad, tensión superficial y densidad son incluidos.
Las más importantes correlaciones para flujo multifásico vertical son aquellas propuestas por:
Duns y Rose, Orkiszewski, Hagederon y Brown así como Beggs y Brill. Otras que han contribuido
en menor proporción son aquellas propuestas por: Poettmann y Carpenter, Gilbert, Fancher y
Brown entre otras.
En esta tesis sólo se verán las correlaciones propuestas por Poettmann y Carpenter, así como la
propuesta por Beggs y Brill.
II.4.1 Correlación de Poettmann y Carpenter
Estos investigadores desarrollaron un método casi empírico en el cual relacionaron la ecuación

general de energía con datos de 34 pozos fluyentes de aceite y 15 pozos con bombeo neumático,
con diámetros de tubería de 2, 2½ y 3 pg. El aceite, gas y agua fueron considerados como una
sola fase, todas las pérdidas incluyendo el resbalamiento fueron combinadas dentro de un tipo de
factor llamado “factor de pérdida de energía”, el cual se asumió constante a lo largo de la longitud
total de la tubería. El término de pérdida de energía fue correlacionado con el numerador del
número de Reynolds.
La ecuación en unidades de campo dada por los autores fue la siguiente:
1 ⎡ ftp (qo M) ⎤
2
dP
= ⎢ρns + ⎥ ---------------2.56
dL 144 ⎢⎣ 7.413 × 1010 ρns ⋅ d5 ⎥⎦
Donde:
qo =Gasto de aceite (bpd)
M =Masa de la mezcla por barril de aceite producido (lbm @ c.s / blso @ c.s)
64
Capítulo II

ρns =Densidad de la mezcla sin resbalamiento (lbm @ c.s / pie3m @ c.s)
Los valores de M y ρns se obtienen de la siguiente manera:
M = 350.5 (γ o + γ w Fw ) + 0.0764 R ⋅ γ g
ρns = ρL λ + ρg (1 − λ )
El factor de fricción (ftp) está en función del numerador del número de Reynolds, el cual se obtiene
mediante la siguiente ecuación:
qo M
dvρns = 1.47 × 10 − 5 ---------------2.57
d
Con este valor se intersecta la curva y se lee el valor de ftp en el eje de las ordenadas. A
continuación se muestra la gráfica general para obtener el factor de fricción de las dos fases
(figura 2.6).
Figura 2.6 Datos de campo del factor de fricción en pozos fluyentes y con BN
Otra forma de obtener el factor de fricción ftp es mediante la siguiente ecuación:
ftp = 5.415 × 10 −3 − 5.723 × 10 −4 a + 1.848 × 10 −4 a2 + 3.5843 × 10 −6 a3 ---------------2.58
65
Capítulo II
Donde:
d × 106
a= ----------------2.59
qo M
Procedimiento de cálculo para obtener el valor de ∆P
1) Se conoce la presión a la entrada de la tubería P1 cuando L=0.

2) Se supone una caída de presión ∆P y se calcula P y P2 de la siguiente forma:
⎛ ΔP ⎞
P = P1 − ⎜ ⎟ Y P2 = P1 − ΔP
⎝ 2 ⎠
3) Se determinan las propiedades de los fluidos (Rs, Bo, Z, Bg, μo, ρo y ρg etc.) a las
condiciones de la presión y temperatura media obtenida P y T .
4) Calcular el colgamiento sin resbalamiento (λ).
5) Calcular la masa de la mezcla (M).
6) Determinar el valor de dvρns
7) Obtener el valor de ftp de la gráfica o de la ecuación.
8) Obtener el valor del gradiente de presión y determinar el ∆L.
ΔPsup uesta
ΔL =
ΔP
ΔL calculada
9) El valor de L ahora será L =Lanterior + ∆L.
10) Sustituir el valor de P1 = P2 y volver a repetir el procedimiento desde el paso 1 hasta
completar la profundidad del pozo.
III.4.2 Correlación de Beggs y Brill
Estos autores establecieron una correlación para calcular la distribución de la presión en tuberías
con flujo multifásico a partir de pruebas de laboratorio, el método es aplicable a flujos horizontales,
direccionales y verticales.
La ecuación general del gradiente de presión es la siguiente:
66
Capítulo II
⎡ g ⋅ senθ ⋅ ρ ftp ⋅ ρns ⋅ v m ⎤

2
⎢ m
+ ⎥
dP 1 ⎢ gc 5.362 d ⎥
=
dL 144 ⎢ v ⋅v ⋅ρ ⎥
⎢ 1 − m sg m ⎥
⎣⎢ 144 gc ⋅ P ⎦⎥
El término del denominador es el término de caída de presión debido a la aceleración que muchas
veces es despreciable por lo que la ecuación más simple es la siguiente:
1 ⎡ g ⋅ senθ ⋅ ρm ftp ⋅ ρns ⋅ v m ⎤

2
dP
= ⎢ + ⎥ ----------------2.60
dL 144 ⎢⎣ gc 5.362 d ⎥⎦
El factor de fricción para las dos fases, se obtiene de la siguiente ecuación:
⎛f ⎞
ftp = ⎜⎜ tp ⎟⎟ fn ----------------2.61
⎝ fn ⎠
Donde:
⎛ ftp ⎞
⎜⎜ ⎟⎟ = eS ----------------2.62
⎝ fn ⎠
ln x
S= ----------------2.63
− 0.0523 + 3.182 (ln x ) − 0.8725 (ln x ) + 0.01853 (ln x )
2 4
λ
x= ----------------2.64
YL (0 )
−2
⎡ ⎛ Nre ⎞⎤
fn = ⎢2 log ⎜⎜ ⎟⎟⎥ ----------------2.65
⎣ ⎝ 4 . 5223 log (Nre ) − 3 . 8215 ⎠⎦
Además:
124 d ⋅ v m ⋅ ρns
Nre = ----------------2.66
μns
El valor de S llega a ser indefinido cuando el valor de x se encuentra en el intervalo (1≤ x ≤1.2);
para este intervalo la función S es calculada de la siguiente manera:
S = ln (2.2 x − 1.2) ----------------2.67
67
Capítulo II
De sus observaciones los autores elaboraron un mapa de patrones de flujo en función de λ

(colgamiento sin resbalamiento) y el número de Froude, el cual es definido de la siguiente manera:
2
wm
NFR = 7734.9 2
----------------2.68
ρns d5
El mapa de patrones de flujo se muestra a continuación.
Patron de Flujo Condiciones

Segregado λ ≤ 0.01 y NFR < L1 o λ ≥ 0.01 y NFR < L2
Transición λ ≥ 0.01 y L2 < NFR ≤ L3
Intermitente 0.01 ≤ λ < 0.4 y L3 < NFR ≤ L1 o λ > 0.4 y L3 < NFR ≤ L4
Distribuido λ > 0.4 y NFR ≥ L1 o λ ≥ 0.4 y NFR > L4
En donde los parámetros de correlación L1, L2, L3 y L4, se obtienen de las siguientes ecuaciones:
L1 = 316 λ 0.302
L 2 = 0.0009252 λ − 2.4684
L3 = 0.10 λ −1.4516
L 4 = 0.5 λ − 6.738
Para flujo vertical se determina el colgamiento que existiría si la tubería fuese horizontal y luego se
corrige por la inclinación real de la tubería que es de ± 90° para flujo vertical.
YL = YL (0 ) ψ ----------------2.69
a λb
YL (0 ) = c
----------------2.70
NFR
ψ = 1 + C [sin (1.8 θ) − 0.333 sin3 (1.8 θ)] ----------------2.71
(
C = (1 − λ ) In d λe NLV NFR
f g
) ----------------2.72
Donde:
YL (0) =Colgamiento para flujo horizontal (θ=0)
NFR =Número de Froude
λ =Colgamiento sin resbalamiento

θ =Angulo de inclinación de la tubería con respecto a la horizontal
68
Capítulo II
Donde los coeficientes están en función del patrón de flujo y se muestran a continuación:
Patron de Flujo a b c d e f g
Segregado 0.98 0.4846 0.0868 - - - -
Intermitente 0.845 0.5351 0.0173 - - - -
Distribuido 1.065 0.5824 0.0609 - - - -
Segregado
- - - 0.011 -3.768 3.539 -1.614
(Cuesta Arriba)
Intermitente
- - - 2.96 0.305 -0.4473 0.0978
(Cuesta Arriba)
Distribuido
- - - C=0 ψ=1
(Cuesta Arriba)
Todos los Patrones
- - - 4.7 -0.3692 0.1244 -0.5056
(Cuesta Abajo)
Para flujo transitorio, el colgamiento se obtiene de la siguiente forma:
YL (0) =A × YL (0) (Segregado) + B1 × YL (0) (Intermitente)
Donde:
L 3 − NFR
A=
L3 − L2
B1 = 1 − A

2) Se supone una caída de presión ∆P y se calcula P y P2 de la siguiente forma:
⎛ ΔP ⎞
P = P1 − ⎜ ⎟ Y P2 = P1 − ΔP
⎝ 2 ⎠
condiciones de la presión y temperatura media obtenidas P y T .
4) Calcular el colgamiento sin resbalamiento (λ), el número de Froude (NFR), los parámetros
de correlación L1, L2, L3, L4 y determinar el patrón de flujo.
5) Para las condiciones dadas calcular el valor de la densidad de la mezcla sin resbalamiento
(ρns).
6) Obtener el valor del colgamiento (YL) con las ecuaciones dadas anteriormente.
69
Capítulo II
7) Calcular el valor de la densidad real de la mezcla (ρm).

8) Determinar el valor de la viscosidad de la mezcla sin resbalamiento (μns).
9) Obtener el valor del Nre usando la correlación propuesta por los autores.
10) Calcular el valor de fn y ftp.
11) En dado caso que el término por aceleración no sea despreciable, calcularlo.
ΔPsup uesta
ΔL =
ΔP
ΔL calculada
completar la profundidad del pozo.
III.4.3 Curvas de gradiente de presión
El análisis del comportamiento del flujo vertical se puede hacer con el auxilio de las gráficas de
gradiente de presión desarrolladas por Gilbert y Kermit Brown. Gilbert da una solución empírica al
problema del flujo bifásico vertical. Efectuó mediciones de la caída de presión en tuberías de
producción bajo distintas condiciones y obtuvo una familia de curvas (figura 2.7).
Los parámetros que midió en un número grande de pozos fluyentes fueron:
Profundidad de la tubería (pie)
Diámetro de la tubería (pg)
Producción de líquidos (bpd)
Relación gas-líquido (pie3 / bl)
Presión en la cabeza del pozo (Ib / pg2)
Presión de fondo fluyendo (Ib / pg2)
70
Capítulo II
Figura 2.7 Curvas de distribución de presión en una tubería vertical
Se considera que la presión de fondo fluyendo depende únicamente de las otras cinco variables.
En la figura 2.7, las curvas a, b, c y d corresponden a diferentes presiones en la cabeza del pozo
(A, B, C y D). Cada una de estas curvas representa la distribución de presión a lo largo de la
tubería de producción para un pozo con: Gasto, relación gas-líquido y diámetro de tubería
determinados.
Del punto B de la curva b, Gilbert trazó una vertical hasta intersectar la curva a y sobreponiendo
éstas obtuvo que la curva b coincidía con una sección de la curva a. Hizo lo mismo con las otras
curvas y concluyó que las curvas a, b, c y d son realmente parte de una misma curva (figura 2.8),
con presiones en la cabeza del pozo A, B, C y D correspondientes a las marcadas en la
figura. 2.7. La curva c por ejemplo, es la curva de la figura 2.8 con el punto x tomando la
profundidad como cero.
Figura 2.8 Curva de gradiente de presión para flujo multifásico
71
Capítulo II
El procedimiento para utilizar las curvas de gradiente de presión es el siguiente:
1) Se escoge la curva de gradiente que corresponda a los valores de qo, RGL, diámetro de
T.P, API, fracción de agua etc.
2) Si queremos conocer la presión de fondo fluyendo (Pwf), partimos de que conocemos el
valor de la presión en la cabeza del pozo (Pwh), este valor es localizado en el eje de las
abscisas (eje de presiones).
3) A partir de este valor se traza una línea recta hasta intersectar la curva y se lee el valor de
la profundidad en el eje de las ordenadas (eje de profundidades).
4) A este valor de profundidad, se le suma la longitud total de la tubería de producción.
5) En este valor de profundidad se traza una línea recta hasta intersectar de nuevo la curva y
se lee el valor de la presión en el eje de las abscisas.
6) Este valor es el valor de la presión de fondo fluyendo (Pwf).
7) En el caso de que lo que se quiera conocer sea la presión en la cabeza del pozo (Pwh), el
procedimiento será el mismo con la diferencia de que la longitud total de la tubería de
producción será restada en lugar de ser sumada.
III.5 Flujo multifásico en tuberías horizontales

Para flujo horizontal, el gradiente de presión debido al cambio de elevación es igual a cero por lo
que la ecuación de balance de energía queda de la siguiente forma:
⎛ ΔP ⎞ ⎛ ΔP ⎞ ⎛ ΔP ⎞
⎜ ⎟ =⎜ ⎟ +⎜ ⎟
⎝ ΔL ⎠T ⎝ ΔL ⎠F ⎝ ΔL ⎠ AC
La ecuación anterior se puede escribir como:
2 2
⎛ ΔP ⎞ ρ v ρ Δv
⎜ ⎟ = ftp m m + m m
⎝ ΔL ⎠T 2gc ⋅ d 2gc ⋅ ΔL
En donde ftp, ρm y vm se refieren a la mezcla y son definidos en forma distinta por los autores de
las diferentes correlaciones.
III.5.1 Correlación de Bertuzzi, Tek y Poettmann
Despreciando el término de caída de presión debido a la aceleración, la correlación propuesta por

Bertuzzi, Tek y Poettmann queda de la siguiente manera:
72
Capítulo II
2
⎛ ΔP ⎞ 174.158 ftp ⋅ w m
⎜ ⎟= ----------------2.73
⎝ ΔL ⎠ ρns d5
Para obtener el factor de fricción para las dos fases se utiliza la siguiente correlación:
φ = (Nreg ) (NreL )
a b
----------------2.74
Donde:
a = ψ / (1 + ψ )
(0.1 ψ )
b = 1/ e
ψ = w g / wL
Para:
0 ≤ φ ≤ 500
log ftp = 1.225 ψ − 0.06561 (log φ) − 0.37 ----------------2.75
φ ≥ 10,000
log ftp = 0.49 ψ − 0.12616 (log φ) − 1.702 ----------------2.76
500 ≤ φ ≤ 10,000
log ftp = F500 − 0.6561 y + (1.1056 + 1.7723 ΔF ) y 2 − (0.46214 + 0.90817 ΔF ) y 3 ----------------2.77
Donde:
y = Log φ − 2.699
∆F = F10,000 − F500
F10,000 = Log (ftp), cuando φ =10,000
F500 = Log (ftp), cuando φ =500
Los números de Reynolds para cada una de las fases se obtienen de las siguientes ecuaciones:
wL
NreL = 22737 ----------------2.78
d μL
wg
Nreg = 22737 ----------------2.79
d μg
73
Capítulo II
El factor de fricción (ftp) también puede ser encontrado a partir de la siguiente gráfica (figura 2.9),
observándose que la gráfica está en función de la relación de los gastos másicos ψ = wg / wL.
Figura 2.9 Gráfica del factor de fricción (Bertuzzi)

2) Se supone una caída de presión ∆P y se calcula P de la siguiente forma:
⎛ ΔP ⎞
P = P1 − ⎜ ⎟ Y P2 = P1 − ΔP
⎝ 2 ⎠
condiciones de la presión media obtenida P y T .
4) Obtener los valores del gasto másico (wm) y el colgamiento sin resbalamiento (λ).
5) Determinar el factor de fricción para las dos fases (ftp).
6) Obtener el valor de la densidad de la mezcla sin resbalamiento (ρns).
ΔPsup uesta
ΔL =
ΔP
ΔL calculada
74
Capítulo II

completar la longitud total de la tubería obteniéndose la presión final P2.
III.5.2 Correlación de Beggs y Brill
Los autores determinaron la siguiente ecuación para calcular el gradiente de presión:
2
ΔP 43.539 ftp ⋅ w m
= ----------------2.80
ΔL ρns d5
Esta ecuación se obtiene al despreciar el término de caída de presión por aceleración en la

ecuación original propuesta por los autores.
Para esta correlación el factor de fricción se obtiene de de la siguiente ecuación:
⎛f ⎞
ftp = ⎜⎜ tp ⎟⎟ fn ----------------2.81
⎝ fn ⎠
Donde:
⎛ ftp ⎞
⎜⎜ ⎟⎟ = eS ----------------2.82
⎝ fn ⎠
ln x
S= ----------------2.83
− 0.0523 + 3.182 (ln x ) − 0.8725 (ln x ) + 0.01853 (ln x )
2 4
λ
x= ----------------2.84
YL (0 )
−2
⎡ ⎛ Nre ⎞⎤
fn = ⎢2 log ⎜⎜ ⎟⎟⎥ ----------------2.85
⎣ ⎝ 4 . 5223 log (Nre ) − 3 . 8215 ⎠⎦
En donde:
124 d ⋅ v m ⋅ ρns
Nre = ----------------2.86
μns
El valor de S llega a ser indefinido cuando el valor de x se encuentra en el intervalo (1≤ x ≤1.2);
para este intervalo la función S es calculada de la siguiente manera:
75
Capítulo II
S = ln (2.2 x − 1.2) ----------------2.87
De sus observaciones los autores elaboraron un mapa de patrones de flujo en función del λ y el
número de Froude el cual es definido de la siguiente manera:
2
wm
NFR = 7734.9 2
----------------2.88
ρns d5
El mapa de patrones de flujo se muestra a continuación.
Patron de Flujo Condiciones

Segregado λ ≤ 0.01 y NFR < L1 o λ ≥ 0.01 y NFR < L2
Transición λ ≥ 0.01 y L2 < NFR ≤ L3
Intermitente 0.01 ≤ λ < 0.4 y L3 < NFR ≤ L1 o λ > 0.4 y L3 < NFR ≤ L4
Distribuido λ > 0.4 y NFR ≥ L1 o λ ≥ 0.4 y NFR > L4
En donde los parámetros de correlación L1, L2, L3 y L4 se obtienen de las siguientes ecuaciones:
L1 = 316 λ 0.302
L 2 = 0.0009252 λ − 2.4684
L3 = 0.10 λ −1.4516
L 4 = 0.5 λ − 6.738
El cálculo del colgamiento real de líquido se obtiene de la siguiente expresión general.
a λb
YL (0 ) = c
----------------2.89
NFR
Donde:
YL (0 ) =Colgamiento para flujo horizontal
λ =Colgamiento sin resbalamiento

NFR =Número de Froude
Donde los coeficientes están en función del régimen de flujo.
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
76
Capítulo II
Para flujo transitorio, el colgamiento se obtiene de la siguiente forma:
YL (0) = A × YL (0) (Segregado) + B × YL (0) (Intermitente)
En donde:
L 3 − NFR
A=
L3 − L2
B = 1− A

2) Se supone una caída de presión ∆P y se calcula P de la siguiente forma:
⎛ ΔP ⎞
P = P1 − ⎜ ⎟ Y P2 = P1 − ΔP
⎝ 2 ⎠
condiciones de la presión media obtenida P y T .
4) Calcular el colgamiento sin resbalamiento (λ), el número de froude (NFR) y los parámetros
de correlación L1, L2, L3, L4 y determinar el patrón de flujo.
5) Determinar el colgamiento YL (0) y la densidad de la mezcla sin resbalamiento (ρns).
6) Obtener el factor de fricción de las dos fases (ftp).
ΔPsup uesta
ΔL =
ΔP
ΔL calculada
completar la longitud total de la tubería obteniéndose la presión final P2.
III.5.3 Curvas de gradiente de presión
Al igual que las curvas de gradiente de presión para flujo multifásico vertical desarrolladas por
Gilbert y Kermit Brown, existen curvas de gradiente de presión para flujo multifásico en tuberías
horizontales, el procedimiento de uso es el mismo para ambos casos.
77
Capítulo II
Algunas compañías petroleras tienen un apartado de curvas de gradiente de presión para flujo
multifásico tanto vertical como horizontal basándose en simuladores que utilizan muchas de las
correlaciones para flujo multifásico existentes en la industria.
II.6 Flujo por espacio anular
Existen muchos pozos en el mundo que producen a gastos muy altos, muchos de los cuales
producen a través del espacio anular. Muchas de las correlaciones pueden ser adaptadas para
flujo por espacio anular, basta con sustituir el diámetro de la tubería de producción (T.P) por el
diámetro hidráulico, recordemos la definición del radio hidráulico.
área de la sección transversa l

rh =
perímetro mojado
Para flujo por espacio anular, tenemos lo siguiente:
rh =
( 2 2
)
π dT.R - dT.P / 4 dT.R - dT.P
= ----------------2.90
π (dT.R + dT.P ) 4
Donde:
dT.R =Diámetro interior de la tubería de revestimiento (T.R).
dT.P =Diámetro exterior de la tubería de producción (T.P).
Para una circunferencia rh =(d / 4), esto implica que el diámetro hidráulico sea igual a 4rh, por
consiguiente obtenemos:
dh = dT.R − dT.P ----------------2.91
Cuando el diámetro aparece a la 1a potencia podemos hacer la sustitución directa de (dT.R − dT.P)
por el diámetro de la tubería, como lo es en el número de Reynolds o en el factor de rugosidad
relativa. En muchas correlaciones el diámetro aparece a la 5ta potencia, por lo que debemos hacer
la siguiente sustitución.
d5 = (dT.R + dT.P ) (dT.R − dT.P )

2 3
----------------2.92
La mayoría de las correlaciones pueden ser utilizadas para flujo por espacio anular, haciendo la
sustitución respectiva para el diámetro de la tubería visto anteriormente. Únicamente unas cuantas
correlaciones han sido propuestas solamente para flujo por espacio anular, como la correlación de
Baxendell.
78
Capítulo II
III.7 Flujo a través de restricciones
Aunque las principales pérdidas de presión en el sistema integral de producción ocurren en el

yacimiento, tubería de producción y línea de descarga. Las pérdidas de presión en restricciones
pueden ser muy importantes en algunos pozos. Los principales tipos de restricciones son:
• Válvulas de seguridad sub-superficiales.

• Estranguladores superficiales o de fondo.
• Válvulas y conexiones.
El control desde la cabeza del pozo con un estrangulador, puede ser necesario por alguna de las
siguientes razones:
• Mantener un gasto de producción deseado.

• Mantener una contrapresión suficiente para prevenir la entrada de arena y en algunas
ocasiones sirve para regular el depósito de parafinas.
• Protección del equipo superficial.
• Prevenir la conificación de agua o gas.
El flujo a través de una restricción puede ser crítico (sónico) o subcrítico (subsónico). Si el flujo es
crítico un disturbio en la presión delante de la restricción (presión de descarga Pd) no causara
ningún efecto sobre el gasto a través de la restricción o la presión antes de esta (presión de
entrega Pe). Como unos de los principales objetivos del estrangulador es controlar el gasto,
usualmente este será dimensionado para que el flujo crítico exista.
Una regla general para distinguir entre flujo crítico y subcrítico es que si la relación de presiones
(Pd / Pe) es menor o igual a 0.528 entonces el flujo será crítico. Este valor es una aproximación
mucho más cercana para una sola fase (gas) que para dos fases (gas y aceite).
“Es de vital importancia saber que en el desarrollo de sus correlaciones, los autores han
supuesto diversas relaciones de presión critica (Pd / Pe), establecer un valor fijo para dicha
relación implica una simplificación que indudablemente se reflejará en la exactitud de las
predicciones que se obtengan al aplicar las correlaciones. Por lo anterior, es recomendable
que al utilizar una correlación se investiguen las fronteras de flujo crítico, esta es una de las
razones por las que el flujo a través de estranguladores es un tema que se debe de tratar
con mucho cuidado”.
II.7.1 Flujo multifásico a través de estranguladores
La siguiente ecuación puede ser usada para determinar la relación entre Pwh, qL, y RGL cuando el
flujo es crítico. Estas es una ecuación empírica, los coeficientes y exponentes pueden variar de
campo a campo o de pozo a pozo.
79
Capítulo II
qL ⋅ RGLc
Pwh = b ----------------2.93
da
Donde:
Pwh = Presión en la cabeza del pozo (lb / pg2)

qL = Gasto de líquido (bpd)
RGL = Relación gas-líquido (pie3 / bls)
d = Diámetro del estrangulador (pg)
Investigador a b c
Ros 2 4.25 x 10-3 0.5
Gilbert 1.89 3.86 x 10-3 0.546
Baxendell 1.93 3.12 x 10-3 0.546
Achong 1.86 1.54 x 10-3 0.65
80
Capítulo III
CAPITULO III
BASES PARA EL ANÁLISIS NODAL
III.1 ¿Qué es el análisis nodal?
Hasta el momento se han tratado los temas del cálculo de las caídas de presión tanto en el
yacimiento como en las tuberías de producción; ¿Pero cómo saber si el yacimiento o las tuberías
instaladas están operando bajo las condiciones idóneas u óptimas? Esta es una pregunta
fundamental para los ingenieros de producción, pues el objetivo de éstos debe ser obtener la
máxima producción de cada uno de los pozos de un campo petrolero, con la menor inversión
económica posible. Para conseguir este objetivo, él ingeniero de producción debe de utilizar
herramientas técnicas para el análisis del sistema integral de producción (SIP), una de estas
herramientas técnicas es el análisis nodal.
La razón fundamental de someter un SIP a la técnica de análisis nodal es porque éste involucra en
sus cálculos todos los elementos del sistema, permite determinar el efecto de su variación de la
capacidad de transporte y tener una imagen en conjunto del comportamiento del pozo. Desde la
perspectiva de evaluación esto es posible; sin embargo en condiciones de diseño (sin restricciones
económicas) es factible dimensionar el SIP en su conjunto para obtener la capacidad de transporte
requerida o bien, la capacidad de transporte idónea, teniendo como límite tan sólo la capacidad de
afluencia del yacimiento.
Las principales pérdidas de presión en el sistema integral de producción son las siguientes:
• Pérdidas de presión en el medio poroso. Representan entre el 10 y el 50% de las pérdidas

totales.
• Pérdidas de presión en la tubería de producción (T.P). Representan entre el 30 y 80% de

las pérdidas totales.
• Pérdidas de presión en la línea de descarga (L.D). Representan entre el 5 y 30% de las

pérdidas totales.
El procedimiento para realizar análisis nodal a un sistema integral de producción, requiere primero
seleccionar un nodo y después calcular la presión en éste, empezando con las presiones
constantes que son usualmente Pr y Psep. El nodo puede ser seleccionado en cualquier punto del
sistema.
En la figura 3.1 se muestra un sistema integral de producción sencillo, en el que se aprecian los
nodos y su respectiva posición. En la figura 3.2 se muestra un sistema integral de producción más
complejo, principalmente instalado en pozos costa afuera.
81
Capítulo III
Figura 3.1 Localización de nodos para un sistema de producción simple (pozos terrestres)
Figura 3.2 Localización de nodos para un sistema de producción complejo (pozos costa afuera)
Las expresiones para el cálculo de presión en el nodo pueden ser expresadas de la siguiente
manera:
82
Capítulo III
Inflow
Pr − ΔP(componente s corriente abajo ) = Pnodo
Outflow
Psep + ΔP(componente s corriente arriba) = Pnodo
Los dos criterios que se deben de satisfacer son los siguientes:
• El gasto de entrada en el nodo debe ser igual al gasto de salida

• Únicamente una presión puede existir en el nodo para un gasto determinado
El gasto y la presión que satisfacen los requerimientos previos pueden ser concluidos de la gráfica
de Presión vs. Gasto en el nodo. La intersección de las curvas de inflow y outflow ocurre en el
gasto que satisface los criterios anteriores. Este gasto será la capacidad de producción del
sistema para un conjunto particular de componentes. Para investigar el efecto en los cambios de
cualquier componente sobre la capacidad de producción de todo el sistema, nuevas curvas de
inflow u outflow pueden ser generadas para cada cambio.
Si el cambio es hecho en un componente corriente abajo del nodo solución, es decir los cálculos
de ∆P están en la ecuación de presión de la curva de inflow, la curva de outflow no sufrirá ningún
cambio y por consiguiente no necesita ser recalculada, inversamente si el cambio es hecho en un
componente corriente arriba del nodo, la curva de inflow permanecerá sin cambios. Esto nos
permite aislar los efectos del cambio de cualquier componente en el SIP. Este método nos permite
determinar si el SIP está funcionando correctamente y también para el diseño de nuevos sistemas.
III.2 Elección del nodo solución
La elección del nodo solución para pozos fluyentes o inyectores depende del componente que se
desee evaluar; esto es, que su análisis muestre convenientemente la respuesta del sistema a las
condiciones dadas y a las que se establezcan como supuestas, de tal forma que se pueda
identificar con certeza el problema y planear la solución técnica para su posterior ejecución.
III.2.1 Fondo del pozo como nodo solución
Probablemente la posición más común como nodo solución es el fondo del pozo, el fondo del pozo
se toma a la profundidad media del intervalo productor o en su caso a la profundidad media del
intervalo disparado.
El sistema entero es dividido en dos componentes, el yacimiento es uno y el otro es el sistema

total de tuberías (T.P + L.D). Suponiendo que no existe ninguna restricción en las tuberías, las
únicas caídas de presión son dadas a continuación:
83
Capítulo III
Las expresiones para las curvas de inflow y outflow son las siguientes:
Pnodo = Pwf
Inflow
Pr − ΔPres = Pnodo
Outflow
Psep + ΔPL.D + ΔPT.P = Pnodo
Procedimiento de cálculo:
1) Partiendo de la presión media del yacimiento (recordando que la Pr es una valor constante, no
varía con el gasto), asumir varios gastos de producción y calcular la presión de fondo fluyendo
(Pwf ) para dichos gastos de producción y construir la curva de IPR.
2) Partiendo de la presión del separador (recordando que Psep es un valor constante, no varía con
el gasto), asumir varios gastos de producción y obtener la presión en la cabeza del pozo (Pwh)
necesaria para mover los fluidos a través de la línea de descarga (L.D) hasta el separador,
utilizando una correlación de flujo multifásico en tubería horizontal o curvas de gradiente de
presión.
3) Usando los mismos gastos asumidos en el punto 2 y las correspondientes presiones en la

cabeza del pozo (Pwh), determinar la presión de fondo fluyendo (Pwf) necesaria para mover los
fluidos a través de la tubería de producción (T.P) hasta la cabeza del pozo, utilizando una
correlación de flujo multifásico en tubería vertical o curvas de gradiente de presión.
4) Graficar la curva de IPR calculada en el paso 1(curva de inflow). Sobre esa misma gráfica, se
grafican los gastos y las presiones de fondo fluyendo (Pwf = Pnodo) obtenidas en el paso 3 (curva de
outflow).La intersección de las curvas mostrará el gasto posible para este sistema, este gasto no
es el gasto máximo u óptimo del sistema, pero es el gasto al cual el pozo produciría para las
condiciones actuales del yacimiento y del sistema total de tuberías que está instalado.
¿Por qué escogemos el fondo del pozo como nodo solución? Se puede observar que el
componente yacimiento ha sido aislado del sistema total de tuberías. Por consiguiente un cambio
en la presión media del yacimiento ocasionará un cambio en el gasto, esto se puede ver en la
siguiente gráfica (figura 3.3).
84
Capítulo III
Figura 3.3 Predicción de gastos futuros
Para casos prácticos de campo, es muy común que las condiciones de la curva de outflow (caídas
de presión en T.P y L.D) también cambien con el depresionamiento. Los principales parámetros
que cambian con el depresionamiento son la relación gas-aceite R y el corte de agua Fw. La
relación de producción gas-aceite R incrementara en cualquier yacimiento en el cual la presión
decline por debajo de la presión de burbuja Pb y el corte de agua incrementara con el tiempo si
existe un acuífero activo o si el agua está siendo inyectada al yacimiento con el objetivo de
mantener la presión constante. Por lo tanto una nueva curva de outflow debe ser construida
cuando existe una caída en la presión media del yacimiento Pr.
Existen otros casos en donde el fondo del pozo como nodo solución es el mejor para ilustrar el
efecto de ciertas variables. El más importante es mostrar el efecto del diámetro de la tubería de
producción (T.P), ya que en algunas ocasiones hasta un 80% de las caídas de presión total se dan
en transportar los fluidos desde el fondo del pozo hasta la cabeza. El procedimiento para aislar el
efecto del diámetro de la T.P será el mismo visto anteriormente, con la diferencia que la presión en
la cabeza del pozo se mantendrá constante (las caídas de presión en L.D no se toman en cuenta).
Otro aspecto puede ser, mostrar el efecto de una buena o mala terminación del pozo, mostrar la
diferencia en el gasto esperado debido a la estimulación del pozo o debido a la remoción del daño
a la formación como se observa en la figura 3.4, la cual muestra un ejemplo de un pozo teniendo
una eficiencia de flujo original de 0.4 (con daño), una eficiencia de 1 (todo el daño ha sido
removido) y una eficiencia de 1.4 obtenida por un proceso de estimulación.
85
Capítulo III
Figura 3.4 Efectos del mejoramiento de la EF sobre los gastos de producción
III.2.2 Cabeza del pozo como nodo solución
El nodo solución más comúnmente usado después del fondo del pozo es la cabeza del pozo. De
nuevo el sistema es dividido en dos componentes para obtener las condiciones óptimas de
producción. El separador y la línea de descarga (L.D) son considerados como un componente y el
yacimiento y la tubería de producción (T.P) son el otro componente.
Pnodo = Pwh
Inflow
Pr − ΔPres − ΔPT.P = Pnodo
Outflow
Psep + ΔPL.D = Pnodo
Al graficar las presiones en la cabeza del pozo (Pwh) vs. Gastos de aceite (qo) para los dos
componentes, se obtiene una gráfica como la siguiente (figura 3.5).
86
Capítulo III
Figura 3.5 Cabeza del pozo como nodo solución.
(Pwf) para dichos gastos de producción utilizando la ecuación que mejor describa el
comportamiento de flujo en el yacimiento.
2) Utilizando los mismos gastos y las presiones de fondo fluyendo (Pwf) obtenidas en el paso 1, se
determinan las presiones en la cabeza del pozo (Pwh) necesarias para mover los fluidos a través
de la tubería de producción (T.P) hasta la cabeza del pozo, utilizando una correlación de flujo
multifásico en tubería vertical o curvas de gradiente de presión.
presión.
4) Graficar los gastos y las presiones en la cabeza del pozo (Pwh = Pnodo), obtenidas en el paso 2
(curva de inflow). Sobre esa misma gráfica, graficar los gastos y las presiones en la cabeza del
pozo (Pwh = Pnodo), obtenidas en el paso 3 (curva de outflow). La intersección de las curvas
mostrará el gasto posible para este sistema, este gasto no es el gasto máximo u óptimo del
sistema, pero es el gasto al cual el pozo produciría para las condiciones dadas del yacimiento y
del sistema total de tuberías que está instalado.
87
Capítulo III
¿Por qué escogemos la cabeza del pozo como nodo solución? Al tomar este nodo como nodo
solución se puede observar que la línea de descarga (L.D) puede ser aislada ya que comúnmente
la presión del separador (Psep) es una presión constante, por lo tanto es fácil mostrar el efecto en
el cambio de los diámetros de la línea de descarga. A continuación se muestra una gráfica que
ejemplifica esto (figura 3.6).
Figura 3.6 Efecto de la variación del diámetro de la L.D sobre los gastos de producción
Para la selección del diámetro óptimo de la línea de descarga (L.D), es necesario tener en mente
que el gasto alcanza un valor máximo a partir del cual el empleo de tuberías de descarga de
mayor diámetro es innecesario. Esto se debe a que otro elemento del sistema (yacimiento, T.P o
estrangulador) limita el ritmo de producción, además de que la elección del diámetro requiere un
análisis económico.
III.2.3 Separador como nodo solución
La selección de la presión de separación es crítica cuando diseñamos un sistema de bombeo

neumático o cuando la presión del gas que sale del separador debe ser incrementada para que el
gas pueda entrar a líneas de alta presión como pueden ser líneas de ventas. La presión de
separación controla la presión de succión del compresor y esta presión está directamente
relacionada con los requerimientos de potencia del compresor.
Sin embrago la presión de separación (Psep) no debe ser indiscriminadamente disminuida o

aumentada sin antes analizar el comportamiento total del sistema, en particular de la línea de
descarga (L.D). Un intuitivo pensamiento es que disminuyendo la presión de separación Psep se
obtendrá un incremento en la producción y esto no es cierto del todo. Existen numerosos casos en
los cuales una disminución en la presión de separación (Psep) no causa efecto en el gasto de
producción aún cuando la productividad de los pozos es muy alta. La razón para esto es que la
T.P o L.D deben ser las causas de la baja productividad del sistema.
88
Capítulo III
En pozos de baja productividad, el yacimiento por si mismo puede ser la restricción y un cambio
en la presión de separación Psep tendrá un efecto muy pequeño sobre la producción. En el bombeo
neumático, una disminución en la presión de separación (Psep) ocasionará una disminución de la
presión en la cabeza del pozo, por lo que el mismo gasto puede ser obtenido con menos gas de
inyección.
Pnodo = Psep
Inflow
Pr − ΔPres − ΔPT.P − ΔPL.D = Pnodo
Outflow
Psep = Pnodo
(Pwf) para dichos gastos de producción, utilizando la ecuación que mejor describa el
de la tubería de producción (T.P) hasta la cabeza del pozo, utilizando correlaciones de flujo
multifásico para tuberías verticales o curvas de gradiente de presión.
3) Utilizando los mismos gastos y las presiones en la cabeza del pozo (Pwh) obtenidas en el paso
2, se determinan las presiones al final de la línea de descarga (L.D), las cuales serán las presiones
a la entrada del separador (Psep), utilizando una correlación de flujo multifásico en tuberías
horizontales o curvas de gradiente de presión.
4) Graficar los gastos y las presiones a la entrada del separador (Psep = Pnodo), obtenidas en el
paso 3 (curva de inflow). Sobre esa misma gráfica, graficar la presión del separador (Psep = Pnodo),
la cual es independiente del gasto y depende de otros factores (curva de outflow). La intersección
de las curvas mostrará el gasto posible para este sistema, este gasto no es el gasto máximo u
óptimo del sistema, pero es el gasto al cual el pozo produciría para las condiciones dadas del
yacimiento y del sistema total de tuberías que está instalado.
¿Por qué escogemos el separador como nodo solución? Si tomamos al separador como nodo
solución es muy fácil visualizar los efectos de la presión de separación (Psep) sobre el gasto de
producción. Veamos la figura 3.7 la cual muestra cuatro diferentes sistemas de pozos.
89
Capítulo III
Es obvio que el pozo A muestra un incremento significante en la producción cuando la presión del
separador es disminuida, pero el pozo D no muestra un cambio esencial en la producción con una
disminución de la presión de separación Psep.
Figura 3.7 Efecto de la presión de separación sobre los gastos de producción
Cada pozo debe ser analizado individualmente para optimizar correctamente la producción. En
todos los casos el criterio final para la selección de la presión de separación es el económico. Una
precaución que siempre debe tomarse en cuenta es que la línea de descarga (L.D) debe ser
siempre analizada previo a la decisión final de la presión de separación. Hay casos de campo en
donde un cambio en la línea de descarga (L.D) muestra un incremento mayor en la producción
que el cambio en la presión de separación (Psep).
IV.3 Análisis nodal para pozos con restricciones
En los puntos vistos anteriormente en este capítulo se tomaba en cuenta que no existe ninguna
restricción en ninguna de las tuberías (T.P y L.D). Muchos pozos son equipados con
estranguladores superficiales y muchos pozos costa afuera son equipados con válvulas de
tormenta. Un estrangulador usualmente será diseñado de tal forma que el flujo a través del
estrangulador sea crítico, mientras que el flujo a través de la válvula de tormenta será subcrítico.
IV.3.1 Estranguladores
Muchos pozos fluyentes y algunos con sistemas artificiales de producción son equipados con
estranguladores para controlar la producción o la presión después de éste. Estos estranguladores
son usualmente colocados en la cabeza del pozo, pero en algunos casos son colocados cerca del
90
Capítulo III
separador, la localización de los estranguladores puede tener un efecto considerable sobre la

capacidad productiva de los pozos, especialmente si los pozos tiene una línea de descarga (L.D)
muy extensa.
Pnodo = Pwh
Inflow
Pr − ΔPres − ΔPT.P =Pnodo
Outflow
a) Para un estrangulador localizado en la cabeza del pozo en el cual exista flujo crítico se tienen
las siguientes relaciones para encontrar la presión en el nodo.
1) Para una relación de presiones Pd / Pe = 0.5
2 (Psep + ΔPL.D ) = Pnodo
b) Para un estrangulador localizado en el separador, la presión a la entrada del estrangulador, es

la presión al final de la línea de descarga L.D por lo tanto se obtienen las siguientes expresiones:
1) Para una relación de presiones Pd / Pe = 0.5
2 Psep + ΔPL.D = Pndo
de la tubería de producción (T.P) hasta la cabeza del pozo, utilizando correlaciones de flujo
multifásico para tubería vertical o curvas de gradiente de presión.
el gasto), asumir varios gastos de producción y dependiendo de la localización del estrangulador y
del tipo de fluido que se transporta (flujo de gas o flujo multifásico) se escoge alguna de las
ecuaciones presentadas anteriormente y se obtiene la presión en la cabeza del pozo (Pwh = Pnodo)
necesaria para mover los fluidos a través del estrangulador y de la línea de descarga (L.D) hasta
91
Capítulo III
el separador, utilizando una correlación de flujo multifásico en tubería horizontal o curvas de

gradiente de presión.
4) Graficar los gastos y las presiones en la cabeza del pozo (Pwh = Pnodo) obtenidas en el paso 2
(curva de inflow). Sobre esa misma gráfica, graficar los gastos y las presiones en la cabeza del
pozo (Pwh = Pnodo) obtenidas en el paso 3 (curva de outflow). La intersección de las curvas
¿Por qué escogemos el estrangulador como nodo solución? Aunque estrictamente el nodo
solución no es el estrangulador si no la cabeza del pozo, este procedimiento nos permite calcular
el diámetro óptimo y la mejor localización del estrangulador para obtener la mayor producción,
asegurando que a través del estrangulador exista flujo crítico.
Para todos los nodos vistos anteriormente se tiene la siguiente nomenclatura:
Pnodo =Presión en el nodo
Pwf =Presión de fondo fluyendo
Pwh =Presión en la cabeza del pozo
Psep =Presión del separador
ΔPres =Caída de presión en el yacimiento
ΔPL.D =Caída de presión en la línea de descarga
ΔPT.P =Caída de presión en la tubería de producción
IV.3.2 Válvulas de tormenta
El análisis del efecto de las válvulas de tormenta en la tubería de producción (T.P) puede ser el
mismo que para los estranguladores, sin embargo las válvulas de tormenta operan en flujo
subcrítico y por consiguiente las caídas de presión a través de las válvulas deben ser calculadas.
SI el nodo 4 (ver figura 3.2) es escogido, la curva de outflow incluye la sección de la tubería de
producción por arriba de la válvula de tormenta. Es conveniente escoger la presión en el nodo
como la presión inmediatamente arriba de la válvula de tormenta, debido a que las ecuaciones
para calcular la caída de presión a través de la válvula de tormenta dependen sobre las
condiciones de presión y temperatura a la entra de la válvula.
Pnodo = PV.T arriba
92
Capítulo III
Inflow.
Pr − ΔPres − ΔPT.P debajo − ΔPV.T = Pnodo
Outflow
Psep + ΔPL.D + ΔPT.P arriba = Pnodo
Donde:
ΔPV .T =Caída de presión a través de la válvula de tormenta
PV.T arriba =Presión corriente arriba de la válvula de tormenta
ΔPT.P debajo =Caída de presión en la tubería de producción desde el fondo del pozo hasta la
válvula de tormenta
ΔPT.P arriba =Caída de presión en la tubería de producción de la válvula de tormenta hasta la
cabeza del pozo
determinan las caídas de presión en la tubería de producción (T.P) hasta la profundidad donde se
encuentra localizada la válvula de tormenta (PT.P debajo), utilizando una apropiada correlación de
flujo multifásico en tubería vertical o curvas de gradiente de presión.
3) Utilizando los mismos gastos y las presiones a la entrada de la válvula de tormenta (PT.P debajo)
obtenidas en el paso 2, se calculan las caídas de presión a través de ésta (∆PV.T), obteniendo la
presión en el nodo (PV.T arriba).
el gasto), asumir varios gastos de producción y calcular la presión en la cabeza del pozo (Pwh)
utilizando una apropiada correlación de flujo multifásico en tubería horizontal o curvas de gradiente
de presión.
4, se determinan las presiones a la profundidad donde se encuentra localizada la válvula de
tormenta (PV.T arriba) necesarias para transportar los fluidos desde este punto hasta la cabeza del
93
Capítulo III
pozo, utilizando una apropiada correlación de flujo multifásico en tubería vertical o curvas de
6) Graficar los gastos y las presiones en el nodo (PV.T arriba) obtenidas en el paso 3 (curva de
inflow). Sobre esa misma gráfica, graficar los gastos y las presiones en el nodo (PV.T arriba)
obtenidas en el paso 5 (curva de outflow). La intersección de las curvas mostrará el gasto posible
para este sistema, este gasto no es el gasto máximo u óptimo del sistema, pero es el gasto al cual
el pozo produciría para las condiciones dadas del yacimiento, del sistema total de tuberías que
está instalado y del diámetro de la válvula de tormenta utilizada.
¿Por qué escogemos la válvula de tormenta como nodo solución? Usando este análisis, la
curva de outflow no cambiara para diferentes diámetros de la válvula de tormenta y únicamente la
caída de presión a través de la válvula de tormenta cambiara en el cálculo de la curva de inflow.
Un análisis para este sistema es mostrado en la figura 3.8 el cual resultaría para diferentes
tamaños de diámetro de la válvula de tormenta, la figura también muestra el caso cuando no se
tiene válvula de tormenta.
Figura 3.8 Efecto de la válvula de tormenta sobre los gastos de producción
IV.4 Análisis nodal para tuberías telescopiadas
Existen numerosas razones por las cuales las tuberías telescopiadas de producción son instaladas
en un pozo, pero la principal involucra el uso de liners (tuberías de revestimiento cortas) en la
terminación de los pozos. Numerosas terminaciones, especialmente en aguas profundas, utilizan
liners en la sección más baja del pozo.
94
Capítulo III
El liner restringe el diámetro de la tubería de producción que puede ser colocada, pero la tubería
de revestimiento larga (T.R) permite que una tubería de producción de mayor diámetro pueda ser
colocada desde la parte superior de liner hasta la superficie, ver figura 3.9. El análisis nodal a este
tipo de sistemas de producción, permite una fácil solución para determinar el efecto de los
diferentes diámetros de la tubería de producción por arriba del liner.
Pnodo = Ptop liner
Inflow.
Pr − ΔPres − ΔPT.P 1 = Pnodo
Outflow
Psep + ΔPL.D + ΔPT.P 2 = Pnodo
Donde:
Ptop liner =Presión en la parte alta del liner (donde se juntan las dos tuberías de producción)
ΔPT.P 1 =Caída de presión en la T.P que corre del fondo del pozo hasta la parte alta del liner
ΔPT.P 2 =Caída de presión en la T.P que corre de la parte alta del liner hasta la cabeza del pozo
Figura 3.9 Esquema de un SIP para un pozo con tuberías telescopiadas
95
Capítulo III
determinan las presiones en la parte alta del liner (donde se juntan las dos tuberías de producción
de diferente diámetro), utilizando una correlación de flujo multifásico en tubería vertical o curvas de
presión.
3, se determinan las presiones en la parte alta del liner (donde se juntan las dos tuberías de
producción de diferente diámetro), utilizando una correlación de flujo multifásico en tubería vertical
o curvas de gradiente de presión.
5) Graficar los gastos y las presiones en la parte alta de liner (Ptop liner = Pnodo) obtenidas en el paso
2 (curva de inflow). Sobre esa misma gráfica, graficar los gastos y las presiones en la parte alta de
liner (Ptop liner = Pnodo) obtenidas en el paso 4 (curva de outflow). La intersección de las curvas
IV.5 Análisis nodal para sistemas multipozos
Hasta este momento sólo hemos visto como utilizar el análisis nodal para un pozo, pero el mismo
procedimiento también puede ser aplicado para el análisis de sistemas con múltiples pozos o
campos enteros. El procedimiento será ilustrado cualitativamente y nos enfocaremos en el sistema
mostrado en la figura 3.10. En este caso un cambio hecho en cualquier componente, afectara la
capacidad total de producción del sistema. Algunos de los cambios que pueden ser considerados
son:
• Colocar algún sistema artificial de producción en algunos pozos.

• Adicionar nuevos pozos al sistema.
• Cerrar algunos pozos existentes.
• Cambios en las características de producción con el tiempo.
96
Capítulo III
• Modificaciones en las líneas de descarga (L.D)

• Instalación de bombas o compresoras.
• Cambio en la presión de salida del sistema (PD).
Figura 3.10 Sistema multipozos
La elección del nodo para el análisis final debe ser aquel punto en el cual no exista mezclado de
flujo sobre la corriente después del nodo, estos nodos son C y D en nuestro sistema
(ver figura 3.10). Como se observa entre los puntos C y D sólo existe una sola corriente de flujo.
El análisis siempre debe empezar con la presión media del yacimiento (Pr) que es independiente
del gasto y terminar con alguna presión de salida (PD) que sea también independiente del gasto.
Para nuestro sistema escogeremos el punto C como el nodo solución.
Empezaremos el análisis con los puntos A y B. La curva de inflow para el punto A será calculada
de la siguiente manera:
Pr − ΔPres − ΔPT.P − ΔPL.D = PA
Esta expresión será evaluada para cada uno de los pozos que alimenta el nodo A (pozos 1, 2 y 3),
para un rango supuesto de gastos de producción. Esto resultará en una gráfica como la mostrada
en la figura 3.11.
Una gráfica similar para el comportamiento de presión en el nodo B puede ser construida
considerando los pozos (4 y 5), esto es ilustrado en la figura 3.12. La relación gas líquido RGL y
fracción de agua Fw usada en los cálculos de las caídas de presión en el sistema de tuberías para
97
Capítulo III
estos puntos (A y B) serán aquellas que correspondan a cada uno de los pozos debido a que no
existe mezclado de corrientes de flujo.
El siguiente paso es calcular la presión en el punto C, para este punto ya existe el mezclado de
corrientes de flujo provenientes de los pozos 1, 2 y 3 para el nodo A y de los pozos 4 y 5 para el
nodo B.
Figura 3.11 Curva de Inflow en el nodo A
Figura 3.12 Curva de Inflow en el nodo B
98
Capítulo III
Las expresiones de la curva de inflow para el punto C que proviene de las corrientes de flujo de
los puntos A y B son:
PA − ΔPAC = PC
PB − ΔPBC = PC
Esto resultará en una relación entre la presión en el punto C y el gasto en el mismo punto
(qC =qA + qB), como se ilustra en la figura 3.13.
Figura 3.13 Curva de Inflow en el nodo C
Los cálculos sobre la caída de presión entre los puntos A y C (∆PAC) y entre los puntos B y C
(∆PBC) es diferente, debido al hecho de que la relación Gas-Líquido (RGL) y la Fracción de agua
(Fw), tienen diferentes valores para cada uno de los pozos del sistema. En este caso, la correcta
RGL y Fw para obtener ∆pAC, se calcula de la siguiente manera:
RGL AC =
∑ q RGL
i i
qA
FwAC =
∑q F
i wi
qA
Expresiones similares son usadas para determinar estos valores a la hora de calcular ∆PBC. La
expresión para calcular la curva de outflow en el punto C es la siguiente:
99
Capítulo III
PD + ΔPCD = PC
Cálculos de ∆PCD para varios gastos requerirán de nuevo determinar la correcta relación
Gas-Líquido (RGL) y Fracción de agua (Fw), correspondiente a cada gasto en el punto C
(qC =qA + qB). Un cambio en la presión de salida (PD) o en la línea de flujo entre los puntos C y D
resultará en diferentes curvas de outflow y además diferentes capacidades de sistema como es
mostrado en la figura 3.14
Figura 3.14 Capacidad del sistema
Esta es la forma en que se aplica el análisis nodal a un sistema con varios pozos, lo más
importante para estos casos es tomar como nodo solución aquel nodo donde no exista mezclado
de corrientes de flujo provenientes de otros puntos y que además las presiones al inicio y al final
del sistema Pr y PD (en nuestro ejemplo) son valores constantes y no dependen del gasto.
100

Teoria - Ejemplo de Aplicación (Pozo Guaricho 12) PDF

Cargado por

Copyright:

Formatos disponibles

Teoria - Ejemplo de Aplicación (Pozo Guaricho 12) PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Teoria - Ejemplo de Aplicación (Pozo Guaricho 12) PDF

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO

“Aplicación del Análisis Nodal para

DIRECTOR DE TESIS: M. en I. Tomás Eduardo Pérez García.

MÉXICO, D.F., CIUDAD UNIVERSITARIA, “2008”

I COMPORTAMIENTO DE FLUJO EN EL YACIMIENTO

I.1 Flujo de aceite

I.3 Flujo de dos fases: Aceite y Gas

I.5 Curvas de IPR para pozos de gas

I.6 IPR de dos fases para pozos direccionales ----------------------------------------- 30

I.9 Curvas de IPR futuras

II.1 Fundamentos de flujo a través de tuberías

II.2 Flujo de una sola fase

II.3 Flujo multifásico en tuberías

II.5 Flujo multifásico en tuberías horizontales

II.6 Flujo por espacio anular -------------------------------------------------------------------- 78

II.7 Flujo a través de restricciones

III BASES PARA EL ANÁLISIS NODAL

III.1 ¿Qué es el análisis nodal? ------------------------------------------------------------------ 81

III.2 Elección del nodo solución

III.3 Análisis nodal para pozos con restricciones

III.4 Análisis nodal para tuberías telescopiadas ----------------------------------------- 94

III.5 Análisis nodal para sistemas multipozos --------------------------------------------- 96

CONCLUSIONES ---------------------------------------------------------------------------------------------- 113

NOMENCLATURA ------------------------------------------------------------------------------------------- 115

BIBLIOGRAFÍA ------------------------------------------------------------------------------------------------ 118

• Cambio del diámetro del estrangulador

• Producir por espacio anular

Del estudio del análisis nodal se concluyó lo siguiente:

La razón fundamental de someter un sistema integral de producción (SIP) a la técnica de análisis

1) Determinar el gasto al cual un pozo productor de aceite o gas, producirá considerando la

COMPORTAMIENTO DE FLUJO EN EL YACIMIENTO

El yacimiento es el espacio físico en donde se encuentran almacenados los fluidos petroleros

Pr =Presión media del yacimiento

Pwf =Presión de fondo fluyendo

En la actualidad se cuenta con buenos simuladores de yacimientos que pueden pronosticar de

I.1 Flujo de aceite

I.1.1 Ley de Darcy

q =Gasto del fluido (cm3 / s)

dP =Diferencial de Presión (atm)

I.1.2 Ecuaciones de flujo

7.08 × 10 −3 (k oh)(Pr − Pwf )

7.08 × 10 −3 (k oh)(Pr − Pwf )

qo =Gasto de aceite (bpd @ c.s)

k o =Permeabilidad efectiva al aceite (md)

h =Espesor de la zona productora (pie)

rw =Radio del pozo (pie)

I.2 ÍNDICE DE PRODUCTIVIDAD Y EL IPR

I.2.1 Definición del índice de productividad (J)

Un parámetro del yacimiento comúnmente usado es el índice de productividad (J), el cual

Figura 1.1 J vs. Np para diferentes mecanismos de producción

Tipos de mecanismos naturales de producción del yacimiento.

Expansión del sistema roca-fluidos.

Expansión por gas disuelto.

Empuje de acuífero activo.

Expansión del casquete de gas.

En el diseño de sistemas artificiales de producción, el índice de productividad es expresado en bpd

qo =Gasto de aceite (bpd @ c.s)

I.2.2 Definición de IPR

I.2.3 Índice de productividad para yacimientos bajo saturados