Ing. de Procutividad de Pozos

UNIV E R SID A D
GR UP O CE DI P
(CENTRO DE EDUCACIÓN Y DESARROLLO EN
INFORMÁTICA PERSONAL)
WWW.CEDIP.EDU.MX
UNIVERSIDAD
GRUPO CEDIP
MAESTRÍA EN
INGENIERÍA PETROLERA
Titulación Oficial Mexicana
TÍTULO
INGENIERÍA DE PRODUCCIÓN
Y LA PRODUCTIVIDAD DE POZOS
OFICIAL
CONOCE LOS PROGRAMAS ACADÉMICOS OFICIALES DE GRUPO CEDIP

/GRUPO CEDIP @GRUPO CEDIP @GRUPO CEDIP TV CEDIP
(833) 228 19 99 Blvd.Loma Real #302
www.grupocedip.edu.mx (833) 228 28 02 Fracc. Lomas del Chairel, Tampico Tamps.
Introducción
La ingeniería de producción es una de las grandes áreas de la ingeniería

petrolera, esta es la encargada del manejo y conducción de los hidrocarburos del
yacimiento a las áreas de comercialización, así como el ritmo de explotación,
método de explotación, diseño de instalaciones superficiales y optimización de los
pozos.
El ingeniero de producción tiene como prioridad administrar la energía del

yacimiento valiéndose de diversas herramientas que ayuden a minimizar la
perdida de energía, manteniendo un gasto de producción óptimo, así como la
planeación de intervenciones que ayudaran a mantener un adecuado desarrollo
de los pozos.
El transporte de hidrocarburos, el diseño de instalaciones superficiales, el gasto de

producción, el diseño de los métodos de explotación del campo, etc. Al ser
realizados de manera apropiada reducirán costos de producción, lo que hará que
el campo tenga una vida útil más duradera. El siguiente documento es una breve
recopilación de las herramientas básicas para el desarrollo de pozos, las cuales
juegan un papel importante en la ingeniería de producción, ya que proveen la
información para poder llevar a cabo tareas que generen condiciones óptimas para
mantener la rentabilidad del campo.
1
2
CAPITULO 1 USO DE CORRELACIONES PARA OBTENER LAS
PROPIEDADES DE LOS FLUIDOS A DIFERENTES
CONDICIONES DE FLUJO
1. USO DE CORRELACIONES PVT PARA OBTENER LAS PROPIEDADES DE

LOS FLUIDOS A DIFERENTES CONDICIONES DE FLUJO
Introducción
La evaluación de propiedades físicas de mezclas de hidrocarburos es un
importante paso para el diseño de diferentes etapas de operación de un campo
petrolero. Las propiedades de los fluidos varían conforme cambia la presión y la
temperatura, la mayoría de estas propiedades deben ser evaluadas por la
ingeniería de yacimientos e ingeniería de producción. Un sistema de tuberías con
dos o más fases fluyendo, requiere de la predicción de propiedades de los fluidos
como el gas disuelto, los factores de formación, compresibilidad del aceite, la
viscosidad en distintos puntos de la tubería. El comportamiento de estas
propiedades cambia con la temperatura y son determinadas mediante
correlaciones numéricas o a través de ecuaciones de estado con la ayuda de
simuladores numéricos.
1.1 Propiedades de los Fluidos Petroleros
1.1.1 Propiedades del Gas Natural
Las propiedades físicas de los componentes puros que se presentan en los gases
evaluadas a condiciones estándar de 14.7 lb/pg2abs y 60 °F, se listan en la Tabla
1.1.Estas propiedades, incluyen la formula química, el peso molecular,
temperatura y presión crítica, densidades del líquido y gas y viscosidad del gas
(para componentes más ligeros que el pentano).
Estas propiedades se utilizan en cálculos que se basan en reglas de mezclado

desarrollando las propiedades pseudo para mezcla de gases, incluyendo el peso
molecular aparente y la densidad relativa del gas. Las propiedades físicas
requeridas para estudios de ingeniería de yacimientos y de producción, se
describen a continuación.
3
Tabla 1.1 Propiedades físicas de los gases a 14.7 lb/pg2abs y 60 °F
Peso Presión Densidad Densidad Viscosi

Temperat
Formula molecular crítica del líquido del gas dad del
Componente ura crítica
química gas
(lbm/ (lb/pg2ab (lbm/ (lbm/
(°R)
lbm-mol) s) pies3) pies3) (cp)
Hidrógeno H2 2.109 59.36 187.5 4.432 0.005312 0.00871
Helio He 4.003 9.34 32.9 7.802 0.010550 0.01927
Agua H2O 18.015 1,164.85 3,200.14 62.336 - ~ 1.122
Monóxido de
CO 28.010 239.26 507.5 49.231 0.073810 0.01725
carbono
Nitrógeno N2 28.013 227.16 493.1 50.479 0.073820 0.01735
oxigeno O2 31.99 278.24 731.4 71.228 0.084320 0.02006
Sulfuro de 0.330898
H2S 34.08 672.35 1,306.0 49.982 0.01240
hidrógeno 10
Dióxido de
CO2 44.010 547.58 1,071.0 51.016 0.116000 0.01439
carbono
Aire - 28.963 238.36 546.9 54.555 0.076320 0.01790
Metano CH4 16.043 343.00 666.4 18.710 0.042280 0.01078
Etano C2H6 30.070 549.59 706.5 22.214 0.079240 0.00901
Propano C3H8 44.097 665.73 616.0 31.619 - 0.00788
i-Butano C4H10 58.123 734.13 527.9 35.104 - 0.00732
n-Butano C4H10 58.123 765.29 550.6 36.422 - 0.00724
i-Pentano C5H12 72.150 828.77 490.4 38.960 - -
n-Pentano C5H12 72.150 845.47 488.6 39.360 - -
n-Hexano C6H14 86.177 913.27 436.9 41.400 - -
n-Heptano C7H16 100.204 972.37 396.8 42.920 - -
4
n-Octano C8H18 114.231 1,023.89 360.7 44.090 - -
n-Nonano C9H20 128.256 1,070.35 331.8 45.020 - -
n-Decano C10H22 142.285 1,111.67 305.2 45.790 - -
Se dan valores de densidad de los líquidos para estos componentes, que puedan existir como
líquidos a 60 °F y 14.7 (lb/pg2abs) , se estima la densidad del líquido para componentes que
son gases naturales a estas condiciones.
Densidad del Gas
La densidad del gas (g) en lb/pie3, pude determinarse aplicando la ley de los
gases reales:
(1. 1)
Donde la constante R tiene un valor de 10.73 si las otras variables tienen las
unidades de lb/pg2, pie3, lb-mol, °R. El numero de moles puede ser expresado
como:
(1. 2)
Combinando esta ecuación y resolviendo la densidad (m/V) :
(1. 3)
Donde:
: Densidad del gas
: Gravedad especifica del gas
: Factor de compresibilidad del gas
: Presión
: Temperatura
5
Como el cambio de volumen esta dado por el factor volumétrico del gas B g , la
densidad del gasa P y T está dada por:
(1. 4)
Donde:
: Factor volumétrico del gas
Densidad específica de un gas (densidad relativa)
La densidad relativa de un gas, γg, es la relación de la densidad del gas a la

densidad del aire seco, ambos medidos a la misma presión y temperatura. La
densidad específica del gas en forma de ecuación se expresa como:
(1. 5)
en donde, g es la densidad de la mezcla de gases en lbm/pies3 y aire es la

densidad del aire en lbm/pies3.
A condiciones estándar, ambos el aire y el gas seco se modelan (comportan) en

forma muy precisa de acuerdo a la ley de los gases ideales (baja presión y
temperatura moderada). Bajo estas condiciones, si se emplea la definición de
número de moles, (n=m/M), y de densidad (=m/V), así como la ecuación de
estado para gases ideales para el aire y el gas, entonces la densidad relativa de
una mezcla de gases se puede expresar como:
(1. 6)
en donde, γg es la densidad relativa del gas (aire =1.0), M es el peso molecular

aparente del gas en lbm/lbm-mol y Maire es el peso molecular del aire e igual a
28.9625 lbm/lbm-mol
Si no se conoce el peso molecular del gas , pero se conoce los pesos

moleculares de las fracciones molares de la mezcla, Mj, se determina el peso
molecular aparente de la misma mediante :
6
(1. 7)
Donde:
: Peso molecular del componente j en la mezcla
: Componente j en la mezcla
: Número de componentes de la mezcla
Factor de desviación de los gases reales (factor de desviación z)
La ley de los gases ideales establece que la presión, temperatura volumen y la

cantidad de un gas ideal se relaciona con
(1. 8)
en donde, p es la presión en lb/pg2abs, V es el volumen en pies3, n es el número

de lbm-mol del gas, T es la temperatura en °R y R es la constante universal de los
gases e igual a 10.732 (lb/pg2abs-pies3)/(lbm-mol-°R).
Para gases reales, el factor de desviación z se define como:
(1. 9)
en donde, Vreal es el volumen real de gas que corresponde a una presión y

temperatura dadas específicamente, y Videal es el volumen de gas calculado con la
ecuación para gases ideales (ecuación 1.5).
Para gases reales, el factor de desviación z se calcula como,
(1.10)
7
El factor de compresibilidad z es adimensional. Comúnmente, z es muy próximo a

la unidad. Para gases cercanos al punto crítico, z se encuentra entre 0.3 y 0.4; si
la temperatura y la presión aumentan el factor de compresibilidad z incrementa
hasta un valor de 2 o mayor.
A partir de la ecuación 1.10 se desarrollan expresiones para definir el factor de

volumen de gas de formación, Bg, la densidad del gas, ρg, y la compresibilidad del
gas, cg.
Factor de volumen del gas de formación, Bg
El factor de volumen del gas de formación se define como la relación del volumen
de una cantidad de gas a las condiciones del yacimiento al volumen de la misma
cantidad del gas a las condiciones estándar,
(1.11)
para una p @c.e.=14.65 lb/pg2abs y una T @c.e.=60°F (519.59 °R), se tiene:
(1.12)
en donde Bg presenta unidades de pies3 @c.y./pies3 @c.e. La Fig. 1.1 presenta el

comportamiento isotérmico común del factor de volumen del gas de formación
respecto a la presión del yacimiento.
8
Bg(ft3@ c.y./ft3 @c.s.)
Presión del yacimiento, py, (lb/pg2 abs)
Fig. 1.1 – Forma común del factor de volumen de gas de formación, Bg, en función
de la presión del yacimiento a temperatura constante del yacimiento.
Coeficiente de compresibilidad isotérmico del gas, cg
El coeficiente de compresibilidad isotérmico del gas (compresibilidad del gas

isotérmico ó compresibilidad del gas), se define como el cambio fraccional del
volumen en función del cambio de presión a una temperatura constante; con signo
(-) para mantenerse positiva:
(1.13)
Relacionando esta ecuación con la de los gases reales:
(1.14)
De donde:
(1.15)
Derivando con respecto a “p” y sustituyendo en la ecuación 1.13:
9
(1.16)
La relación de cg respecto a la presión del yacimiento para un gas seco a

temperatura constante se presenta la Fig. 1.2.
cg(1/lb/pg2 abs.)
Presión del yacimiento, py, (lb/pg2 abs)
Fig. 1.2 – Forma común del coeficiente de compresibilidad isotérmica del gas, cg,
como una función de la presión del yacimiento y a temperatura constante del
yacimiento.
Las unidades del cg son 1/lb/pg2abs ó 1/p. Para un gas ideal la cg es

aproximadamente igual a 1/p (independiente de la temperatura). Esta expresión
(1/p) proporciona una estimación muy útil (aproximación) para el cálculo de cg,
especialmente a presiones bajas.
Se debe de entender que el término compresibilidad del gas, cg, se utiliza para
designar el coeficiente de compresibilidad isotérmica, cg, por lo que, el término de
factor de compresibilidad z se refiere al factor z, el coeficiente en la ecuación de
estado. Aunque ambos términos se relacionan para explicar el efecto de la presión
sobre el volumen de gas, ambos no son equivalentes.
Coeficiente de la viscosidad del gas
El coeficiente de viscosidad es una medida de la resistencia al flujo ejercida por un

fluido.
10
Si se mide el esfuerzo cortante y el gasto cuando un fluido se encuentra en

movimiento entre dos placas paralelas, en donde una placa se mueve con
respecto a la otra placa para un gasto cualquiera, se puede encontrar que el
esfuerzo cortante es directamente proporcional al gasto, es decir.
(1.17)
en donde la constante de proporcionalidad se define como la viscosidad del fluido,

μ, es decir:
(1.18)
luego entonces, la viscosidad es una medida de la resistencia que opone un fluido

al flujo. Al reciproco de la viscosidad se le llama fluidez. En unidades de campo,
la viscosidad se expresa en centipoises, cp, y se le denomina viscosidad dinámica.
La viscosidad dinámica se determina con:
Viscosidad dinámica = (densidad del gas)(viscosidad cinemática) (1.19)
(1.20)
en donde μ es la viscosidad dinámica en cp, g es la densidad del gas en gr/cm3 y

v es la viscosidad cinemática en centistokes. Un centistoke se define como un
centipoise dividido por la densidad; un centistoke es igual a 1 cm2/100 seg y un
centipoise es igual a 1 gr/100 seg-cm.
La viscosidad del gas, g, decrece conforme la presión del yacimiento

decrece. A baja presión las moléculas de los gases están separadas y se mueven
fácilmente una sobre otra.
La Fig. 1.3 muestra la forma de la viscosidad del gas como una función de la
presión del yacimiento para tres temperaturas diferentes de yacimiento. Se
observa que a presiones bajas la viscosidad del gas se incrementa conforme la
11
temperatura se incrementa. Sin embargo, a presiones altas la viscosidad del gas

decrece conforme la temperatura incrementa.
Viscosidad del gas, g
Incremento de
Temperatura, T
Presión del yacimiento, py
Fig. 1.3 – Comportamiento de la viscosidad del gas en función de la presión a tres

temperaturas diferentes.
1.1.2 Propiedades del Aceite
Densidad relativa del aceite, γo
La densidad específica o relativa de un aceite, ɣo, se define como la relación de

densidad del líquido a la densidad del agua, a las mismas condiciones de presión
y temperatura, es decir:
(1.21)
En el sistema Inglés de unidades se tiene, la o se expresa en lbm aceite/pies3

aceite y la w se expresa en lbm agua/pies3 agua.
Asimismo, la densidad relativa del aceite, ɣo, se puede expresar como la

densidad relativa 60°/60°, lo que significa que las densidades del líquido y del
agua se midieron a 60°F a la presión atmosférica. En la industria petrolera se
emplea la densidad en grados API que se define como:
12
(1.22)
En donde ɣo, es la densidad relativa del aceite a 60°/60°.
Factor de volumen de formación del aceite, Bo
El volumen de aceite que se produce en el tanque de almacenamiento a

condiciones estándar, es menor que el volumen de aceite que fluye del yacimiento
hacia el fondo del pozo productor. Este cambio en volumen del aceite se debe a
tres factores:
1. Liberación del gas disuelto en el aceite conforme la presión decrece desde la
presión del yacimiento a la presión de la superficie.
2. La reducción en la presión causa una expansión ligera del aceite remanente.
3. El aceite remanente también se contrae debido a la reducción en la
temperatura.
El factor de volumen de formación del aceite, Bo, se define como el volumen
de aceite del yacimiento que se necesita para producir un barril de aceite a
condiciones atmosféricas. El volumen de aceite del yacimiento incluye el gas
disuelto en el aceite.
(1.23)
o bien,
(1.24)
El volumen de aceite a condiciones de superficie o de tanque se reportan siempre

a 60°F, independiente de la temperatura del tanque; el volumen de líquido del
tanque de almacenamiento, al igual que el volumen de gas en superficie, se
reporta a condiciones estándar. Al Bo también se le llama factor de volumen de la
formación ó factor de volumen del yacimiento.
El factor de volumen de la formación del aceite, Bo, también representa el volumen

del yacimiento que ocupa un barril de aceite a condiciones estándar más el gas en
solución a temperatura y presión de yacimiento.
13
Liberación de gas en el espacio

poroso del yacimiento (el
líquido remanente @c.y. tiene
menos gas en solución).
Expansión del
líquido en el
yacimiento.
Bo
Bls aceite  gas disuleto @ c. y.
Bls aceite @ c.s.
pb
Presión del yacimiento,py, (lb/pg2 abs)
La Fig 1.4 representa el comportamiento típico del factor de volumen del aceite de
formación a la presión del yacimiento para un aceite negro.
Si la presión del yacimiento se pudiera reducir a la presión atmosférica, el valor del

factor de volumen de formación sería muy cercano a 1 [email protected]./[email protected]. Luego,
una reducción en temperatura a 60°F sería requerida para obtener un valor del
factor de volumen de formación igual 1 [email protected]./[email protected].
Por arriba de la presión de burbuja, el factor de volumen de formación disminuye

al tiempo que aumenta la presión (debido a la compresibilidad del aceite). Por
debajo de la presión de burbuja, el factor de volumen de la formación decrece al
disminuir la presión (por ejemplo, se vaporizan los componentes ligeros).
Relación gas en solución-aceite, RGA, o relación de solubilidad, Rs.
A la cantidad de gas disuelto en el aceite a condiciones de yacimiento se le

denomina relación gas en solución-aceite, RGA, o relación de solubilidad, Rs.
La relación de gas en solución-aceite, es la cantidad de gas que se libera del

aceite desde el yacimiento hasta las condiciones de superficie. La relación de gas
en solución-aceite, Rs, es la relación del volumen de gas producido a condiciones
estándar respecto al volumen de aceite producido a condiciones estándar (medido
a condiciones del tanque de almacenamiento), como resultado de la producción de
un volumen original de aceite a condiciones de yacimiento. La relación gas
14
disuelto-aceite, Rs, se define en términos de las cantidades de gas y aceite que se

producen en la superficie:
(1.25)
(1.26)
Es decir, cuando un barril de aceite a condiciones de yacimiento se produce a la

superficie a través de un separador hacia el tanque de almacenamiento, el aceite
podría estar acompañado por una cantidad de gas.
Los volúmenes en superficie del gas y el líquido se referencían a condiciones

estándar, por lo que, las unidades para la relación gas disuelto-aceite son pies
cúbicos a condiciones estándar por barril a condiciones de tanque o condiciones
estándar, pies3 @ c.e./Bls a c.s.
Volumen de gas
liberado @ c.y.
Rsb (cantidad de gas en

solución a pb
(Rsi-Rsb)
Rs(Bls aceite + gas disuelto)
@ c.y./Bls aceite @ c.s.
(Rsb-Rs) = gas liberado en el

yacimiento (gas libre)
pb
Presión, p, (lb/pg2 abs)
15
La Fig. 1.5 muestra el comportamiento de la relación gas en solución-aceite para

un aceite negro respecto al cambio de presión del yacimiento a una temperatura
de yacimiento constante.
A presiones del yacimiento por arriba de la presión de burbuja, se observa que

existe una línea horizontal (relación de solubilidad constante). Esto se explica
debido a que estas presiones el gas no se libera en el espacio poroso y la mezcla
total de líquido se produce dentro del pozo.
A presión de yacimiento por debajo de la presión de burbuja, la relación gas

disuelto-aceite, Rs, decrece conforme decrece la presión del yacimiento. Esto se
explica debido a que más y más gas se libera en el yacimiento, quedando
atrapado en el casquete de gas y no dejando que fluya hacia los pozos
productores, dejando menos cantidad de gas disuelto en el líquido.
Compresibilidad del aceite, Co
En general, la compresibilidad isotérmica de un fluido , Co en (lb/pg2)-1 , se define

como el cambio fraccional en volumen cuando la presión es cambiada a
temperatura constante, esto es:
(1.27)
30
La co (aceite negro) es virtualmente constante con

excepción de presiones cercanas a la pb. Los
co (lb/pg2 abs)-1
valores de co raramente exceden 35 x 10-6 (lb/pg2

abs)-1
0
pb
Presión del yacimiento, p, (lb/pg2 abs)
16
Fig. 1.6 - Comportamiento de la co respecto a la py para un aceite negro a

temperatura constante a p>pb
Coeficiente de viscosidad del aceite, μo.
La viscosidad es una medida de la resistencia al flujo ejercida por un fluido. La

viscosidad del aceite generalmente se expresa en centipoise. La viscosidad del
aceite es afectada por la presión y la temperatura es decir, un incremento en la
temperatura provoca un decremento en la viscosidad, una disminución en la
presión provoca una disminución en la viscosidad, un decremento en la cantidad
de gas en solución en el líquido provoca un incremento en la viscosidad, siendo la
cantidad de gas en solución una función directa de la presión. La Fig. 1.7 presenta
la relación entre la viscosidad de un aceite negro respecto a la presión, a una
temperatura del yacimiento constante.
Por arriba de la presión de burbuja, la viscosidad del aceite en un yacimiento

decrece casi lineal conforme la presión decrece. A presiones por debajo de la
presión de burbuja, la viscosidad del aceite se incrementa conforme la presión
decrece.
En el rango de p≥pb conforme p decrece la μo decrece debido a que las moléculas

se alejan más una de otra y se mueven más fácilmente. Conforme la presión del
yacimiento decrece por debajo de la presión de burbuja, el líquido cambia sus
composición, el gas que se libera toma las moléculas más pequeñas (ligeras) del
líquido, dejando al líquido remanente en el yacimiento con más moléculas con
formas más complejas (más pesadas). Este cambio en la composición del liquido
provoca un gran incremento en la viscosidad del aceite en el yacimiento conforme
la presión decrece por debajo del punto de burbuja.
En un yacimiento de aceite negro conforme la producción de aceite se reduce, la

presión en el yacimiento decrece y se tiene un menor empuje del aceite hacia los
pozos productores, debido a que el gas libre trata de ocupar el espacio para fluir,
asimismo la viscosidad del aceite se incrementa.
17
Existe un cambio de
composición en el líquido
(moléculas más complejas)
Si p decrece la µo decrece (las
moléculas del fluido se
separan y se mueven una
respecto a otra con mayor
facilidad)
o (cp)
pb
Presión del yacimiento, py, (lb/pg2)
Fig. 1.7 – Forma común de la viscosidad del aceite como una función de la presión
a una temperatura de yacimiento constante.
Factor de volumen total de la formación o factor de la fase mixta, Bt
El factor de volumen total de la formación se define como:
(1.28)
En donde cada término de la ecuación 1.28 se expresa como:
(1.29)
(1.30)
(1.31)
18
(1.32)
Bt incrementa si py disminuye (Vt aumenta)

(El factor Bg(Rsi-Rs) incrementa
Bt (Bls aceite y gas en solución + Bls gas) a
Bt (Bls aceite y gas en solución + Bls gas) a

significativamente conforme py disminuye)
Si la ppb, Bt=Bo (todo
el gas se encuentra en
(Bls aceite a c.e.)
solución en el líquido a
c.y.)
Bt  Bo  Bg ( Rsi  Rs )
(Bls aceite a c.e.)

c.y.
c.y.
Bodisminuye si py
disminuye (Vo decrece)
(el espacio poroso en el Bo=Bt aumenta si py
yacimiento se va Bo disminuye (incrementa
acumulando más y más Bo ligeramente debido a la
gas, y el líquido expansión del líquido
remanente en el
pb en el yacimiento)
Presión del yacimiento, p, (lb/pg2 abs)
yacimiento tiene menos
gas en solución)
Fig. 1.8 – Comportamientos típicos del Bt y Bo para un aceite negro como una
función de la presión a temperatura del yacimiento constante.
1.2 Clasificación General de los Yacimientos de Acuerdo a su Diagrama de

Fase
La determinación del tipo de fluido es de suma importancia, ya que este estipula el
tipo y medidas del equipo en superficie, las técnicas para determinación de
reservas, tipo de terminación del pozo, tipo de fluido a inyectar en la etapa de
recuperación secundaria, método de recuperación mejorada a utilizar, etc.
Existen cinco tipos de fluidos de yacimiento. Usualmente llamados: aceite negro,

aceite volátil, gas y condensado, gas húmedo y gas seco. El tipo de fluido de un
yacimiento puede ser confirmado solamente por observación en el laboratorio.
Pero de manera más práctica y mediante información de datos producción,
podemos definir de manera aproximada con solo tres propiedades el tipo de fluido
producido: relación gas-aceite (RGA= ), gravedad del
líquido en el tanque de almacenamiento y color del líquido en el tanque de
19
almacenamiento. Siendo la RGA el indicador más importante para la

determinación del tipo de fluido, otra manera es por medio de su comportamiento
termodinámico de una mezcla natural de hidrocarburos, tomando como base su
diagrama de comportamiento de fases, el cual indicara dependiendo la presión y la
temperatura, la fase en la que se encuentran los fluidos del yacimiento.
En una grafica temperatura-presión, la curva llamada envolvente de fases, que

resulta de unir los puntos de burbujeo y de roció que presenta la muestra a
diferentes temperaturas, las curvas se unen en el punto denominado crítico (Pc).
Fig. 1.9 Diagrama de fase
La envolvente de fase divide el diagrama en tres regiones, la primera, llamada

región de líquidos, está situada fuera de la envolvente de fases y a la izquierda de
la isoterma critica; la segunda, llamada región de gases, también esta situada
fuera de la envolvente de fases y a la derecha de la isoterma critica; la ultima,
encerrada por la envolvente de fases, se conoce como región de dos fases; en
esta región, se encuentran todas las combinaciones de temperatura y presión en
que la muestra de hidrocarburos puede permanecer en dos fases en equilibrio,
existiendo dentro de ella, las llamadas curvas de calidad, que indica el porcentaje
del total de hidrocarburos que se encuentra en estado líquido. Todas estas curvas
inciden en el punto crítico. Se distinguen, además, en el mismo diagrama la
Cricondenterma y la Cricondenbara, que son la temperatura y presión máximas,
respectivamente, a las cuales la mezcla de hidrocarburos puede permanecer en
dos fases en equilibrio.
Yacimiento de Aceite Negro
En este tipo de yacimientos al alcanzar la presión de burbujeo, empieza a variar la

composición de los fluidos producidos y por lo tanto cambiara el diagrama de
fases de los hidrocarburos remanentes. El punto en el que se tiene la presión y la
20
temperatura en el separador, indica que aproximadamente el 85 % del aceite

producido es líquido, esto es un porcentaje promedio alto, de ahí que este aceite
es denominado aceite de bajo encogimiento. Estos tipos de yacimientos producen
generalmente un liquido negro o verde negruzco, con una densidad relativa mayor
de 0.800 y una relación gas-aceite instantánea menor a 200 (m3 de gas/ m3 de
aceite).
En nuestro país este tipo de yacimientos es de los que producen mayor volumen y
se les denomina como de aceite negro o de aceite y gas disuelto de bajo
encogimiento, por que la reducción en el volumen de aceite, al pasar de la presión
del yacimiento a la presión del tanque, es relativamente pequeña, normalmente su
reducción es menor a un 30 %.
Fig.1.10 Diagrama de fase del

aceite negro
Color: negro o verde negrusco
Densidad relativa: mayor de 0.800
RGA: menor 200 m3/m3.
Yacimiento de Aceite Volátil
La división práctica entre el aceite negro y el aceite volátil se da cuando las

ecuaciones de aceite negro dejan de ser adecuadas. La Tc > Ty pero son muy
cercanas. Llamado de alto encogimiento dado que pequeñas caídas de presión
liberan grandes cantidades de gas. La Rsi anda en el rango de 2000 a 3200 [pie3/
bl], el Bo es mayor de 2, la densidad es mayor de 40 °API, color café, verde y
naranja.
Los yacimientos de aceite volátil o también denominado de alto encogimiento,

contienen fluidos que, como su nombre lo indica, tienden a volatilizarse o
evaporarse significativamente con reducciones en la presión, una vez alcanzado el
punto de burbujeo. Como resultado de este fenómeno, el volumen ocupado por el
21
líquido se reduce considerablemente y de ahí el nombre de aceite de alto

encogimiento.
Se caracterizan por tener una RGA relativamente alta (con valores intermedios
entre la RGA de los yacimientos saturados y la RGA de los yacimientos de gas y
condensado) y un factor de volumen del aceite bastante alto (cercano o mayor a
dos).
Fig. 1.11 Diagrama de fase para un

yacimiento de aceite volátil
Color: ligeramente oscuro
. Densidad relativa: entre de 0.800 y
0.75
RGA: entre 200 y 1,000 m3/m3.
Yacimiento de Gas y Condensado

Contiene menos componentes pesados que los aceites, la temperatura del
yacimiento se encuentra entre la Tc y la cricondenterma Tc<Ty.
En este tipo de yacimientos si se encuentra por arriba de la presión de roció, se

tiene gas 100% en el yacimiento, pero al bajar la presión se condensará formando
liquido (este fenómeno va en contra de la naturaleza general de los gases, por lo
cual se le llama condensación retrograda).
Este tipo de yacimientos se encuentra en fase gaseosa a las condiciones iniciales,

la condensación retrograda, en el yacimiento es mayor que la indicada en curvas
presentadas en los reportes de análisis PVT a volumen constante a la temperatura
del yacimiento con decrementos de presión relativamente grandes (500 a 1000
psi). Esto se explica por el hecho de que una vez que empieza la condensación de
líquidos (los cuales permanecen inmóviles dentro de la roca) la composición de la
mezcla de los fluidos remanentes en el yacimiento se va haciendo cada vez más
rica en compuestos pesados y como consecuencia la tendencia a la acumulación
22
de líquidos se incrementa, esto se evidencia por el hecho de que a bajo de la

presión de rocío, el contenido de condensado por unidad de volumen de fluido
producido disminuye considerablemente y generalmente aumenta la RGA.
Fig. 1.12 Diagrama de fase para

un yacimiento de gas y
condensado
Color: café ligero
Densidad relativa: entre 0.74 y
0.78
RGA: entre 1,000 y 10,000
m3/m3.
Yacimiento de Gas Húmedo

No se forma líquido en el yacimiento en ninguna etapa de su explotación, durante
su producción se forma líquido en algún punto de las tuberías de producción
llamado "condensado", cuando estos fluidos son llevados a la superficie entran en
la región de dos fases, generando que la RGA varié entre 10000 y 20000 m3/m3, el
líquido recuperable tiende a ser transparente, con densidades menores a 0.75
gr/cm3 y el contenido de licuables en el gas, generalmente es bajo, menos de 30
Bls/106 pies3.
23
Fig. 1.13 Diagrama de fase de un

yacimiento de gas húmedo
Color: acuoso
Densidad relativa: menor a 0.740
RGA: entre 10,000 y 20,000
m3/m3.
Yacimiento de Gas Seco
Constituido principalmente por metano, posee algunos intermedios, no se forma
líquido en el yacimiento ni en ningún punto del sistema integral de producción. La
palabra seco se refiere a que este gas no contiene los suficientes componentes
pesados para formar líquidos a condiciones superficiales, los yacimientos de gas
seco, contienen solamente gas a lo largo de toda la vida productiva de éste, ni a
las condiciones de yacimiento, ni a las condiciones de superficie se entra en la
región de dos fases, durante su explotación por lo que siempre se está en la
región de estado gaseoso, teóricamente, los yacimientos de gas seco no producen
líquidos en la superficie, sin embargo, la diferencia entre un gas seco y un gas
húmedo es arbitraria generalmente un sistema de hidrocarburos que produzca con
una relación gas-aceite mayor de 20000 m3/m3, se considera gas seco.
24
Fig. 1.14 Diagrama de fase de un

yacimiento de gas seco
Color: acuoso
Densidad relativa: menor a 0.740
RGA: mayor a 20,000 m3/m3.
1.3 Análisis de los Fluidos Petroleros
El objetivo de la caracterización de los fluidos petroleros es el de predecir el

comportamiento termodinámico de los hidrocarburos, en cualquier punto donde
estos se muevan, desde el yacimiento hasta las baterías de separación.
Para conocer las propiedades de los fluidos, a diferentes presiones y

temperaturas, se utilizan generalmente dos métodos, directos (Pruebas PVT,
Análisis Cromatográfico, Análisis Sara) indirectos (Ecuaciones de Estado,
Correlaciones Numéricas PVT) .
1.3.1 Métodos directos
Pruebas PVT
Las pruebas PVT consisten en una serie de procedimientos de análisis en el

laboratorio, diseñado para proveer valores de las propiedades físicas del aceite y
del gas.
Los estudios realizados a partir de estas pruebas dan como resultado la obtención
de las siguientes propiedades:
 Presión en el punto de burbuja
25
 Factor de volumen del aceite

 Factor de solución gas-aceite
 Coeficiente de compresibilidad isotérmica del gas
 Viscosidad del aceite y del gas
Las pruebas se realizan colocando una cantidad de la mezcla de hidrocarburos en

unas celdas llamadas celdas PVT, las cuales contienen un fluido que ejercerá
presión sobre la mezcla analizada.
El mercurio es el fluido que se utiliza frecuentemente en este proceso.
Con este experimento se determina el comportamiento de los fluidos a

determinada temperatura y realizándolo a otras temperaturas se obtienen los
llamados “diagramas de fase”.
Existen varios procedimientos de laboratorio para análisis PVT: separación a

composición constante (este procedimiento también es llamado vaporización flash,
liberación flash o expansión flash), separación diferencial por el método
convencional, separación diferencial por el método a volumen constante y
simulación de condiciones de separación en el campo. Los tres primeros tratan de
simular el comportamiento de los fluidos en el yacimiento y el último, en superficie.
Análisis cromatográfico
La Cromatografía Gaseosa es un procedimiento de análisis para separar,

identificar y cuantificar los diferentes componentes de una mezcla. Las mezclas a
analizar pueden estar inicialmente en estado gaseoso, líquido o sólido, pero en el
momento del análisis la mezcla debe estar vaporizada.
La composición de una muestra de gases es, probablemente, el parámetro más

importante de la misma. Esta característica es el resultado de la facilidad con que
pueden estimarse las distintas propiedades de un gas a partir de su composición,
ya sea mediante correlaciones, ecuaciones de estado o cálculo de propiedades
aditivas.
Cuando se conoce la composición una mezcla de gases es posible obtener los

siguientes parámetros con razonable exactitud
 Densidad.
 Viscosidad.
 Poder Calorífico.
26
 Presión de Rocío.
El análisis cromatográfico de una mezcla gaseosa proporciona los datos

necesarios para calcular, con razonable exactitud, casi todos los parámetros de
interés de dicha mezcla. Esto es así gracias a dos factores concurrentes:
 El cromatograma de una mezcla de gases proporciona la composición

individual de todos los componentes mayoritarios de la mezcla.
 Las propiedades de los gases son básicamente aditivas. En otras palabras,
las propiedades de la mezcla son calculables a partir de las propiedades de
los componentes individuales.
En el caso de mezclas líquidas de hidrocarburos, las dos condiciones

mencionadas no se cumplen o sólo lo hacen parcialmente.
 El sistema cromatográfico convencional retiene componentes "pesados"

tales como parafinas de elevado peso molecular, resinas, asfaltenos, etc.
Estos componentes, presentes en todos los Petróleos Negros, suelen
representar una fracción particularmente importante en los petróleos con
menos de 40 °API.
 Muchas de las propiedades de las mezclas líquidas (y en especial la
viscosidad) no son propiedades aditivas.
Análisis SARA
Es un método de caracterización de aceites pesados basado en fraccionamiento,

por el cual se se separa una muestra de petróleo pesado en cantidades más
pequeñas o fracciones, de modo de que cada fracción tenga una composición
diferente. El fraccionamiento se basa en la solubilidad de los componentes de
hidrocarburo en diversos solventes utilizados en esta prueba cada fracción costa
de una clase de solubilidad que contiene un rango de diferentes especies de peso
molecular. En este método, se fracciona el petróleo crudo en cuatro clases de
solubilidad, llamadas colectivamente “SARA”: saturados, aromáticos, resinas y
asfáltenos. Los saturados son, en general, iso y ciclo parafinas, mientras que los
aromáticos, las resinas y los asfáltenos forman un continuo de moléculas con peso
molecular, aromaticidad y contenido de heteroátomos en aumento. Los asfáltenos
también pueden contener metales como níquel y vanadio. En ocasiones, este
método se denomina análisis de deposición de asfalteno/cera/hidrato.
27
1.3.2 Métodos indirectos
Ecuaciones de estado
Una ecuación de estado es una expresión analítica que relaciona la presión

respecto a la temperatura y el volumen. Una descripción adecuada de una relación
PVT para fluidos hidrocarburos reales es básica para determinar los volúmenes de
los fluidos y el comportamiento de fase de los fluidos petroleros, así como en el
pronóstico del comportamiento de los fluidos para el diseño de las instalaciones de
separación en la superficie.
Durante el desarrollo de la ecuación de estado para gases ideales, se

consideraron dos suposiciones, la primer suposición considera que el volumen de
las moléculas de gas es insignificante en comparación con el volumen del
contenedor (recipiente que lo contiene) y la distancia entre las moléculas, y la
segunda suposición considera la inexistencia de fuerzas de atracción o de
repulsión entre las moléculas del gas o entre el gas y las paredes del recipiente
que las contiene.
Correlaciones Numéricas PVT
Las correlaciones PVT son desarrolladas apartir de datos de laboratorio y/o de

campo y formuladas de manera que puedan ser utilizadas con datos obtenidos sin
mucha pérdida de tiempo y/o inversión de esfuerzo. Estos datos son:
a) Propiedades de los fluidos

 Gravedad específica del aceite
 Gravedad específica del gas
 Relación gas-aceite
b) Propiedades del Yacimiento
 Presión
 Temperatura
28
Antes de determinar una propiedad, se debe asegurar que la aplicación de interés

se encuentra dentro del rango de los datos para la cual la correlación fue
desarrollada. Una vez hecho esto, la correlación utilizada será adecuada y se
obtendrán resultados confiables.
1.4 Ecuaciones de Estado
Ecuación de estado de van der Waals
Una de las primeras ecuaciones de estado para representar el comportamiento de

los gases reales fue desarrollada por van der Waals en 1873. van der Waals se
enfocó en tratar de eliminar las dos suposiciones mencionadas anteriormente
mediante su ecuación de estado, explicando que las moléculas de gas ocupan una
fracción significante de volumen a presiones altas.
Para contrarrestar la primera suposición van der Waals propuso que el volumen de
las moléculas, representado por el parámetro b debe de restarse del volumen
molar real VM en la ecuación de los gases ideales; arreglando esta última ecuación
se tiene:
(1.33)
en donde el parámetro b representa el volumen de las moléculas que se restan y

VM es la variable que representa el volumen molar real, ambos en ft3 por una
unidad de lbm-mol.
Para eliminar la segunda suposición, van der Waals restó un término de

corrección, denominado a/V2M a la ecuación
(1.34)
O bien
(1.35)
29
Los parámetros a y b representan constantes y caracterizan las propiedades

moleculares de cada componente de la mezcla.
Ecuaciones de estado en el punto crítico
Van der Waals al tratar de determinar experimentalmente los valores de las

constantes a y b para cualquier sustancia pura, en un diagrama de fase de presión
contra volumen molar observó que la isoterma crítica, representa una pendiente
horizontal y un punto de inflexión en el punto crítico, por lo que las constantes a y
b se determinan de las siguiente manera (Mcain, the properties of petroleum fluids,
pag. 131-132)
(1.36)
(1.37)
Ecuación de estado cúbica de van der Walls de dos parámetros
30
La ecuación de estado de van der Walls en forma cúbica y en términos de VM se

expresa como:
(1.38)
El término de ecuación de estado cúbica implica una ecuación que expandida

contiene términos de volumen molar, VM, elevados a la primera, segunda y tercera
potencia. Una de las características de mayor importancia de la ecuación 1.38, es
la caracterización del fenómeno condensación–líquido y del cambio de estado de
la fase gas a la fase líquida conforme el gas se comprime.
Expresando la ecuación 1.38 en términos del factor de compresibilidad z, se tiene:
(1.39)
Donde:
(1.40)
(1.41)
La solución de la ecuación de estado 1.38 proporciona una raíz real para z en la

región de una fase (en algunas regiones súper críticas esta ecuación proporcionan
tres raíces reales para z, seleccionando la mayor raíz como el valor con
significado físico para el factor de compresibilidad, z) y tres raíces reales en la
región de dos fases (en donde la presión del sistema es igual a la presión de vapor
de la sustancia). En este último caso, la raíz real con el mayor valor numérico
corresponde al factor de compresibilidad de la fase vapor (gas), zv, mientras que la
raíz real con el menor valor numérico corresponde al factor de compresibilidad de
la fase líquida, zL.
Ecuación de estado de Redlich-Kwong (EdE de RK)
En 1948 Redlich y Kwong modificaron el término de fuerzas de atracción (fuerzas

de presión a/V2) de la ecuación de van der Waals, lo cuál mejora en forma
considerable la predicción de las propiedades físicas y volumétricas de la fase
gas. Redlich-Kwong sustituyeron el término de fuerzas de atracción de presión,
con un término general de dependencia de la temperatura. La ecuación de
Redlich-Kwong se expresa como:
31
(1.42)
Dónde
(1.43)
expresando la ecuación de estado de Redlich-Kwong (ecuación 4.30) a las

condiciones en el punto crítico
(1.44)
Dónde:
(1.45)
(1.46)
en donde a y b son constantes e iguales a 0.427481 y 0.08664,

respectivamente.
Expresando la ecuación de Redlich-Kwong en términos del factor de

compresibilidad se tiene:
(1.47)
De donde:
(1.48)
(1.49)
32
Ecuación de estado de Soave-Redlich-Kwong
Soave en 1972 realizó una modificación en la evaluación del parámetro a en las

presiones de atracción de la ecuación de Redlich-Kwong. Soave reemplazo el
término de a/T1/2 con un término dependiente de la temperatura, aT, es decir:
(1.50)
En donde aT está dado por:
(1.51)
siendo  el término dependiente de la temperatura y es adimensional. Cuando la

temperatura del sistema es igual a la temperatura crítica (T =Tc),  tiene el valor
de una unidad y ac es el valor de aT a la Tc.
A temperatura diferente a la temperatura crítica, el parámetro  se define por:
(1.52)
en donde el parámetro m se correlaciona con el factor acéntrico de Pitzer, ,

como:
(1.53)
siendo el factor acéntrico de Pitzer, , definido a una Tr de 0.7, como:
(1.54)
en donde pvr es la presión de vapor reducida evaluada a una Tr de 0.7. El factor

acéntrico de Pitzer, , es un valor constante para cada sustancia pura.
Expresando la ecuación de Sove-Redlich-Kwong en forma cubica:
33
(1.55)
Introduciendo el factor z dentro de la ecuación 1.55 al reemplazar el volumen

molar, VM , en la ecuación con (zRT/p) y arreglando:
(1.56)
en donde:
(1.57)
(1.58)
Ecuación de Estado de Peng-Robinson (EdE de PR)
Peng y Robinson, PR, en 1975 realizaron un exhaustivo estudio para evaluar el

uso de la ecuación de estado de Soave-Redlich-Kwong, SRK, y predecir el
comportamiento de los hidrocarburos. Peng y Robinson, PR, mejoraron la
ecuación de estado de SRK para predecir las densidades de líquidos y otras
propiedades físicas, principalmente en la vecindad de la región crítica. Ellos
propusieron un ligero cambio en el término de atracción molecular de la ecuación
1.42, es decir:
(1.59)
rearreglando la ecuación 1.59 en la forma del factor de compresibilidad,
(1.60)
en donde A y B se dan para componentes puros mediante:
34
(1.61)
(1.62)
1.5 Correlaciones Numéricas PVT

Las correlaciones son de suma importancia ya que nos darán un valor aproximado
de las propiedades del aceite por medio de una ecuación empírica, enunciaremos
las más importantes o bien las más empleadas en la industria por su credibilidad y
certeza.
Propiedades del aceite saturado
Correlación de M.B. Standing

Establece las relaciones empíricas observadas entre la presión de saturación y el
factor de volumen del aceite, en función de la razón gas disuelto-aceite, las
densidades del gas y del aceite producido, la presión y la temperatura.
La correlación se estableció para aceites, gases producidos en California y

sistemas de crudos de bajo encogimiento, simulando una separación instantánea
en dos etapas a 100° F. La primera etapa se realizo a una presión de 250 a 450
lb/pg2 abs. Y la segunda etapa a la presión atmosférica.
Debe entenderse que la densidad del aceite producido en el tanque de

almacenamiento dependerá de las condiciones de separación (etapas , presiones
y temperatura) .Mientras más etapas de separación sean, el aceite será más ligero
(mayor densidad API).
Así, el Bo fue correlacionado con Rs, la temperatura, la densidad relativa del gas y
la densidad del aceite.
Procedimiento:
La presión del aceite saturado se correlacionó con la relación de solubilidad (Rs)

de la siguiente forma:
35
(1.63)
Por lo que despejando la relación gas disuelto-aceite (Rs) de la ecuación anterior

se tiene:
(1.64)
El factor de volumen del aceite fue relacionado con la relación gas disuelto-aceite,
la temperatura, la densidad relativa del gas y la densidad del aceite. Se obtuvo la
siguiente expresión:
(1.65)
Dónde:
(1.66)
: Relación de solubilidad [pies3/barril de aceite]
: Gravedad especifica del aceite
: Gravedad especifica del gas disuelto
: Presión del sistema [lb/pg2]
: Densidad del petróleo
T: Temperatura en [°F]
: Factor de volumen del aceite [pies3@ c.y. /pies3@ c.s.]
Correlación de Vázquez-Beggs
Para establecer estas correlaciones se usaron más de 6000 datos de Rs, Bo y μo, a
diferentes presiones y temperaturas. Como el valor de la densidad relativa del gas
es un parámetro de correlación importante, se decidió usar un valor de dicha
36
densidad relativa normalizada a una presión de separación de 100 lb/pg2

manométrica.
Procedimiento:
El primer paso para usar estas correlaciones para el aceite saturado consiste en
obtener el valor de la densidad relativa del gas a dicha presión.
(1.67)
La correlación para determinar Rs se afinó dividiendo los datos en dos grupos, de

acuerdo con la densidad del aceite. Se obtuvo la siguiente ecuación:
(1.68)
Los valores de los coeficientes son:
Constante ° API ≤ 30 ° API > 30
C1 0.0362 0.0178
C2 1.0937 1.1870
C3 25.724 23.931
C4 4.677 x 10-4 4.67 x 10-4
C5 1.751 x 10-5 1.1 x 10-5
C6 -1.811 x 10-8 1.337 x10-9
Tabla 1.1 Coeficientes utilizados en la correlación de Vazquez-Beggs
La expresión que se obtuvo para determinar el factor de volumen es:
(1.69)
37
: Densidad relativa del gas normalizada a una presión de separación de 100

lb/pg2 manometricas.
: Densidad relativa del gas a condiciones de Ps y Ts
: Temperatura en la primera etapa de separación [°F]
: Presión de la primera etapa de separación [lb/pg2]
Correlación de Oistein-Glaso
Esta correlación fue establecida utilizando muestras de aceite producido en el Mar
del Norte, donde predominan los aceites de tipo volátil. Por lo que, el cálculo de Rs
y Bo se efectúan dependiendo si es un aceite tipo volátil o aceite negro.
Procedimiento:
Los valores de la relación de solubilidad y del factor de volumen del aceite

saturado se obtienen mediante los pasos siguientes:
Calcule p* con:
(1.70)
Calcule Rs con:
(1.71)
Donde:
Para lograr un mejor ajuste se puede variar el valor del exponente a.
Calcule B*o con:
(1.72)
38
Determine Bo con:
(1.73)
Donde:
(1.74)
Correlación de J. A. Lasater
Para el cálculo de Rs, esta correlación se basa en 158 mediciones experimentales
de la presión en el punto de burbujeo de 137 sistemas independientes, producidos
en Canadá en el centro y Oeste de los estados Unidos y América del Sur. El error
promedio en la representación algebraica es del 3.8% y el máximo error
encontrado es del 14.7%.
El peso molecular del aceite en el tanque de almacenamiento se correlacionó con

los grados API (de 15 a 40 y 40 a 55).
Procedimiento:
Las siguientes ecuaciones corresponden a la correlación de Lasater para un aceite

saturado:
(1.75)
Donde pf es el factor de la presión en el punto de burbujeo, el cual fue relacionado

con la fracción molecular del gas (ɣg) a cuya curva resultante le fue ajustada la
siguiente ecuación:
(1.76)
La fracción molar del gas se calcula con la siguiente expresión:
(1.77)
39
El peso molecular del aceite en el tanque (Mo) se correlacionó con los °API del
aceite en el tanque de almacenamiento, a cuya curva se le ajuntaron las
siguientes expresiones:
(1.78)
(1.79)
La expresión para determinar Rs se obtuvo a partir de la ecuación 1.77:
(1.80)
A la fracción molar del gas en función de pf, se le ajusto la siguiente ecuación:
(1.81)
: Masa de aceite [lbmol @ c.s. / blo @ c.s.]
: Fracción molar del componente parte gaseosa
: Factor de presión de burbuja [lb/pg2]
Densidad del aceite saturado
La densidad del aceite saturado se calcula con la siguiente expresión:
(1.82)
Viscosidad del aceite saturado
La viscosidad del aceite saturado se puede calcular de la manera siguiente:
40
(1.83)
(1.84)
(1.85)
(1.86)
(1.87)
(1.88)
(1.89)
Tensión superficial del aceite saturado
La tensión superficial del aceite saturado se puede determinar con la siguiente

expresión:
(1.90)
Propiedades del aceite bajosaturado
Compresibilidad del aceite bajosaturado
La ecuación siguiente sirve para determinar la compresibilidad del aceite

bajosaturado:
41
(1.91)
Donde:
Densidad del aceite bajosaturado
La densidad del aceite bajosaturado esta dado por la siguiente expresión:
(1.92)
Viscosidad del aceite bajosaturado
La viscosidad del aceite bajo saturado se obtiene de la manera siguiente:
(1.93)
(1.94)
Donde:
Factor de volumen del aceite bajosaturado
Para el aceite bajosaturado se tiene la ecuación:
42
(1.95)
Correlaciones para obtener la presión en el punto de burbujeo
El primer paso para obtener las propiedades del aceite bajosaturado es la

determinación de la presión de saturación del aceite. Las expresiones a aplicar se
obtienen de las ecuaciones previamente establecidas.
Standing:
(1.96)
Vázquez:
(1.97)
Oistein:
(1.98)
Donde:
(1.99)
Lasater:
(1.100)
Donde pf se obtiene con la ecuación 1.76.
Por lo general la presión de saturación del agua se considera igual a la presión de

saturación del aceite.
43
1.6 Ejercicio de Aplicación
El siguiente ejemplo corresponde al procedimiento que se tiene que efectuar para

seleccionar las correlaciones a emplear en el cálculo de Rs, Bo y Pb con sus
respectivos ajustes y el cálculo de las propiedades del aceite a ciertas condiciones
de presión y temperatura.
Suponiendo que los datos recopilados son los siguientes:
44
Solución:
El primer paso es seleccionar las correlaciones que mejor determinan los valores
de Pb, Rs y Bo con sus respectivos factores de corrección.
Para hacer las correlaciones, se tomarán los datos de laboratorio Rsfb, Bofb y Pb.
Aplicando las correlaciones de Standing, Vázquez, Oistein y Lasater

respectivamente se calculan los valores de pb.
Standing:
Vázquez:
Oistein:
Donde:
45
Lasater:
Se observa que con la correlación de Lasater se ajusta de manera mas

aproximada, el valor de pb.
Para la relación de solubilidad se aplican las siguientes ecuaciones:
Standing:
Vázquez:
Oistein:
46
Lasater:
La correlación de Lasater se ajusto de manera más aproximada al valor de Rs.
Para el factor de volumen de aceite se aplican las ecuaciones:
Standing:
Vázquez:
Oistein:
Con la correlación de Standing se calcula de manera mas aproximada el valor de

Bo
A continuación se muestra una tabla de resultados en los cuales se comparan los

valores obtenidos con cada una de las correlaciones y los valores obtenidos del
laboratorio.
47
CORRELACIONES Rs Bo Pb
Standing 736 1.4723 3352
Vázquez 610 1.3729 3626
Oistein 653 1.4324 2971
Lasater 752 - 3299
Laboratorio 750 1.4550 3300
48
CAPITULO 2 COMPORTAMIENTO DE AFLUENCIA
2. COMPORTAMIENTO DE AFLUENCIA
2.1. Flujo del Yacimiento al Pozo
Uno de los componentes más importantes de un Sistema Integral de Producción

es el yacimiento. En esta área de flujo la pérdida de energía se encuentra en un
rango de 10 a 30 % del total (Beggs, 1991). En consecuencia, el flujo hacia el
pozo depende de la caída de presión en el yacimiento hasta el fondo del pozo, es
decir, la presión del yacimiento menos la presión de fondo fluyendo (Pws-Pwf). La
relación entre el gasto y la caída de presión ocurrida en el medio poroso es muy
compleja y depende de los parámetros tales como propiedades de los fluidos,
propiedades de las rocas, saturación de los fluidos contenidos en la roca, daño a
la formación, turbulencia y mecanismos de empuje.
 Ley de Darcy
En ingeniería petrolera se utiliza con mucha

frecuencia la Ley de Darcy para describir el
comportamiento de flujo (laminar) en el
yacimiento, la que fue establecida por Henry
Darcy (1856), para establecer dicha ley Darcy
realizo un experimento relativamente sencillo , un
recipiente lo llenó de arena e hizo fluir agua a
través del empacamiento hasta saturarlo
completamente.
A partir de esto, Darcy encontró que la velocidad

de un fluido a través de un medio poroso es
proporcional al gradiente de presión e
inversamente proporcional a la viscosidad, μ. De
esta forma, Darcy llegó a establecer la siguiente Fig. 2.1 Experimento de
expresión: Darcy , flujo de agua a
través de un
empacamiento de
arena
(2. 1)
49
Wycoff y Muskat (1936), establecieron que la constante de proporcionalidad de la

Ley de Darcy puede ser expresada en términos de dos propiedades
correspondientes al sistema roca-fluidos, es decir, viscosidad y permeabilidad,
quedando finalmente expresada como se conoce en la actualidad a la Ley de
Darcy:
(2. 2)
O bien en términos de gasto volumétrico:
(2. 3)
: Área total transversal del flujo
: Permeabilidad absoluta del medio poro
: Gasto volumétrico a condiciones de flujo
: Velocidad del fluido [
: Viscosidad del fluido [cp]
: Gradiente de presión en dirección al flujo
El signo negativo asociado a la ec.2.2 y ec.2.3, se debe a que si dx es medido en

la dirección de flujo, la presión “P” declina conforme “x” se incrementa. Esta
declinación dará como resultado un valor negativo para . Por consiguiente, el
signo deberá ser incluido en la ecuación para asegurar que la velocidad V, sea
positiva.
La Ley de Darcy, sólo es válida bajo las siguientes consideraciones:
 Medio homogéneo e isotrópico.

 Medio poroso saturado al 100% por un fluido de viscosidad constante.
 Temperatura constante.
50
 Flujo laminar.
2.2 Ecuación de Afluencia
El comportamiento de afluencia de un pozo representa la capacidad de un pozo

para aportar fluidos. Es decir, el comportamiento de flujo indicará la respuesta de
la formación a un abatimiento de presión en el pozo productor.
Para calcular la caída de presión (abatimiento) en un yacimiento, se requiere una

expresión que muestre las pérdidas de energía o presión, debido al esfuerzo
viscoso o fuerzas de fricción como una función de la velocidad o gasto. Por tanto,
para establecer la ecuación de afluencia para un determinado pozo productor,
será necesario aplicar y combinar las siguientes ecuaciones:
a) Ecuación de conservación de la masa.

b) Ecuación de movimiento.
c) Ecuación de estado.
La Ley de Darcy se debe considerar siempre en la predicción de los gastos de

flujo, desde el yacimiento hasta la cara del pozo. Evinger y Muskat (1943),
establecieron la siguiente ecuación, que puede ser aplicada para predecir
cualquier condición de flujo:
(2. 4)
: Constante, la que en unidades de campo es igual a 0.00708
: Alguna función de presión
: Espesor de la zona productora
: Permeabilidad absoluta
: Presión en la frontera exterior
: Presión de flujo en la pared de pozo
: Gasto de aceite
: Radio de drene
51
: Radio del pozo
La ec.2.4 es aplicable a todos los sistemas porosos, sin embargo, la solución o

forma que adquiera, dependerá de las condiciones iniciales y de frontera (daño,
almacenamiento, fracturas, penetración parcial), así como también de la geometría
y tipo de flujos establecidos en el sistema.
2.3 Geometría de Flujo
En el flujo de fluidos del yacimiento hacia el pozo se han observado diferentes

geometrías de flujo, las que se indican a continuación:
Cilindrica/radial Convergente
Lineal Elíptica
Esférica
Hemisférica
52
De acuerdo con Golan y Whitson, los flujos lineal y radial son los más comunes en
los pozos productores de aceite y gas.
Flujo Lineal
Para flujo lineal, considerando área de flujo constante, la Ec.2.3 será integrada
para obtener la caída de presión ocurrida a lo largo de la longitud L.
(2. 5)
Fig. 2.2 Geometría flujo lineal
Si se considera que k, μ, y q son independientes de la presión, o que puedan ser

evaluadas a la presión promedio en el sistema, la Ec.2.5 quedaría como:
(2. 6)
Cuyo resultado seria.
(2. 7)
O bien:
(2. 8)
Donde C es un factor de conversión. El valor correcto para C es 1.0 para unidades

de Darcy y 1.127 x 10-3 para unidades de campo.
53
Considerando, el hecho de que el gasto másico debe ser constante y expresando

la densidad en términos de la presión, la temperatura y la densidad relativa del
gas, se puede mostrar que la Ec.2.7 se transforma en:
(2. 9)
: Área
: Permeabilidad absoluta
: Longitud
: Presión
: Gasto a condiciones estándar
: Temperatura
: Viscosidad del fluido
Este tipo de flujo se puede subdividir específicamente para flujo lineal para aceite
y para gas y se emplean diferentes ecuaciones como las muestra la tabla
siguiente:
Flujo Lineal para Flujo de Aceite Flujo lineal para Flujo de Gas
Ec.(2.10) Ec.(2.11)
54
: Área abierta al flujo : Área abierta al flujo
:Factor de volumen del aceite :Permeabilidad del gas
:Gasto de gas a 14.7 y 60°F

:Permeabilidad del aceite
:Presión corriente arriba

: Temperatura del flujo
:Presión corriente abajo :Factor de compresibilidad del gas
:Gasto de aceite :Coeficiente de velocidad
:Coeficiente de velocidad :Densidad relativa del gas
:Viscosidad del aceite :Viscosidad del gas a T y P ]
: Densidad del aceite
Ec.(2.12)
Tipo de a b
formación
Consolidada 2.329 x 1010 1.2
No 1.470 x 107 0.55

consolidada
Flujo Radial
La Ley de Darcy puede ser usada para calcular el flujo hacia el pozo, donde el
fluido converge radialmente. En este caso, el área abierta al flujo no es constante,
por lo tanto, deberá ser incluida en la integración de la Ec.2.3, Haciendo referencia
a la geometría de flujo ilustrada en la ilustración, el área de la sección transversal
abierta al flujo para cualquier radio sería A = 2 r h.
55
Fig. 2.3 Sistema de Flujo Radial
Del mismo modo, definiendo como negativo el cambio en la presión con respecto
a la dirección de flujo, - . Haciendo estas sustituciones en la Ec.2.3 se obtiene:
(2. 13)
Flujo radial para Aceite Flujo radial para Gas
56
2.4 Regímenes de Flujo
Adicionalmente a los patrones de flujo descritos, Slider (1983), León (1984),

Muskat (1937), clasificaron los diferentes tipos de flujo que se presentan en el
medio poroso, de acuerdo a la dimensión, geometría, tiempo y fase. Tal
clasificación es la siguiente:
REGIMENES DE FLUJO
Dimensión Geometría Tiempo Fase
Estacionario
Unidimensional Lineal Monofásico
Pseudoestacionario (una sola fase)
Bidimensional Radial
Transitorio Multifásico
(bifásico, trifásico)
Tridimensional Esférico
57
Fig. 2.3 Grafica del periodo de flujo estacionario
Fig. 2.4 Grafico del periodo de flujo variable
58
Fig. 2.5 Grafico del periodo de flujo pseudo-

estacionario
2.5 Análisis del Pozo Fluyente
Para llevar a cabo el análisis de un pozo fluyente, es necesario cubrir dos

aspectos fundamentales: en primer lugar tener una concepción muy clara del
mecanismo de flujo que siguen los fluidos producidos, desde la frontera de drene
del yacimiento, hasta la central de recolección o batería de separadores. En
segundo término, disponer de la metodología y herramientas de cálculo, que
permiten predecir el comportamiento del sistema en general. El sistema integral
del flujo está constituido por cuatro partes principales que son:
 Flujo en el yacimiento.
 Flujo en el pozo, a través de tuberías verticales o inclinadas.
 Flujo en el estrangulador.
 Flujo en la línea de descarga.
Flujo en el Yacimiento
Se muestran las curvas típicas que pueden representan el comportamiento de flujo

en el yacimiento de un pozo a través del tiempo de explotación. En la línea A, la
tendencia es una recta que se presenta cuando la presión de fondo fluyendo es
mayor a la presión de saturación. A presiones de fondo fluyendo menores a Pb el
59
comportamiento observa la tendencia de la línea B. Al de presionarse el

yacimiento, puede esperarse un comportamiento como el de las líneas C y D.
Fig. 2.6 Curvas típicas del comportamiento de afluencia del yacimiento al pozo
Cuando la presión de fondo fluyendo es mayor a la presión de saturación la

pendiente “m” de la recta (línea A) es constante y entonces:
(2. 26)
Cuando Pwf < Pb se considera un comportamiento no lineal al que se conoce

como IPR (curvas B, C y D). El gasto teórico que se obtendría de un pozo cuando
Pwf = 0 se conoce como qmax.
IP en Yacimientos Bajosaturados
Suponiendo un Índice de productividad constante, independientemente de la

producción a condiciones superficiales y con producción de aceite y agua, se
puede emplear la siguiente ecuación:
(2. 27)
60
O bien, considerando flujo radial para un yacimiento homogéneo, horizontal,

uniforme y de poca compresibilidad, la ecuación de Darcy:
(2. 28)
: Es el factor total de daño a la formación, el que puede ser determinado

mediante pruebas de presión en los pozos
: Es el término por flujo turbulento, generalmente despreciado cuando se está

produciendo a gastos bajos y para formaciones de baja permeabilidad.
IPR (Inflow Performance Relationship) en yacimientos saturados
Cuando existe flujo en dos fases en el yacimiento la relación de la Ec.2.26 no se

cumple, pues el valor de la pendiente cambia continuamente en función del
abatimiento en la presión
Se concluye que el Fig.2.7

IP paraVariación
cualquierdel IP para
gasto yacimientossiempre
de producción, saturados
que Pwf > Pb,
será la primera derivada del gasto con respecto al abatimiento de presión, esto es:
IP = IPR = dq / dPwf
Para cálculos de IPR en yacimientos saturados se tiene los siguientes métodos de

cálculo:
61
Método de Vogel
(2. 29)
: Gasto de flujo total o líquido [bpd]

: Gasto de aceite [bpd]
: Gasto de Agua [bpd]
: Permeabilidad relativa del aceite [Darcy]
: Permeabilidad relativa del agua [Darcy]
: Viscosidad del aceite [cp]
: Viscosidad del agua [cp]
: Factor de volumen del agua
: Factor de volumen del aceite

: Espesor de la arena productora [pies]
: Radio de drenaje del pozo [pies]
: Radio del pozo [pies]
: Presión estática del yacimiento
: Presión de fondo fluyente
62
Fig. 2.8 Curva de afluencia para pozos sin daño de un

yacimiento con empuje por gas disuelto
Para yacimientos saturados (Pwf < Pb) Vogel ofreció una solución al problema de
flujo bifásico en el yacimiento desarrollando la siguiente ecuación:
(2. 30)
Donde:
: Gasto de líquido [bpd]

: Tasa máxima de líquido [bpd]
: Presión estática o del yacimiento en la vecindad del pozo
: Presión De Fondo Fluyendo
63
Método de Standing
Standing extendió el trabajo de Vogel para tomar en cuenta cuando la eficiencia

de flujo (EF) era distinta a 1. Se define eficiencia de flujo a la relación existente
entre el índice de productividad real y el ideal:
(2. 31)
Donde:
: Presión de fondo fluyendo ideal “libre de daño”

: Presión de fondo fluyendo real o actual
: Presión estática del yacimiento
Fig. 2.9 Curvasde afluencia para pozos con EF

diferente a 1 de yacimientos con empuje por gas
disuelto.
Standing construyó curvas IPR (fig. 2.9) para diferentes EF, así la ecuación de
Vogel se utiliza directamente:
64
(2. 32)
Tomando en cuenta que:
(2. 33)
Donde qmax viene siendo la máxima tasa de producción que se tiene del pozo si
S=0 y EF=1. Standing finalmente llegó a las siguientes ecuaciones:
Para EF < 1:
(2. 34)
Para EF > 1:
(2. 35)
Método de Fetkovich
M.J. Fetkovich en 1973 partió de la ecuación básica de flujo en estado estable

desarrollada por Evinger y Muskat y bajo una serie de artificios y suposiciones,
llegó a la conclusión de que el comportamiento de afluencia de pozos productores
de yacimientos de gas o de petróleo puede determinarse mediante la siguiente
expresión:
(2. 36)
Donde:
: Constante del Flujo

: Presión promedio estática del yacimiento
: Exponente con valor promedio entre 0.5-1

: Presión de fondo fluyendo
65
: Gasto de producción [bpd]
2.6 Curvas IPR Futuras
Del estudio del comportamiento de flujo del yacimiento al pozo y en relación al

gasto, se distinguen dos problemas a resolver. El primero es la forma de la curva
de presión contra gasto de una etapa de explotación en particular. El segundo es
la manera en que el IPR decrece a medida que continúa la explotación.
La determinación de las IPR futuras reside en el pronóstico del gasto del pozo,
tanto para pozos fluyentes como para aquellos que tienen instalado un sistema
artificial de producción. También son relevantes cuando se desea obtener el
memento más conveniente para la conversión a un sistema artificial de
producción.
2.6.1 Método de Fetkovich
Muskat en 1957, presento la ecuación 2.37, que relaciona el índice de

productividad de un tiempo t1 a otro t2
(2. 37)
De un balance de materia para un yacimiento con empuje con gas disuelto,

Fetkovich, encontró el comportamiento de kro es casi lineal con la presión, y se
puede calcular aproximadamente con la siguiente expresión:
(2. 38)
O bien
(2. 39)
Donde kro se toma con respecto a ki, se define en aquel lugar geométrico donde no
existe abatimiento.
De esta forma:
(2. 40)
66
Que define el gasto contra abatimiento de presión; por otra parte la ec.2.39
proporciona la corrección por declinación de la presión. Finalmente la ecuación
propuesta por Fetkovich es:
(2. 41)
Realizando una prueba de flujo de tres o cuatro puntos para el tiempo presente, es
posible construir las curvas IPR para otras presiones estáticas utilizando la ec.
2.41.
2.6.2 Método de Eickemer
Para n = 1 y pwf = 0, el resultado del cociente de la ec. 2.40 a t1, entre la ec. 2.41 a
t2 permite obtener el qmax a t2.
(2. 43)
La solución a la ec. 2.43 requiere una prueba de pozo al tiempo presente t 1, a

partir de la cual se puede obtener qmax1. Conociendo qmax1 y pws1, se puede
resolver a qmax2 a cualquier otra presión estática considerada y, como
complemento, puede construirse una curva completa de IPR para t2 empleando la
ecuación propuesta por Vogel (2.30).
2.6.3 Método de Standing
Standing desarrollo un procedimeinto que requiere cálculos de balance de materia

para determinar saturaciones futuras, a partir se pueden obtener valores de kro. De
su estudio, la ec. 2.40 toma la siguiente forma:
(2. 44)
Y de la ec. 2.27 se obtiene:
67
(2. 45)
Si el valor J se toma para abatimientos bajos, es decir, cuando la

ec. 2.45 se convierte en:
(2. 46)
Dividiendo la ec.2.45 entre la ec. 2.46, se tiene que:
(2. 47)
Que permite calcular a partir de un valor medido de J. de lo anterior se infiere

*
que un valor futuro de J ( ) se puede calcular a partir del valor presente para
cualquier pws futura. Este valor se determina con la siguiente ecuación:
(2. 48)
Finalmente si q, se elimina de la ec. 2.30 y 2.47, la curva del IPR futura se puede
construir a partir de la ec. 2.49.
(2. 49)
Opara la aplicación de la ecu. 2.48, es necesario contar con los valores presentes
y futuros de kro, µo y Bo pueden usarse correlaciones numéricas (capitulo 1), para
obtener la kro puede obtenerse de la correlación de Corey.
(2. 50)
68
Donde So y Sw representan las saturaciones del agua y del aceite

respectivamente.
Ejemplo 2.1
Un pozo drena un área circular de 80 acres (re=932 pies) de un yacimiento de

aceite bajo saturado. Si la presión del yacimiento (Pe) es 1000 [lb/pg 2], la
permeabilidad es 50 [mD], el espesor neto de la arena es 20 pies, la viscosidad del
aceite es 3 [cp], el factor de volumen del aceite es 1.25 y el pozo está terminado y
disparado en una TR de 7 [pg] de diámetro:
a) ¿Cuál es el índice de productividad del pozo?
b) ¿Cuál será el gasto de producción del pozo para una presión de fondo fluyendo
de 500 [lb/pg2]?.
Solución:
Datos:
re  932 pies   o  3[cp ]

Py  Pb  bl @ c. y 
Bo  1.25  
 lb   bl @ c.s 
Py  1000  2   Pws
 pg 
 TR  7 [ pg ]
k o  50 [mD] rw  3.25 [ pg ]  0.2917 [ pies ]
h  20  pies   lb 
Pwf  500  2 
 pg 
0.00708 K o h
J
r 
μo Bo ln e 
 rw 
(0.00708)(50)(20) 7.08  bl / día 

a) J    0.234 
 932  30.2601  lb / pg 2 
(3)(1.25) Ln 
0.2917 
b) q  J ( Pws  Pwf )  0.234 (1000  500)  117 [ BPD ]
69
Ejercicio 2.2
Un pozo fluyente con 300 pies de tubería de producción tiene una presión de
fondo fluyendo de 580 lb/pg2 cuando la producción es de 42 bl/día, y de 328 lb/pg2
cuando se producen 66 bl/día. Si su presión de burbuja es de 225 lb/pg2 ;
Determinar:
a) El índice de productividad del pozo.
b) Su presión estática
c) Su gasto máximo o potencial.
Solución:
Datos:
h  300 [ pies ]
 lb 
Pwf1  580  2  qo 1 42 [ BPD ]
 pg 
 lb 
Pwf 2  328  2  qo2  66 [ BPD ]
 pg 
Como pwf > pb J es lineal por lo que:
1
J 
m
pwf 2  pwf 1
m
qo2  qo1
q
J
Pws  Pwf
 
328  580 1  BPD 
a) m  10.5 J   0.0952  
66  42  10.5  lb 
 pg 2 
q 66  lb 
b) Pws  ( )  Pwf   328  1021.27  2 
J 0.0952  pg 
c) qo máx  J Pws  (0.0952)(1021.27)  97.22 BPD 
70
Ejemplo 2.3
Un pozo fluyente, al que se realizaron pruebas de incremento y decremento de

presión, así como un análisis PVT obtuvo la siguiente información:
2
Presión [lb/pg ] Bo Qo [BPD]
6280.17 1.428 0
3519.99 1.472
2665.03 1.488
2016.73 1.408
1516.78 1.349
1015.69 1.285 200
 bl / día 
Considerando un índice de productividad lineal de 0.098  2 
y una
 lb / pg 
pws =2016.73 lb / pg 2 ; Determine:
a) El gasto máximo o potencial del pozo
b) Elabore la curva IPR de dicho pozo para: p>pb y p<pb
Solución:
a)
Del comportamiento del factor volumétrico del aceite a la presión de burbuja B o es

máximo, por lo que el pozo está en un yacimiento saturado (p<pb).
Por el método de Voguel:
71
= = 287.20 barriles
b)
Para pwf > pb el comportamiento es lineal
 bl / día 
Con J= 0.098  2 
y suponiendo presiones por debajo de la presión inicial.
 lb / pg 
2
pwf [lb/pg ] Qo=J(pws-pwf) [barriles]
6000 27.45
4500 174.45
3000 321.45
2700 350
Curva IPR para pwf > pb
7000
6000
5000
pwf [lb/pg2]
4000
3000
2000
1000
0
0 100 200 300 400
qo [barriles]
 P  P  
2
qo  qomáx 1  0.2 wf   0.8 wf   [barriles]

2
pwf [lb/pg ]   Pws   Pws  
2016.73 0
1516.78 114.0351
1200 171.6749
72
1015.69 199.9938
Curva IPR para pwf < pb
250
200
pwf [lb/pg2]
150
100
50
0
0 500 1000 1500 2000 2500
qo [barriles]
Ejemplo 2.4
Se tiene un pozo de aceite con las siguientes caracteristicas:
Pws=4200 lb/pg2
J= 2.0
Pb= 3000 lb/pg2
Calcular:
a) qb
b) qmax
c) q para pwf=1500 lb/pg2
Solución:
1) qb= j(pws-pb)=2.0(4200-3000)=2400 B/D

2) qmax= qb+ =2400+ =5733.33 B/D
3) qo= qb+
=2400+ (5733.33-2400)
=4733.331 B/D
73
Ejemplo 2.5
Se hace pasar un fluido incompresible de manera horizontal a gasto constante de

0.5cm3/seg a través de un medio homogéneo e isotrópico cuya área transversal
es de 4.5 cm2, generando una caída de presión de 4 atmosferas. Si dicho fluido
tiene una viscosidad de 2cp y atraviesa una sección de 15cm determine la
permeabilidad en milidarcys.
Solución:
De la ecuación de Darcy para flujo lineal despejamos para la permeabilidad.
Ejemplo 2.6
Para un pozo de aceite con las siguientes características:
k= 120mD
h= 30 pies
µo= 0.8 cp
Bo=1.15
Determine su índice de productividad:
Solución:
De la ecu.
74
Ejemplo 2.7
Pozo de aceite:
k= 20 mD
h= 40 pies
Pws= 2600 [lb/pg2]
Haga una evaluación del pozo.
Solución:
J= kh= (20)(40)=800 md-pie=0.8 Darcy-pie
Para J= cte
qmax=0.8(2600-0)=2080 B/D
Por Voguel
qmax=
Este es un buen pozo, del potencial obtenido por Voguel de 1155 B/D, que es un
valor mínimo dependiendo de la presión de burbuja.
Ejemplo 2.8
Pozo de aceite:
k= 5 mD
h= 10 pies
Pws= 1000 [lb/pg2]
Solución:
J=(5)(10)=50 mD-pie
J=
qmax= 0.05(1000-0)= 50 B/D
75
Este no es un buen pozo y probablemente la presión del yacimiento se encuentre

debajo de su presión de burbuja, si calculamos qmax por el método de Voguel,
50/1.8= 28 B/D. La producción por un sistema de levantamiento artificial
probablemente sería de 15 a 20 B/D.
Ejemplo 2.9
Para un pozo de aceite con las siguientes características:
k= 300mD
h= 40 pies
µo= 800 cp
Bo=1.15
Pws= 3000 [lb/pg2]
Solución:
kh= (300)(40) =12000 mD-pies
J=12
kh> 5000 indicando un excelente pozo. Sin embargo, la alta viscosidad reduce el
índice de productividad como a continuación se muestra:
J=
qmax=0.013(3000-0)= 39 B/D
Para pozos que producen aceite pesado la productividad es reducida

drásticamente por la alta viscosidad.
Ejemplo 2.10
Para un pozo de aceite con los siguientes datos:
Pws= 2600 [lb/pg2 a]
qo= 500 B/D para pwf=1800 a
E.F= 0.6
76
Calcular
a) qmax para E.F=1

b) qmax para E.F=0.6
c) encontrar qo para pwf=1300 y E.F=0.6
Solución:
a)
E.F=
Encontramos Pwf para cuando E.F=1
= 2600-(2600-0)0.6 = 1040 lb/pg2
= 2600-(2600-0)0.6 = 1040 lb/pg2
b)
Pws=2600
P´wf=2600-(2600-0)(0.6)=1040 lb/pg2
c)
Pws=2600
P´wf=2600-(2600-1300)(1.3)=1298 lb/pg2
77
Ejemplo 2.11
Para Pws = 408.2 lb/pg2 analizar los datos de la prueba de potencial mostrados en
la tabla y determinar la ecuación de afluencia particular del pozo de gas y el gasto
máximo de este en MPCD:
Prueba Pwf [lb/pg2 ] q g [MPCD]
1 403.1 4.288
2 394.0 9.265
3 378.5 15.552
4 362.6 20.177
Patm 14.7 Qgmáx
Solución:
qg  C ( Pws2  Pwf2 ) n
Pwf [lb/pg2] qg [MPCD] PWs2 - PWf2 [lb/pg2]
408.2 ----------
403.1 4.288 4137.63
394.0 9.265 11391.24
78
378.5 15.552 23364.99
362.6 20.177 35148.48
14.7 Pot. Abs. 166411.15
Determinamos las constantes C y n
log( Pws2  Pwf2 ) 2  log( Pws2  Pwf2 )1

m
log( q g ) 2  log( q g )1
log( 35148.48) 2  log( 4137.63)1

m
log( 20.177) 2  log( 4.288)1
0.9292
m  1.3815
0.6726
=0.010329
qg  C ( Pws2  Pwf2 ) n  0.010329( Pws2  Pwf2 ) 0.7239
79
Ejercicios propuestos
2. A
Un pozo drena un área circular de un yacimiento de aceite bajo saturado

(re=575 mts). Si la permeabilidad es 10-4[m2], el espesor neto de la arena es 10
mts, la viscosidad del aceite es 2X10-3 [Pa-s], el factor de volumen del aceite es
1.1 y el pozo está terminado y disparado en una TR de 0.2 [mts] de diámetro.
Determine la caída de presión (pws-pwf) en pascales. Considere S=0
Solución: (pws-pwf) = 3.043x107
2. B
Desarrolle una curva IPR y calcule el gasto máximo para el ejercicio 2.2,
considerando una pb= 1200 lb/pg2.
Solución: qomáx= 66.83 barriles, para una pwf = 582 lb/pg2 y qo= 42 barriles.
1200
1000
800
pwf [lb/pg2]
600
400
200
0
0 20 40 60 80
qo [barriles]
80
CAPITULO 3 FACTOR DE DAÑO Y SU RELACIÓN CON
COMPORTAMIENTO DE AFLUENCIA
3. FACTOR DE DAÑO Y SU RELACIÓN CON COMPORTAMIENTO DE

AFLUENCIA
3.1 Factor de Daño
Durante la perforación, terminación o producción de un pozo, es posible que una

zona de permeabilidad alterada pueda desarrollarse alrededor de las paredes del
pozo. La zona con la permeabilidad alterada es llamada “zona dañada” y su efecto
sobre la presión o comportamiento de flujo del pozo, es denominado como efecto
de daño.
Asimismo, el daño se define como un factor que causa, en o alrededor del pozo,
una caída de presión adicional a la que ocurre cuando el yacimiento es
homogéneo y el pozo penetra totalmente a la formación. Si el lodo de perforación
provoca reducción a la permeabilidad, o si el pozo penetra parcialmente a la
formación o si existe flujo No - Darciano (pozo de gas), se tiene un daño. El factor
de daño (en su definición tradicional) representa una caída de presión adicional, la
que ocurre en las cercanías o en la entrada al pozo (Van Everdingen, 1953).
El factor de daño es una medida cuantitativa empleada para evaluar el

comportamiento de un pozo, relativa a la producción ideal de un pozo a partir de
una formación completamente abierta y sin restricciones. La magnitud del daño
indica la necesidad de estimular un pozo, o bien establecer un programa de
reacondicionamiento del pozo.
El concepto de daño fue propuesto originalmente por Hurst y Van Everdingen

(1953). Estos investigadores propusieron el factor de daño como una forma de
cuantificar el flujo no – ideal, señalaron que las presiones medidas en un pozo
frecuentemente no se ajustaban a las soluciones teóricas. Propusieron que la
diferencia era una caída de presión adicional causada por restricciones al flujo
cercanas al pozo.
En relación al IPR de un pozo, el factor de daño considera la diferencia entre el

abatimiento de presión ideal y el abatimiento de presión actual o real.
Generalmente, sólo interesa el factor de daño durante el periodo pseudo-

estacionario y se desprecia el efecto de daño en el periodo de transición debido a
los tiempos cortos.
Para condiciones de periodo pseudo-estacionario de un pozo, el factor de daño es

incluido en el cálculo de la caída de presión total, esto es, p y – pwf. Lo cual se
puede establecer de la siguiente manera:
81
(3. 1)
Donde:
: Presión del yacimiento
: Presión de fondo fluyendo a condiciones reales
: Presión de fondo fluyendo considerando un caso ideal
La diferencia de presión de fondo fluyendo entre la ideal y la actual, P’wf - Pwf,

representa una pérdida de presión adicional debido al daño de la formación,
estimulación de la zona alrededor del pozo y otras restricciones al flujo a la
entrada del pozo.
Usualmente se representa la caída de presión adicional como . El factor de

daño adimensional S, proporcional a , es definido como:
(3. 2)
O bien:
(3. 3)
El daño bien es la alteración, natural o inducida de las propiedades petrofísicas de

la roca en la vecindad del pozo que caracterizan el flujo de fluidos a través de la
misma.
82
Fig. 3.1 Distribución de presión de un pozo con una zona

alterada en la cercanía del pozo (Golan y Whitson, 1991)
83
Ahora bien, considerando flujo pseudoestacionario para un caso ideal, se tiene

que:
(3. 4)
Para expresar la caída de presión actual o real en términos de factor de

daño se tiene que:
(3. 5)
Y para el gasto nos da como resultado la ec:
(3. 6)
Para cuantificar el efecto actual del factor de daño, es necesario calcular la

constante ( ) S, la que proporciona las pérdidas de presión debido al
daño.
En algunas ocasiones resulta más útil expresar el efecto de daño y estimulación

en términos de la eficiencia de flujo, EF.
La eficiencia de flujo EF, se define como la relación entre el gasto actual y el gasto
ideal para un abatimiento de presión dado.
(3. 7)
En términos de la presión de fondo fluyendo ideal y actual, requeridas para

producir a un gasto determinado, la eficiencia de flujo tiene la siguiente forma:
(3. 8)
84
En términos del factor de daño, se puede mostrar fácilmente que la eficiencia de

flujo esta expresada como:
(3. 9)
Para la mayoría de pozos el término tiene un rango de valores entre 6.5 y
8.5. Utilizando un promedio de , se puede escribir una expresión

aproximada para la eficiencia de flujo en términos del factor de daño:
(3. 10)
La cual se puede utilizar como una regla de dedo.
En la industria se reportan a menudo condiciones no ideales en términos de una

relación de daño la que no es más que el recíproco de la eficiencia de flujo.
(3. 11)
Otra ecuación o expresión utilizada para cuantificar las condiciones de flujo ideales
el radio aparente del pozo
(3. 12)
El que es sustituido simplemente en la ecuación 76 por el radio del pozo.
EFECTO CUANTIFICABLE DEL FLUJO NO LINEAL
CONDICON DEL
POZO
85
DAÑADO
INALTERADO(SIN
DAÑO)
ESTIMULADO
Por otra parte, el factor de daño considerado hasta este momento se puede tratar
desde otro punto de vista. Hawkins (1956) estableció la siguiente forma para
calcular
(3. 13)
Multiplicando y dividiendo la Ec.83 por 2.303 para cambiar el valor logarítmico

obtenemos que;
(3. 14)
Donde:
(3. 15)
Finalmente obtenemos:
(3. 16)
86
El valor del factor de daño, usualmente se calcula a partir del análisis de datos de
pruebas de incremento o decremento de presión.
En relación a las mediciones del daño, es importante reconocer que el factor de

daño calculado a partir del análisis de una prueba de incremento de presión, por
ejemplo, un gráfico de Horner (1951), representa el efecto de daño experimentado
por el pozo en el instante en que se cierra.
El efecto de daño total S, para un pozo se puede establecer de la siguiente forma:
(3. 17)
3.2 Flujo Restringido
El daño a la formación o flujo restringido, está asociado a diferentes factores que

reducen la permeabilidad efectiva alrededor del pozo. El flujo restringido puede ser
causado por el taponamiento de los espacios porosos mediante partículas sólidas
generadas por trituramiento mecánico de la formación durante la perforación del
pozo o disgregación del medio poroso, o bien, por efecto del fluido producido tales
como creación de emulsiones o cambios en la permeabilidad efectiva.
El taponamiento de los poros mediante partículas es uno de los mecanismos más

comunes en el flujo restringido, el que puede resultar de diferentes causas,
incluyendo la inyección de sólidos en la formación (para fracturamiento de la
formación), dispersión de arcillas presentes en la formación, precipitación y
crecimiento de bacterias.
Mc Leod (1984) establece que las causas que originan el flujo restringido pueden
clasificarse en tres categorías:
 Físicas.
 Químicas.
 Biológicas.
87
Causas Físicas
Un medio poroso es un complejo ensamble de granos de mineral con espacios

vacíos (poros) de forma y distribución irregular, que proporciona una trayectoria
para el transporte del fluido. Esta complicada estructura puede ser idealizada
como una colección de largas cámaras (los cuerpos de los poros) conectados por
estrechas aberturas. La permeabilidad del medio está controlada ampliamente por
el número y conductividad de las cámaras de los poros.
Los fluidos producidos de la formación hacia el pozo, contienen diminutas

partículas o finos, que se pueden adherir a los cuerpos de los poros, o bien,
aglutinarse de tal forma que pueden obturar las aperturas de los poros.
Fig. 3.2 Modelo de entrampamiento de

la partícula
Calculo del pseudo- Daño por el efecto de los disparos
Karakas y Tarquiq (1988) presentaron una solución semianalítica para el cálculo

del pseudodaño debido a los disparos, el cual dividieron en 3 componentes: El
efecto en el plano de flujo (Sh), el efecto de la convergencia vertical (Sv) y el efecto
debido al pozo Swb. El pseudodaño total causado por los disparos está dado por:
Sp = Sh + Swb + Sv (3. 18)
88
La fig.3.3 muestra los componentes para el cálculo del pseudodaño. Esta incluye
el radio del pozo rw, el radio de la perforación rperf, la longitud de la perforación lperf,
el ángulo de la fase de la perforación, la distancia entre cada perforación h perf.
Fig. 3.3 Daños por disparos
El procedimiento para el cálculo de cada componente del pseudodaño se muestra

a continuación:
lp (3. 19)
rw´ ( )  ; si θ0
4
(3. 20)
rw´ ( )  αθ rw  lp ; si θ0
 r  (3. 21)
Sh  ln ´ w 
 rw ( ) 
89
hk h /k v  (3. 22)

0.5
hD 
lp
rpD  rperf
1 k /k  
v h
0.5 (3. 23)
2h
(3. 24)
a  a1logrpD   a2
(3. 25)
b  b1rpD  b2
(3. 26)
Sv  10a hDb-1 rpD
b
rw (3. 27)
rwD 
lp  rw
(3. 28)
Swb  c1 ec 2rwD
Procedimiento de cálculo:
1. de la tabla 1 - 2 : α θ
lp
2. rw´ ( )  ; si θ0
4
rw´ ( )  α θ rw  lp  ; si θ0
90
 r 
3. S h  ln ´ w 
 rw ( ) 
hk h /k v 
0.5
4. hD 
lp
5. rpD  rperf
1  k /kv h 0.5 
2h
6. de la tabla 1 - 3 : a1, a 2 , b1 y b 2
7. a  a1logrpD   a 2
b  b1rpD  b 2
8. Sv  10 a hDb-1 rpD
b
9. de la tabla 1 - 4 : c 1 y c 2
rw
rwD  S wb  c 1 e c 2rwD
l p  rw
10. Sp = Sh + Swb + Sv
Pseudodaño debido al ángulo de desviación del pozo
Cinco-Ley et al. (1975) presento una solución semianalítica para determinar el

valor del pseudodaño causado por la terminación parcial, desviación del pozo y
elevación de la terminación ( Sθc ). La fig.3.4 muestra las variables relevantes para
dicho cálculo. Donde hw es el espesor perforado, Zw es la elevación del punto
medio del espesor perforado a la base del yacimiento, h el espesor del yacimiento,
θ es el ángulo de desviación del pozo con respecto a la vertical y r w el radio del
pozo.
91
Fig. 3.4 Daño debido al ángulo de

desviación del pozo (Economide,1994)
Las siguientes ecuaciones son utilizadas para el cálculo del pseudodaño debido a
la terminación:
h (3. 29)
hD 
rw
zw (3. 30)
z wD 
rw
hw (3. 31)
h wD 
rw
hwDcosθ (3. 32)

 h'wD cosθ
hD
92
Procedimiento de cálculo:
Cálculo de Sθ c Sθ  Sc 
h
1. hD 
rw
zw
2. z wD 
rw
z wD
3.
hD
hw
4. h wD 
rw
h wD cosθ
5.  h' wD cosθ
hD
6. de la tabla 1 - 1 : S θc
Causas Químicas
Los finos, responsables del obturamiento de los poros pueden provenir de

orígenes externos u originarse en el mismo espacio poroso. Los finos en el medio
poroso pueden ser movilizados por un cambio en la composición química del agua
congénita o mecánicamente debido a las fuerzas cortantes aplicadas al
movimiento del fluido. El flujo restringido, a menudo es causado por la dispersión
de partículas finas de arcilla cuando la salinidad del agua intersticial es reducida o
la composición iónica es alterada. De esta forma, cualquier fluido que esté en
contacto con la formación productora (fluido de perforación, fluidos de terminación,
fluidos de estimulación), deberá tener una composición iónica que no dañe la
formación.
La precipitación de sólidos de la salmuera o lodo de perforación y del aceite en la

formación puede causar severo daño a la formación cuando estos sólidos obturan
los espacios porosos. Los precipitados pueden ser compuestos inorgánicos de la
salmuera o especies orgánicas del aceite producido. En cualquier caso, la
precipitación ocurre por cambios en la temperatura o presión en las cercanías del
pozo, o a partir de alteraciones en la composición de la fase de los fluidos
inyectados.
93
Otra causa común que origina el flujo restringido son las parafinas y los asfaltenos.
Las parafinas son cadenas largas de hidrocarburos que se precipitan de ciertos
hidrocarburos cuando la temperatura se reduce, o la composición del aceite
cambia debido a la liberación de gas conforme la presión se reduce. Los
asfaltenos son compuestos aromáticos con alto peso molecular que pueden tener
forma coloidal, dispersos en el aceite. Este estado coloidal es estabilizado por la
presencia de resinas en el aceite; cuando estas resinas son removidas, los
asfaltenos pueden flocular, creando partículas sumamente grandes que causan el
flujo restringido. El cambio químico en el aceite puede reducir la concentración de
resinas y de esta forma se lleva a cabo el depósito de los asfaltenos en la
formación.
Causas Biológicas
Las bacterias inyectadas en la formación, básicamente bacterias anaeróbicas,

pueden crecer rápidamente en la formación, obturando los espacios de los poros,
o bien, con precipitados resultado de la actividad biológica de los organismos. La
reducción de la permeabilidad causada por la bacteria puede ser significante, por
tanto, se debe considerar la posibilidad de inyección de agua con bactericidas
para reducir al mínimo este problema.
Por otra parte, mientras que el factor de daño es adimensional, la zona dañada
asociada no lo es. La figura 3.5 es una representación típica de la condición en las
cercanías del pozo, con y , siendo la profundidad de daño y la permeabilidad
alterada respectivamente. Fuera de esta zona el yacimiento permanece sin
disturbio con permeabilidad k.
Fig. 3.5 Zona alrededor del pozo con

permeabilidad alterada (Economide,1994)
94
Hawkins (1956) establece una expresión que relaciona el efecto de daño con las
variables y .
Si la permeabilidad en la cercanía del pozo es la permeabilidad del yacimiento (sin

daño), entonces la caída de presión entre la presión en la frontera externa ( ) y la
presión en el pozo, dará como resultado una presión ideal, , expresada
como:
(3. 33)
Si la permeabilidad en la cercanía del pozo es alterada a , entonces la presión

de fondo fluyendo real está relacionado como:
(3. 34)
La diferencia entre y es exactamente la caída de presión debido

al efecto de daño Por lo tanto se tiene que:
(3. 35)
Simplificando dicha ecuación se tiene que:
(3. 36)
Un problema inherente con el uso práctico de la Ec. 3.33, es que la zona alterada
es muy difícil de cuantificar en términos de y . Esto es, la permeabilidad
alterada y el radio no pueden ser medidos directamente. En el mejor de los casos,
si el daño de la zona alterada puede ser aislado del daño total (calculado a partir
de análisis de pruebas de incremento o decremento de presión), una estimación
del radio de la zona alterada permitirá el cálculo de la permeabilidad de la zona
alterada:
95
(3. 37)
Resolviendo para el radio de la zona alterada se obtiene:
(3. 38)
Simplificando las expresiones para una situación típica de un pozo estimulado

cuando >>> k:
(3. 39)
(3. 40)
Ambas ecuaciones permitirán saber la situación o estado actual el pozo en

producción.
3.3 Factor de Daño Compuesto
El factor de daño total o compuesto, S, que se evalúa a través del análisis de las
pruebas de incremento y decremento de presión; ha sido visualizado en términos
de permeabilidad reducida o mejorada. Sin embargo, en los últimos años, este
concepto ha sido extendido para incluir una variedad de efectos tales como,
penetración parcial, desviación del pozo, perforaciones, etc.
El efecto total de todas las condiciones no ideales cerca del pozo, es una
composición de todos los efectos individuales. Por esta razón, el factor de daño
total ha sido definido incluyendo todos estos efectos de la siguiente manera:
(3. 41)
96
Donde:
: Factor de daño compuesto
: Factor de daño por flujo restringido
: Factor de daño asociado a alta velocidad de flujo
: Factor de daño causado por disparos
: Factor de daño causado por fracturas
: Factor de daño por empacamiento de grava
: Factor de daño debido a pozos horizontales
: Factor de daño por formación de condensados
3.4 Diagnostico del Daño
Existen diversas causas por las cuales un pozo no aporta la producción que se
espera, o declina su producción con el tiempo. Por ello es importante realizar un
estudio, ya que puede atribuirse la baja productividad de un pozo de un daño en el
mismo.
3.4.1 Determinación del daño a la formación
Si un pozo tiene daño es necesario evaluarlo y diagnosticarlo para así reducir o

minimizar sus efectos nocivos en la producción. Existen varios métodos para
diagnosticar el daño a la formación entre ellos están:
•Pruebas de restauración de presión con tubería de perforación
•Registros de resistividad
•Revisión histórica de perforación, terminación, y reparación del pozo
•Análisis de estimulaciones previas
•Comparación con pozos vecinos
•Análisis de pruebas de presión
97
•Análisis nodal
•Registros de producción
•Pruebas y análisis de núcleos
Prueba DST
Una prueba DST puede dar indicación de la presencia de daño cuando restaura
rápidamente la presión en un periodo de cierre, y hay gran diferencia entre la
presión de flujo inicial y la final en poco tiempo, indicando que hay transmisibilidad
de la presión pero muy poca al flujo, lo cual puede deberse a la obstrucción del
gasto
Registros de resistividad.
La combinación de un registro dual de inducción y el latero log pueden dar idea de

la profundidad de la invasión de los fluidos hacia la formación. Tomando en cuenta
el registro de calibre del agujero para hallar el espesor de revoque en las zonas
permeables.
Revisión histórica de perforación, terminación, y reparación del pozo
Son un paso crítico del proceso de diagnóstico del daño a la formación, ya que es
durante estas operaciones que la mayoría de los daños son detectados. Donde se
debe tener en cuenta:
•El historial del día a día de la perforación detectando los aditivos del lodo
• La detección de pegamientos de tubería, perdida de circulación, fallas de

funcionamiento del equipo de control de sólidos y cambios en las tasas de
penetración
•El tipo de fluido presente durante el cañoneo
•El método de cañoneo (penetración de los disparos)
•El niel de filtración de los fluidos
98
•Es importante ver los cambios en la curva de producción que estén asociados a
eventos de reparación o estimulación del pozo ya que es el mejor indicio del
proceso que género el daño.
Comparación con pozos vecinos
La superposición de curvas del historial de producción entre pozos vecinos del

mismo yacimiento podrá dar una indicación del comportamiento irregular de
alguno de ellos.
Análisis nodal
El estudio del análisis nodal de pozo sirve para determinar el diseño correcto de
tuberías, estranguladores, equipo de levantamiento artificial por gas y líneas de
flujo. Para lograr esto lo que interesa es que la caída de presión en el yacimiento
sea mínima para lograr la máxima presión de fondo fluyente para levantar los
fluidos producidos.
Registros de Producción
El objetivo principal de los registros de producción es el análisis del sistema de

movimiento de los fluidos, definiendo sistema como el tipo de régimen de flujo en
el pozo. El Registro de Producción es un método que mide y registra el flujo de
fluidos, o el efecto del flujo de fluidos, al colocar los instrumentos de medición a
varias profundidades en una producción o inyección de pozo.
Pruebas y análisis de núcleos
Se realizan pruebas en núcleos extraídos del yacimiento con el fin de reproducir

los fenómenos que han ocurrido en el yacimiento, permitiendo determinar la
existencia del daño, los mecanismos que lo provocan y las posibles soluciones al
mismo. Las técnicas empleadas son:
 Análisis petrográfico
 Análisis mineralógico
 Pruebas de flujo a través de los núcleos
 Análisis químico de los fluidos
 Compatibilidad de los fluidos
99
Pruebas de incremento de presión (build up test)
Estas pruebas consisten en una serie de mediciones de presión de fondo durante

un periodo de tiempo, luego de cerrar el pozo después de haber estado fluyendo a
un gasto constante y estabilizado.
Son utilizadas para encontrar:
 Presión estática promedio en el área de drene o yacimiento (p ws)

 Permeabilidad promedio en el área de drene (k)
 Daño total
 Presencia de limites o heterogeneidades (fallas, contactos, barreras
estratigráficas)
 Interferencia o comunicación entre pozos/fallas
Fig. 3. 6 Incremento de la presión respecto del tiempo en una

prueba de incremento presión
Fig. 3.7 Decremento del gasto respecto del tiempo en una prueba de
Incremento de presión
100
La prueba del pozo requiere que se produzca a un gasto estabilizado durante un

cierto tiempo, denominado tiempo de flujo (tp), para lograr una distribución
homogénea en la presión antes del cierre.
Al cerrar el pozo se mide la presión de fondo fluyendo (p wf @ t=0) y se empieza a

medir en función del tiempo de cierre ( t)
El tiempo de flujo (tp) se define como sigue:
(3. 42)
Dónde:
Np= petróleo acumulado por el pozo antes del cierre [bls]
Qo= gasto de aceite con flujo estabilizado antes del cierre [bpd]
Ecuación de Horner
Horner en 1951 establece que para flujo no estable se cumple la siguiente
ecuación:
(3. 43)
Donde:
pi = presión inicial del yacimiento
La ecuación de Horner sugiere que la relación entre pws y es una línea recta
en escala semi-log.
(3. 44)
Donde:
a = ordenada al origen (pi)
101
m = pendiente
La presión inicial se puede estimar al extrapolar el tiempo de cierre t para valores

muy grandes ( t infinito).
Para un t muy grande tiende a 1, estas suposición solo es válida cuando el

yacimiento posee poco tiempo de producción.
Para la determinación del daño se estima mediante la siguiente formula:
(3. 45)
μ
Dónde:
pw f(δt=0): Presión de fondo inmediatamente antes del cierre del pozo
p1hr: Presión leída de la recta en el gráfico de Horner a δt=1 hr
m: Pendiente de la recta en el gráfico de Horner
Ejemplo 3.1
Asuma que un pozo tiene un radio rw igual a 0.3048, con una penetración de 3
pies de daño y un rs= 3.328 pies. Con una relación de k/ks =5 y k/ks = 10;
Determine:
a) El daño para cada una de las relaciones de permeabilidad alterada.
b) ¿Cuál sería el espesor de la zona dañada para para proveer el mismo daño que
la relación k/ks=5?
102
Solución:
a)
Para la rs y rw dados y k/ks=5 obtenemos
S= (5-1) =9.3
Para k/ks=10
S= (10-1) =20.9
b)
Para obtener el mismo daño que con la relación k/ks=5
Si S=20.9 tenemos:
rs= rw = (0.3048 ) = 61 pies
Se necesitaría una profundidad de daño de 61 pies para producir el mismo

daño.
Ejemplo 3.2
Asuma un pozo productor de aceite tiene las siguientes características:
re  2980 pies  o  1.7 [cp ]

 lb   bl @ c. y 
Py  5661  2   Pws Bo  1.1 
 pg   bl @ c.s 
k  8.2 [mD] rw  0.328 [ pies ]
h  53  pies   lb 
Pwf  4500  2 
 pg 
Calcular el gasto de aceite si la presión de fondo fluyendo es de 4500 lb/pg 2.
Considere S= + 10. Describa dos mecanismos que aumente el gasto de aceite en
un 50 % muestre los cálculos efectuados
Solución:
Para incrementar en un 50% el gasto de aceite:

103
1) Incrementar la caída de presión (pws-pwf) en un 50%
 lb 
(5651-pwf) = 1.5 = 3925  2 
 pg 
 lb 
pwf= 5651-(1.5 = 3924.5  2 
 pg 
2) Disminuir el daño. En este caso
S 2= 3.628
Ejemplo 3.3
 lb 
Un pozo con flujo radial tuvo una caída de presión de 700  2  cuando S=0,
 pg 
si el pozo en la actualidad tiene una eficiencia de flujo de 0.75. Determine el
daño actual del pozo, si este tiene las mismas características del ejercicio 2.2.
Solución:
De la definición de eficiencia de flujo
EF=
= (0.75)(700)=525 lb/pg2
Sustituyendo en la ecuación de Darcy para flujo radial
104
525 lb/pg2=
Despejando y sustituyendo
S=
S=-0.473
Ejemplo 3.4
Un pozo de aceite tiene un rw=0.3048 pies y un área de drene circular de 16

acres, si tiene una eficiencia de flujo de 0.6. Determine el daño que causa esta
eficiencia de flujo.
Solución:
Utilizando la ecuación…… donde E.F esta en términos de re y rw
S=5.644
105
Ejemplo 3.5
Determine el daño por ángulo de inclinación de los disparos de un pozo de

aceite con las siguientes características
rw  0.1 [m] θ  45

h w  79.5 [m] h  225 [m]
z w  197 [m]
Solución:
h 225
1. hD    2250
rw 0.1
z w 197
2. z wD    1970
rw 0.1
z wD 1970
3.   0.875
hD 2250
h w 79.5
4. h wD    795
rw 0.1
h wD cos45 (795)(0.71)
5.   0.25
hD 2250
6. S θc hD100  5.968

S θc hD1000  10.366
S θc hD2250  ?
Interpolamos
Se calcula la pendiente
y 2  y1 10.366  5.968
m   4.8867 103
x 2  x1 1000  100
Calculo de la ordenada al origen
b  5.968  (4.8867 103 )(100)  5.4793

106
y  mx  b  (4.8867 103 )(2250)  5.4793  16.4743

Sθc hD2250  16.4743
Ejemplo 3.6
Determine el daño por ángulo de inclinación de los disparos de un pozo de aceite

con las siguientes características:
rw  0.5 [ft] θ  45

h w  65.75 [ft] h  50 [ft]
z w  40 [ft]
Solución:
h 50
1. hD    100
rw 0.5
zw 25
2. z wD    50
rw 0.5
z wD 50
3.   0.5
hD 100
h w 65.75
4. h wD    131.5
rw 0.5
h wD cos45 131.5(0.71)
5.   0.99  1
hD 100
107
z wD
Para una  0.5
hD
y
h wDcos45
1
hD
Sθc hD100  -1.178

Ejemplo 3.7
Calcule el daño debido al efecto de los disparos de un pozo de aceite con las
siguientes características:
rw  0.3 [pies] lp  9 [pg] θ  120

h  0.5 [pg] rperf  0.25 [pg] k h /k v  10
Solución:
1. αθ  0.648
9 
2. r' w ( )  αθ rw  lp   0.648 0.3    0.68

 12 
 r   0.3 
3. Sp1  ln w   ln   0.818
 r' w ( )   0.68 
hk h /k v  0.510
0.5 0.5
4. hD    0.176
lp 9

r 1 k v /k h 
0.5

0.25
12

1 1/10 
0.5

5. rpD  perf  0.32906
2h (2)(0.5/12)
108
6. a1  2.018 a 2  0.0634
b1  1.6136 b 2  1.7770
7. a  a1 logrpD   a 2
 2.018 log(0.329) 0.0634
 1.03
b  b1 rpD  b 2
 (1.6136)(0.329)  1.777
 2.3079
8. Sp2  10a hDb-1 rpD
b
 (101.03 )(1.1761.3079 )(0.3291)2.3079

 0.0862
9. c 1  6.6x10 -3 c 2  5.320
rw 0.3
rwD    0.2857
lp  rw 9/12  0.3
Sp3  c 1 ec 2rwD  6.6x10 3 e5.320*0.2857  0.0302
10. Sp  Sp1  Sp2  Sp3  0.818  0.0862  0.0302

Sp  0.7024
Ejemplo 3.8
Calcule el daño debido al efecto de los disparos de un pozo de aceite con las
siguientes características:
Cálculo de Sp :
rw  0.328 [ft]
lp  8 [pg]
θ  120
h  0.5 [ft]
rperf  0.25 [pg]
k h /k v  5
109
Solución:
1. αθ  0.648
 8 
2. r'w ( )  αθ rw  lp   0.648 0.328    0.644
 12 
 r   0.328 
3. Sp1  ln w   ln   0.675
 r'w ( )   0.644 
hk h /k v  0.55  (12)

0.5 0.5
4. hD    1.677
lp 8

rperf 1  k v /k h 
0.5

0.25
12

1  1/5 
0.5
5. rpD   0.0301
2h (2)0.5
6. a1  2.018 a 2  0.0634
b1  1.6136 b 2  1.7770
7. a  a1 logrpD   a 2  2.018 log(0.0301)  0.0634  3.132

1
b  b1 rpD  b 2  (1.6136)(0.0301)  1.777  1.8255
8. S p2  10 a h Db-1 rpD
b
 (10 3.1321 )(1.6771.6136 )(0.03011.8255)  5.212
9. c 1  6.6x10 -3 c 2  5.320
rw 0.328
rwD    0.329
l p  rw 8/12  0.328
S p3  c 1 e c 2rwD  6.6x10 3 e 5.320*0.329  0.0379
110
Ejemplo 3.9
La tabla siguiente muestra los datos de una prueba de incremento de presión de

un pozo con un radio de drene estimado de 2640 pies. Antes del cierre el pozo
produjo 4900 bpd por 310 horas. Los datos del yacimiento son:
profundidad = 10476 pies
rw = 4.25 pg
Ct = 22.6x10-6 (lb/pg2)-1
qo = 4900 bpd
h = 482 pies
pwf (Dt=0) = 2761 lb/pg2
mo = 0.20 cp
f = 0.09
Bo = 1.55
dTR = (6.276/12) pies
tp = 310 horas
Dt (horas) pw (lb/pg2)
Dt (horas) pw (lb/pg2) 3.04 3283
0 2761 3.46 3286
0.1 3057 4.08 3289
0.21 3153 5.03 3293
0.31 3234 5.97 3297
0.52 3249 6.07 3297
0.63 3256 7.01 3300
0.73 3260 8.06 3303
0.84 3263 9 3305
0.94 3266 10.05 3306
1.05 3267 13.09 3310
1.15 3268 16.02 3313
1.36 3271 20 3317
1.68 3274 26.07 3320
1.99 3276 31.03 3322
2.51 3280 34.98 3323
37.54 3323
111
Solución:
Graficando t vs pw se tiene:
Gráfico semi-log de Horner
Dt (horas) pw (lb/pg2) (tp+Dt)/Dt

0 2761
0.1 3057 3,101.0
0.21 3153 1,477.2
0.31 3234 1,001.0
0.52 3249 597.2
0.63 3256 493.1
0.73 3260 425.7
0.84 3263 370.0  tp  Δt 
0.94 3266 330.8  Δt  Δt1hr
t=1 1.05 3267 296.2
hr 1.15 3268 270.6
1.36 3271 228.9
1.68 3274 185.5
1.99 3276 156.8
2.51 3280 124.5
3.04 3283 103.0
3.46 3286 90.6
4.08 3289 77.0
5.03 3293 62.6
5.97 3297 52.9
6.07 3297 52.1
7.01 3300 45.2
8.06 3303 39.5
9 3305 35.4
10.05 3306 31.8
13.09 3310 24.7
16.02 3313 20.4
20 3317 16.5
26.07 3320 12.9
31.03 3322 11.0
34.98 3323 9.9
37.54 3323 9.3
112
Despejando k de la pendiente
162.6 49001.550.20
k  12.8 md
40482
Sustituyendo valores en la ecu. 3.45
 3266 - 2761
s  1.1513 

 log
12.8122  
  3.2275  8.6
 40  
6 2 
 0.090.20 22.6 x 10 4.25  
113
114
4. CURVAS DE DECLINACIÓN.
4.1 Declinación Transitoria
La Declinación Transitoria se considera una declinación natural causada por la

expansión del aceite, gas y agua en una región de drene con un incremento
continuo del radio de drene. Al abrir un pozo se altera el estado de equilibrio del
yacimiento y se crea una respuesta de presión del pozo.
Fig.4.1 Disturbio de presión

propagado lejos del pozo
El disturbio de presión se propaga gradualmente lejos del pozo, incrementando

el área de drene del pozo. Conforme el disturbio de presión se propaga hacia las
fronteras externas del yacimiento, las condiciones de producción cambian
rápidamente en función del tiempo.
Existen dos modos de producción transitoria, los cuales son considerados y

tratados matemáticamente:
a) Gasto de producción constante.
b) Presión de fondo fluyendo constante.
Gasto de producción constante. Presión de fondo fluyendo constante
: radio aparente y es definido como rwa = rw e-s
115
4.2 Declinación en Estado Pseudoestacionario
En este tipo de declinación, será necesario considerar la producción de aceite

como un conjunto o serie de producciones en estado estacionario para describir el
comportamiento de un pozo. El inicio del abatimiento (de presión) está
determinado por el tiempo en el cual, el radio de drene ha alcanzado las fronteras
externas de no- flujo. De allí en adelante como resultado de la producción, la
región total drenada por el pozo comienza a ser declinada y de este modo, la
caída de presión a lo largo del área total de drene, tal como se muestra en la Fig
4.1. El rango de declinación de la presión depende de los siguientes factores:
 Rapidez con la que los fluidos son producidos.
 Expansión de los fluidos del yacimiento.
 Compactación del volumen de poros.
Cuantificar la declinación de presión será labor del ingeniero de yacimientos, y

para realizarlo, usualmente tendrá que calcularla a partir de un balance de materia
volumétrica. El efecto más importante de la declinación es el deterioro del
comportamiento de afluencia, reflejado mediante la declinación de la presión
media del yacimiento y el incremento de la resistencia al flujo.
La Fig.4.2 se ilustra dos casos de declinación:
1. Declinación a gasto constante.
2. Declinación a presión constante
.
Fig. 4.2 Declinación de presión. Gasto constante. Presión
constante 116
Declinación a Gasto Constante
El balance de materia que relaciona la declinación de la presión media del

yacimiento py con la producción acumulada de aceite Np es:
(4. 1)
(4. 4)
Dónde:
A = Área [pie2]
BO=Factor de volumen del aceite [BY/BS]
Ct= Compresibilidad total de la formación [(lb/pg2)-1]
h= Espesor de la formación [pies]
= Presión estática [lb/pg2]
= presión de fondo fluyendo [lb/pg2]
Np= Producción acumulada de aceite [bls]
Declinación a Presión Constante
En caso de una declinación a Presión constante, la expresión para la declinación

de la Presión de un yacimiento bajosaturado es más complicado. El gasto de
producción en yacimientos declinados está expresada por la ecuación de flujo
radial:
117
(4. 5)
Donde la presión en la frontera externa Pe(t) declina en función del tiempo. Una
ecuación de balance de materia relaciona la Np con la presión Pe(t). Esta expresa
la producción acumulada en función de la Compresibilidad total aparente del
sistema Cta, el volumen de poros Vp(1-Sw), y la caída de Presión en el yacimiento
(Pi - Pe(t)).
(4. 7)
La base para calcular la Declinación del Gasto de Producción es un conjunto de

curvas características, definidas como “curvas tipo”, las cuales fueron
desarrolladas por Fetkovich (1980).
Estas curvas son el resultado de investigaciones empíricas y desarrollos

matemáticos aplicados para soluciones rápidas en una gran variedad de
problemas relacionados con la declinación de la producción.
Los principales períodos de declinación de un pozo productor son:
1.- Declinación transitoria.
2.- Declinación en estado pseudoestacionario.
Así mismo, dentro de la declinación en estado pseudoestacionario se encuentran

otros tres tipos de declinación las cuales son:
• Declinación exponencial.
• Declinación hiperbólica.
• Declinación Armónica.
Son métodos para estimar reservas y regímenes de producción futuros, están

basados en la suposición de que los pozos permanecerán sin la influencia de
controles artificiales; que se les permite constantemente mantener la producción
máxima de que son capaces, y de que siempre producen bajo las mismas
condiciones de operación.
118
1. Declinación Exponencial (geométrica, semilog o de porcentaje constante):
Los cambios en la producción por unidad de tiempo son constantes.
2. Declinación Hiperbólica (log - log): La caída en la producción por unidad de

tiempo, expresada como una fracción del gasto de producción, es una
fracción elevada a una potencia entre cero y uno.
3. Declinación Armónica: Es un caso particular de la declinación hiperbólica,

donde la potencia de la fracción del gasto producido es la unidad.
Declinación exponencial
El gráfico de producción de aceite contra tiempo para un pozo, podría ser

extrapolado en el futuro para proporcionar una estimación de los gastos futuros de
producción. Conociendo los gastos futuros de producción es posible determinar la
producción futura total o reserva del yacimiento en cuestión.
Cuando el gasto de producción se grafica contra el tiempo, se puede observar que

el gasto declina con el tiempo tal como se ilustra en la fig 4.3 .
Este tipo de declinación se presenta cuando se tiene un yacimiento cerrado que

produce a una presión de fondo fluyendo constante y cuando se sienten los
efectos de frontera.
Cuando el logaritmo de gastos de producción es trazado contra tiempo lineal, a

menudo resulta una línea recta. Este fenómeno se refiere como Declinación
Exponencial, se expresa:
q  qi e ( t / a ) (4. 8)
donde “a” es una constante positiva.
119
Fig.4.3 Gasto de producción vs tiempo
Pero esta también se puede expresar como el gasto contra la producción

acumulada, así:
N p  a(qi  q) (4. 9)
En función de la rapidez de declinación, D=1/a, las expresiones anteriores quedan:
q  q i e (  Dt ) (4. 10)
Al transformar la ecuación 105 de logaritmo natural a logaritmo base 10,

obtenemos:
D
log( q)  log( q i )  t (4. 11)
2.3
Que gráficamente queda representada como una línea recta de pendiente (-D/2.3)
y ordenada al origen qi.
Al extrapolar esta línea hasta el LE (Limite Económico) puede conocerse la vida

futura del pozo.
120
Porcentaje de declinación
El cual queda representado como sigue:
100
% de Declinacio n  ( qi  q) (4. 12)
Np
Tiempo de vida del yacimiento
El tiempo de vida útil del yacimiento se determina sustituyendo el valor del LE:
(4. 13)
Limite económico
El Límite Económico se obtiene mediante una extrapolación de la curva de

declinación, debe llevarse hasta un punto en el que el valor de la producción sea
equivalente a los gastos de producción, el cual nos indica hasta cuando es
rentable la explotación del pozo en estudio. El valor de la producción mínima que
121
sufraga los costos de operación, mantenimiento de equipo, personal empleado,
etc., es conocido como Límite Económico.
La expresión matemática que permite determinar este valor es la siguiente:
(4. 14)
dónde:
 $ 
C: Costo estimado de operación al límite económico,  Año 
 bl 
L.E : Límite económico,  Año 
$
O: Precio de aceite,  bl 
$
S: Monto de regalías, impuestos, etc.,  bl 
Declinación hiperbólica
No todos los pozos productores exhiben declinación exponencial durante la

declinación. En muchos casos se puede observar una declinación hiperbólica
gradual donde el comportamiento de gasto contra tiempo es estimada de mejor
forma que a partir de la solución analítica.
La declinación hiperbólica es el resultado de energías (mecanismos de empuje)

naturales o artificiales que disminuyen el abatimiento de presión comparado con el
abatimiento causado por la expansión de un aceite ligeramente compresible.
• Gas en solución
• Expansión del casquete de gas
• Empuje de agua
Su expresión matemática es:
122
(4. 15)
Donde b representa la constante de declinación o ritmo de declinación (constante

positiva) 0<b<1.
Al integrar en dos ocasiones la ecuación anterior tenemos:
(4. 16)
Donde Di es la rapidez de declinación cuando el gasto qi prevalece, y t es el

tiempo que tarda en declinar el gasto de qi a q.
Esta ecuación, se ajusta a la ecuación de una línea recta en papel log-log cuando
cambia horizontalmente sobre la distancia (1/Di b), donde 1/b es la pendiente de la
recta.
Para determinar la ecuación de gasto-producción acumulada se debe integrar

respecto al tiempo la ecuación anterior por lo que nos queda:
(4. 17)
El porcentaje de declinación mensual se obtiene mediante:
(4. 18)
(4. 19)
Por lo que finalmente se expresa como:
123
(4. 20)
Por otra parte, tanto las ecuación para la declinación exponencial como para la
hiperbólica nos permiten realizar una extrapolación matemática de datos, difiriendo
una respecto de la otra en lo que corresponde al valor de b, ya que para la
declinación exponencial b= 0, mientras que para la declinación hiperbólica, el valor
de b oscila entre 0.25 y 0.6, sin llegar a exceder este último valor.
El comportamiento tanto de la declinación exponencial como el de la hiperbólica al

ser graficados en papel semilogarítmico es el siguiente:
Fig.4.4 Grafica de la declinación exponencial
Declinación Armónica
Algunas veces cuando la producción es controlada predominantemente por

segregación gravitacional, la rapidez de declinación D es proporcional al gasto q.
124
Este tipo de declinación es un caso particular de la declinación hiperbólica y ocurre
cuando el valor de b=1.
qi
q ………ec.4.22
(1  Di)
Fig.4.5 Representación gráfica de la

declinación armónica en papel cartesiano.
Al integrar esta última ecuación obtenemos la expresión de gasto producido-

producción acumulada:
(4. 24)
La ecuación anterior puede ser representada por una línea recta en papel
semilogarítmico, graficando gasto producido en la escala logarítmica.
En lo que respecta al porcentaje de declinación mensual, se tiene lo siguiente:
(4. 25)
dq qi D

dt (1  Dt) 2
Por lo que finalmente nos queda:
(4. 26)
100 D
Porcentaje de declinació n  
(1  Dt)
El tiempo de vida útil para declinación armónica:
(4. 27)
125
4.3 Ejercicios de Aplicación
Ejemplo 4.1
El gasto de producción de un yacimiento declina de 11300 a 9500 BPD en un periodo

de un año. Calcular la producción acumulada (Np) para un gasto de abandono de 850
BPD. Utilizando los métodos de declinación exponencial (b=0), hiperbólica (b=0.5) y
armónica (b=1).
Solución:
Calculando la tasa de declinación
1) Declinación exponencial
D=
2) Declinación hiperbólica
3) Declinación armónica
126
Ejemplo 4.2
Elabore un perfil de producción de crudo de un campo y calcule el ritmo de

producción para los meses 1-20 por el método de declinación exponencial, con los
siguientes datos:
qi=1500 BPD
D=0.020 (mensual)
Solución:
Utilizando la ecuación
qo t
14970.03 1
14910.2695 2
14821.0757 3
14702.9801 4
14556.683 5
14383.0467 6
14183.087 7
13957.9634 8
13708.9678 9
13437.512 10
13145.1149 11
12833.3879 12
12504.0201 13
12158.7637 14
11799.4179 15
11427.8139 16
11045.7993 17
10655.2231 18
10257.9211 19
9855.7023 20
127
128
CAPITULO 5 REGISTROS DE PRODUCCIÓN
5. REGISTROS DE PRODUCCIÓN
Los registros de producción nos dan un conocimiento de la naturaleza y

comportamiento de los fluidos en el pozo durante periodos de producción o
inyección. Estos registros se pueden tomar después que se han cementado las
tuberías de revestimiento, permitiendo conocer con más detalle no solo el
comportamiento de los pozos, sino también de las formaciones. Mediante los
registros de producción se pueden determinar características tales como:
 Evaluación temprana para lograr el diseño de una terminación eficiente.

 Detección de cambios en el comportamiento de los fluidos.
 Zonas que toman fluidos (zonas ladronas).
 Canalización a través el cemento.
 Perforaciones taponadas.
 Determinación de zonas productoras o receptoras de fluidos.
 Evaluación de la eficiencia en el proceso de inyección.
 Control sobre los procesos de producción.
 Guía esencial para diseño de programas de reparación de pozos,
proyectos de recuperación secundaria y terciaria.
Su aplicación más común es la medición del perfil del flujo del pozo, es decir, la
distribución del fluido dentro y fuera del pozo, y con ello determinar el potencial
productor de las distintas zonas presentes en el/los intervalo(s) productore(s).
Dentro de la evaluación del rendimiento de un pozo se pueden hacer las

siguientes determinaciones:
 Contribución de cada zona.

 Tipos y porcentajes de fluidos por zonas.
 Medidas de temperatura y presión.
 Flujos cruzados o perdida de fluidos es zonas de bajas presiones.
129
Dentro de los análisis de problemas mecánicos se pueden detectar:
 Fisuras en la Tubería de Revestimiento.

 Fisuras en la Tubería de Producción.
 Escapes a través de empaquetaduras.
 Comunicación entre zonas por el espacio anular debido a mala
cementación.
Los registros de producción más comúnmente usados son:
 Medidores de flujo.
 Medidores de densidad.
 Medidores de cortes de agua.
 Medidores de temperatura.
Los registros de producción convencionales se basan en mediciones hechas con

sensores centrados en el pozo en los pozos verticales o casi verticales, el
comportamiento de fluidos es relativamente simple y lo sensores convencionales
de los registros de producción son más precisos para medir los parámetros del
flujo de fluido zona por zona.
Los ingenieros de producción aplican la toma de registros de producción como una

herramienta para diagnosticar las causas de las bajas productividades de los
pozos. Como tales, los registros de producción pueden indicar acciones de ajustes
a ser tomados en cuenta para incrementar la productividad del pozo.
5.1 Productividad Anormalmente Baja
La causa de baja productividad en un pozo varía de un problema fundamental de

yacimiento a restricciones de flujo en las vecindades del pozo o dentro del mismo.
Causas potenciales incluyen la baja permeabilidad relativa a la fase de
hidrocarburos, daño en la formación, pobre penetración o sello de las
perforaciones (u otras restricciones en la terminación, tales como empacador de
grava parcialmente sellado) y restricciones del pozo.
Aquí se define baja productividad de un pozo como uno, que tiene anormalmente
un bajo Índice de Productividad (J); esto es diferente de un pozo con baja tasa de
130
productividad, como uno de baja tasa de producción, puede deberse a insuficiente

cierre debido a fallas de mecanismos de bombeo o excesiva caída de presión en
la tubería.
El primer paso, en la evaluación de la baja productividad del pozo es medir el

índice de productividad. Si se encuentra que es anormalmente bajo, (comparado
con la etapa inicial de la vida del pozo o con pozo similares localizados en las
vecindades, por ejemplo). Es necesario distinguir entre la baja capacidad de la
formación y las restricciones de flujo en la vecindad del pozo o en la terminación.
Este es el intento para efectuar pruebas de presión transitoria para medir el kh del
yacimiento y factor de daño (s).
5.2 Producción Excesiva de Gas o Agua
La producción excesiva de gas o agua es un problema comúnmente visto en los

pozos petroleros y las causas pueden ser por roturas de la tubería de
revestimiento, canalizaciones detrás de la tubería de revestimiento, flujo
preferencial a través de la zona de alta permeabilidad en el yacimiento, o la
conificación. El registro de producción, se puede emplear para localizar la fuente
de la producción de gas o de agua y son apropiados para determinar la causa de
la baja producción no deseada.
Canalización
La canalización entre la tubería de revestimiento y la formación causada por

condiciones de pobre cementación, son algunas veces las causas de producción
de altas tasas de agua o gas. El registro de cementación o registros de pulso eco-
ultrasónico pueden indicar la posibilidad de canalización al medir las propiedades
del cemento detrás de la tubería de revestimiento. Para identificar una
canalización, es necesario un registro de producción que puede responder al flujo
detrás de la tubería de revestimiento. Los registros que pueden apoyar a este
propósito son temperatura, trazador radioactivo y registros de ruido (frecuencia o
sónico de cementación). El tratamiento que se aplica para eliminar la canalización
es la cementación forzada; para diseñar el tratamiento de la cementación, se debe
conocer toda la zona canalizada que aporta la producción no deseada.
131
Fig. 5.1 Canalización de gas
Fig. 5.2 Canalización de agua
Flujo preferencial de gas o agua a través de capas de alta permeabilidad

El flujo preferencial de agua o gas a través de zonas de alta permeabilidad (a
menudo referidas como zonas ladronas), como se ilustra en las Figuras (Clark y
Schultz, 1956), es una causa común de alta producción de gas o agua en pozos
de aceite. Las entradas no deseadas de gas o de agua pueden ser localizadas con
registros de producción.
Una excesiva producción de agua puede ser resultado de inyección de agua en un

acuífero o en una invasión del acuífero. Un registro detallado del pozo productor
puede identificar la localización de la zona de alta permeabilidad o zonas que
contribuyen a la alta producción. Sin embargo, la localización de la entrada de
agua no es generalmente información suficiente para identificar la causa de
132
producción de agua, como siendo un flujo a través de la zona ladrona.

Particularmente, si la zona de entrada de agua está al fondo del intervalo
terminado, la fuente de agua puede ser canalizada o conificada de zonas
inferiores.
Debido a que la interpretación y cálculo de registros en la zona de flujo multifásico

en pozos productores son generalmente menos resolutivos (y más caros) que
aquéllos en un flujo de una sola fase, en operaciones de desplazamiento de agua,
la distribución del agua en un yacimiento es a menudo monitoreada al calcular los
perfiles de inyección y asumiendo la continuidad de las capas del yacimiento entre
los pozos inyectores y los productores.
Excesiva producción de gas, puede resultar de flujo de gas inyectado o de una

zona de casquete de gas. De nuevo, un registro calculado en un pozo productor,
identificará las zonas de entrada, o las zonas de alta permeabilidad causando alta
producción de gas, puede ser inferida de cálculo de registros en los pozos de
inyección de gas, cuando el gas es inyectado en el yacimiento.
Fig. 5.3 Irrupción de agua temprana en capas

altamente permeables
133
Fig. 5.4 Irrupción temprana de gas en capas altamente

permeables
Conificación de gas y agua
La conificación de gas resulta cuando un pozo es terminado cerca del contacto

gas/aceite y existe suficiente permeabilidad vertical para que migre el gas hacia la
parte inferior en el pozo a medida que la presión disminuye alrededor del pozo.
Similarmente, el agua puede ser conificada hacia arriba, proveniente de un
acuífero inferior, si la permeabilidad vertical es suficientemente alta.
Fig.5.5 Conificación de agua
134
Fig. 5.6 Conificación de gas
La conificación es un fenómeno difícil de identificar conclusivamente con los

registros de producción. Considera un pozo que está experimentando conificación
de agua. Un registro de flujo, indicará producción de agua de la parte inferior del
intervalo productor.
Esta agua podría resultar de la canalización localizada debajo del intervalo

perforado, ya sea a través de la zona de alta permeabilidad del flujo en la parte
inferior del intervalo o bien conificada. Un registro que responde al flujo externo de
la tubería de revestimiento, como el registro sónico, se puede usar para identificar
y eliminar la canalización como la fuente de agua (note que la temperatura del
registro responderá similarmente a la conificación y canalización y no será
distinguible entre ellas).
5.3 Uso de Registros de Producción para la Evaluación de Trabajos de

Reparación
Los registros de producción proveen información útil para planear los trabajos de
reparación, primero porque dan información de la distribución de los flujos de las
diversas fases en el pozo. En forma similar, pueden ayudar a evaluar el éxito o
fracaso de reparación de pozos. La mayoría de los trabajos de reparación que se
pueden beneficiar con los registros de producción, son la cementación,
perforaciones adicionales, acidificación, fracturamiento y eliminación de agua o
registros de tratamientos de modificación.
135
La mayoría de las aplicaciones avanzadas en las reparaciones, es antes y

después del tratamiento al pozo. Para la evaluación de los tratamientos al pozo, el
uso de los registros de producción es similar a su aplicación en la planeación de
las reparaciones: indican, que regiones del pozo fueron afectadas por la
reparación y su impacto. Adicionalmente, al cálculo al perfil de producción
después del tratamiento, algunos registros pueden usarse para evaluar
directamente del tratamiento mismo. Los ejemplos más comunes son los registros
de temperatura, trazadores radiactivos para medir el tamaño de la fractura en las
cercanías del pozo.
5.4 Diagnostico de Pozos Inyectores
Los registros de producción se usan en inyección de pozos para monitorear el

desempeño del yacimiento y para evaluar los problemas observados en los pozos
de inyección o del yacimiento.
La mayoría de los problemas que pueden surgir son anormalmente bajos o de alta
inyectabilidad, presión anormal o nivel de fluido en el espacio anular y la baja o
alta productividad de agua en los pozos productores. Los registros de producción
se usan para evaluar estos problemas en los pozos inyectores de forma similar a
la descrita para los pozos productores, que es por la medición del flujo de cada
intervalo del yacimiento para verificar un intervalo aislado por la localización de
una zona de alta permeabilidad y para definir debilidades en el equipamiento del
pozo.
La información fundamental que se solicita con un registro de producción en un

pozo inyector es el flujo, la cantidad de fluido que está siendo inyectado en cada
intervalo. Perfiles de flujo son calculados en pozos inyectores con temperatura,
trazador radioactivo, y registro de medición de flujo (PLT) El registro de
temperatura permitirá obtener información cualitativa de los intervalos de inyección
en la formación, mientras que el medidor de flujo o trazador radioactivo definen
más precisamente la distribución de flujo existente en el pozo.
La causa de un cambio de tasa en un pozo, es a menudo fácil de diagnosticar si

los registros de producción se han corrido periódicamente a través de la vida del
pozo. Por ejemplo, El bajo perfil de inyección de agua en el pozo, puede cambiar
gradualmente en la vida del pozo, en la medida que la distribución de la saturación
cambia en el yacimiento. Los registros obtenidos ocasionalmente deberían mostrar
esta situación como una progresión natural en el barrido de agua. Sin el
136
conocimiento de este cambio gradual, una corrida obtenida años después de que
se inició la inyección, puede aparecer suficientemente diferente del perfil y
ocasionar una conclusión errónea que la canalización ha desarrollado o que otros
Cambios Drásticos Han Sucedido.
EJEMPLO 5.1
Uso del perfil de flujo (PLT) para evaluar un pozo dañado.

La tasa de producción del pozo A-1 en el yacimiento Alfa rápidamente ha
declinado a menos de la mitad de la producción inicial en un periodo de 6 meses.
Estimaciones de la presión del yacimiento y mediciones de la presión de fondo
fluyendo (Pwf), se mostró que el índice de productividad del pozo estaba 50%
debajo de los pozos productores vecinos. Se realizó una prueba de presión-
producción y el factor de daño se calculó de 20, mientras que la kh es cerca del
valor esperado.
Solución:
Estrategia y análisis de registros de producción. De la rápida declinación y

el alto factor de daño, el daño a la formación en las inmediaciones del pozo es la
causa de la baja productividad del pozo. Para apoyar el diseño de un tratamiento
de acidificación para eliminar el daño, se corrieron registros de temperatura y de
medición de flujo (PLT) con el fin de medir el perfil de flujo del yacimiento. Los
resultados interpretados se ilustran en la figura.
137
Registro de Medición de flujo y de temperatura
La pista de medición de flujo, muestra que la zona A produce menos del 10% del
flujo total, la zona B produce cerca de 70% y la zona C contribuye con el 25% de
la producción. El registro de temperatura, cualitativamente confirma la
interpretación del registro de medición de flujo.
Aparentemente, la zona A ha sido significativamente dañada durante la

producción, quizás por la migración de finos a la vecindad del pozo. El registro de
producción, muestra la necesidad de un tratamiento para este pozo en forma
selectiva para la zona A, quizás una menor cantidad de estimulación en la zona C.
La zona B no requiere tratamiento; de hecho la pista de flujo muestra que una
buena desviación y eliminación del tratamiento de estimulación se requiere para
minimizar la inyección en la zona B.
138
Ejemplo 5.2
Localización de una zona canalizada de gas con registros de temperatura y de

ruido.
Los registros de temperatura y de ruido (frecuencia, o sónico de cementación)

mostrados en la Figura 5.4 se obtuvieron en un pozo, produciendo un alto
contenido de Relación Gas-Aceite (RGA). Ambos registros claramente indican que
el gas está siendo producido de una capa de arena con gas y existe canalización
hacia abajo a las perforaciones superiores de la zona de aceite.
Ambos registros responden a la expansión de gas a través de las restricciones; el

registro de temperatura muestra anomalías de zonas frías causadas por el
enfriamiento Joule-Thompson en la zona de expansión de gas, mientras que el
registro de ruidos (frecuencia, o sónico de cementación) mide la amplitud de
frecuencia incrementado en las misma zona. De esta forma, ambos registros
responden a la fuente de flujo de gas, a la restricción en el canal detrás de la
tubería de revestimiento y la zona localizada de entrada de gas al pozo.
Registros de temperatura y de ruido
139
Solución:
Para eliminar la excesiva producción de gas, se requiere una cementación que

debe bloquear el flujo de gas en la zona canalizada. Esto puede ser mejor
completado al perforar cerca de la zona de gas y circular el cemento a través del
canal (Nelson, 1990).
Note que la medición de la pista del registro en este pozo podría no ser
particularmente útil para localizar la causa de la alta producción de gas o en la
planeación de las acciones correctivas. Una pista de perfil de flujo podría mostrar
producción de gas en el pozo en la parte superior de la zona de aceite. Esto
podría ser debido a la canalización (como fue el caso) o la alta saturación en la
parte superior de la zona de aceite, como puede ocurrir en un desarrollo de
casquete de gas. Únicamente al correr los registros que claramente identifican la
canalización pueden ser apropiados para la planeación de los trabajos de
reparación.
Ejemplo 5.3
Excesiva producción de gas de una zona ladrona.
Un pozo en el yacimiento Beta está produciendo a una inusual cantidad de gas,

con una producción baja de aceite, comparado con pozos similares en el campo.
¿Qué registros de producción o pruebas se deberían correr para determinar mejor
la migración de gas del casquete a través de una zona ladrona?
Solución:
Una prudente investigación sería primero, correr los registro de temperatura y

registros de densidades de fluidos. Esos registros deberían localizar
cualitativamente la entrada o entradas de gas; en adición, el registro de
temperatura, ayudará a diferenciar entre la producción de una zona ladrona y la
producción de gas resultado de una canalización.
140
Localización de entrada de gas en los registros de temperatura y densidad de

fluidos.
La Figura muestra los registros de temperatura y gradiomanómetros (densidad de

fluidos) que claramente indican producción de gas de una zona ladrona en dicho
pozo. De la anomalía de enfriamiento en el registro de temperatura y la
disminución de la densidad del fluido, la zona B es identificada como una zona
ladrona. A partir de la producción de aceite de la zona A arriba de esta zona, se
muestra por el ligero incremento de la densidad del fluido a través de la zona A, la
producción de la zona de gas no es canalizada o conificada hacia abajo a este
nivel. El registro de temperatura también indica que no ocurre la canalización.
Ejemplo 5.4
Determinar las causas de excesiva producción de agua en el fondo del pozo.
Una serie de registros de producción (temperatura, canasta de medición de flujo y

densidad del fluido), se corrieron en un pozo productor de aceite que está
produciendo una excesiva cantidad de agua, muestra que el agua está siendo
producida del fondo a 20ft de las perforaciones. El campo está siendo barrido por
agua en un patrón de cinco-spot, con la inyección de pozos de agua cercanos
aproximadamente a 800ft.
¿Qué otra información (pruebas de pozos, registros de producción, etc.) sería útil
en la planeación correctiva para este pozo y más generalmente, hacer más
eficiente el manejo del barrido de agua?
141
Solución:
La información primaria requerida para definir el exceso de producción de agua en

este pozo en particular (canalización o conificación) o puede ser debida al
excedente de inyección de agua en la parte inferior del yacimiento. Primero, los
registros de producción se deben enfocar para cualquier evidencia de canalización
o conificación en el pozo. El registro de temperatura, en particular, puede indicar si
la ocurrencia del flujo proviene de la zona que esté debajo de las perforaciones.
Enseguida, el funcionamiento de la inyección de agua debe ser investigada en los
pozos que están alrededor. Se deben registrar perfiles de inyección si no se han
obtenido recientemente.
Si el patrón de inyección de agua en uno o más pozos inyectores muestra altos

volúmenes de agua en la zona inferior, con un flujo preferencial a través de la
zona de alta permeabilidad, ésta puede ser la causa de la elevada producción de
agua.
Finalmente, el volumen total de producción de agua en el pozo podría ser

disminuido si se encuentra que la producción de agua ha cesado, y la conificación
de agua ha sido identificada como la causa principal de excesiva producción de
agua.
Ejemplo 5.5
Medición de la fractura
El pozo 2, uno de los primeros pozos del campo D, es candidato a ser fracturado
hidráulicamente. Para apoyar en el diseño de la fractura, se tomaron los registros
de temperatura antes y después de la misma. Y las últimas 10,000 lb de
sustentante se etiquetaron radioactivamente para medir el tamaño de la fractura.
Determinar el tamaño de la fractura a partir de los registros de temperatura y de la
investigación post fractura de los rayos gamma.
Solución:
Registro de Temperatura de tamaño de la fractura.
Debido a que el fluido de inyectado por lo general es significativamente más frío

que la formación que está siendo fracturada a típicos rangos de inyección, la
temperatura del fluido en la fractura del pozo, será cercana a la temperatura de la
superficie. Como el tratamiento procede, la formación no fracturada alrededor del
pozo es enfriada por el fenómeno de conducción radial de calor, como ocurre en
cualquier pozo inyector, mientras que el fluido frío es introducido en la fractura.
Cuando el pozo es cerrado, opuestamente al pozo, la formación no fracturada
142
inicia a regresar a la temperatura geotérmica por la inestable conducción radial de

calor radial, mientras en la región fracturada la temperatura del pozo es afectada
por la conducción linear de calor de la formación a la fractura. Debido a que la
transferencia radial de calor en la región no fracturada es más rápida que la
conducción linear de calor en la fractura, la región fracturada se calentará más
lentamente, dando lugar a una anomalía fría en un registro de temperatura. De
esta forma, el tamaño de la fractura puede identificarse por la localización de la
anomalía fría en una corrida de registro de temperatura después de un breve
cierre (unas horas) posterior al fracturamiento.
Las anomalías calientes en el registro post fractura que corresponden a anomalías

calientes en el registro prefractura resultan de las variaciones propiamente
térmicas; esas regiones no se deberían incluir en la interpretación de la zona
fracturada. Cuando una anomalía caliente aparece en un registro postfractura y no
corresponde a la anomalía existente en el registro prefractura, la anomalía caliente
es aparentemente originada por movimiento del fluido en la fractura después del
cierre o la derivación de la fractura de la intersección con el pozo. La anomalía
caliente de la región, donde se ubica el pozo, debería ser incluida como parte
interpretación de la zona fracturada.
En la Figura se muestran los registros de temperatura que se corrieron en el pozo

D-2 después de circular el pozo con fluido frío y antes del fracturamiento y
posteriormente al proceso de fracturamiento. La extensión vertical de la fractura
está indicada por la región donde existe divergencia de las dos curvas de
registros, mostrando que la fractura está localizada en este caso entre los 10,100
y los 10,300ft. Las anomalías de temperatura que aparecen en el registro de
temperatura, aparentemente se deben a variaciones en difusibilidad térmica de la
143
formación.
Registros de temperatura pre y postfractura
Ejemplo 5.6
En una operación de desplazamiento de agua en el yacimiento A, el agua está

siendo distribuida en diferentes pozos inyectores en un sistema común de
inyección; el agua está distribuida aproximadamente a la misma presión de
cabeza de pozo. Mediciones rutinarias de la inyección individual de pozo muestran
que un pozo estuvo recibiendo aproximadamente 40% más de agua que los pozos
vecinos. La suma de kh de todos los pozos inyectores fue aproximadamente la
144
misma y fueron terminados casi a la misma profundidad. ¿Cuáles son las posibles
causas de la anormalmente alta cantidad de inyección en este pozo y cuáles
registros de producción u otras pruebas que se pueden realizar para diagnosticar
el problema y las acciones de solución?
Solución:
La mayoría de las posibles causas de la alta inyección son roturas en la tubería de

producción, tubería de cementación, o empacador o canalización a otra zona.
Fracturamiento no es la misma causa, porque la similaridad de terminación en los
pozos inyectores que está alrededor, tiene la misma presión de cabeza, aún no
exhiben anormalmente alto volumen de inyectabilidad. Otra menor pero posible
causa es que todos los pozos inyectores vecinos estén dañados en extensión
similar, mientras que el alto volumen de inyección en el pozo es relativamente no
dañado, en comparación con el resto de los pozos.
Para este escenario, los registros de producción pueden detectar debilidades o

canalizaciones y deberían ser corridos en un rango elevado de inyección, una
combinación de un registro de temperatura y de ruidos sería una buena selección
para localizar la rotura o canalización.
145
146
CAPITULO 6 FUNDAMENTOS DEL FLUJO MULTIFÁSICO; FLUJO
MULTIFÁSICO HORIZONTAL Y VERTICAL
6. FUNDAMENTOS DEL FLUJO MULTIFÁSICO, FLUJO MULTIFÁSICO

HORIZONTAL, FLUJO MULTIFÁSICO VERTICAL
6.1 Flujo Multifásico en Tuberías
La ecuación general que gobierna el flujo de fluidos a través de una tubería, se

obtiene a partir de un balance macroscópico de la energía asociada a la unidad de
masa de un fluido, que pasa a través de un elemento aislado del sistema, tal como
se muestra en la Figura 6.1.
Δx
ΔH
α
1
Fig. 6.1 diagrama de flujo en

un conducto aislado
La ecuación general de energía expresa un balance de energía entre dos puntos

en un sistema de flujo. De acuerdo con el principio de conservación de la energía,
se establece que la energía de un fluido que entra en la sección 1 de una tubería,
más el trabajo adicional realizado sobre el fluido entre las secciones 1 y 2, menos
cualquier pérdida de energía en el sistema entre esas secciones, 1 y 2, es igual a
la energía del fluido que sale de la sección 2.
La ecuación general de energía se utiliza para resolver muchos problemas que

involucran flujo multifásico en dirección vertical, horizontal o inclinada. A partir del
principio de conservación de la energía se tiene que:
(6. 1)
Dónde:
E1 = Energía por unidad de masa en la sección 1
E2 = Energía por unidad de masa en la sección 2
147
W s = Perdida de energía pro fricción.
W f = Perdida o adición de trabajo externo.
Para E1 y E2 contempla las siguientes energías
A. Energía de expansión (Ee)
(6. 2)
Ee =
Donde:
p = presión
v = volumen especifico
B. Energía de potencial (Ep):
Ep = (6. 3)
Donde:
g= Aceleración gravitacional
gc= Constante gravitacional
C. Energía cinetica (Ec)
(6. 4)
Al sustituir las energías correspondientes a las secciones 1 y 2,

suponiendo que el volumen específico no cambia, reordenando términos
igualando a cero, considerando despreciables las pérdidas de energía
148
por trabajo externo y multiplicando por ( w s = 0) (Transporte de

Hidrocarburos, Garaicochea, pag.31-32), se tiene:
(6. 5)
A la ec. 6.5 se acostumbra escribirse de la siguiente forma:
(6. 6)
Donde:
6.1.1 Pérdidas de presión por fricción
Ecuación de Darcy
Darcy, Weisbach y otros, dedujeron en 1857 una la siguiente ecuación para
determinar las pérdidas de presión por fricción en unidades consistentes:
(6. 7)
Donde:
f = Factor de fricción
ρ = Densidad del fluido
gc = Constante gravitacional (32.174)
d = Diámetro interior de la tubería [pg]

149
Ecuación de Faning
Una ecuación similar a la de Darcy fue establecida posteriormente por Fanning,

quien obtuvo valores de f cuatro veces menores que los de Darcy, esta diferencia
es por el uso del radio hidráulico en lugar del diámetro de la tubería en su
correlación. La ecuación establecida por Fanning es:
(6. 8)
Dónde:
Rh = Radio hidráulico
Factor de fricción
El valor del factor de fricción (f) es función del número de Reynolds (NRe) y la
rugosidad (ξ) de la tubería, esto es
f = f (NRe, ξ)
El número de Reynolds se define como:
(6. 9)
Donde
µ= viscosidad del fluido [cp]
En el caso de la rugosidad los valores mas comúnmente empleados en la industria

son:
Tubería Estriada 6x10-5 [pg]

Tubería de Producción 6x10-4 [pg]
Tubería de Escurrimiento 7x10-4 [pg]
Tubería Galvanizada 6x10-3 [pg]
Tabla 6.1
Valores más comunes para ξ
[pg] en la industria
150
Para conocer el valor de f es necesario determinar el régimen de flujo (laminar,

critico o turbulento). Para determinarlo se utiliza comúnmente el número de
Reynolds, el cual determinara la ecuación a utilizar para la obtención de f, de la
siguiente manera:
NRe<2300 Flujo laminar

NRe=3100 Zona de transición
2300<NRe=3100 Flujo critico
NRe>3100 Flujo turbulento
Flujo laminar
(6. 10)
Zona de transición
(6. 11)
Flujo turbulento
(6. 12)
Flujo critico
(6. 13)
6.1.2 Perdidas de presión por elevación
Para calcular las pérdidas de presión por elevación, es necesario predecir con
precisión el colgamiento (HL) considerando el resbalamiento entre las fases.
151
El término “resbalamiento” se usa para describir el fenómeno natural del flujo,

cuando una de las dos fases fluye a mayor velocidad que la otra.
6.1.3 Colgamiento (HL)
Se define como la relación entre volumen de líquido existente en una sección de

tubería a las condiciones de flujo entre el volumen de la sección aludida.
La correlación más general para obtener el colgamiento del líquido es la de

Mukherjee y Brill (1983). La ecuación que establecieron es la siguiente:
 N Cgv5 
HL   2 2
 exp C1  C 2 sen θ  C 3 sen θ  C 4 N L

 

(6. 14)
 N CLv6 
Dónde:
C1 a C6 = Coeficientes dependientes del tipo de flujo [adimensional]
HL = Colgamiento del líquido, [adimensional]
NL = Número de la viscosidad del líquido, [adimensional]
 = Ángulo de inclinación de la tubería con respecto a la horizontal [grados]
Donde NL, Ngv, NLv están dados por:

0.25
 1 3
NL  0.15726 μ L  σ  (6. 15)
 ρL 
Número de la viscosidad del gas

0.25
 ρL 
N gv  1.938 vsg   (6. 16)
σ
Número de la velocidad del líquido

0.25
 ρL  (6. 17)
N Lv  1.938 VsL  
σ
152
Las unidades utilizadas en los términos de estos números adimensionales son:
vsg = Velocidad superficial de gas  pie 

 seg 
 
vsL = Velocidad superficial del líquido  pie 

 seg 
 
L = Viscosidad del líquido [cp]
L = Densidad del líquido  lb f 

 pie3 
 
 = Tensión superficial  dinas 

 cm 
 
Dirección
Del Flujo Tipo De Flujo C1 C2 C3 C4 C5 C6
Horizontal
- 0.38011 0.12988 - 0.11979 2.34323 0.47569 0.28866
Y Ascendente Todos
Estratificado - 0.33028 4.80814 4.17158 56.26227 0.07995 0.50489
Descendente
Cualquier otro - 0.51664 0.78981 0.55163 15.51921 0.37177 0.39395
Tabla 6.2 Coeficientes de la Ec. 6.14
De acuerdo con Mukherjee y Brill (1983), sí:
(0.321 - 0.017 N gv - 4.267sen θ - 2.972N L - 0.033(log N gv ) 2 - 3.925 sen 2 θ)

N Lv  10 (6. 18)
entonces, el flujo es descendente estratificado.
Velocidades superficiales. Es la velocidad que tendría cualquiera de las dos fases

si ocupara toda la tubería, y se definen con las siguientes expresiones:
153
vsL 
0.01191 q 0 Bo  q w Bw 
2 (6. 19)
d
vsg 
0.002122 q R
o - R s  Bg  (6. 20)
2
d
donde:
 pie3g @ c.y.
 3 
Bg = Factor de volumen del gas,  pieg @ c.s. 
 pie3o @ c.y.
 3 
Bo = Factor de volumen del aceite,  pieo @ c.s. 
 pie3w @ c.y.
 3 
Bw = Factor de volumen del agua,  piew @ c.s. 
d = Diámetro interno de la tubería, [pg]
 bl 
 día 
qo = Gasto de aceite,
 bl 
 día 
qw = Gasto del agua,
 pie3 
 bl 
R = Relación gas- aceite,  
 pie3 
 bl 
Rs = Relación de solubilidad,  
La densidad real de la mezcla de fluidos se obtiene a partir del colgamiento de

líquido, HL, de la siguiente manera:
ρ m  ρ L H L  ρg (1 - H L ) (6. 21)
154
Dónde:
g = Densidad del gas  lb 

 pie3 
 
L = Densidad del líquido  lb 

 pie3 
 
ns = Densidad de la mezcla  lb 

 pie3 
 
Otro concepto que se emplea con frecuencia en los cálculos de gradientes de

presión para flujo multifásico, es el colgamiento sin resbalamiento, . Este se define
en la misma forma que HL y puede ser calculado directamente como:
gasto de líquido
λL  (6. 22)
gasto total
O bien, a partir de las condiciones de presión y temperatura existentes

considerando las producciones obtenidas en la superficie (qo y R), esta es:
1
λL 
q o (R - R s )Bg (6. 23)
1
5.615 (q o Bo  q w Bw )
Para este caso, la densidad de la mezcla se calcula:
 ns  ρ L λ L  ρ g (1 - λ L ) (6. 24)
6.2 Flujo Multifásico Vertical
El trayecto de los fluidos a través de la tubería, consume la mayor parte de

presión disponible para llevarlos del yacimiento a las baterías de separación. Se
ha establecido que la caída de presión en esta área de flujo es alrededor de 40 a
80% del total, y depende de variables tales como diámetro de la tubería,
profundidad del pozo, gasto de producción y relación gas-líquido (RGL).
Para describir el comportamiento de flujo a través de la tubería vertical, es decir el

pozo, se requiere un estudio de pérdidas de presión en tuberías verticales que
155
conducen mezclas en dos fases (Donohue y Lang, 1986). La siguiente tabla

muestra los principales patrones de flujo en la tubería de producción y sus
características.
Tipo de flujo Descripción

Flujo de burbujas de gas dispersas en
Flujo burbuja
un medio líquido continuo
Baches de aceite que contienen
pequeñas burbujas de gas están
separados entre sí por bolsas de gas
que ocupan toda la sección transversal
Flujo bache
de la tubería, excepto por una película
de aceite que se mueve relativamente
despacio a lo largo de la pared de la
tubería
Bolsas de gas atravesando baches del
aceite más viscoso, con el resultado de
Flujo anular que el gas forma una fase continua
cerca del centro de la tubería, llevando
hacia arriba gotitas de aceite en ella
Una fase continua de gas en donde las
Flujo niebla gotitas de aceite son transportadas
conjuntamente con el gas
Tabla 6.3 principales patrones de flujo en la tubería de producción.
6.2.1 Correlaciones para flujo multifásico
Las diversas correlaciones existentes para el cálculo de distribuciones de presión

con flujo multifásico en tuberías, pueden clasificarse en tres grupos basándose en
el criterio utilizado para su desarrollo (Brown y Beggs, 1977, vol 1; Garaicochea,
1991; Donohue, 1986).
Grupo I
No se considera resbalamiento entre las fases. La densidad de la mezcla se

obtiene en función de las propiedades de los fluidos, corregidas por presión y
temperatura. Las pérdidas por fricción y los efectos de colgamiento se expresan
por medio de un factor de fricción correlacionando empíricamente. No se
distinguen patrones de flujo. En este grupo se emplean los siguientes métodos:
Poettman y Carpenter (1952)
156
Baxenderll y Thomas (1961)
Fancher y Brown (1963)
Grupo II
Se toma en cuenta el resbalamiento entre las fases. La densidad de la mezcla se

calcula utilizando el efecto del colgamiento. El factor de fricción se correlaciona
con las propiedades combinadas del gas y el líquido. No se distinguen regímenes
de flujo. En este grupo se utiliza el método de:
Hagendorn y Brown (1965)
Fig. 6.2 Patrones de flujo

durante el flujo vertical
Grupo III
Se considera resbalamiento entre las fases. La densidad de la mezcla se calcula

utilizando el efecto de colgamiento. El factor de fricción se correlaciona con las
157
propiedades del fluido en la fase continua. Se distinguen diferentes patrones de

flujo. Las principales correlaciones en este grupo son:
Duns y Ros (1963)
Orkiszewski (1967)
Beggs y Brill (1973)
Gould y Tek (1974)
La precisión de las correlaciones de flujo multifásico es bastante buena, lo cual

permite las siguientes aplicaciones:
 Seleccionar el diámetro de la tubería de producción.
 Predecir cuándo un pozo dejará de fluir y por lo tanto predecir el tiempo para
instalar un sistema de producción artificial.
 Diseñar instalaciones de producción artificial.
 Determinar las presiones de fondo fluyendo.
 Determinar los índices de productividad de los pozos.
 Predecir el gasto máximo de producción.
 Otros.
El gradiente de presión total (o cambio en la presión con respecto a la longitud de

flujo) para flujo multifásico vertical es la suma de tres factores: gradiente de
presión por elevación, gradiente de presión por fricción y gradiente de presión por
aceleración, pero debido a que las caídas de presión por aceleración son muy
pequeñas en comparación con las otras dos, se pueden considerar despreciables
por lo que:
O bien:
(6. 25)
158
2
 Δp  g fv ρ
   ρ h sen θ  (6. 26)
 ΔL  T gc 2 gc d
Donde:
θ = Angulo de flujo, para flujo vertical θ=90°
Por lo que la ec. 6.26 se reduce a:

2
 Δp  g fv ρ
   ρh  (6. 27)
 ΔL  T gc 2 gc d
6.2.2 Correlaciones para flujo multifásico vertical
Poettmann and Carpenter
(6. 28)
Orkiszewski
Régimen de burbuja
Gradiente por elevación
(6. 29)
De donde:
(6. 30)
159
(6. 31)
(6. 32)
Gradiente por fricción
(6. 33)
Régimen de bache
Gradiente por elevacion
(6. 34)
De donde:
(6. 35)
Donde:
(6. 36)
160
Régimen de transición bache-niebla
(6. 37)
Y el término por fricción, por:
(6. 38)
Donde a y b se refieren a la ponderación lineal y están dadas por:
(6. 39)
(6. 40)
De donde:
(6. 41)
(6. 42)
Donde:
Régimen de niebla
161
(6. 43)
(6. 44)
6.3 Flujo Multifásico Horizontal
El siguiente componente superficial es la línea de flujo o descarga, conocida

comúnmente como línea de escurrimiento, cuya función es conducir la producción
del pozo hacia el separador. Aquí se presenta una caída de presión adicional,
cuyo rango varía entre 10 y 15% del total.
Asimismo, dependiendo del tipo de fluido producido por el pozo, los patrones de
flujo presentes en la tubería horizontal son los siguientes (Gómez, 1984; Brown y
Beggs, 1977, vol 1): figura 6.3.
A. Flujo Estratificado. En éste flujo, casi la mayor parte del líquido se desplaza en
el fondo de la tubería y el gas en la parte superior de la tubería, formando así, una
interfase entre el líquido y el gas, Figura 6.3 (a).
B. Flujo Ondulado. Este flujo es similar al flujo estratificado con la diferencia de

que el gas se mueve a una velocidad más alta y la interfase (gas – líquido) es
distorsionada por ondas en la dirección de flujo Figura 6.3 (b).
C. Flujo Anular. En este flujo, el líquido forma una película alrededor de la pared
interna de la tubería y el gas fluye a una velocidad mayor a lo largo de la tubería
Figura 6.3 (c).
162
D. Flujo Enchufe. En este flujo, burbujas de gas y líquido en forma alternada se

desplazan en la parte superior de la tubería Figura 6.3 (d).
E. Flujo Bache. En éste flujo, una onda es levantada periódicamente debido al

rápido desplazamiento del gas formando así, una ola espumosa la que se
desplaza a lo largo de la tubería, a una velocidad mucho mayor que la velocidad
promedio del líquido Figura 6.3 (e).
F. Flujo Burbuja. En este tipo de flujo, las burbujas de gas se desplazan a la parte
superior de la tubería aproximadamente a la misma velocidad que el líquido Figura
6.3 (f).
G. Flujo Niebla.
En éste flujo, casi la mayor parte del líquido se pulveriza formando gotas diminutas
de gas, desplazándose a la misma velocidad que el líquido restante Figura 6.3
(g).
163
Fig. 6. 3 Patrones de flujo observados en

tuberías horizontales
Los regímenes de flujo horizontal se pueden predecir con diagramas de régimen

de flujo. Baker (1953) propuso el siguiente diagrama, el cual fue modificado por
Scott (1963), mostrado en la Figura 6.4. Los ejes de esta gráfica son Gg/(λ) y Gl
(λ) Fi/Gg donde Gl y Gg son las masas de flujo de líquido y gas respectivamente
(lbm/hr-pie2) y los parámetros λ y ф son:
(6. 28)
(6. 46)
Donde las densidades están en lbm/pie3, µ en cp, y σ1 en dinas/cm. La región

sombreada en el diagrama indica que las transiciones de un régimen de flujo a
otro no son abruptas pero ocurren en esos rangos de condiciones de flujo.
164
Fig. 6.4 Diagrama de régimen de flujo
Otro diagrama de régimen de flujo comúnmente usado es el de Mandhane et.al.,

(1974), Figura 6.5. Como muchos diagramas de régimen de flujo vertical, éste usa
las velocidades superficiales de gas y líquido como las coordenadas.
La correlación de Beggs y Brill, está basada en un diagrama de régimen de flujo

horizontal que divide el dominio en tres categorías de régimen de flujo,
segregadas, intermitente y distribuido, este diagrama es mostrado en la Figura 6.5.
El número de Froude está definido por: NFr=um2/gD.
Fig. 6. 5 Diagrama de régimen de flujo
Taitel y Dukler (1976) desarrollaron un modelo teórico de transición del flujo de

régimen en un flujo horizontal Gas-Líquido; su modelo puede ser usado para
generar diagramas de flujo para fluidos particulares y tamaños de tubo. La figura
165
6.6, muestra una comparación de sus predicciones de régimen con las de

Mandhane et.al., para un flujo Aire-Agua en un tubo de 2.5cm.
Para flujo horizontal, el gradiente de presión debido al cambio de elevación es

igual a cero, por lo que:
(6. 29)
O bien:
(6. 48)
La mayoría de los investigadores han adoptado la ecuación anterior para evaluar

las características del flujo de dos fases y posteriormente determinar el gradiente
de presión total. El problema de la variación de las características de flujo se
elimina al suponer que la mezcla Gas - Líquido es homogénea en un intervalo
pequeño de la tubería. Así, la E. 6.31 se puede escribir como:
(6.49)
En donde ftp, y vm se refieren a la mezcla y son definidos en forma distinta por

los autores de las diferentes correlaciones.
6.3.1 Correlaciones para flujo horizontal
Poetmann, Bertuzzi y Tek
La ecuación para obtener el gradiente de presión por fricción es:
(6. 50)
De donde:
166
(6. 51)
Dónde:
[Adimensional]
Correlación de Eaton, Andrews, Knowels y Brown
(6. 52)
De donde:
(6. 53)
Correlación de Beggs y Brill
(6. 54)
Correlación de dukler
167
(6. 55)
Dónde:
(6. 56)
Dónde:
Ejemplo 6.1
Calcular la perdida de fricción mediante la ecuación de darcy para un fluido con

densidad de 50 °API en tuberías de diámetros de ½” , ¾” y 3 ½” con una
rugosidad de 0.0003, una velocidad de 5 m/s2 y determine en cuál de estas
tuberías existe la menor perdida por fricción.
Solución:
Con la ec. De Darcy
Para ½”
Para ¾”
168
Para 3 ½”
0179
Por lo tanto lo que causa la menor perdida por fricción es la tubería de 3 ½”
Ejemplo 6.2
Calcular la caída de presión por fricción en una tubería de 3000 pies de largo,
3.937 pulgadas de diámetro interno y con una rugosidad (ε) de 0.0006 pulgadas,
donde fluye aceite de densidad relativa de 0.9, μ = 46 (cp), si el gasto es:
a)2,560 bl/día
b)12,800 bl/día
Solución:
a)
Ya que el , el flujo es laminar y f se calcula con la ecuación:
b)
Ya que el flujo es turbulento .El factor de fricción ,según la

ecuación 5 se puede obtener con el siguiente procedimiento iterativo.
Iniciando con se obtienen los resultados que se muestran a

continuación:
ITERACION
1 0.03700
169
2 0.03576
3 0.03596
4 0.03592
5 0.03593
Por lo tanto se tiene como resultado :
Ejemplo 6.3
Determinar con el método de Poettman y Carpenter el gradiente de presión para

un pozo con flujo vertical que presenta las siguientes condiciones:
: 943.5 [bl/día] 153 °F
: 1122 pie3/bl 0
: 1.995 pg 0.85
: 960 lb/pg2 0.85
: 297.59 0.91229
: 1.1353 : 0.016784
: 0.55532 2.5278
: 0.93043 49.976
: 17.497 lbm/pie3
Solución:
De las ec. 7.5 y 7.4
Aplicando la ec. 7.2
170
Ejemplo 6.4
Con los datos del problema anterior, determinar el gradiente de presión con el
método de Orkizsewski.
: 943.5 [bl/día] 153 °F
: 1122 pie3/bl 0
: 1.995 pg 0.85
: 960 lb/pg2 0.85
: 297.59 0.91229
: 1.1353 : 0.016784
: 0.55532 2.5278
: 0.93043 49.976
: 49.976 1.0619
: 3.2054 : 8.6639
: 6.9605 18.8136
10.1659
Solucion:
De las ec. 7.10 , 7.16 y 7.42
Dado que : 0.6846 > -136.859 y 18.8136< 361.9004 se deduce que es flujo bache,
de las ec.
171
Con las ec. 7.18 y 7.17
De la ec. 7.40
Aplicando la ec. 7.1 obtenemos como resultado final
172
Capítulo 7 FLUJO A TRAVÉS DE ESTRANGULADORES
7. FLUJO A TRAVÉS DE ESTRANGULADORES
7.1 Estranguladores
Los estranguladores son restricciones en una línea de flujo que provoca una caída
de presión con el objetivo de controlar las condiciones de producción del pozo
.
Fig.7.1 Estrangulador
Clasificación de los Estranguladores
Diámetro fijo
Manual
Por el tipo de restricción
Diámetro variable Eléctrico
Hidráulico
Superficial
Por su ubicación en el pozo
Fondo
* Estranguladores submarinos
173
Las principales razones para instalar un estrangulador superficial en el pozo son

(Gilbert, 1954; Economides, 1994):

 Conservar la energía del yacimiento, asegurando una declinación
más lenta de su presión.
 Mantener una producción razonable.
 Proteger el equipo superficial.
 Mantener suficiente contrapresión para prevenir entrada de arena.
 Prevenir conificación de gas.
 Prevenir conificación de agua.
 Obtener el gasto de producción deseado.
7.2 Flujo en Estranguladores
Una vez que los fluidos producidos por el pozo llegan a la superficie, éstos pueden
o no pasar por alguna restricción denominada "estrangulador". La Figura 7.2
muestra esquemáticamente un pozo fluyente con un estrangulador instalado.
Figura 7.2. Estrangulador superficial
174
Las presiones presentes en un estrangulador son:
: Presión corriente abajo (presión en la línea de descarga)
: Presión corriente abajo (presión en la cabeza del pozo)
Cuando un gas o una mezcla de gas-líquido fluyen a través de un estrangulador,

el fluido es acelerado de tal manera que alcanza la velocidad del sonido en el
interior del estrangulador. Al ocurrir esto, el flujo es llamado “flujo crítico”, en
términos generales esta condición se alcanza cuando la presión de entrada al
estrangulador es aproximadamente el doble de la presión de salida.
Cuando se tiene flujo crítico (supersónico) en el estrangulador, las perturbaciones

de presión corriente abajo del estrangulador no afectan a los componentes que
están corriente arriba. Esto es, un cambio en la presión corriente abajo, por
ejemplo, la presión en el separador no afectará el gasto o presión corriente arriba
(presión en la cabeza del pozo).
El flujo crítico ocurre cuando:
(7.1)
Varios autores han desarrollado expresiones para calcular la caída de presión en

el estrangulador. Las expresiones propuestas por Gilbert (1954), Ros (1960),
Baxendell (1961) y Achong (1974) tienen la misma forma, sólo difieren en los
valores de las constantes empleadas por cada uno de ellos.
Varios métodos han sido desarrollados para describir el comportamiento del flujo
multifásico a través de estranguladores. Entre ellos podemos mencionar los
siguientes:
1 Gilbert.
2 Ros.
3 Baxendell.
4 Achong.
5 Poettmann y Beck.
6 Ashford.
175
7 Ashford-Pierce.
8 Omaña.
9 Pilehvari etc.
Gilbert, Ros, Baxendell y Achong.
A partir de datos de producción, Gilbert (1954) desarrolló una expresión aplicable

al flujo simultáneo Gas - Líquido a través de estranguladores. En su trabajo
describe en forma detallada el papel del estrangulador en un pozo y analiza cuál
es el efecto sobre la producción de cambios bruscos en el diámetro del orificio.
Tomando como base la relación entre las presiones antes y después de un orificio
para flujo sónico de una fase, Gilbert recomendó para tener flujo critico (sónico)
una relación de 0.588 o menor, entre la presión promedio en el sistema de
recolección (después del estrangulador) y la presión en la boca del pozo (antes del
estrangulador).
Utilizando datos adicionales, Baxandell (1961) actualizó la ecuación de Gilbert,

modificando los coeficientes de su ecuación.
Ros (1960) orientó su trabajo al flujo de mezclas con alta relación gas - aceite en
las que el gas fue la fase continua. En su desarrollo llegó a una expresión similar a
Gilbert pero con coeficientes diferentes.
Achong (1974) también revisó la ecuación de Gilbert y estableció una expresión

que validó comparándola con más de cien pruebas de campo.
Las ecuaciones propuestas por Gilbert (1954), Ros (1960), Baxandell (1961), y
Achong (1974) tienen la misma forma, sólo difieren en los valores de las
constantes empleadas. La forma general es la siguiente:
(7.2)
Donde:
: Coeficientes empíricos [adimensionales]
: Diámetro interior del estrangulador [ 64 avos de pulgada]
176
: Presión corriente arriba
: Gasto del liquido total
: Relación gas – liquido
La tabla 7.1 muestra los valores obtenidos por los diferentes autores para dicha
ecuación, coeficientes empíricos obtenidos de diversos experimentos
INVESTIGADOR A (adimensional) B (adimensional) C (adimensional)

Ros (1954) 17.40 0.500 1.89
Gilbert (1960) 10.00 0.546 2.00
Baxandell (1961) 9.56 0.546 1.93
Achong (1974) 3.82 0.650 1.88
Tabla 7.1 Valores de los coeficientes para la ecuación 7.2
La ecuación solo describe el flujo a través del estrangulador si se tiene un flujo

supersónico, como es inversamente proporcional a la pendiente de la recta, ésta
se hace menor al aumentar el diámetro Figura 7.3.
Fig. 7.3 Variación de la pendiente de la línea recta con el diámetro del

estrangulador
177
Poettmann y Beck. (P y B)
Este modelo fue establecido a partir del trabajo presentado por Ros. La precisión
de los resultados obtenidos se comprobó comparándolos con 108 datos medidos.
El método fue establecido a partir de un análisis teórico del flujo simultáneo Gas -
Líquido a velocidad sónica a través de orificios y una correlación para el
comportamiento PVT de los fluidos. No se consideró producción de agua.
Para que exista flujo crítico se supuso que la presión corriente abajo, debe ser al
menos 0.55 de la presión en la boca del pozo. Bajo estas condiciones el gasto
(E.8) en el estrangulador es sólo función de la presión corriente arriba y de la
relación gas-aceite a condiciones de flujo:
(7.3)
Donde:
(7.4)
(7.5)
(7.6)
Siendo:
: Relación gas libre – aceite a condiciones de flujo
: Volumen especifico del líquido
178
: Masa de líquido por unidad de masa de mezcla
Ashford
A partir de un balance de energía y considerando que el fluido se expande

politrópicamente al pasar por el estrangulador, Ashford derivó una ecuación que
describe el flujo multifásico, bajo condiciones sónicas, a través de un orificio.
Para compensar la ecuación por las suposiciones incluidas en su desarrollo, se

introdujo en ella un coeficiente de descarga. Sin embargo, al evaluarla,
comparando sus resultados con datos medidos en 14 pozos, se encontró que el
coeficiente de descarga resultaba muy cercano a la unidad.
(7.7)
Donde:
(7.8)
(7.9)
En su derivación Ashford supuso una relación de calores específicos k = 1.04 y

una relación de presiones, para obtener flujo sónico en el orificio de 0.544.
Ashford y Pierce
Establecieron una ecuación que describe la dinámica de las caídas de presión y

capacidades de flujo en condiciones de flujo multifásico. Este modelo relaciona el
comportamiento del estrangulador en ambos regímenes de flujo: crítico y no
crítico.
La capacidad y caídas de presión que se presentan en la restricción se han

relacionado con sus dimensiones y las propiedades de los fluidos manejados.
Para validación del modelo, se diseñó una prueba de campo en un pozo fluyente.
179
Tanto las caídas de presión como el gasto se midieron directamente y luego se

compararon con datos análogos obtenidos del modelo. Esta información se usó
para determinar “el coeficiente de descarga del orificio (c)”, definido por la relación
de gasto medido entre el gasto calculado.
Los resultados obtenidos en las pruebas, para diámetros de estrangulador de

4/64, 16/64 y 20/64 de pg son:
DIAMETRO DE ESTRANGULADOR
C
(1/64 pg)
14 1.1511
16 1.0564
20 0.976
Para diámetros (d) menores de 20.81/64 pg puede aproximarse el coeficiente de

descarga con la siguiente ecuación, que es el resultado del ajuste de la relación
entre el diámetro del estrangulador y “c”:
(7.10)
Para valores mayores, el valor de “c” es constante e igual a 0.95.
Omaña
Desarrolló una correlación (para flujo crítico) entre el gasto, la presión corriente
arriba del estrangulador, la relación Gas - Líquido, la densidad de los líquidos y el
diámetro del orificio. Dicha correlación se obtuvo a partir de datos experimentales.
En vista de que estos datos estuvieron dentro de rangos muy limitados, su
aplicación sólo se recomienda para orificios hasta de 14/64 pg y gastos máximos
de 800 bl/d. Las condiciones de flujo crítico se fijaron para una relación de
presiones igual o menor de 0.546 y una relación Gas - Líquido mayor de 1.0.
7.3 Flujo de Gas a Través de Estranguladores
Para la adecuada interpretación del comportamiento de un pozo fluyente, se hace

necesario tener un buen conocimiento del efecto que sobre éste tiene la
colocación de estranguladores superficiales. Prácticamente todos los pozos
fluyentes utilizan algún tipo de restricción, lográndose por este medio regular el
gasto de producción. Muy pocos pozos producen sin ningún tipo de restricción,
180
encontrándose por lo tanto, produciendo al gasto máximo que sus condiciones le

permiten.
El hecho de regular la presión en la cabeza del pozo puede requerirse por alguna
de las siguientes razones:

 Ejercer la contrapresión suficiente para evitar la entrada de
arena
 Conificación de agua y/o gas hacia el pozo.
 Conservar la energía del yacimiento, asegurando una
declinación más lenta de su presión.
 Proteger el equipo superficial y subsuperficial.
 Producir el yacimiento a un gasto adecuado.
Los estranguladores que se instalan en la boca del pozo, para controlar la

producción, están basados en el principio de flujo crítico.
El flujo crítico es un fenómeno de flujo definido por el flujo de gases compresibles

en la sección de estrangulamiento de una restricción, cuando su velocidad es
sónica (velocidad del sonido en el fluido) o el número de Mach es uno. El número
de Mach (M) es una relación adimensional dada por el cociente de la velocidad del
fluido entre la velocidad del sonido:
(7.11)
Para:
ó í
Existe flujo subcrítico sí:
(7.12)
Y se tiene flujo crítico cuando:
181
(7.13)
Donde:
(7.14)
: Calor específico a presión constante
: Calor específico a volumen constante
El diámetro del estrangulador se puede obtener con:
(7.15)
Donde:
(7.16)
para p = 14.7 psi y To = 60 ªF, el coeficiente de descarga Cd = 0.00019083.
Cook y Dotterweich por su parte, plantearon la siguiente ecuación para calcular el

diámetro del estrangulador:
(7.17)
Donde :
Para el coeficiente de descarga se obtiene con la siguiente

ecuación:
(7.18)
182
para valores de do superiores a 32/64 pg el coeficiente de descarga es constante e

igual a 0.828; finalmente el diámetro del estrangulador en 1/64 de pg, se calcula
mediante la siguiente ecuación:
(7.19)
Cuando se tiene flujo crítico a través del estrangulador, la presión corriente arriba
es independiente de la presión que prevalece en el sistema de recolección (línea
de descarga, separadores, etc.). Se infiere que el control de la producción se
logrará cuando las variaciones de presión en el sistema de recolección no se
reflejan en la formación productora, provocando fluctuaciones en la producción.
Esta situación prevalecerá al usar un estrangulador que permita obtener la
producción deseada bajo condiciones de flujo crítico.
Ejemplo 7.1
Dada la siguiente información, calcular el diámetro del estrangulador en 64 avos

de pulgada empleando las correlaciones de Gilbert y Ros
qo = 1800 bpd RSL = 335 pcn/bn

qg = 1.116 MMpcd Bo = 1.10 by/bn
dt= 2.994 pg = o.249 pies Z1 =0.78
Ap= 0.0487 pies2 gl = 6.856 lbm/pie3
C=1 0l = 52.59 lbm/pie3
API = 30, Υo = 0.876 T = 128 °F
Υg = 0.70 P1 = 1663 lb/pg2
σL1 = 9 dinas/cm fo=62.4
Solución:
Despejando la ec. 8.2 despejamos d
183
De la tabla 8.1
Para Gilbert
A = 10.0
B = 0.546
C = 1.89
Para Ros
A = 17.40
B = 0.500
C = 2.00
Ejemplo 7.2
Determinar el diámetro del estrangulador que se requiere para calcular la

producción máxima de un pozo mediante la correlación de Achong y Baxandell,
cuya información es:
ql = 480 bpd
p1= 120 lb/pg2
Rp= 300 pcn/bn
184
Solución:
Despejando la ec. 8.2 despejamos d
Para Baxandell
A = 9.56
B = 0.546
C = 1.93
Para Baxandell
A = 3.82
B = 0.650
C = 1.88
Ejemplo 7.3
Con los datos del ejemplo 8.1 determine el diámetro del estrangulador mediante
la correlación de Poettmann y Beck.
Solución:
Calculando los para parámetros para la correlación de Poettmann
185
Despejando Ao de la correlación de Poettmann.
d= 0.324 pg = 20.73 64 avos de pg
186
Ejemplo 7.4
Dada la siguiente información, calcular el diámetro del estrangulador en 64 avos

de pulgada empleando las correlaciones de Ashford
qo = 1100bpd RSL = 400 pcn/bn

C=1 Bo = 1.05 by/bn
API = 30, Υo = 0.876 Z1 =0.7
Υg = 0.70 ρgl = 6.3 lbm/pie3
σL1 = 9 dinas/cm ρ0l = 59.6 lbm/pie3
T = 100 °F P1 = 1500 lb/pg2
Solución:
Despejando d de la correlación de Ashford tenemos
γ γ
Calculando los parámetros
Sustituyendo los parámetros
d=14.23/64 pg
187
Ejemplo 7.5
Dada la siguiente información determine el diámetro optimo del estrangulador ,

mediante el método de Omaña
qo = 300 bpd
σL1 = 9 dinas/cm
vsg=1.2 pie/seg
vsl=2.5 pie/seg
ρgl = 7 lbm/pie3
ρol = 60 lbm/pie3
P1 = 1000lb/pg2
Solución:
Calculando los parámetros para Omaña
Sustituyendo valores en la ec…… y despejando para do tenemos
188
189
190
SISTEMA INTEGRAL DE PRODUCCIÓN,
Capítulo 8 ANÁLISIS NODAL
8. SISTEMA INTEGRAL DE PRODUCCIÓN, ANÁLISIS NODAL
8.1 Sistema Integral de Producción
Básicamente, un sistema integral de producción es un conjunto de elementos que

transporta los fluidos del yacimiento hacia la superficie, los separa en aceite, gas y
agua, y finalmente los envía a instalaciones para su almacenamiento y/o
comercialización. Asimismo, un sistema integral de producción puede ser
relativamente simple o puede incluir muchos componentes.
Los componentes básicos de un sistema integral de producción son (Golan y

Whitson, 1991):
 Yacimiento
 Pozo
 Tubería de descarga
 Estrangulador
 Separadores y equipo de procesamiento
 Tanque de almacenamiento
191
Fig.8.1 Sistema integral de producción.
Yacimiento
Se entiende por yacimiento, la porción de una trampa geológica que contiene

hidrocarburos, la que se comporta como un sistema intercomunicado
hidráulicamente. Los hidrocarburos que ocupan los poros o huecos de la roca
almacén, se encuentran a alta presión y temperatura, debido a la profundidad que
se encuentra la zona productora.
Pozo
Es un agujero que se hace a través de la roca hasta llegar al yacimiento; en este

agujero se instalan sistemas de tuberías y otros elementos, con el fin de
establecer un flujo de fluidos controlados entre la formación productora y la
superficie.
Tubería de descarga
Las tuberías son estructuras de acero, cuya finalidad son transportar el gas, aceite
y en algunos casos, agua desde la cabeza del pozo hasta el tanque de
almacenamiento. Los costos específicos en el transporte, tanto de aceite como de
gas, disminuyen cuando la capacidad de manejo aumenta; esto se logra si el
aceite, gas y agua se transportan en tuberías de diámetro óptimo, para una
capacidad dada.
Estrangulador
Es un aditamento que se instala en los pozos productores, con el fin de establecer

una restricción al flujo de fluidos. Es decir, permite obtener un gasto deseado,
además de prevenir la conificación de agua, producción de arena y sobre todo,
ofrecer seguridad a las instalaciones superficiales.
Separadores
Los separadores como su nombre lo indica, son equipos utilizados para separar la
mezcla de aceite y gas, y en algunos casos aceite, gas y agua que proviene
directamente de los pozos. Los separadores pueden clasificarse por su forma o
192
geometría en horizontales, verticales y esféricos; y por su finalidad, separar dos

fases (gas y líquido) o tres (gas, aceite y agua).
Tanques de almacenamiento
Son recipientes de gran capacidad para almacenar la producción de fluidos de uno

o varios pozos. Los tanques de almacenamiento pueden ser estructuras cilíndricas
de acero instalados en tierra firme, o bien, buque - tanques, usualmente utilizados
en pozos localizados costa afuera. En la industria petrolera, los tanques pueden
tener una capacidad de almacenamiento que va desde 100,000 hasta 500,000
barriles.
Para analizar el comportamiento de un pozo fluyente terminado, es necesario

analizar las tres áreas de flujo, las que se tienen que estudiar en forma separada y
unirlas después, antes de obtener una idea precisa del comportamiento de flujo
del pozo productor.
Estas áreas de flujo son (Nind, 1964):
 Flujo del yacimiento al pozo

 Flujo en tuberías
 Flujo en estrangulador
8.2 Análisis Nodal
La razón fundamental de someter un sistema de producción a la técnica de

Análisis Nodal, es simplemente porque ésta involucra en sus cálculos a todos los
elementos del sistema, permite determinar el efecto de su variación en la
capacidad de transporte y tener una imagen de 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 sistema de
producción en su conjunto para obtener la capacidad de transporte requerida, o
bien, la capacidad de transporte idónea, teniendo como limitante tan sólo la
capacidad de afluencia del yacimiento al pozo.
193
El Análisis Nodal se puede aplicar a pozos fluyentes, inyectores o productores ya

sea fluyentes o con algún sistema artificial de producción. Por otra parte, del
análisis del comportamiento de los elementos del sistema se pueden identificar las
posibles restricciones que modifiquen negativamente la capacidad de transporte
del mismo. Asimismo, es posible estudiar y comprender, con relativa facilidad, el
comportamiento de cada uno de los componentes del sistema integral de
producción (desde el yacimiento hasta la batería de recolección).
Análisis cualitativo del sistema
En la Figura 8.2 se muestra un sistema de producción, en el que se aprecian las

posibles pérdidas de presión desde el yacimiento hasta los separadores.
Fig.8.2 Pérdidas de presión en un sistema completo de producción.
p 1: Pws – Pwsf Pérdidas de presión en el yacimiento.

p 2: Pwsf – Pwf Pérdidas de presión en el radio de daño.
p 3: Pur – Pdr Pérdidas de presión por restricciones en la TP.
194
p 4: Pusv – Pdsv Pérdidas de presión en la válvula de seguridad.

p 5: Pth – Pe Pérdidas de presión en el estrangulador superficial.
p 6: Pe – Ps Pérdidas de presión en la línea de flujo.
p 7: Pwf – Pth Pérdidas de presión totales en la TP.
p 8: Pth – Ps Pérdidas de presión en el estrangulador y la LD.
En la figura 8.3 se muestran en el mismo sistema de producción los nodos

principales para su estudio.
Fig. 8.3 Pérdidas de presión y nodos principales en un sistema

básico de producción.
p 1: pérdidas de presión en el medio poroso. Representan entre el 10 y el 50% de

las pérdidas totales.
p 2: pérdidas de presión en la tubería vertical. Representan entre el 30 y el 80%

de las pérdidas totales.
195
p 3: pérdidas de presión en la línea de descarga. Generalmente, constituyen

entre el 5 y el 30% de las pérdidas totales.
En la Figura 8.4 se muestran las pérdidas de presión asociadas a cada elemento

de la Figura 8.3.
Fig. 8.4 Distribución de presión en un sistema de producción
Flujo en el yacimiento
Partiendo del análisis de la Figura 8.4 se puede distinguir un principio, que a

medida que el gasto se incrementa, la diferencia entre la presión estática y la
presión de fondo fluyendo se acentúa. Esta diferencia depende de las
características del sistema roca-fluidos y de la eficiencia de la terminación.
Comportamiento del flujo por la TP
Nind, para relacionar el gasto de producción con la RGL supone constante esta
última, mientras varía el gasto de líquido y muestra que:”para cualquier diámetro y
196
profundidad de tubería dados, hay un gasto de producción que ocasiona la mínima

pérdida de presión en la tubería de producción”.
Es de esperarse, sin embargo, que la declinación de la presión del yacimiento

permitirá un incremento de la RGL, que en principio beneficiará la producción de
fluidos, pero después su continuo aumento, podría llegar a producir pérdidas por
fricción paulatinamente mayores.
Fig. 8.5 Pérdidas de presión en función del gasto para varias RGL
(tomado de Gillbert, W.E.)
La Figura 8.6 muestra los gastos que limitan el flujo estable. Un pozo que
produzca con un gasto menor generalmente estará operando con “cabeceo” o flujo
inestable. Las condiciones para obtener flujo estable deben ser tales, que al
agregar a la curva anterior la curva IPR, se obtenga un resultado semejante al de
la Figura 8.7.
197
Fig. 8.6 Curvas típicas de gasto vs Pwf para diferentes diámetros de TP
Fig. 8.7 Gráfica típica de condiciones de flujo estable.
Terminación del flujo natural
La determinación de la presión estática Pws a la que el pozo dejará de fluir es una

aplicación importante del flujo multifásico vertical. El procedimiento consiste en
graficar los valores de la Pwf obtenidos a partir del comportamiento del flujo en el
yacimiento y del flujo por la TP, Figura 9.7. Las curvas de IPR corresponden a
198
presiones estáticas de 1,200 y 1,300 lb/pg2. El pozo tiene una TP de 3 1/2 pg,
Pwh = 100 lb/pg2. El pozo no fluirá a una Pws < 120 lb/pg2. A una Pws = 1,150 lb
/pg2 el pozo estará “muerto”.
Fig. 8.8 Curvas de comportamiento de flujo

correspondiente a un pozo “muerto”
Efecto del diámetro de la TP
Nind, demostró el efecto del cambio del diámetro de la TP, sobre el gasto de
producción y la presión de fondo fluyendo. En sus cálculos consideró un pozo de
10,000 pies de profundidad y Pwh = 0 lb/pg2.
Sus resultados se muestran en las Figuras 8.9 y 8.10. En ellas se distingue que a
gastos bajos se reducen las pérdidas de presión al utilizar diámetros menores de
TP (Figura 8.9). En relación a la Figura 8.10 se observa que para diámetros de TP
pequeños, aumentan las caídas de presión al incrementar el gasto, pero esta
situación se invierte para diámetros de TP mayores.
Fig.8.9 Determinación de la presión estática a la que el pozo deja 199

de fluir.
Fig. 8.10 Efecto del gasto sobre las pérdidas de presión por flujo
vertical: diferentes diámetros de tubería de producción.
Se puede observar de la ecuación desarrollada por Gilbert (Figura 8.11), para un

diámetro dado de estrangulador y una R constante, es una línea recta que pasa
por el origen, Figura 8.4. Suponiendo un gasto muy pequeño, la presión en la
cabeza y la presión corriente abajo tenderían a igualarse a la presión en el
separador. Al fluir el pozo, el comportamiento del estrangulador sería semejante al
mostrado en la Figura 8.12.
Fig.8.11 Curvas de afluencia para pozos con E.F. diferente a 1 de yacimientos

con empuje por gas disuelto.
200
Fig. 8.12 Efecto del diámetro de la TP sobre las pérdidas de presión

en flujo vertical con diferentes gastos
Comportamiento del flujo por la LDD
La caída de presión en la LD, se debe fundamentalmente a la fricción del fluido

con las paredes de la tubería y sus características más significativas se reflejan en
las figuras 8.13, 8.14 y 8.15.
Fig.8.13 Gráfica de Pth vs q para diferentes diámetros de estrangulador.
201
Fig. 8.14 Variación de las DPLD vs Gasto para diferentes

diámetros de la LDD
Fig. 8.15 Variación de las DPLD vs Gasto para diferentes

relaciones gas-aceite.
Para la selección del diámetro óptimo de la línea de descarga, 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 (el yacimiento, la TP, el estrangulador o bien la presión del
separador) limita el ritmo de producción.
202
8.3 Optimización de un Sistema de Producción
En conclusión, puede afirmarse que la aplicación de la técnica nodal a un sistema

de producción, permite identificar los elementos que limitan la capacidad de flujo
del sistema; es decir, que el Análisis Nodal se emplea para diagnosticar la
variación del gasto de producción al realizar alguna de las modificaciones
siguientes:
a) Presión de separación.
b) Eliminar o sustituir válvulas o conexiones inapropiadas.
c) Colocar separadores a boca del pozo. En este caso se pueden analizar
dos opciones:
c1) Separar con la presión necesaria para transportar el líquido
(aceite + agua) hasta la central de recolección para continuar con su
proceso.
c2) Separar a baja presión y bombear el aceite hasta la central de
recolección para continuar con su proceso.
d) Cambiar diámetro de la TP.
e) Cambiar diámetro de la LD o construir una adicional.
f) Instalar un sistema artificial de producción.
Es evidente que la selección de las modificaciones a un sistema y el orden de su

aplicación deben basarse en un análisis económico, en el que se comparan los
incrementos en la producción, al efectuar algún cambio, con la inversión adicional
que sea necesario realizar.
El Análisis Nodal involucra en sus cálculos a todos los elementos del sistema,
permite determinar el efecto de su variación en la capacidad de transporte y tener
una imagen de conjunto del comportamiento del pozo. El Análisis Nodal se puede
aplicar a pozos fluyentes, inyectores o productores ya sea fluyentes o con algún
sistema artificial de producción. Por otra parte, del análisis del comportamiento de
los elementos del sistema se pueden identificar las posibles restricciones que
modifiquen negativamente la capacidad de transporte del mismo. Asimismo, es
posible estudiar y comprender, con relativa facilidad, el comportamiento de cada
uno de los componentes del sistema integral de producción (desde el yacimiento
hasta la batería de recolección).
203
La elección del nodo de 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, tomando en cuenta una justificación
económica, para su posterior ejecución.
EJEMPLO 8.1 (Primer método Pwf):
Calcular el gasto máximo que puede fluir por una tubería de producción de 2´´ D.I.
de un pozo del cual se tienen los siguientes datos:
Profundidad de la T.P. = 5000 pies
¿Cuál sería el gasto para Pth = 120 [lb/pg2]?
Pws = 2000 lb/pg2
qo = 350 bl/dia
Pwf = 1250 lb/pg2
R = 300 pie3/bl
Pwf < Pb
Solución:
Se traza la curva de IPR, utilizando la curva de referencia de Vogel.
Pwf lb/pg2 2000 1750 1500 1250 1000 750 500 250 0
qo bl/día 0 130 248 380 434 505 558 601 620
204
qo Prof. Equivalente Prof. Equivalente Pwf

a Pth a Pwf
(bpd) (ft) (ft) (lb/pg2)
100 1900 6900 935
200 1700 6700 930
300 1600 6600 980
400 1400 6400 1000
500 1250 6250 1020
600 1150 6150 1050
700 1000 6000 1140
Utilizando las curvas de gradiente de presión, con gasto y diámetro de tubería

conocidos, se determina la Pwf correspondiente para Pth = 120 [lb/pg2].
El proceso se presenta en la siguiente tabla:
205
Prof. Equivalente a Pwf=Prof. Equivalente a Pth + Prof de la TP
La intersección de la curva de Pwf con la curva de IPR, determina el gasto

máximo que puede fluir hasta la superficie y la Pwf necesaria para vencer una
contrapresión de 120 lb/pg2 en la boca del pozo.
Ejemplo 8.2
Se calcula Pth para gastos que coincidan con las curvas de gradiente conocidas a
partir de la Pwf correspondientes a cada gasto obtenido de la curva de IPR. Los
pasos a seguir se presentan en la siguiente tabla:
a Pwf a Pth
(bpd) (lb/pg2) (ft) (ft) (lb/pg2)
100 1800 9400 4400 400
200 1600 8700 3700 330
300 1370 7800 2800 250
400 110 6750 1750 160
500 770 5300 300 25
* Pwf se obtiene de la curva de IPR

** Prof. equivalente a Pth = Prof. Equivalente a Pwf - Prof. de la TP
206
a intersección de la curva de Pwf con la curva de IPR, determina el gasto máximo

que puede fluir hasta la superficie y la Pwf necesaria para vencer una
contrapresión de 120 lb/pg2 en la boca del pozo.
Considerando los siguientes datos:
q0 = 1000 bl/día (100% aceite)
R = 1000 pie³/bl
Psep = 100 lb/pg²
Longitud de la tubería de descarga = 5000 piesCalcular la presión de flujo en la

boca del pozo Pwh para un diámetro de tubería de descarga de 3 pg.
Solución:
Diámetro de la tubería de descarga = 3 pg. Entrando en la abscisa con un valor de

presión de 100 lb/pg2 y bajando verticalmente hasta intersectar con la curva de
R = 1000 pie³/bl.
Se lee en el eje de las ordenadas una longitud equivalente de 2000 pies, la que
sumada a la longitud de la tubería de descarga es ahora 7000 pies.
Con esta profundidad y moviéndose ahora horizontalmente hasta intersectar la

curva de R = 1000 pie³/bl y subiendo hasta el eje de las presiones se lee un valor
de 220 psi que es la presión solicitada.
207
Ejemplo 8. 4
Considerando los siguientes datos:
diámetro de la línea de descarga = 3 pg
Psep = 160 lb/pg2
Pwh = 600 lb/pg2
Longitud de la tubería de descarga = 4500 pies
R = 5000 pies³/bl
Encontrar el gasto posible de flujo a través de la tubería de descarga.
Solución:
1.- Suponer varios gastos.
2.- Determinar la presión en la cabeza del pozo para cada gasto, y preparar una
tabla como la siguiente:
Gastos supuestos (bpd) Pwh (psi)
3000 415
4000 565
5000 665
3.- Elaborar una gráfica de q vs Pwh.
4.- Con una presión en la cabeza de 600 lb/pg2 se intersecta verticalmente la

curva construida.
5.- Se lee horizontalmente hacia la izquierda hasta intersectar el eje de las

abscisas y se lee el gasto de producción real de 4350 bl/día.
208
700
600
500
Pwh lb/pg2
400
300
200
100
0
0 1000 2000 3000 4000 5000 6000
qo bpd
209
210
Capítulo 9 DISEÑO DEL APAREJO DE PRODUCCIÓN
9 DISEÑO DEL APAREJO DE PRODUCCIÓN
El principal objetivo de un aparejo de producción es conducir los fluidos

producidos por el yacimiento (gas, aceite, agua o ambos) hacia la superficie, y los
de inyección hacia la formación en forma controlada.
Además, debe soportar los esfuerzos a los que estará sometido durante la
terminación del pozo (inducciones, pruebas de admisión, estimulaciones,
fracturamientos, fluidos corrosivos, etc.) y su vida productiva
9.1 Opciones de Terminación de un Pozo
La zona del yacimiento alrededor del pozo no debería perturbarse por las
operaciones de la perforación. Muchos pozos son cementados con tubería de
revestimiento para mantener la misma en seguridad del pozo y también para aislar
el yacimiento. En caso de no haber aislamiento total, se puede hacer una
terminación en agüero descubierto, siempre y cuando la formación sea lo
suficientemente consolidada. Se puede usar la terminación de tubería ranurada,
siempre y cuando la estabilidad del pozo lo permita, accede a no afectar al
yacimiento con la cementación y tiene mayor área de comunicación del yacimiento
con el pozo.
Para combatir el arenamiento del pozo y dañar al sistema y sus componentes, se

usan los empacadores de grava que evitarán la incorporación de arena y granos
finos a las tuberías ocasionando erosiones, sobre todo si las arenas son
cuarcíferas; se usan mucho en sedimentos poco consolidados y tienen la función
de preservar la permeabilidad e inclusive reducir la caída de presión por efecto de
los empacadores.
Pozo fluyendo con T.P franca.- En la figura 9.1se muestra la T.P que va colgada y
situada a determinada profundidad sobre el intervalo productor. Los fluidos pueden
ser explotados por dentro y fuera de la T.P, aunque no se recomienda que se
produzca por espacio anular debido a que la T.R seria expuesta a danos por
fricción y corrosión.
211
Fig 9.1 pozo fluyendo con TP

franca
Después de que la zona productora sea perforada, una tubería de revestimiento o

liner es introducida y cementada. Posteriormente es perforada a través de pistolas
las cuales crean un conducto del yacimiento al pozo. Estas perforaciones deben
atravesar la tubería de revestimiento, el cemento y la zona invadida del fluido de
perforación, para que el flujo de hidrocarburos no pase por una zona dañada.
Fig 9.2 Terminación en

agujero revestido
212
Después de haber perforado el intervalo productor se introduce una tubería

ranurada o liner ranurado que se ancla por medio de un empacador cerca de la
zapata de la tubería de revestimiento que por lo general se encuentra en la cima
del intervalo productor. Esta tubería no es cementada.
Fig 9.3 Terminación en

tubería ranurada sin cementar
9.2 Clasificación de los Aparejos de Producción

De acuerdo a la función que realizan, los aparejos de producción, se clasifican en:
Bombeo Neumático Continuo (BNC)

En la figura 9.4 se observa el BNC, donde se introduce un volumen de gas a alta
presión de forma constante por el espacio anular, donde mediante una válvula de
fondo el gas pasara a la
T.P. La columna de fluido por encima del punto de inyección es aligerada por
aireación causada por la baja densidad del gas inyectado. Esta disminución del
peso de la columna de fluido permita la presión de fondo un incremento en la
producción.
213
Fig 9.4 Bombeo Neumático

Continuo.
Bombeo Neumático Intermitente (BNI)

En la figura 1.15 se observa el BNI, donde la producción del aceite es periódica,
esto es debido a que el gas que es inyectado periódicamente en múltiples puntos
o bien en un solo punto.
Un regulador en la superficie controla los tiempos de cada ciclo de inyeccion-

producciòn.
Usualmente se utiliza en pozos con bajos volúmenes de fluido, un alto IP y bajas

presiones de fondo, o bien, bajo IP y altas presiones de fondo. Debido a esto las
válvulas serán abiertas para inyectar cierta cantidad de gas a distintos tiempos,
donde la primera válvula será la que se encuentre por debajo de la columna o
bache fluido, el fluido aligerado llegara a la siguiente válvula ubicada más arriba de
la primera válvula, donde también será inyectado gas, que levantara aún más este
fluido; esto se repite con las demás válvulas hasta llevar el fluido a la superficie.
214
Fig 9.5 Bombeo Neumático

Intermitente.
Sistema de Bombeo Hidráulico (BH)

El principio de funcionamiento del sistema de BH se basa en transmitir potencia
mediante el uso de un fluido presurizado a través de la tubería, este fluido es
conocido como fluido de potencia o fluido motriz, este fluido es utilizado por una
bomba subsuperficial, la cual transforma la energía del fluido en energía potencial
o de presión para el fluido producido, el cual será enviado a la superficie. En la
figura 9.6 se ilustra el sistema de BH.
Fig 9.6 Sistema de Bombeo

Hidráulico (BH)
215
Sistema de Bombeo Electrocentrífugo (BEC)

Un sistema BEC estándar consiste de instalaciones subsuperficiales como una
bomba electrocentrífuga de etapas múltiples, intake y/o separador de gas,
protector, motor eléctrico y cable de potencia. En sus instalaciones superficiales
tiene un transformador, variador de frecuencia, caja de venteo y conexiones
superficiales. También van incluidos todos aquellos accesorios que aseguran una
buena operación como lo son: flejes de cable, extensión de la mufa, válvula de
drene, válvula de contra presión, centradores, sensor de presión y temperatura de
fondo, dispositivos electrónicos para control del motor, caja de unión y controlador
de velocidad variable. En la figura 9.7 se muestra el sistema BEC.

Electrocentrífugo (BEC)
Sistema de Bombeo Mecánico

La figura 9.8 muestra el sistema de bombeo mecánico, que básicamente consiste
en instalar en el fondo de la T.P una bomba subsuperficial, la cual succiona aceite
debido al movimiento reciprocante de un embolo, el cual se desplaza en forma
ascendente y descendente en el interior de la bomba al ser puesto en operación
desde la superficie por medio de un mecanismo conocido como unidad de bombeo
mecánico, siendo accionado por la energía proporcionada de un motor eléctrico o
de combustión interna, transmitiendo esta energía hasta el embolo a través de una
sarta de varillas metálicas, las cuales van a unir a la unidad de bombeo mecánico
con la bomba subsuperficial, siendo indispensable que la bomba se encuentre
completamente sumergida en el fluido del pozo.
216

Mecánico (BEC)
9.3 Selección del Diámetro del Aparejo de Producción
El diámetro del aparejo de producción se realiza mediante un análisis nodal, que

estudia simultáneamente el comportamiento de flujo del yacimiento al pozo y el
IPR. El análisis nodal, relaciona las diferentes pérdidas de presión en el sistema, y
en particular para este caso en el pozo, con la finalidad de optimizar la producción,
para su análisis se selecciona el subsistema del pozo con el del flujo del
yacimiento al pozo, en donde se encuentra la unión de estos subsistemas en el
punto denominado nodo funcional.
La Figura 9.3, muestra un sistema de producción y los nodos más representativos

y utilizados para el análisis del sistema mencionado:
217
Fig. 9.9 Caídas de presión en el aparejo de producción.
Las caídas de presión en el aparejo de producción es resultado de los siguientes

efectos:
Elevación
Siendo el componente más importante, representa el 90% de la caída de presión

total y es una función de la densidad de la mezcla (agua, aceite y gas).
Fricción
Depende del tamaño y rugosidad de la tubería y es, principalmente, una función

de la viscosidad del fluido.
218
Aceleración
Es el menos importante, en tuberías verticales, y es función del cambio de

velocidad en la tubería.
El diámetro del aparejo debe ser tal que permita transportar los fluidos producidos
con los gastos esperados, ya que si es muy pequeño restringirá la producción y si
es muy grande se tendrían flujos inestables e intermitentes.
Fig. 9.10 Curvas típicas de IPR y comportamiento de flujo en

aparejos de distintos diámetros.
La Figura 9.4 muestra un análisis de sensibilidad, con el que se determina el

diámetro óptimo del aparejo de producción, observando que con un diámetro de
31/2 pg, se obtiene el máximo gasto. Las curvas tienen tres características
importantes:
 La sección de pendiente negativa a gastos bajos es una indicación de flujo

inestable en la tubería.
 El punto de inflexión de la curva indica el gasto crítico por debajo del cual el
gas se separa del líquido y el pozo se cargará, cualquier gasto por abajo del
gasto crítico mataría el pozo.
 La sección de pendiente positiva de la curva, muestra los gastos y

presiones en los que la velocidad es lo suficientemente alta para llevar los
fluidos a superficie.
219
9.4 Empacadores de Grava
El control de arena es posible con empacadores de arena que se emplean en las

terminaciones del pozo, Suman, et al (1983), hacen un resumen de las mejores
prácticas de control de arena.
Un empacador de grava en la terminación, tiene arena de mayor tamaño que la

arena de la formación y ésta es colocada dentro del empacador, sin embargo hay
partículas más finas que se incorporan al flujo de fluidos en el sistema. Dos de los
más comunes empacadores son los que están dentro de la tubería de
revestimiento y los que se encuentran en agujero descubierto o denominados
también como empacadores de grava anillados a la tubería de revestimiento
(figura 9.5), éste último es exitoso ya que retiene la arena y ofrece la mínima
resistencia al flujo dentro del empacador de grava, sin embargo, está limitado a
terminaciones sencillas.
Fig. 9.11 Tipos comunes de empacadores en

terminaciones.
220
9.5 Diseño de Grava y de la Malla
Un elemento crítico en el diseño de los empacadores de grava es el apropiado

tamaño de la grava y la malla. Para mejorar el control de arena y maximizar la
permeabilidad del empacador, la grava debe ser tan pequeña que los finos de la
formación sean retenidos en el empacador y lo suficientemente grandes para que
las arcillas y otras partículas pasen a través de empacador.
El primer paso para determinar el tamaño de la grava del empacador, es hacer la

medición de la distribución del tamaño de las partículas de arena de la formación
que constituye al yacimiento. Una muestra representativa del material de la
formación debe ser obtenida, en orden de preferencia, fragmentos caídos, de
núcleos o de núcleos de pared. Asimismo, las muestras de arena producida en
pozos perforados no deben ser incluidas para este análisis del pozo, sin embargo,
sí a manera de hacer una base de datos de la formación de tal forma que sea
representativa en la distribución areal, siempre y cuando represente a la misma
formación del yacimiento. En el pozo, la arena producida tenderá a tener
partículas muy finas, mientras que la arena tendrá mayor proporción de granos de
partículas más grandes.
El tamaño de la arena de la formación es obtenido con un análisis “Sieve”, usando

una serie del tamaño estándar de Sieve para los Estados Unidos, estandarizan el
tamaño de la malla que está dado en la tabla 9.1, Perry, (1963). Los resultados del
análisis de Sieve son reportados comúnmente en una escala semi-logarítmica de
pesos acumulativos del material de la formación retenidos vs el tamaño del grano.
Típicas distribuciones del tamaño de la arena en California y de la Costa del Golfo
de los Estados Unidos, en arenas no consolidadas se muestran en la figura 9.6
Suman, et al (1983).
221
Tamaño estándar de la malla U.S. TAMAÑO SIEVES

(pg) (mm)
2½ 0.315 8.00
3 0.265 6.73
3½ 0.223 6.68
4 0.187 4.76
5 0.157 4.00
6 0.132 3.36
7 0.111 2.83
8 0.0937 2.38
10 0.0787 2.00
12 0.0661 1.68
14 0.0555 1.41
16 0.0469 1.19
18 0.0394 1.00
20 0.0331 0.840
25 0.0280 0.710
30 0.0232 0.589
35 0.0197 0.500
40 0.0165 0.420
45 0.0138 0.351
50 0.0117 0.297
60 0.0098 0.250
70 0.0083 0.210
80 0.0070 0.177
100 0.0059 0.149
120 0.0049 0.124
140 0.0041 0.104
170 0.0035 0.088
200 0.0029 0.074
230 0.0024 0.062
270 0.0021 0.053
325 0.0017 0.044
Tabla 9. 1 Tamaños estándar Sieve.
222
FIG.9.12 Distribución del tamaño de las arenas para

California y Costa del Golfo, EUA (Suman et al, 1983).
Schwartz (1969) y Saucier (1974), presentaron correlaciones un tanto diferentes

de tamaños de grava óptimos basados en el tamaño de la distribución del grano
de formación. Las correlaciones de Schwartz dependen de la uniformidad de la
formación y de la velocidad de la arena a través de la malla, pero para mayores
condiciones de la arena no uniforme es:
(9. 1)
Donde Dg40 es el tamaño recomendado del tamaño de la grava y Df40 es el

diámetro de la formación para la cual 40 wt % de los granos son de gran tamaño.
Para fijar la distribución del tamaño de la grava, Schwartz recomienda que la
distribución del tamaño de la grava se debería graficar como una línea recta en la
gráfica estándar semi-logarítmica y un coeficiente de uniformidad, Uc, definido
como:
223
(9. 2)
Debería ser 1.5 o menos. De esos requerimientos encontramos:
(9. 3)
(9. 4)
Donde Dg, min y Dg, máx son los tamaños mínimos y máximo de la grava a ser
usados, respectivamente. Las ecuaciones (9.3 y 9.4) definen el rango del tamaño
de la grava recomendado.
Saucier recomienda que la geometría del tamaño de la grava sea cinco o seis
veces que la media del tamaño de la formación, o:
(9. 5)
Saucier no dio recomendaciones acerca del tamaño de distribución de la grava. Si

aplicamos el criterio de Schwartz, entonces:
(9. 6)
(9. 7)
224
9.6 Productividad de Pozos con Empacadores de Grava
La productividad de un pozo con empacador de grava, está afectada por caída de

presión en el empacador, si la caída de presión en el empacador es significativa al
compararla con la caída de presión de la formación. En un empacador de agujero
descubierto, la caída de presión debería ser muy pequeña comparada con la
formación, a menos que la permeabilidad del empacador haya sido severamente
reducida por las partículas de la formación. Si la productividad está expresada con
base en el radio de la tubería, el empacador de grava debería contribuir a un
efecto de daño negativo, así que funcionará como un pozo de mayor diámetro.
Para terminaciones con empacador dentro de la tubería, la caída e presión a

través de las perforaciones pueden contribuir significativamente a una caída de
presión. También desde la productividad de una arena no consolidada que
requiere empacador de grava, es alta, la turbulencia en las perforaciones pueden
agregar al flujo laminar, caída de presión a través de las perforaciones. Las
contribuciones de caídas de presión a través de las perforaciones, pueden ser
representadas por un factor de daño del empacador, sg y el coeficiente para el
relleno de grava en la perforación, Dg.
Golan y Whitson (1991), presentaron las siguientes ecuaciones para empacadores

dentro de la tubería:
(9. 8)
Y para pozos de gas
(9. 9)
Mientras que para pozos de aceite:
225
(9. 10)
Dónde
kh = es el producto permeabilidad – espesor (md-pie)
lperf = es la longitud de la perforación del empacador (pg)
kg = es la permeabilidad de la grava (md)
Dperf = es el diámetro de perforación (pg)
γg = es la gravedad del gas, μ es la viscosidad (cp)
ρ = es la densidad del fluido (lbm/pie3)
n = es el número de perforaciones y (adimensional)
g = es el factor de turbulencia de la grava (adimensional)
El factor de turbulencia es correlacionado con la permeabilidad de la grava

(Cooke, 1973), como:
(9. 11)
Los valores a y b son constantes del tamaño de la grava común, están en la tabla
9.2 (Golan y Whitson, 1991).
Tamaño de la malla, Diámetro significativo Permeabilidad

E.U.A (pg) (mD)
a b
5
40/60 0.014 1.2x10 1.6 2.12x1012
20/40 0.025 1.2x105 1.54 3.37x1012
10/20 0.056 5x105 1.34 8.4x1011
8/12 0.080 1.7x106 1.24 5.31x1011
Tabla 9.2 Propiedades de la grava.
226
Las ecuaciones (9.8) a la (9.11) están basadas en el uso del radio de la tubería de
revestimiento, no de la malla del radio como rw. El término, sg + Dgq, cuando se
agrega al término [ln(re/rw) + s] se toma en cuenta para el aparente efecto de
daños del empacador de grava. Otros efectos de daño están presentes, en
particular, el efecto de daño de las perforaciones, ocurriendo también en el
empacador que está dentro de la tubería y se agregará al factor de daño del
empacador.
Finalmente, usando las ecuaciones (9.8) a la (9.11), las longitudes de las

perforaciones únicamente deberían ser las longitudes a través de la tubería y del
cemento, que es la diferencia entre el radio de la perforación y el radio interno de
la tubería de revestimiento. Esto es porque parece poco probable que un túnel
perforado sea mantenido en un medio no consolidado.
Ejemplo 9.1
Selección del óptimo tamaño de la grava y de la malla.

Usando la correlación de Schwartz, determine el tamaño óptimo de la grava y de
la malla para la arena no consolidada, cuya distribución del tamaño está dada en
la figura 9.6.
Solución.
Correlación Schwartz. La distribución del tamaño de la arena no consolidada está

graficada en la figura 9.7. Leyendo de la gráfica para una fracción de peso
acumulativo de 40%, encontramos que D f40 = 0.0135 pg. El 40% del tamaño del
grano de la grava es entonces (6)(0.0135pg) = 0.081 pg. El 90% del tamaño de
grano es Dg40/1.5 = 0.054 pg. La distribución del tamaño del grano recomendado
se muestra gráficamente como la línea achurada en la figura 9.6; las
intersecciones con la línea del peso acumulativo % 0 10, y % = 0, define el tamaño
mínimo y máximo de la grava, respectivamente, calculado con las Ecs. (9.3) y
(9.4) y son 0.05 y 0.11.
227
De la tabla 9.1, el tamaño Sieves más próximo corresponde al máximo y mínimo

del tamaño de la grava son malla7 y malla16; una malla de arena 8/16 podría ser
seleccionada, ya que la malla 7 es raramente usada. El tamaño de la malla
debería ser menor que 0.0469 pg, así que todas gravas de malla 16 serán
retenidas.
Ejemplo 9.2
Selección del óptimo tamaño de la grava y de la malla.

Usando la correlación de Saucier, determine el tamaño óptimo de la grava y de la
malla para la arena no consolidada, cuya distribución del tamaño está dada en la
figura 9.6.
Solución:
Correlación Saucier. De la distribución del tamaño de la arena, la media (Df50) es

encontrada que es 0.0117 pg. El tamaño de grava recomendado es (5 ó 6)(0.0117
pg) = 0.059 ó 0.070 pg, y de las Ecs. (9.6) y (9.7), el tamaño mínimo de grava es
0.039 a 0.047 pg, mientras que el máximo tamaño de la grava es 0.088 a 0.105
pg. Este rango se muestra como la zona sombreada en la figura 9.6. De la tabla
9.1, el tamaño de granos corresponde al tamaño Sieves de malla 8 y 16 ó 18. Una
malla 8/16 podría ser seleccionada, con una malla menor a 0.0469 pg.
FIG.9 13 Predicciones de la distribución del tamaño de la grava

por las correlaciones de Saucier y Schwartz.
228
Una vez que el tamaño de la grava ha sido seleccionado, es importante verificar

que la grava usada esté conforme a este tamaño. La API (1986), recomienda que
un mínimo de 96% del empacador de grava, la arena debiera pasar el curso
designado por Sieve y ser retenida en el tamaño fino de Sieve
Ejemplo 9.3
Determine el gasto de un pozo de aceite con flujo radial con las siguientes
características:
Pws = 6200 lb/pg2
Pwf = 4900
µ= 2 [cp]
ρo= 39.27 [lb/pie3]
Bo= 1.63
ρdisparos= 8 cargas/pie
Dperf = 0.5 [pg]
Lperf= 6 [pg]
h= 23 [pies]
kg= 120,000 [mD]
re= 835 [pies]
rw= 3 ½ [pg]=0.2917 [pies]
tamaño de la malla= 20-40
Solución:
Calculando el daño Sg de la ec. 9.8
229
De la ecu. Para flujo radial de aceite vista en el capítulo 2
Tenemos
Ejemplo 9.4
Con la información del pozo del ejemplo 9.3 determine el coeficiente para relleno
de grava en la perforación.
Solución:
Calculando el factor de turbulencia de la grava con la ec. 9.11 y obteniendo a y b

de la tabla 9.2 se tiene
Sustituyendo en la ecu. 9.10 tenemos
230
Capítulo 10 DISEÑO DE ESTIMULACIONES Y
FRACTURAMIENTO
10. DISEÑO DE ESTIMULACIONES Y FRACTURAMIENTO
El sistema óptimo de producción es aquél que tiene las menores caídas de presión
en sus componentes, uno de los aspectos que se pueden trabajar y optimizar es la
permeabilidad de la formación que desde un principio es apretada y dada por las
condiciones del medio ambiente de baja energía que ha ocasionado que la misma
sea baja. También es necesario remover el daño que se ha producido durante la
explotación del yacimiento con el fin de que el yacimiento provea de los
hidrocarburos al sistema yacimiento–pozo–instalaciones superficiales.
Para realizar la estimulación o el fracturamiento, es necesario que se tenga el
conocimiento de los ambientes de depósito con el fin de conocer las mineralogías
que constituyen las formaciones de los yacimientos, el factor de daño kh y tener
un equipo de trabajo integrado alrededor del yacimiento en las diferentes
especialidades y los análisis de laboratorio.
10.1.1 Diseño de Acidificación en Areniscas
La acidificación que se emplea normalmente para remover el daño en las

inmediaciones del pozo, es la técnica de estimulación matricial. El análisis
contempla estudios de sensibilización de las rocas y su mineralogía con los ácidos
y concentraciones a nivel de laboratorio, para definir las concentraciones y
cantidades, así como volúmenes y ritmos de inyección. Se diseña un tratamiento
para las areniscas que consiste en etapa de prelavado con la inyección de HCL,
con 50 gal/ft, seguido de la inyección de 50 a 200 gal/ft de una mezcla de HF/HCl.
Un post-lavado de diesel, agua salobre, o HCl que desplaza al HF/HCl del tubo del
pozo. Una vez realizado el tratamiento, el gasto de ácido debe ser inmediatamente
recuperado con el fin de evitar un daño por la reacción de los productos.
Selección del ácido

El tipo de ácido y las concentraciones que generalmente se seleccionan, han sido
bajo la experiencia con las formaciones. Se han realizado tratamientos
estandarizados que consisten en 15 % HCl para formaciones carbonatadas y una
mezcla de 3% HF, 12% HCl. Para las areniscas se utiliza una etapa prelavado con
15% HCl como etapa. De hecho los 3/12 de la mezcla HF/HCl ha sido tan común
que se ha referido genéricamente como ácido de lodo. En años recientes, se ha
dado la tendencia a reducir la concentración de HF (Brannon et al, 1987) cuyo
beneficio de reducción de concentración de HF es que disminuye el daño por los
231
FRACTURAMIENTO
precipitados de la reacción y evitar el riesgo de generar no consolidación de la

formación en inmediaciones del pozo.
Mc Leod (1984), presentó algunas guías para la selección del ácido basado en
experiencia intensa de campo. Sus recomendaciones para yacimientos de
areniscas se muestran en la Tabla 10.1. Generalmente, las pruebas se realizan en
núcleos de 1 pg de diámetro y 3 pg de longitud.
En la Figura 10.1 (Smith y Hendrickson, 1965), muestran la permeabilidad del
núcleo como una función del ácido a través del volumen de poros y denominada
como “curva de respuesta de cido”. Las curvas muestran que las bajas
concentraciones de HF disminuyen el daño en las etapas iniciales de inyección.
Debe tomarse en cuenta que una prueba de laboratorio en un núcleo pequeño no
es representativo a la escala de pozo, pero es una guía para el tratamiento de
campo y tener el mejor modelo de acidificación. Experimentos de núcleos más
grandes, como los reportados por Cheung y Van Arsdale (1992), reflejan mejor las
condiciones de escala del campo.
Tabla 10.1 Acidificación de Arenisca

HCI solubilidad > 20% Uso solo de HCI
Alta permeabilidad (+100 md)
Alto cuarzo (80%), baja arcilla (<5%) 10% HCI-3% HFa
Alto feldespato (>20%) 13.5% HCI-1.5% HFa
Alta arcilla (>10%) 6.5 HCI-1% HFb
Alta arcilla de cloruro de hierro 3% HCI-0.5% HFb
Baja permeabilidad (10 md o menos)

Baja arcilla (<5%) 6% HCI-1.5% HFc
Alta clorita 3% de HCI-0.5% HFd
a
Prelavado con 15% HCI
bPrelavado con 5% HCI secuestrado
CPrelavado con 7.5% de HCI o 10% de ácido acético
dPrelavado con 5% de ácido acético
232
FRACTURAMIENTO
Figura 10.1. Curvas de respuesta de ácido (Smith y Hendrickson, 1965).
10.1.2 Volumen de Ácido y Gasto de Inyección

Las dos metas fundamentales de la acidificación son: remover el daño y
disminuirlo por el proceso mismo de acidificación. Por ejemplo, si la zona dañada
existe cerca de las perforaciones (i.e. 2 pg o menos de la región dañada), el daño
puede ser removido con un mínimo gasto de ácido con una inyección a bajo gasto,
así que el mayor volumen del ácido es reaccionado dentro de las 2 pg de la región
dañada. El gasto bajo puede permitir precipitados de la reacción, reduciendo la
efectividad del ácido. Así que el gasto óptimo está relacionado a la disolución de
los minerales en la zona dañada.
La selección del volumen de ácido necesario está en función de: la profundidad de
la zona dañada y de la curva de o respuesta del ácido en pruebas de laboratorio o
análisis de sensibilidad de la roca al modelo de acidificación.
Modelos de Acidificación de Areniscas
Modelo de dos minerales.

El modelo común empleado es el de ¨dos minerales¨ (Hill et al, 1981; Heckim et al,
1982; Taha et al, 1989) que ubican a los minerales en una a dos categorías:
especies de reacción rápida y de reacción lenta. Schechter (1982) categoriza los
feldespatos, arcillas autigénicas y sílice amorfo como reacción rápida, mientras
que los detritos de partículas minerales y granos de cuarzo son los minerales
233
FRACTURAMIENTO
primarios de reacción lenta. El modelo consiste de balance de materia aplicado al

ácido HF y el reactivo de minerales, los que para flujo lineal, como en un núcleo
impregnado por ácido, puede ser escrito como:
(10. 1)
(10. 2)
(10. 3)
Cuando son dimensionales, asumiendo que la porosidad permanece constante,

esas ecuaciones son:
(10. 4)
(10. 5)
(10. 6)
Donde las variables dimensionales están definidas como:
(10. 7)
234
FRACTURAMIENTO
(10. 8)
(10. 9)
(10. 10)
(10. 11)
Para el núcleo impregnado: L es la longitud del núcleo.

De las Ecs. (10.4) a la (10.6), dos grupos adimensionales aparecen para cada
mineral, NDa el número Damkohler y NAc el número de capacidad del ácido. Esos
dos grupos describen la cinética y la Stokiometría de la reacción del HF con los
minerales.
El número Damkohler es la relación entre el gasto de ácido consumido y el gasto
de ácido convectivo, para la que la reacción rápida de mineral es:
(10. 12)
El número de capacidad del ácido es la relación de la cantidad de mineral disuelto

por el ácido por cada unidad de espacio poroso de la roca a la cantidad de mineral
presente en la unidad de volumen de roca, para la que la reacción rápida de
mineral es:
(10. 13)
Donde la concentración de ácido está en fracción de peso.
Como el ácido es inyectado en la arenisca, se establece un frente de reacción

entre el HF y los minerales de reacción rápida. La forma de este frente depende
de . Para los números Damkohler, el gasto de convección es alto con
235
FRACTURAMIENTO
respecto al gasto de reacción y el frente será difuso. Con un elevado número

Damkohler, el frente de la reacción será relativamente de cuña porque el gasto de
reacción es alto comparado con el gasto de concentración. La figura 10.2 (da
Motta et al, 1992ª) muestra perfiles de concentración para valores altos y bajos de
.
Figura 10.2. Perfiles de concentración de ácido y mineral de reacción rápida.
Schechter (1992) presento una solución aproximada que es válida para números
altos de Damkohler y es útil para propósitos de diseño. Esta solución
aproxima el frente de reacción rápida de minerales como un frente de cuña, de tal
forma que detrás del frente todos los reactivos han sido removidos.
Contrariamente, en el frente no ha ocurrido la disolución. La reacción entre los
minerales de lenta reacción y el HF detrás del frente, sirve para eliminar la
concentración de HF alcanzada en el frente. La localización del frente está dada
por:
(10. 14)
236
FRACTURAMIENTO
La que relaciona el tiempo adimensional a la posición adimensional del frente

dividido por la longitud del núcleo para flujo linear. La concentración adimensional
detrás del frente es:
(10. 15)
Una característica particularmente conveniente de esta aproximación es que es

aplicable a flujo laminar, radial y elipsoidal. El flujo radial representa el flujo de
ácido de un agujero abierto de la terminación y puede también ser una
aproximación razonable a la del flujo de un pozo perforado con suficiente densidad
de perforaciones. La geometría del flujo elipsoidal aproxima al flujo alrededor de
una perforación (fig. 10.3). Las variables apropiadas adimensionales y grupos de
esos tres campos de flujo (tabla 10.2.). Para la geometría de perforación, la
posición del frente, Єf, depende sobre la posición a lo largo de la perforación.
Figura 10.3. Flujo elipsoidal alrededor de una perforación.
237
FRACTURAMIENTO
Tabla 10.2 Grupos adimensionales en modelo de acidificación en areniscas.

Geometría de
Flujo
Linear
Radial
Penetración desde el punto de la

perforación
Penetración adyacente al pozo

Elipsoidal
El número Damkohler para la reacción lenta de minerales y el número de

capacidad del ácido para la reacción rápida de minerales, son únicamente los
grupos adimensionales que aparecen en esta solución. regula que tanto el
HF alcanza el frente; si la reacción lenta de minerales reacciona relativamente
rápido al gasto de concentración, poco ácido será disponible para propagar el
frente de mineral rápido.
El número de capacidad del ácido para la reacción lenta de minerales no es
importante, debido a que el suministro de fluido de reacción lenta de minerales es
en la mayoría constante detrás del frente. afecta directamente al gasto del
frente de propagación: a mayor fluido de reacción rápida de mineral está presente,
será lento el frente. no aparece porque se ha asumido que el frente tiene
forma de cuña, implicando que es infinito. Esta solución puede ser usada
para estimar el volumen de ácido necesario para remover el fluido de reacción
rápida de minerales de una región dada alrededor del pozo o perforación.
y pueden ser calculados con la E.10.3 y la Tabla 10.2 basado en la
mineralogía de la roca o puede ser obtenida de experimentos.
238
FRACTURAMIENTO
Modelo de tres minerales, dos ácidos

Bryant (1991) y da Motta et al (1992b), presentaron evidencias que el proceso de
acidificación de las areniscas no está descrito adecuadamente por el modelo de
dos minerales, particularmente a elevadas temperaturas. Esos estudios sugieren
que la reacción del ácido fluorsilíceo ( ) con los minerales aluminosilicatos
(reacción rápida) puede ser significativo. Esto agrega la siguiente reacción al
modelo de dos minerales:
(10. 16)
Las implicaciones prácticas de esta reacción son significativas, porque se requiere

menos cantidad de HF para consumir el fluido de reacción rápida de minerales con
un volumen dado de ácido porque el ácido fluorsilíceo, también reaccionará con
esos minerales y el producto de la reacción, (gel de sílice) será
precipitado. Esta reacción permite al HF penetrar más dentro de la formación
.
El modelo de tres minerales – dos ácidos, no ha sido probado extensamente, se
sugiere que usando este modelo para predecir el volumen de ácido requerido será
una investigación conservativa, principalmente para aplicaciones de alta
temperatura.
Modelos de precipitación
Walsh et al (1982) propuso un modelo que asume que todas las reacciones están
en equilibrio; todos los gastos de reacciones son infinitamente rápidos. En la figura
10.3 se muestra el resultado de este modelo, que es un diagrama de tiempo –
distancia para la inyección de 4% de concentración de HF/11, concentración de
HCl dentro de la formación que contiene calcita, kaolinita y cuarzo. La gráfica
muestra regiones donde la sílice amorfa y el fluoruro de aluminio tenderán a
precipitar. Una línea vertical, representa la especie de mineral presente como una
función de la distancia si todas las reacciones están en equilibrio.
239
FRACTURAMIENTO
Figura10.3. Diagrama tiempo–distancia mostrando regiones de posible

precipitación.
Sevougian et al (1992) presentó un modelo geoquímico que incluye la cinética

para las reacciones de disolución y precipitación. Este modelo muestra que el
daño de la precipitación será cada vez menos, si tanto la disolución como la
precipitación no son instantáneas. La figura 10.4 muestra predicciones de regiones
y concentraciones de precipitados para cuatro diferentes gastos de reacciones de
precipitados (minerales AC y DB son precipitados). En la medida que el gasto
disminuye la cantidad de daño del precipitado formado disminuye.
240
FRACTURAMIENTO
Figura 10.4. Efecto de precipitación cinética sobre regiones de precipitados de

formación.
Modelos de permeabilidad
Para predecir la respuesta de acidificación de la formación, es necesario predecir
los cambios en la permeabilidad, en la medida que el ácido disuelve algunos
minerales de la formación y otros minerales precipitados. La permeabilidad se
incrementa a medida que el poro y garganta de poro son agrandados por la
disolución del mineral. Al mismo tiempo, pequeñas partículas son removidas del
material cementante y algunas veces esas partículas se alojarán (quizás
temporalmente) en las gargantas de poro, reduciendo la permeabilidad.
Cualquier precipitado formado, tenderá a disminuir la permeabilidad. La formación
de como minerales de carbonato son disueltos y pueden ser causa temporal
de reducción relativa de la permeabilidad a los líquidos.
La compleja naturaleza de la permeabilidad ha hecho impráctica su predicción
como resultado se han generado correlaciones empíricas relacionando el
incremento de la porosidad durante la acidificación. Las correlaciones más
comunes son las de Labrid (1975), Lund y Fogler (1976) y Lambert (1981):
241
FRACTURAMIENTO
(10. 17)
Donde :
son la permeabilidad y porosidad inicial y son la permeabilidad y

porosidad después de la acidificación. M y n son constantes empíricas, reportadas
a ser 1 y 3 para la arenisca Fontainbleu.
(10. 18)
Donde mejor ajustada para la arenisca Phacoides
(10. 19)
La expresión Lambert es idéntica a Lund y Fogler cuando:

Usando los valores de las constantes sugeridas, la correlación de Labrid predice el
más pequeño incremento de permeabilidad, seguido de Lambert y luego Lund y
Fogler. El mejor acercamiento en usar esas correlaciones es seleccionar
constantes empíricas basadas en respuestas de núcleos del yacimiento, si están
disponibles. Careciendo de datos de la formación, la ecuación de Labrid permitirá
el diseño más conservativo.
10.2 Fracturamiento Hidráulico

El objetivo del fracturamiento hidráulico es incrementar la productividad del pozo
creando una ruta o camino altamente conductivo (en comparación con la
permeabilidad del yacimiento) a partir de cierta distancia del pozo y hacia el
interior de la formación. Generalmente la conductividad es mantenida por un
sustentante o apuntalante de arena para mantener la fractura abierta logrando así
que los hidrocarburos tengan una ruta o canal por donde fluir hacia el pozo.
El fracturamiento hidráulico consiste en mezclar una serie de componentes

químicos para crear un fluido fracturante, dicho fluido es bombeado hacia la
formación productora a presiones y gastos lo suficientemente altos para crear y
propagar la fractura a través de la formación.
En primer lugar, se bombea un “colchón” o “pad” de fluido sin apuntalante, es

decir, fluido o gel lineal para iniciar y establecer la propagación de la fractura a
242
FRACTURAMIENTO
través de la formación productora. Esto, es seguido por el gel mezclado con un

sustentante o apuntalante. Este gel continúa siendo bombeado hasta extender la
fractura y simultáneamente transportar el sustentante a través de la formación.
Después de que el fluido es bombeado y se han alcanzado los gastos y presiones
deseadas para establecer la geometría de la fractura, el gel químicamente se
rompe, es decir, baja su viscosidad logrando así que este fluya hacia fuera del
pozo, dejando así una fractura altamente conductiva para que el aceite y/o gas
fluyan fácilmente hacia el pozo, dejando así una fractura altamente conductiva
para que el aceite y/o gas fluyan fácilmente hacia el pozo.
10.2.1 Diseño de un Fracturamiento Hidráulico

Muchos factores influyen en la efectividad y costo de un tratamiento de
fracturamiento hidráulico. En esencia, se tiene muy poco control en cómo y dónde
se propagaran las fracturas. Nuestros esfuerzos se han concentrado y están
limitados en la selección del tipo de apuntalante, los volúmenes apropiados de
materiales, los gastos de inyección y los volúmenes para ser bombeados, así
como un programa de bombeo. Se han tenido algunos desarrollos exitosos en el
crecimiento vertical de la fractura controlando la viscosidad y/o usando diferentes
tipos de aditivos durante el fracturamiento hidráulico.
Con la tecnología actual, el proceso completo de un diseño puede usar un
conjunto de datos y/o características como las mencionadas a continuación, para
evaluar el potencial de producción del yacimiento y para especificar la información
apropiada para el diseño del fracturamiento hidráulico.
1. Radio de drene y configuración
2. Distribución de la zona productora
3. Permeabilidad de la formación, porosidad, saturación de hidrocarburos y perfil
de distribuciones de estos parámetros
4. Propiedades de los fluidos de formación
5. Presión estática del yacimiento
6. Temperatura de la formación
7. Altura de la fractura y crecimiento de la extensión de la misma que ocurrirá
durante el tratamiento de fracturamiento
8. Extensión de la fractura y perfiles de esfuerzos de cierre
9. Presión neta critica de fractura
10. Relación de Poisson y perfiles de densidad.
243
FRACTURAMIENTO
11. Características reológicas del fluido de fracturación, puede ser posible que sea
necesario especificar estos valores como función de la velocidad de corte, tiempo
y temperatura.
12. Características de flujo y pérdidas de presión por los disparos
13. Perdida del fluido de fracturamiento y si es necesario, la dependencia de este
con la temperatura.
14. Coeficiente de filtrado del fluido de fracturamiento, y si es necesario, el
comportamiento como una función de la presión diferencial y la temperatura.
15. Extensión vertical y altura neta del filtrado.
16. Propiedades térmicas del fluido.
17. Tamaño del apuntalante.
18. Densidad del apuntalante.
19. Conductividad del apuntalante como función de los esfuerzos de cierre, tipo de
apuntalante, tamaño de la distribución del apuntalante, concentración del
apuntalante en la fractura, y empacamiento dentro de la formación.
20. Presión de empacamiento y/o integración de la formación
21. Configuración de los disparos (intervalos, disparos por pie, y tamaño de los
agujeros)
22. Configuración y características del árbol de válvulas, tamaños y presiones que
soportan.
Los puntos 1 al 4 pertenecen principalmente a las características de
comportamiento del yacimiento. Los puntos 5 y 6 pertenecen tanto a aspectos del
yacimiento como del fracturamiento. Los puntos 7 a 22 se refieren principalmente
al diseño del fracturamiento hidráulico.
Aunque la lista parece ser un poco sencilla, no refleja fielmente los muchos
factores que intervienen en el diseño de fracturamiento.
La sensibilidad y veracidad de la predicción de resultados depende de la calidad
de la información y el costo del fracturamiento; de la naturaleza y experiencia de
los equipos de trabajo. En algunas áreas, puede ser típico y normal para algunos
operadores progresar a través de un extensivo proceso y probando con diferentes
tipos de fluidos, tamaños del tratamiento y procedimientos que los llevan a
conseguir resultados aceptables. En esas áreas donde los tratamientos de
fracturamiento constituyen una pequeña parte del total de los costos de
perforación y terminación por ejemplo, formaciones de alta permeabilidad donde
las fracturas son muy cortas, estas aproximaciones generalmente son usadas para
244
FRACTURAMIENTO
obtener resultados efectivos. De cualquier manera, en formaciones de baja

permeabilidad donde se requieren que los tratamientos sean muy penetrantes, la
calidad de los parámetros antes mencionados tiene que ser muy buenos. En áreas
donde existen grandes tratamientos y las cantidades son cercanas a la mitad del
total de los costos del pozo, la importancia del fracturamiento es igual o mayor que
el desarrollo de pozos para incrementar las reservas.
Existen métodos para cuantificar adecuadamente algunos de los parámetros de
fracturamiento- como son longitud de la fractura, anchura, conductividad, altura,
azimut, forma, y/o geometría cercana al pozo- aunque algunos de esos están en
etapa de experimentación. Esto hace extremadamente difícil poder evaluar
adecuadamente como podemos predecir el comportamiento y efectividad para un
conjunto de datos dados en el diseño.
Ejemplo 10.1
Se realizó una prueba de inyección de ácido en el laboratorio, a un núcleo de 0.87
[pg] de diámetro con una longitud de 1.57 [pg]. El flujo de ácido fue de 0.0114
ft/min. Determine el volumen de ácido necesario para remover todos los minerales
que reaccionen a una distancia de 3[pg] y 6[pg] de un pozo con un rw= 0.328
[pies]. Asuma que el gasto de ácido inyectado es de 0.1 bl/min-pie de espesor y la
porosidad es de 0.2 y con
Solución:
La posición adimensional del frente para flujo radial es
De la ec. 10.14
De la definición de definición de para flujo radial
245
FRACTURAMIENTO
Ejemplo 10.2
Calcule el volumen de acido necesario en gal/pie para remover todos los

minerales que reaccionan a una distancia de 3 [pg]. En un pozo que fue disparado
con 4 cargas/pie, basado en las pruebas de laboratorio del ejemplo 10.1. la
longitud de las perforaciones es de 6 [pg], la porosidad es de 0.2 y el gasto de
inyección es de 0.1 bl/min-pie.
Solución:
De la tabla para flujo elipsoidal tenemos.
Utilizando la expresión de la tabala relacionamos la posición adimensional del

frente para la penetración del acido
Para el volumen inyectado por perforación
Para el gasto de inyección por pie
246
FRACTURAMIENTO
Ejemplo 10.3
Una arenisca con una porosidad inicial de 0.2 y una permeabilidad inicial de 20
[mD] contiene 10 vol% de carbonato y minerales de reacción rápida en la región
dañada. Calcular la permeabilidad después de remover todos estos minerales
utilizando la correlación de Labrid.
Solución:
El cambio de la porosidad debido a la acidificación es la fracción inciial de volumen
total(1-Φ) multiplicado por la fracción de solido disuelto esto es:
Por lo que la porosidad después de la la acidificación es 0.28. la permeabilidad

para después de la acidificación por la correlación de Labrid es:
Ejemplo 10.4
utilizando la correlación de Lun and Floger
Solución:
total (1-Φ) multiplicado por la fracción de solido disuelto esto es:

para después de la acidificación por la correlación de Lun and Floger es:
Ejemplo 10.5
247
FRACTURAMIENTO

utilizando la correlación de Lun and Floger
Solución:
total (1-Φ) multiplicado por la fracción de solido disuelto esto es:

para después de la acidificación por la correlación de Lambert es:
1704.21 mD
248
Conclusiones.
Las propiedades de los fluidos petroleros se determinan más eficientemente con

una caracterización directa del fluido, ya que a pesar de que las correlaciones
numéricas generan resultados buenos, estos no representan el valor de una
prueba de laboratorio, esto debido a la diferencia que existe en los fluidos a
caracterizar y los fluidos en los cuales se basaron las correlaciones numéricas
propuestas por los autores mencionados en esta tesis. Las correlaciones
numéricas sirven como instrumento alternativo para la caracterización de los
fluidos, siempre y cuando el fluido a ser caracterizado tenga las mismas
condiciones, que el fluido para el cual fue desarrollado la correlación a utilizar.
No se puede tener un pozo que no tenga algún tipo de daño, sin embargo el daño
puede ser despreciado, si por alguna razón este no representa una caída de
presión significativa.
La medición del potencial de un pozo es indispensable, para obtener un estimado

de lo que es capaz de producir el pozo, lo cual servirá como referencia para poder
establecer el método de explotación más apropiado, tanto técnica como
económicamente.
El análisis del sistema integral de producción genera una discretización, de los

principales puntos en los que se tiene las caídas de presión más significativas en
el transporte de hidrocarburos, lo cual ayuda para la implementación de líneas de
acción para la disminución de estas caídas.
249
Bibliografía
1. Michael J. Economides; A. Daniel Hill; Christine Ehlig-Economides.

Petroleum Production System. New Jersey: Precentice Hall, 1944
2. Francisco Garaicochea Petrirena; Cesar Bernal Huicochea: Oscar López

Ortiz. Transporte de Hidrocarburos por Ductos, México: Colegio de
Ingenieros Petroleros de México A.C.
3. Michael Golan; Curtis H. Whitson. Well Performance. 2da Edicion. Noruega:

Tapir, 1996
4. Kermit E. Brown; H. Dale Beggs.The Technology Artificial Methods.

Volumen 1. Tulsa, Oklahoma: Penn Well Books, 1984
5. Kermit E. Brown; H. Dale Beggs.The Technology Artificial Methods.

Volumen 4. Tulsa, Oklahoma: Penn Well Books, 1984
6. H. Dale Beggs. Production Optimizaction, Using Nodal Analysis. OGCI

Publication. Tulsa, 1991.
7. William D. McCain. The Properties of Petroleum Flluids. 2da Edición. Tulsa,

Oklahoma: Penn Well Books, 1990.
8. Cuatli Hernández M.E (2005).Propiedades de los Fluidos Petroleros y

Aplicaciones. Tesis de Licenciatura no publicada, Universidad Nacional
Autónoma de México, Distrito Federal, México.
9. Segura Flores J.E (2010). Trabajos de Reparación y Mantenimiento de

Pozos para los Activos de la Región Norte. Tesis de licenciatura no
publicada. Universidad Nacional Autónoma de México, Distrito Federal,
México.
250
Contenido
Introducción ............................................................................................................. 1
1. USO DE CORRELACIONES PVT PARA OBTENER LAS PROPIEDADES DE
LOS FLUIDOS A DIFERENTES CONDICIONES DE FLUJO ................................. 3
Introducción ......................................................................................................... 3
1.1 Propiedades de los Fluidos Petroleros ....................................................... 3
1.1.1 Propiedades del Gas Natural ............................................................... 3
1.1.2 Propiedades del Aceite .......................................................................... 12
1.2 Clasificación General de los Yacimientos de Acuerdo a su Diagrama de
Fase 19
Yacimiento de Aceite Negro ........................................................................... 20
Yacimiento de Aceite Volátil ........................................................................... 21
Yacimiento de Gas y Condensado ................................................................. 22
Yacimiento de Gas Húmedo ........................................................................... 23
Yacimiento de Gas Seco ................................................................................ 24
1.3 Análisis de los Fluidos Petroleros ............................................................ 25
1.3.1 Métodos directos.................................................................................... 25
1.3.2 Métodos indirectos ................................................................................. 28
1.4 Ecuaciones de Estado ................................................................................. 29
1.5 Correlaciones Numéricas PVT ..................................................................... 35
Propiedades del aceite saturado .................................................................... 35
Propiedades del aceite bajosaturado ............................................................. 41
Correlaciones para obtener la presión en el punto de burbujeo ..................... 43
1.6 Ejercicio de Aplicación ................................................................................. 44
2. COMPORTAMIENTO DE AFLUENCIA ............................................................. 49
2.1. Flujo del Yacimiento al Pozo ....................................................................... 49
2.2 Ecuación de Afluencia .................................................................................. 51
2.3 Geometría de Flujo ...................................................................................... 52
Flujo Lineal ..................................................................................................... 53
Flujo Radial ..................................................................................................... 55
2.4 Regímenes de Flujo ..................................................................................... 57
2.5 Análisis del Pozo Fluyente ........................................................................... 59
Flujo en el Yacimiento..................................................................................... 59
IP en Yacimientos Bajosaturados ................................................................... 60
IPR (Inflow Performance Relationship) en yacimientos saturados ................. 61
Método de Vogel ............................................................................................. 62
Método de Standing ........................................................................................ 64
Método de Fetkovich....................................................................................... 65
2.6 Curvas IPR Futuras...................................................................................... 66
2.6.1 Método de Fetkovich .............................................................................. 66
2.6.2 Método de Eickemer .............................................................................. 67
2.6.3 Método de Standing ............................................................................... 67
2.7 Ejercicio de Aplicación ................................................................................. 69
3. FACTOR DE DAÑO Y SU RELACIÓN CON COMPORTAMIENTO DE
AFLUENCIA .......................................................................................................... 81
3.1 Factor de Daño ............................................................................................ 81
3.2 Flujo Restringido .......................................................................................... 87
Causas Físicas ............................................................................................... 88
Causas Químicas ........................................................................................... 93
Causas Biológicas .......................................................................................... 94
3.3 Factor de Daño Compuesto ......................................................................... 96
3.4 Diagnostico del Daño ................................................................................... 97
3.4.1 Determinación del daño a la formación.................................................. 97
3.5 Ejercicio de Aplicación ........................................................................... 102
4. CURVAS DE DECLINACIÓN. ......................................................................... 115
4.1 Declinación Transitoria ........................................................................... 115
4.2 Declinación en Estado Pseudoestacionario ........................................... 116
Declinación a Gasto Constante .................................................................... 117
Declinación a Presión Constante .................................................................. 117
Declinación exponencial .................................................................................. 119
Declinación hiperbólica ................................................................................. 122
Declinación Armónica ................................................................................... 124
4.3 Ejercicios de Aplicación ............................................................................. 126
5. REGISTROS DE PRODUCCIÓN .................................................................... 129
5.1 Productividad Anormalmente Baja ............................................................. 130
5.2 Producción Excesiva de Gas o Agua ......................................................... 131
Canalización ................................................................................................. 131
Flujo preferencial de gas o agua a través de capas de alta permeabilidad .. 132
Conificación de gas y agua ........................................................................... 134
5.3 Uso de Registros de Producción para la Evaluación de Trabajos de
Reparación ....................................................................................................... 135
5.4 Diagnostico de Pozos Inyectores ............................................................... 136
6. FUNDAMENTOS DEL FLUJO MULTIFÁSICO, FLUJO MULTIFÁSICO
HORIZONTAL, FLUJO MULTIFÁSICO VERTICAL ............................................ 147
6.1 Flujo Multifásico en Tuberías ..................................................................... 147
6.1.1 Pérdidas de presión por fricción........................................................... 149
6.1.2 Perdidas de presión por elevación ....................................................... 151
6.1.3 Colgamiento (HL) ................................................................................. 152
6.2 Flujo Multifásico Vertical ............................................................................ 155
6.2.1 Correlaciones para flujo multifásico ..................................................... 156
6.2.2 Correlaciones para flujo multifásico vertical ......................................... 159
6.3 Flujo Multifásico Horizontal ........................................................................ 162
6.3.1 Correlaciones para flujo horizontal....................................................... 166
7. FLUJO A TRAVÉS DE ESTRANGULADORES .............................................. 173
7.1 Estranguladores ......................................................................................... 173
Clasificación de los Estranguladores ............................................................ 173
7.2 Flujo en Estranguladores ........................................................................... 174
Gilbert, Ros, Baxendell y Achong. ................................................................ 176
Poettmann y Beck. (P y B) ............................................................................ 178
Ashford ......................................................................................................... 179
Ashford y Pierce ........................................................................................... 179
Omaña .......................................................................................................... 180
7.3 Flujo de Gas a Través de Estranguladores ................................................ 180
8. SISTEMA INTEGRAL DE PRODUCCIÓN, ANÁLISIS NODAL ....................... 191
8.1 Sistema Integral de Producción ................................................................. 191
Yacimiento .................................................................................................... 192
Pozo.............................................................................................................. 192
Tubería de descarga ..................................................................................... 192
Estrangulador ............................................................................................... 192
Separadores ................................................................................................. 192
Tanques de almacenamiento........................................................................ 193
8.2 Análisis Nodal ............................................................................................ 193
Análisis cualitativo del sistema ..................................................................... 194
Flujo en el yacimiento ................................................................................... 196
Comportamiento del flujo por la TP .............................................................. 196
Terminación del flujo natural ......................................................................... 198
Efecto del diámetro de la TP......................................................................... 199
Comportamiento del flujo por la LDD ............................................................ 201
8.3 Optimización de un Sistema de Producción ............................................... 203
9 DISEÑO DEL APAREJO DE PRODUCCIÓN ................................................... 211
9.1 Opciones de Terminación de un Pozo ....................................................... 211
9.2 Clasificación de los Aparejos de Producción ............................................. 213
Bombeo Neumático Continuo (BNC) ............................................................ 213
Bombeo Neumático Intermitente (BNI) ......................................................... 214
Sistema de Bombeo Hidráulico (BH) ............................................................ 215
Sistema de Bombeo Electrocentrífugo (BEC) ............................................... 216
Sistema de Bombeo Mecánico ..................................................................... 216
9.3 Selección del Diámetro del Aparejo de Producción ................................... 217
Elevación ...................................................................................................... 218
Fricción ......................................................................................................... 218
Aceleración ................................................................................................... 219
9.4 Empacadores de Grava ............................................................................. 220
9.5 Diseño de Grava y de la Malla ................................................................... 221
9.6 Productividad de Pozos con Empacadores de Grava ................................ 225
10. DISEÑO DE ESTIMULACIONES Y FRACTURAMIENTO ............................ 231
10.1.1 Diseño de Acidificación en Areniscas ................................................... 231
Selección del ácido ....................................................................................... 231
10.1.2 Volumen de Ácido y Gasto de Inyección ............................................... 233
Modelos de Acidificación de Areniscas ......................................................... 233
10.2 Fracturamiento Hidráulico ........................................................................ 242
10.2.1 Diseño de un Fracturamiento Hidráulico ............................................ 243
10.3 Ejercicios de Aplicación ........................................................................... 245
Conclusiones. ...................................................................................................... 249
Bibliografía .......................................................................................................... 250

Ing. de Procutividad de Pozos

Cargado por

Copyright:

Formatos disponibles

Ing. de Procutividad de Pozos

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ing. de Procutividad de Pozos

Cargado por

Copyright:

Formatos disponibles

UNIV E R SID A D

CONOCE LOS PROGRAMAS ACADÉMICOS OFICIALES DE GRUPO CEDIP

La ingeniería de producción es una de las grandes áreas de la ingeniería

El ingeniero de producción tiene como prioridad administrar la energía del

El transporte de hidrocarburos, el diseño de instalaciones superficiales, el gasto de

1. USO DE CORRELACIONES PVT PARA OBTENER LAS PROPIEDADES DE

1.1 Propiedades de los Fluidos Petroleros

1.1.1 Propiedades del Gas Natural

Estas propiedades se utilizan en cálculos que se basan en reglas de mezclado

Tabla 1.1 Propiedades físicas de los gases a 14.7 lb/pg2abs y 60 °F

Peso Presión Densidad Densidad Viscosi

Hidrógeno H2 2.109 59.36 187.5 4.432 0.005312 0.00871

Helio He 4.003 9.34 32.9 7.802 0.010550 0.01927

Agua H2O 18.015 1,164.85 3,200.14 62.336 - ~ 1.122

Nitrógeno N2 28.013 227.16 493.1 50.479 0.073820 0.01735

oxigeno O2 31.99 278.24 731.4 71.228 0.084320 0.02006

Aire - 28.963 238.36 546.9 54.555 0.076320 0.01790

Metano CH4 16.043 343.00 666.4 18.710 0.042280 0.01078

Etano C2H6 30.070 549.59 706.5 22.214 0.079240 0.00901

Propano C3H8 44.097 665.73 616.0 31.619 - 0.00788

i-Butano C4H10 58.123 734.13 527.9 35.104 - 0.00732

n-Butano C4H10 58.123 765.29 550.6 36.422 - 0.00724

i-Pentano C5H12 72.150 828.77 490.4 38.960 - -

n-Pentano C5H12 72.150 845.47 488.6 39.360 - -

n-Hexano C6H14 86.177 913.27 436.9 41.400 - -

n-Heptano C7H16 100.204 972.37 396.8 42.920 - -

n-Octano C8H18 114.231 1,023.89 360.7 44.090 - -

n-Nonano C9H20 128.256 1,070.35 331.8 45.020 - -

n-Decano C10H22 142.285 1,111.67 305.2 45.790 - -

Densidad del Gas

Combinando esta ecuación y resolviendo la densidad (m/V) :

: Densidad del gas

: Gravedad especifica del gas

: Factor de compresibilidad del gas

: Factor volumétrico del gas

Densidad específica de un gas (densidad relativa)

La densidad relativa de un gas, γg, es la relación de la densidad del gas a la

en donde, g es la densidad de la mezcla de gases en lbm/pies3 y aire es la

A condiciones estándar, ambos el aire y el gas seco se modelan (comportan) en

en donde, γg es la densidad relativa del gas (aire =1.0), M es el peso molecular

Si no se conoce el peso molecular del gas , pero se conoce los pesos

: Peso molecular del componente j en la mezcla

: Número de componentes de la mezcla

Factor de desviación de los gases reales (factor de desviación z)

La ley de los gases ideales establece que la presión, temperatura volumen y la

en donde, p es la presión en lb/pg2abs, V es el volumen en pies3, n es el número

Para gases reales, el factor de desviación z se define como:

en donde, Vreal es el volumen real de gas que corresponde a una presión y

Para gases reales, el factor de desviación z se calcula como,

El factor de compresibilidad z es adimensional. Comúnmente, z es muy próximo a

A partir de la ecuación 1.10 se desarrollan expresiones para definir el factor de

Factor de volumen del gas de formación, Bg

para una p @c.e.=14.65 lb/pg2abs y una T @c.e.=60°F (519.59 °R), se tiene:

en donde Bg presenta unidades de pies3 @c.y./pies3 @c.e. La Fig. 1.1 presenta el

Bg(ft3@ c.y./ft3 @c.s.)

Presión del yacimiento, py, (lb/pg2 abs)

Coeficiente de compresibilidad isotérmico del gas, cg

El coeficiente de compresibilidad isotérmico del gas (compresibilidad del gas