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:

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 cual nos indica
el comportamiento del índice de productividad (J) con respecto a la producción acumulada para
diferentes mecanismos de producción.
FIGURA Nº1 MECANISMOS DE PRODUCCION

Fuente:(http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1)
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:

Dónde:
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.

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).
FIGURA Nº2 CURVA IPR

Fuente(http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1)

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:

I.2.4 Índice de productividad para yacimientos saturado

En el caso de que Pwf < Pb, la curva de IPR se comporta

FIGURA Nº3 CURVA DE IPR PARA YACIMIENTOS SATURADOS

Fuente:http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1

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 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
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:
El término es una función de la presión y a su vez el valor de la permeabilidad relativa al
aceite (rok ) es una función de la saturación de aceite ( o S ). Una gráfica típica de esto se
muestra a continuación
FIGURA Nº 4 PRESION

Fuente:http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1
La integral puede ser en dividida en dos partes, como se muestra a continuación:

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 del Pb.

FIGURA Nº5 Curva de IPR para yacimientos con empuje de gas disuelto (Vogel)

Fuente:http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1

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:
  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.
 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.
 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:

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.
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:

FIGURA Nº5 SOLUCION ANALITICA

Fuente:http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1
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.

Dónde:
wf P =Presión de fondo fluyendo si el pozo no tuviera daño
wf P =Presión de fondo fluyendo del pozo
R P =Presión media del yacimiento
Despejando el valor de Pwf

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
FIGURA Nº6 CORRELACION DE POZOS

Fuente;http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1
Solución:

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

FIGURA Nº7 ERROR EN LA EXTRAPOLACION

Fuente;http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1
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:

Dónde:
J =Í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 – Pwf 2)
vs. Qo dará una línea recta en una gráfica Log-Log como se muestra en la figura 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).
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
FIGUTA Nº9 AMPLIFICACION DE TRABAJO

Fuente:http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1

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.
I.3.4 Ecuación de Jones, Blount y Glaze
Jones, Blount y Glaze propusieron la siguiente relación entre presión y gasto.

En un plano cartesiano graficamos lo cual nos dará una línea recta con ordenada al
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.

I.3.5 Ecuación de Klins y Majcher

Klins y Majcher(13) desarrollaron un IPR incorporando el valor de Pb.
I.3.6 Ecuación de Sukarno
Sukarno 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:

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.
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:

FIGURA Nº10 COMPORTAMIENTO DE VOGEL

Fuente:http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1
Sumando las dos ecuaciones obtenemos el gasto máximo para dicho yacimiento.

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:

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:

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.

Dónde:
r P =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. podemos conocer el índice de
productividad J cuando se tiene flujo de dos fases y con este valor recurriendo a la ecuación
general de índice de productividad podemos obtener 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:

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.
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:

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:

Dónde:
ρ =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:

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:
Dónde:
g q =Gasto de Gas (Mpcd)
k =Permeabilidad en (md)
Capítulo I
25
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
e r =Radio de drene (pie)
w r =Radio del pozo (pie)
r P =Presión en el radio de drene en (lb / pg2)
wf P =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:
Dónde:
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.
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
pseudoestacionario.
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.

Dónde:
g q =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:

En la ecuación 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 – Pwf 2) vs. qg, 1 / n es la pendiente de la línea graficada y
el valor de C =qg cuando (Pr 2 – Pwf 2) =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:

Si graficamos los datos en una gráfica Log-Log se obtiene la siguiente gráfica.

FIGURA Nº11 PRUEBAS DE POZO

Fuente:http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1

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:

Utilizando la prueba 4 para calcular el valor de C:

Por lo tanto la ecuación que representa el comportamiento de este pozo es la siguiente:

Dónde:
g q =Gasto de gas (Mpcd)
r P =Presión media del yacimiento (lb / pg2)
wf P =Presión de fondo fluyendo (lb / pg2)
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.

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:

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.
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.
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.

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:

Wang y Wiggins 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:

La segunda IPR fue desarrollada basada en diferentes etapas de recuperación. La ecuación a la

cual los autores llegaron es:

Donde d está en función del factor de recuperación. La siguiente tabla muestra los valores de d
para las diferentes etapas.
Por ultimo desarrollaron una ecuación que permite el cálculo de IPR futuras:

Para utilizar esta ecuación, primero se debe de encontrar el gasto máximo presente con
cualquiera de las dos ecuaciones dadas anteriormente y con este dato, obtener el gasto máximo
futuro de la ecuación Con el gasto máximo futuro regresar a las ecuaciones para poder graficar
la curva de IPR futura.
I.8 IPR para flujo de tres fases (aceite, gas y agua)
M.L Wiggins 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.

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:
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
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.
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 cuándo 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 Nº 12 CURVAS IPR

Fuente;http://www.ptolomeo.unam.mx:8080/xmlui/bitstream/handle/132.248.52.100/8610/Tesis_Completa.pdf?
sequence=1
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:

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:

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.

Dónde:
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.
40
I.9.2 Método de Fetkovich
Fetkovich propuso la siguiente ecuación para el comportamiento de la curva de IPR futura.
( 2 )n
wf
2
r,F
r,P
' r,F
o o,P P P
P
P
Jq−⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
= ---------------1.46
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騁odo 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
r,F
r,P
omax,F
omax,P
P
P
q
q
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
= ---------------1.47
I.9.4 M騁odo 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:
( ) ( ) (EF)(1 P')[1.8 0.8 (EF)(1 P')]
B
Pk
In 0.472 r / r
q 3.49 A kh
oo
ro
r
ew
o−−−⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
μ
= ---------------1.48
Donde:
A =Factor de conversión =0.001127, para (pie, lb / pg2, cp, md y bpd)
P' =
r
wf
P
P
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.
Cap咜ulo I
41
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.
( 2 )n
wf
2
g F r,F q = C P −P ---------------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.
Debido a que
⎥ ⎥⎦
⎤
⎢ ⎢⎣
⎡
+⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
⋅μ⋅
×
=
−
S
r
T Z In 0.472r
C 703 10 kh
w
e
g
6
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:
()
( )g 1
g2
2
1
Z
Z
C
C
μ
μ
= , por lo tanto para nuestro caso tenemos que:
()
( )g F
gp
FpZ
Z
CC
μ
μ
= ---------------1.50
Donde:
g q =Gasto futuro de gas
P C =Coeficiente al tiempo presente consistente con las unidades de gasto y presión
F C =Coeficiente al tiempo futuro consistente con las unidades de gasto y presión
r P =Presión media del yacimiento a un tiempo dado (lb / pg2)
wf P =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.
Capítulo II
42
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.
Capítulo II
43
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.
E1 =Energía por unidad de masa en la posición uno ⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛−
m
f
lb
lb pie
E2 =Energía por unidad de masa en la posición dos ⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛−
m
f
lb
lb pie
Energía de expansión
La energía de expansión está dada por la siguiente expresión:
Capítulo II
44
PV
lb
V pie
pie
P lb
lb
E lb pie
m
3
2
f
m
f
e=⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
=⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛−
-----------------2.2
Donde:
V =Volumen especifico (1 / ρ).
Energía potencial
La energía potencial está dada por la siguiente expresión:
()h
g
h pie g
lb pie
lb s
g
1
s
g pie
lb
E lb pie
mc
2
f
c
2
m
f
p=⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
−
−
⎟⎠
⎞
⎜⎝
=⎛⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛−
-----------------2.3
Energía cinética
La energía cinética está dada por la siguiente expresión:
c
2
m
2
f
c
2
22
m
f
c 2g
v
lb pie
lb s
g
1
s
pie
2
v
lb
E lb pie = ⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
−
−
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
=⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛−
-----------------2.4
Al sustituir las energías correspondientes a las posiciones 1 y 2 en la ecuación 2.1 se obtiene:
c
22
2
c
fS22
c
2
1
1
c
1 1 2g
hv
g
WWPVg
2g
hv
g
PV + g + + Δ + Δ = + + -----------------2.5
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:
L
W
2g L
v
gL
gh
L
Pf
c
2
cΔ
Δ
+ρ
Δ
Δ
+ρ
Δ
Δ
=ρ
Δ
Δ
A esta ecuación se le acostumbra escribir de la siguiente forma:
T E AC F L
P
L
P
L
P
L
P ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
+ ⎟⎠
⎞
⎜⎝⎛
Δ
Δ
+ ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
= ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
-----------------2.6
Capítulo II
45
Donde:
ΔL T
ΔP⎟⎠
⎞
⎜⎝
⎛ =Gradiente de presión total
EL
P ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
=Gradiente de presión debido a la elevación
AC L
P ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
=Gradiente de presión debido a la aceleración
FL
P ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
=Gradiente de presión debido a la fricción
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:
2g d
fv
L
P
c
2
F
ρ
= ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
-----------------2.7
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:
μ
ρ
Nre = dv -----------------2.8
Donde:
Capítulo II
46
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:
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:
Nre
f = 64 -----------------2.9
Para flujo turbulento (Nre > 3100), (f) está dado por la ecuación de Colebrook y White:
2
f Nre
2.514
3.715 d
f 2log
−
⎥ ⎥⎦
⎤
⎢ ⎢⎣
⎡
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
+
ε
= − ---------------2.10
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:
Tipo de Tubería ε (pg)
Tubería estriada 0.00006
Tubería de producción o
perforación 0.0006
Tubería de escurrimiento 0.0007
Tubería galvanizada 0.006
Capítulo II
47
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:
( ( )) 2 f 2log / 3.715 d − = − ε ---------------2.11
d) Cuando el flujo es crítico (2300< Nre <3100), (f) se puede aproximar con la siguiente
expresión:
0.032
3100 f
2.514
3.715 d
2.3026 log
1.3521
2300
f Nre 2300 2 +
⎥⎥⎥⎥⎥
⎦
⎤
⎢⎢⎢⎢⎢
⎣
⎡
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
+
ε
×
−
= ---------------2.12
Capítulo II
48
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:
T E F ΔP = ΔP + ΔP
La caída de presión por elevación se calcula de la siguiente manera:
P 0.433 h E L Δ = λ Δ ---------------2.13
Donde:
E ΔP =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:
5
2
L
Fd
P 0.06056 f q L
⋅λ⋅⋅
Δ = ---------------2.14
Donde:
q =Gasto de líquido (bpd)
L =Longitud de la tubería (millas)
d =Diámetro de la tubería (pg)
F ΔP =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:
Capítulo II
49
5
2
L
TLd
P 0.433 h 0.06056 f q L
⋅λ⋅⋅
Δ = λ Δ + ---------------2.15
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:
μ
λ
=
d
Nre 92.2 q L ---------------2.16
Donde:
q =Gasto de líquido (bpd)
d =Diámetro de la tubería (pg)
μ =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:
()
dP
q Z T 460 L
T 460
P
519.98
Pf5
g
2
0
0
F⋅
⋅λ⋅⋅+⋅
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
+
⎟⎠
⎞
⎜⎝
Δ=⎛
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 mediaP 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:
()
5
g
2
0
0
2
2
2
2
1d
q Z T 460 L
T 460
P
461.346
PPf
⋅λ⋅⋅+⋅
⎟ ⎟⎠ ⎞
⎜ ⎜⎝
⎛
+
⎟⎠
⎞
⎜⎝
− = ⎛ ---------------2.17
Donde:
q =Gasto de gas (pcd)
L =Longitud de la tubería (millas)
d =Diámetro de la tubería (pg)
Capítulo II
50
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
1 2 P y P =Presiones al inicio y al final P1 > P2 (lb / pg2)
0 0 P y T =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:
PhEgΔ=ρ
Por lo que las pérdidas por elevación en unidades prácticas queda de la siguiente forma:
Z (T 460)
Ph
144
P 2.7044 g
E⋅+
⋅λ⋅
⎟⎠
⎞
⎜⎝
Δ=⎛
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:
Z (T 460)
Ph
P P 0.03756 g
2
2
´2
2
1⋅+
⋅λ⋅
− = ---------------2.18
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 unaP 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.
Capítulo II
51
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.
()
Z (T 460)
Ph
0.03756
d
q Z T 460 L
T 460
P
461.346
PPfg
2
5
g
2
0
0
2
2
2
2
1⋅+
⋅λ⋅
+
⋅λ⋅⋅+
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
+
⎟⎠
⎞
⎜⎝
− = ⎛ ---------------2.19
Como se puede observar en la ecuación anterior, se debe suponer un valor de P2 para poder
obtener el valor deP 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.
()
2.5
0.5
g
2
2
2
2
g 0.5 1 d
Z T 460
0.03756 P h
PP
L
qC
⎥ ⎥⎦
⎤
⎢ ⎢⎣
⎡
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅+
⋅λ⋅
= − − ---------------2.20
Donde:
()
0.5
0g
0
2
Z T 460
1
P
T 460
f
C 461.346
⎥ ⎥⎦
⎤
⎢ ⎢⎣
⎡
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+⋅⋅λ⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛+
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
= ---------------2.21
O bien el diámetro dadas unas condiciones de presión y gasto:
( ( ))
0.2
2g
g
22
2
2
1
q
P P 0.03756 P h / Z T 460
1
C
dL
⎥ ⎥⎦
⎤
⎢ ⎢⎣
⎡
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
−−⋅λ⋅⋅+
⎟⎠
⎞
⎜⎝
= ⎛ ---------------2.22
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:
g
gg
d
q
Nre 0.0201056
μ
λ
= ---------------2.23
Donde:
g q =Gasto de gas (pcd)
Capítulo II
52
g λ =Densidad relativa del gas
d =Diámetro de la tubería (pg)
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
2
2
1
3
2
3
1
PP
PP
3
P 2 ---------------2.24
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:
Linea
Contenido
de liquido
gal /
MMpie3
EF
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
Capítulo II
53
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.
Capítulo II
54
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.
( 1.175 )
o B = 0.972 + 0.000147 F
Donde:
F R 1.25 T
o
gd
s + ⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
γ
γ
=
Capítulo II
55
Vazquez.
o 1 s ( )( gs )( 2 3 s ) B = 1+ C R + T − 60 API / γ C + C R
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:
(1 5.912 10 (API)(T )log (P / 114.7)) s s
5
gs gp
γ=γ+×−
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:
()
P
B 0.02825 Z T 460 g
⋅+
=
Para estos casos se tiene la siguiente nomenclatura:
γ = gd Densidad relativa del gas disuelto
γ = o Densidad relativa del aceite
= o B Factor de volumen del aceite
= g B Factor de volumen del gas
= s R 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
0.25 0.02 API R 10 6 (0.6874 3.5864 API)
gd s γ = + + × − −
Capítulo II
56
Densidad del aceite
Saturado.
o
o s gd
oB
62.4 γ + 0.01362R ⋅ γ
ρ=
Bajo saturado.
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛−
ρ=ρ⋅b
Co P P
e o ob
Donde:
o C =Compresibilidad del aceite
b P =Presión de burbuja
o ρ =Densidad del aceite
ob ρ =Densidad del aceite a la presión de burbuja
Densidad del gas
Z (T 460)
2.7044P gf
g⋅+
⋅γ
ρ=
Donde:
gf γ =Densidad relativa del gas libre
g ρ =Densidad del gas
Densidad del agua saturada
w
wB
ρ = 62.43
Donde:
w ρ =Densidad del agua
w B =Factor del volumen del agua
Capítulo II
57
Viscosidad del aceite
Saturado
b
o om μ = a μ
( ) 0.515
s a 10.715 R 100 − = +
( ) 0.388
s b 5.44 R 150 − = +
10D 1
om μ = −
D = E (T −1.163 )
E = 10F
F = 3.0324 − 0.02023 API
Bajo saturado
m
b
o ob P
P
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
μ=μ
⎥⎥
⎦
⎤
⎢⎢
⎣
⎡
⎟⎠ ⎞
⎜⎝ ⎛
=⋅
−11.513+ − 8.98 × 10− 5 P
m 2.6P1.187 e
Donde:
o μ =Viscosidad del aceite
Relación de solubilidad
Standing
( ) 0.83
1
10 0.0125 API 0.00091T
18
R P gd s ⎥⎦ ⎤
⎢⎣
⎡
⎟⎠
⎞
⎜⎝
=γ⎛−
Vazquez
R C P exp [C API / (T 460)] 3
C
s 1 gs
=γ⋅2⋅+
Donde:
[1 5.912 10 (API) (T )log (P / 114.7)] s s
5
gs g
γ=γ+×−
Capítulo II
58
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
Capítulo II
59
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:
Capítulo II
60
()
⎥ ⎥⎦
⎤
⎢ ⎢⎣
⎡
=+++
6
5
C
Lv
C
2 gv
4 Lμ
2
L123N
N
Y exp C C senθ C sen θ C N ---------------2.25
Donde:
0.25
3
L
L
2
L
N 0.15726 1 ⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
ρσ
= μ μ ---------------2.26
0.25
L
Lv sL v 938 . 1 N ⎟⎠
⎞
⎜⎝
⎛
σ
ρ
= ---------------2.27
0.25
L
sg gv v 938 . 1 N ⎟⎠
⎞
⎜⎝
⎛
σ
ρ
= ---------------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
Flujo Tipo de Flujo C1 C2 C3 C4 C5 C6
Horizontal o
ascendente Todos -0.38011 0.12988 -0.11979 2.34323 0.47569 0.28866
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):
()
()ooww
Lgosg
L
5.615 q B q B
qRRB
1
1
qq
q
+
⋅−⋅
+
=
+
λ = ---------------2.29
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:
()
2
ooww
p
L
sL d
0.01191 q B q B
A
vq
+
= = ---------------2.30
Capítulo II
61
()
2
osg
p
g
sg d
0.002122 q R R B
A
q
v
⋅−⋅
= = ---------------2.31
sL sg
p
Lg
mvv
A
qq
v=+
+
= ---------------2.32
Donde:
Ap es el área de la sección transversal de la tubería. De estas ecuaciones se observa que:
m
sL
v
λ = v ---------------2.33
II.3.5 Velocidad real
Aplicando el concepto de colgamiento, se puede obtener la velocidad real correspondiente a cada
fase:
L
sL
pL
L
L
L
LY
v
AY
q
A
v = q = = ---------------2.34
()()L
sg
pL
g
g
g
g1Y
v
A1Y
q
A
q
v
−
=
−
= = ---------------2.35
Donde:
L g q y q =Gastos de líquido y gas
L g A y A =Áreas por donde pasa el líquido y el gas respectivamente
p A =Á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 L g L ρ = ρ Y + ρ 1− Y ---------------2.36
Algunos autores calculan la densidad de la mezcla sin considerar el resbalamiento entre las fases,
esto es:
Capítulo II
62
ρ = ρ λ + ρ (1− λ) ns L g ---------------2.37
También la densidad puede obtenerse a partir de la siguiente expresión:
m
ns V
ρ = M ---------------2.38
Donde:
M =Masa de la mezcla @ c.s por barril de aceite producido @ c.s
m V =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:
o g w M = M + M + M ---------------2.39
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
=λ
o
o
o o bls
M 350.5 lb ---------------2.40
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
=λ⋅
o
g
g g bls
lb
M 0.0764 R ---------------2.41
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
=λ⋅
o
w
w w w bls
M 350.5 F lb ---------------2.42
Sustituyendo las ecuaciones 2.40, 2.41 y 2.42 en la ecuación 2.39 se obtiene:
( ) ⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
=λ+λ+⋅λ
o
m
o w w g bls
M 350.5 F 0.0764R lb ---------------2.43
Cálculo de Vm (pie3
o+g+w @ c.s / blso @ c.s)
mo o V = 5.615B ---------------2.44
( ) mg s g V = R −R B ---------------2.45
mw w w V = F ×5.615B ---------------2.46
( ) ( ) m o w w s g V = 5.615 B +B F + R −R B ---------------2.47
Sustituyendo las ecuaciones 2.43 y 2.47 en la ecuación 2.38 se obtiene:
Capítulo II
63
()
()()⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
++−
λ+λ+⋅λ
ρ=3
m
owwsg
owwg
ns pie
lb
5.615 B B F R R B
350.5 F 0.0764R
---------------2.48
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:
86400
w qo M
m = ---------------2.49
m o g w w = w + w + w ---------------2.50
Donde:
w q B / 15388 o o o o = ρ
w q B / 15388 w w w w = ρ
w q (R R )B / 86400 g g o s g = ρ −
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
( ) L 1 YL
g
Y
mL
− μ = μ μ ---------------2.52
Donde:
ns μ =Viscosidad de la mezcla sin resbalamiento
m μ =Viscosidad de la mezcla con resbalamiento
Además:
L o o w w μ = μ F + μ F ---------------2.53
ow
o
wwoo
oo
o B WOR B
B
qBqB
FqB
+
=
+
= ---------------2.54
w o F = 1−F ---------------2.55
Capítulo II
64
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:
()
⎥ ⎥⎦
⎤
⎢ ⎢⎣
⎡
×ρ⋅
=ρ+5
ns
10
2
tp o
ns 7.413 10 d
fqM
144
1
dL
dP ---------------2.56
Donde:
o q =Gasto de aceite (bpd)
M =Masa de la mezcla por barril de aceite producido (lbm @ c.s / blso @ c.s)
Capítulo II
65
d =Diámetro de la tubería (pg)
ns ρ =Densidad de la mezcla sin resbalamiento (lbm @ c.s / pie3
m @ c.s)
Los valores de M y ρns se obtienen de la siguiente manera:
( o w w ) g M = 350.5 γ + γ F + 0.0764R⋅ γ
ρ = ρ λ + ρ (1− λ) ns L g
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:
d
dv 1.47 10 5 qo M
ns
ρ = × − ---------------2.57
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:
344263
tp f = 5.415 ×10− − 5.723 ×10− a +1.848 ×10− a + 3.5843 ×10− a ---------------2.58
Capítulo II
66
Donde:
qM
a d 10
o
×6
= ----------------2.59
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:
⎟⎠
⎞
⎜⎝
⎛Δ
=−
2
PPP1YPPP21=−Δ
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.
calculada
sup uesta
L
P
P
L
Δ
Δ
Δ
Δ=
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:
Capítulo II
67
⎥⎥⎥⎥⎥
⎦
⎤
⎢⎢⎢⎢⎢
⎣
⎡
⋅
⋅⋅ρ
−
⋅ρ ⋅
+
⋅ θ ⋅ρ
=
144 g P
vv
1
5.362 d
fv
g
g sen
144
1
dL
dP
c
m sg m
2
tp ns m
c
m
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:
⎥ ⎥⎦
⎤
⎢ ⎢⎣
⎡ ⋅ρ ⋅
+
⋅ θ ⋅ρ
=
5.362 d
fv
g
g sen
144
1
dL
dP 2
tp ns m
c
m ----------------2.60
El factor de fricción para las dos fases, se obtiene de la siguiente ecuación:
n
n
tp
tp f
f
f
f⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
= ----------------2.61
Donde:
S
n
tp e
f
f
=⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
----------------2.62
( ) ( )2 ( )4 0.0523 3.182 ln x 0.8725 ln x 0.01853 ln x
S ln x
−+−+
= ----------------2.63
Y (0)
x
L
λ
= ----------------2.64
()
2
n 4.5223 log Nre 3.8215
f 2log Nre
−
⎥⎦
⎤
⎢⎣
⎡
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
−
= ----------------2.65
Además:
ns
m ns Nre 124 d v
μ
⋅ ⋅ρ
= ----------------2.66
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
Capítulo II
68
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:
25
ns
2
m
FR d
N 7734.9 w
ρ
= ----------------2.68
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:
0.302
1 L = 316 λ
2.4684
2 L = 0.0009252 λ−
1.4516
3 L = 0.10 λ−
6.738
4 L = 0.5 λ−
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.
Y = Y (0)ψ L L ----------------2.69
()c
FR
b
LN
Y0a
λ
= ----------------2.70
ψ = 1+ C[sin(1.8 θ) − 0.333 sin3(1.8 θ)] ----------------2.71
()(g)
FR
f
LV
C = 1− λ In d λe N N ----------------2.72
Donde:
Y (0) L =Colgamiento para flujo horizontal (θ=0)
FR N =Número de Froude
λ =Colgamiento sin resbalamiento
θ =Angulo de inclinación de la tubería con respecto a la horizontal
Capítulo II
69
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
(Cuesta Arriba) - - - 0.011 -3.768 3.539 -1.614
Intermitente
(Cuesta Arriba) - - - 2.96 0.305 -0.4473 0.0978
Distribuido
(Cuesta Arriba) - - - C = 0 ψ = 1
Todos los Patrones
(Cuesta Abajo) - - - 4.7 -0.3692 0.1244 -0.5056
Para flujo transitorio, el colgamiento se obtiene de la siguiente forma:
Y (0) L =A×YL (0) (Segregado) + B1×YL (0) (Intermitente)
Donde:
32
3 FR
LL
ALN
−
−
=
B1A1=−
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:
⎟⎠
⎞
⎜⎝
⎛Δ
=−
2
PPP1YPPP21=−Δ
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 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.
Capítulo II
70
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.
12) Obtener el valor del gradiente de presión y determinar el ΔL.
calculada
sup uesta
L
P
P
L
Δ
Δ
Δ
Δ=
13) El valor de L ahora será L =Lanterior + ΔL.
14) Sustituir el valor de P1 = P2 y volver a repetir el procedimiento desde el paso 1 hasta
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)
Capítulo II
71
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
Capítulo II
72
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:
T F AC L
P
L
P
L
P ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
+ ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
= ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
La ecuación anterior se puede escribir como:
2g L
v
2g d
fv
L
P
c
2
mm
c
2
mm
tp
T⋅Δ
ρΔ
+
⋅
ρ
= ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
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:
Capítulo II
73
5
ns
2
tp m
d
174.158 f w
L
P
ρ
⋅
= ⎟⎠
⎞
⎜⎝
⎛
Δ
Δ
----------------2.73
Para obtener el factor de fricción para las dos fases se utiliza la siguiente correlación:
( ) ( )b
L
a
g φ = Nre Nre ----------------2.74
Donde:
a = ψ /(1+ ψ)
b = 1/ e (0.1ψ)
gLψ=w/w
Para:
0 ≤ φ ≤ 500
log f 1.225 0.06561(log ) 0.37 tp = ψ − φ − ----------------2.75
φ ≥ 10,000
log f 0.49 0.12616 (log ) 1.702 tp = ψ − φ − ----------------2.76
500 ≤ φ ≤ 10,000
()2()3
tp 500 log f = F − 0.6561y + 1.1056 +1.7723 ΔF y − 0.46214 + 0.90817 ΔF y ----------------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:
L
L
Ld
Nre 22737 w
μ
= ----------------2.78
g
g
gd
w
Nre 22737
μ
= ----------------2.79
Capítulo II
74
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)
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 de la siguiente forma:
⎟⎠
⎞
⎜⎝
⎛Δ
=−
2
PPP1YPPP21=−Δ
3) Se determinan las propiedades de los fluidos (Rs, Bo, Z, Bg, μo, ρo y ρg etc.) a las
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).
7) Obtener el valor del gradiente de presión y determinar el ΔL.
calculada
sup uesta
L
P
P
L
Δ
Δ
Δ
Δ=
Capítulo II
75
8) El valor de L ahora será L =Lanterior + ΔL.
9) Sustituir el valor de P1 = P2 y volver a repetir el procedimiento desde el paso 1 hasta
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:
5
ns
2
tp m
d
43.539 f w
L
P
ρ
⋅
=
Δ
Δ
----------------2.80
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:
n
n
tp
tp f
f
f
f⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
= ----------------2.81
Donde:
S
n
tp e
f
f
=⎟
⎟⎠
⎞
⎜ ⎜⎝
⎛
----------------2.82
( ) ( )2 ( )4 0.0523 3.182 ln x 0.8725 ln x 0.01853 ln x
S ln x
−+−+
= ----------------2.83
Y (0)
x
L
λ
= ----------------2.84
()
2
n 4.5223 log Nre 3.8215
f 2log Nre
−
⎥⎦ ⎤
⎢⎣
⎡
⎟ ⎟⎠
⎞
⎜ ⎜⎝
⎛
−
= ----------------2.85
En donde:
ns
m ns Nre 124 d v
μ
⋅ ⋅ρ
= ----------------2.86
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:
Capítulo II
76
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:
25
ns
2
m
FR d
N 7734.9 w
ρ
= ----------------2.88
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:
0.302
1 L = 316 λ
2.4684
2 L = 0.0009252 λ−
1.4516
3 L = 0.10 λ−
6.738
4 L = 0.5 λ−
El cálculo del colgamiento real de líquido se obtiene de la siguiente expresión general.
()c
FR
b
LN
Y0a
λ
= ----------------2.89
Donde:
Y (0) L =Colgamiento para flujo horizontal
λ =Colgamiento sin resbalamiento
FR N =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
Capítulo II
77
Para flujo transitorio, el colgamiento se obtiene de la siguiente forma:
YL (0) = A ×YL (0) (Segregado) + B × YL (0) (Intermitente)
En donde:
32
3 FR
LL
ALN
−
−
=
B = 1− A
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 de la siguiente forma:
⎟⎠
⎞
⎜⎝
⎛Δ
=−
2
PPP1YPPP21=−Δ
3) Se determinan las propiedades de los fluidos (Rs, Bo, Z, Bg, μo, ρo y ρg etc.) a las
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).
7) Obtener el valor del gradiente de presión y determinar el ΔL.
calculada
sup uesta
L
P
P
L
Δ
Δ
Δ
Δ=
8) El valor de L ahora será L =Lanterior + ΔL.
9) Sustituir el valor de P1 = P2 y volver a repetir el procedimiento desde el paso 1 hasta
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.
Capítulo II
78
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.
perímetro mojado
r área de la sección transversal h =
Para flujo por espacio anular, tenemos lo siguiente:
()
()4
dd
dd
r d d / 4 T.R T.P
T.R T.P
2
T.P
2
T.R
h
-=-
π+
π
= ----------------2.90
Donde:
T.R d =Diámetro interior de la tubería de revestimiento (T.R).
T.P d =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:
h T.R T.P d = d − d ----------------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.
( ) ( )3
T.R T.P
2
T.R T.P
d5 = d + d d − d ----------------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.
Capítulo II
79
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.
Capítulo II
80
a
c
L
wh d
P b q RGL
⋅
= ----------------2.93
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
Capítulo III
81
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.
Capítulo III
82
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:
Capítulo III
83
Inflow
r ( ) nodo P − ΔP componentes corriente abajo = P
Outflow
( ) sep nodo P + ΔP componentes corriente arriba = P
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:
Capítulo III
84
Las expresiones para las curvas de inflow y outflow son las siguientes:
nodo wf P = P
Inflow
r res nodo P − ΔP = P
Outflow
sep L.D T.P nodo P + ΔP + ΔP = P
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).
Capítulo III
85
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.
Capítulo III
86
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.
Las expresiones para las curvas de inflow y outflow son las siguientes:
nodo wh P = P
Inflow
r res T.P nodo P − ΔP − ΔP = P
Outflow
sep L.D nodo P + ΔP = P
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).
Capítulo III
87
Figura 3.5 Cabeza del pozo como nodo solución.
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 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.
3) 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.
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.
Capítulo III
88
¿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.
Capítulo III
89
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.
Las expresiones para las curvas de inflow y outflow son las siguientes:
nodo sep P = P
Inflow
r res T.P L.D nodo P − ΔP − ΔP − ΔP = P
Outflow
sep nodo P = P
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, 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 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.
Capítulo III
90
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
Capítulo III
91
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.
Las expresiones para las curvas de inflow y outflow son las siguientes:
nodo wh P = P
Inflow
r res T.P nodo P − ΔP − ΔP =P
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
( ) sep L.D nodo 2 P + ΔP = P
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
sep L.D ndo 2P + ΔP = P
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, 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 correlaciones de flujo
multifásico para tubería vertical o curvas de gradiente de presión.
3) 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 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
Capítulo III
92
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
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 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:
nodo P =Presión en el nodo
r P =Presión media del yacimiento
wf P =Presión de fondo fluyendo
wh P =Presión en la cabeza del pozo
sep P =Presión del separador
ΔPres =Caída de presión en el yacimiento
L.D ΔP =Caída de presión en la línea de descarga
T.P Δ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.
Las expresiones para las curvas de inflow y outflow son las siguientes:
nodo V.T arriba P = P
Capítulo III
93
Inflow.
r res T.P debajo V.T nodo P − ΔP − ΔP − ΔP = P
Outflow
sep L.D T.P arriba nodo P + ΔP + ΔP = P
Donde:
V.T ΔP =Caída de presión a través de la válvula de tormenta
V.T arriba P =Presión corriente arriba de la válvula de tormenta
T.P debajo ΔP =Caída de presión en la tubería de producción desde el fondo del pozo hasta la
válvula de tormenta
T.P arriba ΔP =Caída de presión en la tubería de producción de la válvula de tormenta hasta la
cabeza del pozo
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, 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 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).
4) 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 calcular 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 apropiada correlación de flujo multifásico en tubería horizontal o curvas de
gradiente
de presión.
5) Utilizando los mismos gastos y las presiones en la cabeza del pozo (Pwh) obtenidas en el paso
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
Capítulo III
94
pozo, utilizando una apropiada correlación de flujo multifásico en tubería vertical o curvas de
gradiente de presión.
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.
Capítulo III
95
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.
Las expresiones para las curvas de inflow y outflow son las siguientes:
nodo top liner P = P
Inflow.
r res T.P 1 nodo P − ΔP − ΔP = P
Outflow
sep L.D T.P 2 nodo P + ΔP + ΔP = P
Donde:
top liner P =Presión en la parte alta del liner (donde se juntan las dos tuberías de producción)
T.P 1 ΔP =Caída de presión en la T.P que corre del fondo del pozo hasta la parte alta del liner
T.P 2 ΔP =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
Capítulo III
96
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, 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 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.
3) 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.
4) Utilizando los mismos gastos y las presiones en la cabeza del pozo (Pwh) obtenidas en el paso
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
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.
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.
Capítulo III
97
• 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:
r res T.P L.D A P − ΔP − ΔP − ΔP = P
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
Capítulo III
98
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
Capítulo III
99
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:
A AC C P − ΔP = P
B BC C P − ΔP = P
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:
A
ii
AC q
q RGL
RGL Σ =
A
i wi
wAC q
qF
FΣ=
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:
Capítulo III
100
D CD C P + ΔP = P
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.__
CAPÍTULO 3. CONCLUSIONES
Como resultado de la investigación de las disposiciones de San Telmo, Coipasa y Corregidores
presentadas en el proceso de extracción de petróleo en Bolivia ha sido a lo largo de varios años
motivo de grandes conflictos, puesto que ha generado una destrucción masiva al medio
ambiente. Lo que se ha pretendido en la investigación es realizar una comparación entre dos
escenarios propuestos, por un lado las reservas de los pozos, la cual propone mantener el crudo
pesado del campo en la contribución internacional. Bolivia percibiría en caso de que estas
reservas fuesen explotadas, contrastándola con la explotación petrolera en campos abiertos de
Bolivia en los tres campos descubiertos de San Telmo, Coipasa y Corregidores. Al analizar los
resultados de los valores presentes de ambos escenarios, tanto el Realista, Optimista y
Moderadamente Pesimista, podemos apreciar la realidad de Bolivia. Es cierto que la aplicación
de la explotación genera beneficios económicos en la economía boliviana, pero estos son
inferiores en relación al escenario de explotación. La propuesta de la investigación en
explotación es más sostenible en el tiempo, pues no afecta al medioambiente y los ingresos que
se obtuviese se destinarían a proyectos que tengan un impacto positivo para el país.

La optimización de la producción de petróleo mediante el análisis nodal con los diversos

software que existen en el mercado de la industria petrolera, es una herramienta para
poder lograr incrementar la producción de petróleo en los campos petroleros del Ecuador.
ANEXOS