Introducción

La curva IPR es una representación gráfica de la tasa de producción del pozo contra la presión del fondo
fluyendo. De la cual, la forma de la curva es gobernada por la composición del fluido del yacimiento y el
comportamiento de las fases fluidas en condiciones de flujo. Sirve para determinar el potencial de producción
de un yacimiento con respecto al sistema de producción y(o) a las instalaciones del Sistema Artificial de
Producción. Además de utilizarse para optimizar los parámetros de producción y para determinar el IPR para un
tiempo dado, para lo cual se realizan procedimientos iterativos.

Las curvas analíticas de IPR se pueden desarrollar para cualquier estado de agotamiento, siempre y cuando las
permeabilidades relativas y propiedades PVT de los fluidos se conozcan.

CLASIFICACIÓN DE LOS YACIMIENTOS DE PETRÓLEO

Los yacimientos de hidrocarburos se clasifican generalmente en gasíferos o petrolíferos. Normalmente, ciertos

yacimientos de hidrocarburos contienen fracciones líquidas y/o gaseosas dependiendo de las condiciones
ambientales. Las determinaciones matemáticas o experimentales de estas condiciones se expresan en
diferentes tipos de diagramas llamados “Diagramas de Fases”. El Diagrama de Fases más usado es el “Diagrama
Presión – Temperatura”. La Figura 4 muestra un Diagrama Presión – Temperatura típico con composiciones de
fases específicas.

Comportamiento de Afluencia (IPR) para pozos productores verticales.

Una gráfica de Comportamiento de Afluencia (IPR) representa la capacidad de entrega de fluidos de un

yacimiento petrolífero a un pozo productor. La presión de fondo fluyente es función del caudal de fluido. La
capacidad de entrega de fluidos de un yacimiento depende de las propiedades petrofísicas, las propiedades del
fluido y la presión estática. La caída de presión en gráficas IPR se define como la diferencia entre la presión
estática del yacimiento y la presión de fondo fluyente. Las tasas de producción son directamente proporcionales
a la caída de presión. La
Ecuación 3 presenta la caída de presión entre el yacimiento y el pozo productor.

La Figura 8 muestra la estructura típica de una curva IPR para flujo de petróleo y gas. Se
observa que el corte con la abscisa corresponde al caudal máximo ofrecido por el yacimiento
a un pozo productor. El corte con la ordenada representa la presión del yacimiento.

En 1935 Rawlins y Schellhardt hacen la primera conceptualización de la gráfica

IPR cuando analizaban el efecto del caudal de líquido en el rendimiento de la producción de
gas. Relacionaron por primera vez la presión en cabeza de pozo con la tasa de flujo de gas.

En 1954 Gilbert define el concepto de “IPR” (Comportamiento de Afluencia), enfatizando que

el potencial de producción de un pozo abarca todas las condiciones que afectan el flujo de
fluidos del yacimiento hacia la tubería.

Gilbert propuso un método empírico para determinar la curva de Comportamiento de Afluencia

en pozos con flujo natural basándose en datos medidos de: tasa de producción, relación gas-
líquido (RGL), presión de fondo fluyente (Pwf), presión en cabeza del pozo (Pwh) y presión
estática del yacimiento (Pr). El método empírico de Gilbert permite determinar la curva IPR
para un diámetro de tubería específico a una profundidad determinada. Gilbert evidencia que
para relaciones gas-líquido (RGL) altas, la gráfica IPR tiende a comportarse como una
curvatura. En la Gráfica 4 se observa la tendencia curva cuando aumenta la relación gas –
líquido (RGL).

Después de la definición de IPR por Gilbert, se han creado diferentes correlaciones para
realizar la curva IPR de un pozo. Cada correlación asume diferentes condiciones estáticas e
incorpora diferentes parámetros petrofísicos y del fluido. Fetkovich, Vogel, Wiggins, Del
Castillo y Elias et al son algunos de los autores que han postulado las correlaciones más
usadas actualmente. La correlación de Vogel es ampliamente usada debido a su bajo
requerimiento de datos y sencillez. Esta correlación no modela cambios por las saturaciones
de agua.

Comportamiento de Afluencia (IPR) para yacimientos subsaturados – Método de línea

recta de Darcy.
Este método genera la gráfica IPR para pozos produciendo de yacimientos subsaturados (𝑷𝒓
> Pb) debido a que los hidrocarburos están en una sola fase líquida (flujo monofásico), por
consiguiente, la Ley de Darcy describe su comportamiento. El tipo de flujo más similar al que
se presenta en un yacimiento real es el flujo pseudoestable, por lo tanto, se expondrá la
correspondiente forma de la Ley de Darcy. El coeficiente de turbulencia se considera nulo por
ser producción líquida y se considera que la permeabilidad es constante en todo el
yacimiento.

Matemáticamente el índice de productividad es la pendiente en un punto de una gráfica IPR.

Para los pozos produciendo de yacimientos subsaturados (𝑃𝑟 > Pb) la gráfica IPR se
comporta como una línea recta tal como lo muestra la Figura 9. Debido a esto el índice de
productividad es constante a cualquier caída de presión.
El caudal de petróleo se puede determinar con un medidor de flujo, la permeabilidad efectiva
mediante un análisis de núcleos, el espesor del reservorio mediante registros eléctricos, la
presión media y el factor de daño mediante el análisis de una prueba de presión (Horner), la
viscosidad mediante un análisis PVT del fluido. Odeh basado en el planteamiento de Mathews
y Russell, postuló los Factores de Forma como se muestra en la Figura 10, para poder
modelar pozos con diferentes geometrías de flujo mediante la ecuación derivada de la Ley de
Darcy.
Comportamiento de Afluencia (IPR) para flujo bifásico – Modelo de
Vogel – Yacimientos saturados. En 1968 Vogel desarrolló la correlación IPR más usada en
la actualidad. Dicha correlación empírica fue producto de la programación computacional de
las aproximaciones de Weller para yacimientos por debajo del punto de burbuja con flujo
bifásico (petróleo y gas). Las asunciones de dicha aproximación son:

 Yacimiento circular y completamente cerrado con un pozo en el centro.

 Medio poroso isotrópico con saturación de agua constante en todos los puntos.
 Efectos de la gravedad nulos.
 Compresibilidad de la roca y el agua nulas.
 La composición y el equilibrio son constantes para el petróleo y el gas.
 Existe la misma presión en las fases (petróleo y gas).
  GOR constante en cualquier tiempo.

En la programación computacional, Vogel varió la curva de permeabilidades relativas, las

propiedades PVT, estados de agotamiento, espaciamiento y permeabilidad absoluta del medio
para concluir que el aumento de la viscosidad del crudo, presencia de daño de formación y
yacimientos por encima del punto de burbuja generaban una linealidad en la gráfica IPR. Los
caudales de crudo son inversamente proporcionales a la viscosidad de este a unas
condiciones dadas. A medida que la depletación ocurre, la gráfica IPR se deteriora como se
muestra en la Figura 11, debido a que la presión del yacimiento disminuye y por ello el gas en
solución con el petróleo se libera y ocupa un volumen en las gargantas porales produciendo
una reducción de la permeabilidad relativa al petróleo. Vogel graficó los resultados de la
Figura 11 en variables adimensionales para observar algún comportamiento. En la Figura 12,
se evidencia que para bajos estados de agotamiento la curva IPR muestra un comportamiento
propio de una ecuación cuadrática. Altos estados de agotamiento generan curvas IPR
diferenciadas de la media. La correlación de Vogel se representa en la Ecuación 5.
Para realizar la curva de Comportamiento de Afluencia según el modelo de Vogel, es
necesario disponer de una prueba de flujo (medición de caudal de petróleo a una presión de
fondo fluyente dada). Con estos valores es posible determinar el caudal máximo de petróleo
(se presenta cuando la caída de presión es máxima; Pwf = 0).
Acto seguido, se asumen presiones de fondo fluyente menores a la presión de yacimiento
para obtener los respectivos caudales. El daño de formación según este modelo es nulo. Este
método genera un 20% de error promedio (Vogel, 1968), si se respeta el uso para condiciones
similares a las planteadas por Weller (1966).
Modelo de Vogel para yacimientos subsaturados
El método de Vogel puede ser extendido para yacimientos subsaturados. El comportamiento
rectilíneo de la gráfica IPR finaliza cuando se llega al punto de burbuja. Una vez el yacimiento
se satura, el gas en solución se libera y la gráfica IPR se torna curvilínea. Para realizar el
procedimiento es necesario conocer la presión del yacimiento (Pr), la presión de burbuja (Pb)
y una prueba de flujo.
El primer paso consiste en calcular el Índice de Productividad (J), que es constante cuando el
yacimiento es subsaturado (Pwf>Pb), con la utilización de la prueba de flujo como se muestra
en la Ecuación 6.

La gráfica IPR para este procedimiento simplemente consiste en una sección rectilínea (desde
Pr hasta Pb) y una sección curvilínea (desde Pb hasta Pwf=0).
Una vez se conoce el Índice de Productividad es posible calcular el potencial de producción a
la presión de burbuja (Qb). Acto seguido, se realiza el mismo procedimiento para la construir
la gráfica IPR con el método de Vogel, asumiendo presiones de fondo fluyente (Pwf) menores
a la presión de burbuja (Pb), mediante la Ecuación 7.

Comportamiento de Afluencia (IPR) para pozos con daño de formación –

Método de Standing.
En 1970 Standing basado en la ecuación de Vogel introduce el concepto de Eficiencia de
Flujo para observar las variaciones en los caudales de petróleo por efecto del daño de
formación. La eficiencia de flujo se expresa como la razón entre el índice de productividad real
e índice de productividad ideal. Para un pozo drenando un yacimiento cerrado y en flujo
pseudoestable definió la Eficiencia de Flujo como se muestra en la Ecuación 8.

El método es similar al de Vogel. Se debe conocer una prueba de flujo para poder calcular el
caudal máximo. Después se asumen presiones de fondo fluyentes con FE = 1 (Pwfi) para
después calcular las presiones de fondo fluyente con daño (Pwf).
Seguido a esto se determinan los caudales de petróleo correspondientes para lograr graficar
la curva IPR con Factores de daño incluidos.
En la Gráfica 6, se muestran los efectos presentados por Standing en variables
adimensionales cuando la eficiencia de flujo varía en una gráfica IPR.
Es posible que el Factor de Daño calculado mediante el análisis de una prueba de subida de
presión (PBU) de un pozo no corresponda al Factor de Daño propuesto por Standing cuando
se dispone de una prueba de flujo. Esto ocurre debido a que las consideraciones del modelo
de Vogel se alejan de las condiciones reales de los yacimientos.

Comportamiento de Afluencia (IPR) para flujo multifásico – Modelo de

Wiggins – Yacimientos saturados.
En 1992 Wiggins, Russell y Jennings crearon un método para determinar la gráfica de comportamiento
de afluencia para flujo de tres fases. El método consiste en una ecuación de cuarto grado que requiere
del conocimiento de las gráficas de permeabilidad relativa del yacimiento y de las propiedades del fluido
con su comportamiento con el cambio de presión. No siempre se dispone de esta información por lo
tanto el método se restringió a uso académico.
En 1993 Wiggins determinó una nueva correlación de Comportamiento de Afluencia
(IPR) para flujo multifásico asumiendo los efectos gravitacionales y capilares nulos.
Dicha correlación solo requiere de una prueba de flujo. Wiggins realizó la simulación numérica para
yacimientos completamente cerrados sin entrada externa de agua. Utilizó un amplio rango de
propiedades petrofísicas y del fluido dando paso al análisis del flujo de agua y petróleo de manera
separada. Por ello asignó una correlación IPR de comportamiento cuadrático para cada una de las fases
mencionadas, definiéndolas en la Ecuación 9.
Las pruebas de flujo disponibles generalmente se realizan en pozos con Daño de
Formación presente. Para poder determinar el caudal máximo de oferta de líquidos sin efecto de Daño
de Formación es necesario calcular la eficiencia de flujo (FE) con el Factor de Daño de Formación (S)
mediante la Ecuación 10. Después es necesario calcular el caudal máximo cuando la eficiencia de flujo
es diferente de 100% y de esta manera se hace posible determinar el caudal máximo cuando el Daño de
Formación es nulo.

Se realizó una prueba de flujo a un pozo produciendo de un yacimiento que posee 4500 psi de
presión promedio del yacimiento (Pr). El resultado de la prueba arrojó un caudal (Qo) de 550
BOPD con una presión de fondo fluyente (Pwf) de 1500 psi. Se cotejaron las gráficas IPR
propuestas por Darcy, Vogel y Wiggins en la Figura 13.
Método IPR de Fetkovich.
En 1973 Fetkovich demuestra que los pozos de petróleo pueden ser analizados con las
mismas ecuaciones de flujo introducidas por Rawlins y Schellhardt cuando analizaron el
potencial de producción de gas de un pozo. Para la realización de este método es necesario
disponer de una prueba “Multitasa Isócrona” o una prueba “Flujo tras flujo”. Para yacimientos
de alta permeabilidad ambas pruebas son válidas, en caso contrario, la prueba Isócrona
funciona debido a que se eliminan los efectos transientes que una prueba “Flujo tras flujo”
acarrea.
La correlación de Fetkovich asume que la movilidad del petróleo es función lineal con la
presión de reservorio tal como se muestra en la Figura 14. Fetkovich también asumió que el
mecanismo de producción es el de gas en solución. Para determinar la curva IPR según se
expresó la Ecuación 11.

Este método ha demostrado un porcentaje de error promedio de 7%. El método es confiable pero su
principal restricción es que es necesario disponer de una prueba Multitasa isócrona o Flujo tras flujo
para poder realizar la IPR, el método de Vogel y Wiggins solo requieren de una prueba de flujo.

Método IPR de Klins y Majcher.

En 1992 Klins y Majcher basados en el método de Vogel desarrollan una correlación de
Comportamiento de Afluencia (IPR) que tiene en cuenta el cambio en el punto de burbuja de
yacimientos con empuje por gas en solución. Determinaron la correlación expresada en la Ecuación 12.
Método IPR de Sukarno y Wisnogroho.
En 1995 Sukarno y Wisnogroho desarrollaron una correlación IPR basados en resultados de simulación
numérica de yacimientos. Esta correlación tiene en cuenta la Eficiencia de Flujo que es calculada
mediante el Factor de Daño y las constantes adimensionales propuestas por ellos. La correlación se
define en la Ecuación 13 y requiere del conocimiento de la Tabla 6.

Método IPR de Elias - Fattah et al.

En 2009 Elias et al. Propusieron una nueva correlación de Comportamiento de Afluencia (IPR)
para yacimientos con mecanismo de producción de gas en solución. Consideraron que la
movilidad del petróleo no era función lineal de la presión de reservorio. Para deducir esto
realizaron seis corridas de simulación numérica de yacimientos por debajo del punto de
burbuja encontrando una relación recíproca entre la movilidad del petróleo y la presión de
reservorio tal como se muestra en la Gráfica 7.

La metodología para realizar la gráfica IPR consiste en identificar si la presión del yacimiento
es menor o mayor a 1600 psia para hacer el respectivo cálculo del parámetro IPR del petróleo
(∝). Para utilizar este modelo se requiere una prueba de flujo. Con dicha prueba se calcula el
caudal máximo de petróleo (𝑄𝑜 𝑚𝑎𝑥). Después se determinan los caudales de petróleo a
diferentes presiones de fondo fluyente (𝑃𝑤𝑓) menores a la presión del yacimiento (𝑃 ). La
correlación desarrollada se muestra en la Ecuación 14.

Según Elias et al, este método ofrece el menor porcentaje de error entre los métodos
propuestos de Vogel, Fetkovich, Wiggins y Sukarno, después de analizar 12 casos como se
muestra en la Tabla 7.
Modelos IPR para producción de petróleo.