Explotación Del Gas y Optimización de La Producción

Explotacin del Gas y
Optimizacin de la
Produccin
Usando Anlisis Nodal
Balance de Materiales con entrada de agua
5000
4500
4000
3500
P /Z (P S I)
3000
2500
2000
1500
1000
500
0
0 200000 400000 600000 800000 1000000 1200000 1400000
(Psc*Gp*Tr/Tsc)-Pf/Zf(Wi-WP*Bw+Giny)
Jos Luis Rivero S.

Explotacin del Gas y
Optimizacin de la Produccin
Software y Mtodos Analticos
Por: Jos Luis Rivero Sandoval

Presentacin
Debido al crecimiento de la demanda de la energa, en los distintos campo de la

industria. El gas se convirti en una fuente de energa de rpido crecimiento. Y por lo
tanto su explotacin y su optimizacin de la produccin deben ser analizadas mas
profundamente para obtener un buen desarrollo. El presente libro sobre la explotacin
del gas y optimizacin de la produccin nos da una herramienta para el desarrollo de un
campo gasifero, que depende mucho del la demanda y su mercado.
La publicacin de este libro es muy importante en el mercado del gas para

personas de habla hispana, debido a que la mayora de los libros y folletos estn
publicados en ingles. En el presente libro estamos presentando tres softwares para
facilitar el desarrollo de un campo, primeramente se estudia el reservorio, luego
optimizamos los pozos para tener un entrega de potencial que este de acuerdo a los
requerimientos.
Los tres programas presentados WELLGAS, BALANGAS, PRONOSGAS ,

fueron programados en visual basic, siguiendo paso a paso los clculos realizados,
donde la aplicacin de las formulas se muestran en los mismos programas, los cuales
tratan de cubrir algunos concepto y mtodos de clculos que otros programa no lo
tienen.
Ing. Jos Luis Rivero

[email protected]
Santa Cruz de la Sierra Bolivia

Septiembre 2004
Prohibida su reproduccin
total o parcial. Resolucin
Ministerial No 1-416/2004
Agradecimientos
Primeramente quiero agradecer a la Universidad Gabriel Rene Moreno y YPFB

por el apoyo brindado para la cristalizacin de este libro, principalmente al Ing.Walter
Yabeta Decano de la Facultad de Tecnologa por el apoyo en la publicacin de este
libro.
Al Dr. Hugo Araujo por la revisin criteriosa del libro, comentarios y

sugerencias, sobre el contenido y forma de presentacin de los diferentes captulos.
Al Ing. Luis Kin por las contribuciones recibidas durante el proceso de

confeccin del libro y la revisin de algunos captulos y sugerencias.
Quiero tambin agradecer a los Ings. Gerardo Corcos, Javier Velarde, Joaqun
Texeira, Ming. Jos Escobar por la colaboracin y sugerencias recibidas para la
confeccin de este libro.
A mis colegas y amigos Ings. Esteban Cabrera, Oscar Jalil por la transcripcin
y ajustes de los programas realizados.
Finalmente quiero agradecer a mi esposa Amaly por la comprensin, apoyo e

incentivo para llegar a terminar este proyecto. Tambin quiero agradecer a mis hijos.
Explotacin del Gas y Optimizacin de la Produccin ------------------------------------ 1
SUMARIO
Pag.
Introduccin.7
Captulo 1: Anlisis Nodal Pozos de Gas
1.1.- Introduccin al Anlisis Nodal10

1.2 Anlisis del Sistema de Produccin ..12
1.3 Nodo ...14
1.3.1 Nodo Fijo 14
1.3.2 Nodo Comn ..14
1.4 Elementos usados en el Sistema del Anlisis Nodal....15
1.4.1 Ubicacin de los Nodos Componentes15
1.4.2 Componentes que Intervienen en el Anlisis Nodal.16
1.4.2.1 Separador..16
1.4.2.2 Lnea de Flujo Horizontal.16
1.4.2.3 Lnea de Flujo Vertical.16
1.4.2.4 Choque Superficial16
1.4.2.5 Cabeza de Pozo.16
1.4.2.6 Vlvula de Seguridad16
1.4.2.7 Choque de Fondo..17
1.4.2.8 Presin Fluyente. 17
1.4.2.9 Completacin o perforacin en el fondo...17
1.4.3 Presin Constante 17
1.4.4 Anlisis del Sistema en Fondo 17
1.4.5 Optimizacin de la Tubera de Produccin ..19
1.4.6 Efecto de Agotamiento del Reservorio ....20
1.5 Anlisis del Sistema Nodo en Cabeza de Pozo21
1.6 Anlisis del Sistema en el Separador22
1.7 Seleccin del Compresor.24
1.8 Anlisis del Sistema para pozos con restricciones Superficiales.25
Captulo 2: Clasificacin de los Fluidos en el Reservorio
2.1 Introduccin.29
2.2 Diagrama de Fases (Presin Temperatura)30
2.2.1. Propiedades Intensivas..31
2.2.2 Punto Crtico.31
2.2.3 Curva de Burbujeo .31
2.2.4 Curva de Roco ..31
2.2.5 Regin de dos fases.32
2.2.6 Cricondenbar32
2.2.7 Cricondenterma32
2.2.8 Zona de Condensacin Retrograda...32
2.2.9 Petrleo Saturado..32
2.2.10 Petrleo Bajo Saturado.32
2.2.11 Petrleo Supersaturado.32
2.2.12 Saturacin Crtica del Fluido32

2.3 Clasificacin de los Reservorios ...33
2.3.1 Reservorio de Petrleo..33
2.3.1.1 Reservorio de Petrleo Subsaturado.33
2.3.1.2 Reservorio de Petrleo Saturado...34
2.3.1.3 Reservorio con Capa de Gas 34
2.3.2 Petrleo Negro..35
2.3.3 Petrleo Negro de bajo Rendimiento 36
2.3.4 Petrleo Voltil ..37
2.3.5 Petrleo Cerca al Punto Critico38
2.4 Reservorio de Gas.39
2.4.1 Reservorio de Condensacin Retrograda.39
2.4.2 Reservorio de Gas Condensado cerca al punto Crtico40
2.4.3 Reservorio de Gas Hmedo..41
2.4.4 Reservorio de Gas Seco41
2.5 Correlaciones para determinar el Punto de Roco42
2.5.1 Determinacin del Punto de Roco con la Composicin..42
2.5.2 Determinacin del Punto de Roco con datos de Produccin...45
2.6 Pruebas PVT.46
2.6.1 Tipos de Pruebas PVT ...47
2.6.1.1 Proceso a Composicin Constante...47
2.6.1.2 Proceso a Volumen Constante..48
2.6.1.3 Proceso de Liberacin Diferencial .48
Captulo 3 : Propiedades del Gas Natural
3.1 Introduccin a las Propiedades del Gas Natural.50

3.2 Gas Natural.50
3.2.1.- Composicin del Gas Natural 51
3.2.2.- Comportamiento de los Gases Ideales.51
3.2.3.- Ley de los Gases Ideales..51
3.2.3.1 Ley de Boyle52
3.2.3.2 Ley de Charles..52
3.2.3.3 Ley de Charles y Boyle53
3.2.3.4 Ley de Avogadro...53
3.2.4 Ecuacin para los Gases Ideales...53
3.2.5 Mezclas de Gases Ideales.54
3.2.5.1 Ley de Dalton55
3.2.5.2 Ley de Amagar..56
3.2.5.3 Fraccin Volumtrica57
3.2.5.4 Peso Molecular Aparente..58
3.2.5.5 Fraccin Molar..58
3.2.6 Densidad del Gas..59
3.3.1 Gases Reales.63
3.3.2 Mtodo de obtencin del factor de compresibilidad Z 63
3.3.3 Correlacin de Standing y Kats...64
3.3.4 Correlacin de Brill & Beggs...69
3.3.5 Correlacin de Drank, Purvis y Robinson69
3.3.6 Correlacin de Hall Yarborough71

3.4.1 Factor Volumtrico del Gas..73
3.5.1 Compresibilidad Isotrmica del Gas Natural74
3.5.1.1 Compresibilidad para un Gas Ideal...74
3.5.1.2 Compresibilidad para un Gas Real ..75
3.6.1 Viscosidad del Gas Natural..76
3.6.1.1 Determinacin de la viscosidad: Mtodo Carr, Kobayashi..... 76
3.6.1.2 Determinacin de la viscosidad: Mtodo Lee, Gonzlez y Eakin77
3.7.1 Factor de Compresibilidad para un Sistema Bifsico Gas Retrogrado.83
3.7.1.1 Coordenadas seudo crticas del C+7..84
3.7.1.2 Coordenadas seudo crticas de la Mezcla.84
Captulo 4 : Anlisis de Reservorio
4.1 Introduccin87
4.2 Ley de Darcy..88
4.2.1 Flujo Lineal.89
4.2.2 Flujo Radial.92
4.2.3 Flujo de Gas.93
4.3 Flujo en Estado de Equilibrio..95
4.3.1 Condiciones Semiestable.95
4.3.2 Ecuacin para Flujo Radial en Funcin al Pseudo Potencial .98
4.3.3 Capacidad de Entrega de un Pozo de Gas .99
4.3.4 Flujo Transiente en pozo de Gas100
4.4 Tipo de Prueba 104
4.4.1 Prueba de Flujo tras Flujo 104
4.4.2 Prueba Isocronal..105
4.4.3 Prueba Isocronal Modificada ...105
4.4.4 Prueba de Produccin..106
4.5 Mtodo de Interpretacin de Prueba 106
4.5.1 Mtodo Simplificado...107
4.5.2 Mtodo de Blount Jones and Glaze 109
4.5.3 Mtodo de Brar y Aziz112
4.5.4 Mtodo Inercial LIT 115
4.6 Produccin Pozos Horizontales .120
4.6.1 Impacto del Dao en el Comportamiento de un pozo Horizontal122
4.6.2 Efectos de produccin de Agua y Permeabilidades relativas..123
4.6.3 Relacin del ndice de Productividad para un pozo Horizontal 124
4.6.4 Efecto del Dao de un pozo Horizontal..125
4.7 Factores que afectan la curva del comportamiento del pozo con el tiempo...127
4.7.1 Coeficiente C y la exponente n 127
4.7.2 Permeabilidad del Gas127
4.7.3 Espesor de la Formacin128
4.7.4 Viscosidad del Gas y el Factor de Compresibilidad..128
4.7.5 Radio de Pozo y Drene..128
4.7.6 Factor de Dao...129
4.8 Cada de la Presin a travs de las Perforaciones .129
4.9 Perforaciones de Pozo y Efecto de Dao..........................................................132

4.9.1 Clculo del efecto de Dao en la Perforacin............................134
4.9.2 Dao cerca al pozo y las Perforaciones......................................136
Captulo 5: Perdida de Presin en Lneas de Produccin
5.1 Introduccin..139
5.2 Ecuacin de Energa..140
5.3 Gradiente de Presin Total ...144
5.4 Nmero de Reynolds.145
5.5 Rugosidad Relativa146
5.6 Determinacin del Factor de Friccin149
5.6.1 Flujo Laminar de Fase Simple149
5.6.2 Flujo Turbulento de Fase Simple150
5.6.2.1 Tuberas Lisas 150
5.6.2.2 Tuberas Rugosas150
5.7 Flujo de Fase Simple...151
5.8 Flujo de dos Fases...152
5.8.1 Variable de Flujo de dos Fases ..152
5.8.1.1 Escurrimiento de Lquido HL..152
5.8.1.2 Suspensin de Lquido L153
5.8.1.3 Densidad..154
5.8.1.4 Velocidad156
5.8.1.5 Viscosidad ..156
5.8.1.6 Tensin Superficial 156
5.8.2 Modificacin de la ecuacin de la Gradiente de Presin para flujo
Bifsico...156
5.8.3 Modelo Simplificado para Predecir la Velocidad Mnima para Remover
el Lquido en el Fondo157
5.8.4 Componente de Elevacin..159
5.8.5 Componente de Friccin.....159
5.8.6 Componente de Aceleracin...160
5.9 Modelo de Flujo de dos fases 161
5.9.1 Clculo de la Presin Transversal ..161
5.9.2 Determinacin de la Distribucin de Temperatura 164
5.10 Flujo en Pozos de Gas ...165
5.10.1 Presin de Fondo 165
5.10.1.1 Mtodo de la Presin y Temperatura Media (Esttica) ..165
5.10.1.2 Mtodo de Cullender y Smith (Esttica) 166
5.10.2 Presin Dinmica de Fondo168
5.10.2.1 Mtodo de la Presin y Temperatura Media (Dinmica)169
5.10.2.2 Mtodo de Cullender y Smith (Dinmica)..171
5.10.2.3 Mtodo de Grey (Dinmica)...175
5.11 Flujo de Gas en lneas de Surgencia...179
Captulo 6 : Anlisis de Flujo de Gas a travs de los Choques
6.1 Introduccin183
6.2 Clasificacin de los Choques..184
6.2.1 Choques Superficiales.184
6.2.1.1 Tipo Positivo ..184
6.2.1.2 Tipo Ajustable.....184
6.2.2 Choques de Fondo...184
6.2.2.1 Tipo Fijo.. ...186
6.2.2.2 Tipo Removible .....186
6.3 Factores que Influyen en el Choque. .....187
6.4 Modelos de Flujo...188
6.4.1 Flujo Subcrtico..189
6.4.2 Flujo Crtico189
6.4.2.1 Pruebas de Flujo Crtico..............................................................190
6.5 Flujo de Gas191
6.5.1 Flujo simple Fase191
6.5.2 Flujo dos Fases ...193
Captulo 7: Balance de Materiales Reservorio de Gas
7.1 Introduccin197
7.2 Ecuacin de Balance de Materia.198
7.3 Reservorio Volumtrico de Gas .200
7.3.1 Reservorio Volumtrico Anormalmente Presurizado.201
7.4 Reservorio de Gas con Entrada de Agua205
7.5 Linealizacin de la Ecuacin de Balance de Materiales 206
7.6 Balance de Materiales Propuesto208
7.7 Reservorio de Gas Condensado no Retrogrado..213
7.7.1 Clculo del Volumen de Gas Equivalente al Condensado..213
7.7.2 Clculo del Volumen de vapor de Agua Equivalente al Agua Producida .214
7.8 Reservorio de Gas Condensado Retrogrado ..215
7.9 Factor de Recuperacin .216
Captulo 8 : Intrusin de Agua
8.1 Introduccin218
8.2 Modelo de Van Everdingent & Hurst 218
8.2.1 Acufero Radial ..219
8.2.2 Acufero Lineal ..221
8.3 Efecto de Superposicin 223
8.3.1 Teora del Ajuste de la Entrada de Agua224
8.4 Modelo de Fetkovich..225
8.5 Modelo de Carter Tracy..229
8.6 Modelo de Leung231
8.6.1 Modelo Pseudo Permanente (PSS).231
8.6.2 Modelo Pseudo Permanente Modificado (MPSS).236

8.7 Comparacin de Modelos .245
Captulo 9: Desarrollo de un Campo de Gas
9.1 Introduccin259
9.2 Reservas .260
9.2.1 Comportamiento del Reservorio.260
9.2.2 Desarrollo del Campo 262
9.3 Entrega de Potencial...263
9.4 Espaciamiento de Pozo...266
9.5 Capacidad de los Equipos de Produccin ..268
9.5.1 Capacidad de la Lnea de Flujo ...269
9.5.2 Capacidad de Compresin .270
9.5.3 Capacidad de lnea..270
9.6 Prediccin y Comportamiento del Reservorio ..271
9.7 Optimizacin del Desarrollo del Campo276
Captulo 10: Manual y Aplicacin de los software
10.1 Introduccin279
10.2 Aplicacin del software WELLGAS
10.2.1 Prueba isocronal para pozo de gas seco.........305
10.2.2 Analisis nodal para un pozo de seco con datos de reservorio.................306
10.2.3 Analisis Nodal Pozo de Gas Condensado.......319
10.2.4 Pozo de Gas Condensado Horizontal.....................................333
10.3 Aplicacin del software BALANGAS
10.3.1 Reservorio Gas de Condensado Con entrada de Agua...349
10.3.2 Reservorio de Gas Condensado Con Entrada de Agua e inyeccin de Gas......363
10.4 Aplicacin de software PRONOSGAS
10.4.1 Entrega de potencial Reservorio de Gas para una presin de
planta...................................................................................................................375
10.4.2 Entrega de potencia Reservorio Gas para distintas presiones de
planta..383
10.4.3 Entrega de potencial Reservorio Gas sin produccin Inicial...................392
Introduccin
Sobre la base de la creciente demanda de energa en el mundo, el gas natural se
convirti en una fuente de energa de rpido crecimiento en los ltimos aos, y se
estima una proyeccin exponencial para las prximas dcadas, en cuanto a su
participacin en el consumo total de energa.
El mayor crecimiento registrado fue dado en el campo de la generacin de

electricidad. Las turbinas a gas de ciclo combinado con las plantas de generacin,
constituyen un ejemplo de la ms alta eficiencia econmica disponible y adems, el gas
natural ambientalmente es ms atractivo, debido a que emite menos dixido de sulfuro y
dixido de carbono que el petrleo o el carbn.
En los pases industrializados, dadas estas ventajas econmicas ambientales, la

proyeccin del crecimiento en su consumo es mayor que entre los combustibles
tradicionales.
Para los pases en desarrollo se espera que el gas natural, se utilice no solamente
para la generacin de electricidad, si no tambin para otros usos tales como gas
domiciliario y como combustible en el sector industrial y automotriz.
En base a la experiencia obtenida con el transcursos de los aos en

interpretaciones de pruebas de pozo, tanto para restituciones de cierre, como para
pruebas de flujo, isocronal, isocronal modificada y flujo tras flujo, he observado que
existe un mediano porcentaje de pruebas no representativas, debido a que durante la
prueba no se analiz el caudal mnimo de gas para remover el lquido en el fondo ni se
analiz las condiciones de flujo subcrticos con los cuales esta produciendo el pozo,
dando resultados no representativos creando criterios errneos o equivocados para
desarrollar el reservorio. Por la importancia de este problema hemos planteado el
desarrollo de un software llamado WELLGAS, que analiza la calidad de los datos de
entrada evitando una mala interpretacin de la capacidad real del pozo, y as no sobre
dimensionar los nmeros de pozos a perforar. Si obtenemos datos representativos de
pozo vamos a tener pronsticos ms reales para cualquier desarrollo. En el captulo 10
mostramos el manual de aplicaciones para distintos reservorios de gas.
El clculo de reserva es muy importante ya que nos da una idea general de la

reserva disponible en el yacimiento, conociendo esta reserva podemos saber hasta
cuando un reservorio es econmicamente rentable para explotarlo, de ah la importancia
de tener valores precisos. En los mtodos de balance de materia comnmente usado P/z,
mtodo de la lnea recta, el volumen de gas In-Situ es calculado en base a la linealidad
de los puntos del balance que puede variar de acuerdo a la los datos; algunas veces
estos datos, de los historiales de presin vs. el acumulado de gas, no se encuentran sobre
una lnea recta, debiendo contrapesar los datos, los cuales nos lleva a obtener distintos
resultados del volumen In-Situ de acuerdo al criterio del analista. En el presente libro,
captulo 7 (Balance de Materia para reservorio de gas) proponemos un mtodo basado
en la determinacin del Volumen In-Situ en funcin al tiempo, el cual no debe variar
para las distintas producciones acumuladas, ahorrando la incertidumbre de tomar otros
valores que puedan ser o no los correctos, los que pueden influenciar en nuestro anlisis
econmico.
Otro factor importante que se debe tomar en cuenta en el balance de materiales

es el factor de compresibilidad bifsico para los reservorios de gas condensado, debido a
la condensacin del gas por debajo de la presin de roco, en estos casos si tomamos p/z
monofsico estamos sobreestimando nuestras reservas. Para evitar todos estos
problemas plantados estamos presentando un programa para el balance de Materia
llamado BALANGAS estimando nuestras reservas en forma ms confiable.
En el captulo 9, mostramos el proceso de un desarrollo ptimo de campo de gas

dependiendo de su mercado y su demanda en funcin a su reserva recuperable,
capacidad de entrega de cada pozo, y facilidades superficiales, como ser las presiones
de planta o lnea. Unos de los grandes problemas planteados es la prediccin de los
caudales de entrega de gas, por pozo o por reservorio, en funcin al agotamiento natural
del reservorio. Considerando que se lo debe entregar en un punto de venta con una
presin especifica de lnea o planta y un poder calorfico definido a una presin de
roco. En este libro, tambin estamos presentando un programa para el desarrollo de un
campo llamado PRONOSGAS, que nos da la pauta del desarrollo de un campo
gasifero, los nmeros de pozos que podemos perforar o intervenir de acuerdo a la
reserva recuperable en funcin al agotamiento del reservorio.
En el captulo 1, sealamos los conceptos bsicos del anlisis nodal para la

optimizacin de la produccin, definiendo el sistema total y subdividiendo los distintos
componentes para analizar la cada de presin en cada uno de ellos, determinando los
cuellos de botella del sistema donde existe la mayor perdida de presin.
En el captulo 2-3, clasificamos los fluidos en el reservorio en funcin al

diagrama de fase (Presin- Temperatura) y a su relacin RGC, tambin se determina
presin de Roco en funcin de la composicin del gas y datos de produccin,
definiendo as los sistemas monofsico, bifsicos y multifsicos. En estos captulos
mostramos las propiedades del gas natural y algunas de las correlaciones ms usadas
dentro de la industria, como as tambin, se muestra las correlaciones del factor de
compresibilidad para un sistema bifsico de gas condensado, con las coordenadas seudo
criticas C7+ , que es muy importante en los clculos de balance de materia.
En los captulos 4-5-6, se explica el anlisis de reservorio, los distintas tipos de

pruebas y mtodos de interpretacin, determinando la relacin del ndice de
productividad tanto para pozo vertical como horizontal, definiendo los factores que
afectan al comportamiento del pozo con el tiempo. En estos captulos se analiza la cada
de presin en las lneas de produccin para sistemas monofsico como tambin
multifsicos, teniendo dos correlaciones importantes para la cada de presin en tubera
vertical para un sistema monofsico. y una correlacin para el sistema multifsico. Con
respecto a la lnea horizontal, estamos proponiendo las correlaciones ms utilizadas en
el mercado al igual que para las restricciones o choques.
En el capitulo 7-8, sugerimos un nuevo mtodo de Balance de materiales para

determinar el volumen In-Situ en funcin al tiempo, ya que el volumen Inicial no varia
con el tiempo de agotamiento, adems mostramos y comparamos 5 mtodos de entrada
de agua para poder definir un buen comportamiento.
En general, el libro Explotacin del Gas y Optimizacin de la Produccin nos da

el concepto esencial de Ingeniera de reservorio para la optimizacin y desarrollo de los
reservorios de gas natural; tambin se presenta la optimizacin de la produccin en
funcin al anlisis nodal para evaluar efectivamente un sistema completo de extraccin,
considerando todos los componentes del pozo, comenzando por la presin de reservorio,
incluyendo el flujo a travs del medio poroso, flujo a travs de las perforaciones, flujo a
travs de las tuberas tanto vertical, horizontal, terminando en el separador. Siendo su
objetivo principal el de optimizar la produccin variando los distintos componentes del
sistema para un mejor rendimiento econmico.
El desarrollo de un campo de gas no solo depende del reservorio y la

caracterizacin de los pozos, si no tambin del equipamiento superficial (separadores,
planta) con este libro sealamos la sensibilidad del reservorio con respecto al sistema de
produccin, debido a que si no se tiene un sistema de produccin superficial adecuado
no podemos obtener una mayor recuperacin del volumen In-Situ.
Explotacin del Gas y Optimizacin de la Produccin________________________________________10
1 ANLISIS NODAL
1.1 Introduccin
El anlisis nodal se define como la segmentacin de un sistema de produccin en
puntos o nodos, donde se producen cambios de presin. Los nodos estn definidos por
diferentes ecuaciones o correlaciones.
El anlisis nodal es presentado para evaluar efectivamente un sistema completo de

produccin, considerando todos los componentes del sistema comenzando por la presin de
reservorio Pr y terminando en el separador incluyendo el flujo a travs del medio poroso,
flujo a travs de las perforaciones de terminacin, flujo a travs de la tubera de produccin
con posibles restricciones de fondo, flujo de la lnea horizontal pasando a travs del
estrangulador en superficie hacia el separador.
El objetivo principal del anlisis nodal, es el de diagnosticar el comportamiento de

un pozo, optimizando la produccin, variando los distintos componentes manejables del
sistema para un mejor rendimiento econmico.
Para que ocurra el flujo de fluidos en un sistema de produccin, es necesario que la

energa de los fluidos en el reservorio sea capaz de superar las prdidas de carga en los
diversos componentes del sistema. Los fluidos tienen que ir desde el reservorio hacia los
separadores en superficie, pasando por las tuberas de produccin, equipos superficiales en
cabeza de pozo y las lneas de surgencia. La Figura 1.1 muestra un sistema de produccin
simple, con tres fases:
1. Flujo a travs del medio poroso.

2. Flujo a travs de la tubera vertical o direccional.
3. Flujo a travs de tubera horizontal
La Figura 1.1 muestra todos los componentes del sistema en los cuales ocurren las
prdidas de presin, que va desde el reservorio hacia el separador.
p1 = pr pwfs = Prdidas de presin en medios porosos.

p2 = pwfs pwf = Prdidas de presin a travs de la completacin.
p3 = pUR pDR = Prdidas de presin a travs de las restricciones.
p4 = pUSV pDSV = Prdidas de presin a travs de la vlvula de
seguridad.
p5 = pwh pDSC = Prdidas de presin a travs de choques
superficiales.
p6 = pDSC pSEP = Prdidas de presin en lneas de flujo.
p7 = pwf pwh = Prdida de presin total en la tubera de produccin

p8 = pwh pSEP = Prdida de presin total en la lnea de flujo.
Figura 1.1(Sistema de Produccin)
La prdida total de presin en un sistema de produccin es el punto inicial Pr

menos la presin final del fluido, (PR Psep ). El anlisis de las figuras mencionadas,
indican que esta presin es la suma de las prdidas de presin en cada componente que
conforma el sistema.
La presin en cada componente es dependiente del caudal de produccin, el caudal

puede ser controlado por los componentes seleccionados, siendo por lo tanto muy
importante la seleccin y el dimensionamiento de los componentes individuales en el
estudio de un pozo especfico.
El diseo final de un sistema de produccin, debe ser analizado como una unidad,
puesto que, la cantidad de gas fluyente desde el reservorio hasta superficie en un pozo
depende de la cada de presin en el sistema.
El caudal de produccin de un pozo puede muchas veces estar muy restringido por
el comportamiento de uno de los componentes del sistema. El comportamiento total del
sistema puede ser aislado, y optimizado de manera ms econmica. Experiencias pasadas
han mostrado que se gast una gran cantidad de dinero en estimular la formacin, cuando la
capacidad de produccin del pozo es restringido, porque la tubera o lnea de flujo eran
extremadamente pequeas. Otro ejemplo de error en el diseo de terminacin es

sobredimensionar las tuberas. Esto ocurre frecuentemente en pozos que se espera un
caudal de produccin muy alto y cuyo resultado no es el esperado. En esta prctica, no slo
se invierte dinero en equipamiento innecesario, si no tambin se obtiene un escurrimiento
en el pozo lo cual nos reduce la produccin.
La Interrelacin entre caudal y presin es aprovechada por el Anlisis Nodal para

resolver muchos problemas que se presentan con la excesiva resistencia al flujo y las
variaciones en el caudal durante la vida productiva del pozo, en la etapa de surgencia
natural o en la del levantamiento artificial. En este trabajo, estudiaremos la etapa de
surgencia natural.
1.2.-Anlisis del Sistema de Produccin
La optimizacin de la produccin de pozos de gas para un Sistema de Produccin,

llamado tambin Anlisis Nodal, tiene como objetivo el mejorar las tcnicas de
terminacin, produccin y rendimiento para muchos pozos. Este tipo de anlisis fue
propuesto por Gilbert en 1954, discutido por Nind en 1964 y Brown en 1978.
Se darn algunas definiciones sobre el Anlisis del Sistema de Produccin de

autores reconocidos en la industria del petrleo y gas:
i
En 1978, Kermit E. Brow , El Sistema de Anlisis Nodal esta diseado para
combinar varios componentes de pozos de petrleo y gas para predecir los caudales y para
optimizar los tantos componentes de un sistema
ii
En 1979, Joe Mach, Eduardo Proao, Kermit E. Brown , Un Sistema de Anlisis
Nodal fu presentado para evaluar efectivamente un sistema completo de produccin. Son
considerados todos los componentes, comenzando desde la presin de reservorio y
finalizando en el separador.
iii
En 1985, Kermit E. Brown, James F. Lea , Anlisis Nodal fu definido como un
sistema de optimizacin de pozos de petrleo y gas, y es utilizado como una herramienta
excelente para evaluar a fondo un sistema completo de produccin y optimizar el caudal de
flujo objetivo.
iv
En 1987, R.M Frear Jr., and J.P. Yu and J.R. Blair , Anlisis Nodal es un anlisis
de sistema que puede ser utilizado para optimizar un sistema de produccin para pozos de
petrleo y gas. Cada componente en un pozo de un sistema de produccin es incluido y
analizado para aumentar el caudal para ver una mejor rentabilidad econmica.
v
En 1991, H. Dale Beggs , Un Sistema de Anlisis Nodal, es un mtodo muy
flexible que puede ser utilizado para mejorar el desempeo de muchos sistemas de pozos.
Para aplicar un procedimiento de anlisis de un sistema para un pozo, es necesario calcular
la cada de presin que ocurrir en todos los componentes del sistema mencionados en la
figura 1.1.
Desde que el Anlisis Nodal TM fue propuesto la esencia del contenido se

mantiene. Podemos observar en los conceptos presentados anteriormente, por los distintos
autores en diferentes pocas, que el avance y desenvolmiento en la tecnologa de la
computacin, permite hacer clculos exactos y rpidos de algoritmos complejos y
proporciona resultados fcilmente entendidos, adems que este tipo de anlisis se vuelve
popular en todo tipo de pozos de petrleo y gas. El Anlisis Nodal, es el procedimiento de
anlisis que requiere un sistema.
Primero, la colocacin de los nodos, que se pueden encontrar en diferentes partes

del pozo.
Segundo, la generacin del grfico nodal, presin en el nodo versus el caudal como
una herramienta visual para determinar los parmetros dominantes del pozo. Estas curvas
generadas independientemente para cada segmento, son interceptadas en un punto que
indica la capacidad productiva del sistema para un caso particular de componentes. (ver
Figura 1.2). El anlisis de esta figura muestra que la curva de flujo de entrada (inflow)
representa las presiones (aguas arriba) del nodo y la curva de flujo de salida (outflow)
representa las presiones (aguas abajo) del nodo.
Optimizacin del Dimetro de Tubera

Inflow Tuberia 2.445 Ouflow Linea de 3
9000
8000
7000
Presin psi
6000
5000
Curva de Entrada
4000
3000
2000
Curva de Salida Caudal de mximo
1000
de Operacin
0
0,00 2,00 4,00 6,00 8,00 10,00 12,00
Caudal de Gas MMpcd
Figura 1.2 (Determinacin de la capacidad de flujo)

1.3- Nodo
Un nodo es el punto donde existe un cambio en el rgimen o direccin de flujo. Los

cuales se pueden dividir en nodo Comn y nodo fijo.
1.3.1.1 Nodo Fijo
Son los puntos terminales e inicial del sistema de produccin, donde no existe una cada de
Presin.
1.3.2.1 Nodo Comn
Este es el nombre que recibe una seccin determinada de un sistema de produccin

donde se produce una cada de presin, las cadas de presin estn expresadas por
ecuaciones fsicas o matemticas que relacionan la presin y caudal. La Figura 1.3 muestra
los nodos comn y fijos que se utilizan con ms frecuencia.
Todos los componentes aguas arriba del nodo, comprenden la seccin de flujo de
entrada (inflow), en cuanto a la seccin de flujo de salida (outflow) agrupa todos los
componentes aguas abajo.
Es importante notar que para cada restriccin localizada en el sistema, el clculo de la

cada de presin a travs del nodo, como una funcin del caudal, esta representado por la
misma ecuacin general:
P= Q n
Una vez el nodo es seleccionado, se realiza un balance de presiones que representan al

nodo:
Entrada (Inflow) al nodo:
PR P (componentes aguas.arriba ) = Pnodo EC.(1.1)
Salida (Outflow) del nodo
Psep + P (componentes aguas.abajo ) = Pnodo EC.(1.2)
Estas relaciones deben cumplir los siguientes requisitos:
1) El caudal que ingresa al nodo debe ser igual al de salida.

2) Solamente existe una presin en el nodo.
1.4 .- Elementos usados en el Sistema del Anlisis Nodal
Considerando las variadas configuraciones de pozos de un sistema de produccin,

estos elementos, tambin llamados componentes, pueden ser muchos debido a que existen
sistemas muy complejos de terminacin. Los ms comunes estn representados en la Figura
1.3.
1.4.1.- Ubicacin de los Nodos componentes
Observando la Figura 1.3, podemos determinar las posiciones de los nodos

componentes ms comunes, siendo estos modificados de acuerdo a necesidades y
requerimientos del sistema de produccin o polticas de produccin adoptadas.
Nodo Posicin Tipo
10 Lnea de Petrleo al Tanque Fijo

9 Lnea de venta de gas Fijo
8 Separador Fijo
7 Lnea de Flujo Horizontal Comn
6 Choque Superficial Comn
5 Cabeza de Pozo Fijo
4 Restricciones o choque de fondo Comn
3 Tubera Vertical o Inclinada Comn
2 Vlvula de Seguridad Comn
1 Presin Fluyente de Reservorio Fijo
Fig. 1.3 (Componentes del sistema de produccin)

1.4.2 .- Componentes que intervienen en el Anlisis Nodal
En funcin a la necesidad que se tiene de cada uno de los elementos que intervienen
como componente de un sistema de produccin, definiremos la funcionalidad de los ms
importantes.
1.4.2.1.- Separador. En el proceso de separacin de petrleo y gas en los campos, no existe

un criterio nico para establecer las condiciones ms adecuadas de produccin optima de
los equipos. El anlisis nodal TM, esta orientado a obtener ciertos objetivos puntuales que
nos den condiciones de mxima eficiencia en el proceso de separacin; obteniendo de esta
manera:
Alta eficiencia en el proceso de separacin de gas Petrleo
Mayor incremento en los volmenes de produccin
Incremento en la recuperacin de lquido
Disminucin de costos por compresin
Estabilizacin de gas-condensado
1.4.2.2.- Lnea de Flujo Horizontal. Este componente, es el que comunica la cabeza del
pozo con el separador y donde el fluido presenta un comportamiento que obedece a las
condiciones adoptadas para el sistema de produccin de los pozos.
El tratamiento del componente para flujo horizontal, puede ser analizado usando las
diversas ecuaciones y correlaciones presentadas por investigadores que han estudiado la
incidencia, que puede tener este componente, sobre el conjunto del sistema en su
interrelacin con los dems nodos.
1.4.2.3.- Lnea de Flujo Vertical. Este componente es el que comunica el fondo del pozo
con la superficie, donde el fluido presenta un comportamiento que obedece a las
condiciones de presin y temperatura, que estn de acuerdo a la profundidad. En este
componente existe la mayor prdida de energa del sistema, que va desde el 20 al 50 % de
acuerdo a su relacin gas / condensado y corte de agua.
1.4.2.4.- Choque Superficial. Es el que controla la produccin del pozo con el cual se
puede aumentar o disminuir el caudal de produccin, siendo que en este componente se
produce una presin diferencial que puede ser calculada con una de las muchas ecuaciones
para choques o estranguladores.
1.4.2.5.- Cabeza de Pozo. Es un punto del sistema en el que se produce el cambio de

direccin, de flujo vertical a flujo horizontal, y de donde se toma el dato de la presin de
surgencia para conocer la energa de produccin del pozo, siendo tambin un punto crtico
que es tomado en cuenta para su anlisis dentro del sistema.
1.4.2.6.- Vlvula de Seguridad. Este componente, es un elemento que se instala en la

tubera vertical y que opera en cualquier anormalidad del flujo que puede ocurrir en el
transcurso de la produccin, siendo vital para la seguridad operativa del pozo.
1.4.2.7.- Choque de fondo. De acuerdo a la necesidad de elevar la presin o controlar la

energa en el flujo de la lnea vertical, as como tambin, tener una presin de aporte y
elevacin controlada, se procede a la bajada de este tipo de restriccin, por lo que se va
producir una presin diferencial en la que se tendr una cada de presin que a su ves puede
ser calculada.
1.4.2.8.- Presin fluyente. Esta es muy importante para el sistema, ya que de ella depende
toda la capacidad de la instalacin que se desea conectar al reservorio a travs del pozo y
as producir todo el campo.
Esta presin, es medida en el fondo del pozo y tomada en el punto medio del nivel
productor. Su determinacin se la hace en una forma indirecta utilizando herramienta
mecnica o electrnica de toma de presin, o tambin se la puede calcular utilizando
correlaciones.
1.4.2.9.- Completacin o Perforaciones en el Fondo. Este nodo es muy importante en el

sistema de produccin debido a que comunica el reservorio con el pozo, y de l depende
mucho el potencial de entrega de pozo, debido a la disminucin del rea por donde debe
pasar el fluido, la cual puede ser expresada por correlaciones.
1.4.3.- Presin Constante

El nodo 8, ubicado en un sistema de produccin en el separador, establece que existen
dos presiones que no estn en funcin del caudal de produccin del reservorio. La presin
de separacin es usualmente regulada a una presin de entrega de gas, planta o la presin
de succin del compresor nodo 8. Por lo tanto, la presin del separador Psep ser ( )
( )
constante para cualquier caudal de flujo. La presin del reservorio PR , nombrada por el
nodo 1, ser tambin considerada constante en el momento de la prueba o anlisis. El
balance de presin para el nodo en el choque se puede definir como:
(Psep = Pr es Pcomplet. Ptub.vert. Pchoque Ptub.horz.) EC.(1.3)
1.4.4.- Anlisis del sistema en el fondo de pozo
Si colocamos el nodo solucin en el fondo de pozo, esto nos permite aislar el

reservorio de las tuberas tanto vertical como horizontal; dando la posibilidad de estudiar
varios efectos, podemos estudiar la sensibilidad al dimetro de tubera manteniendo los
parmetros de reservorio constante y la sensibilidad de los parmetros de reservorio como
la permeabilidad, dao, conductividad. Ver Figuras 1.5 y 1.6.
La ecuacin de flujo de entrada y salida respectivamente son:

Pr eserv. = (Psep + Ptub.horz + Pchoque + Ptub.vert + Pperf .) EC.(1.4)
Entrada = Salida
Sensibilidad al Dao de Formacin
5000
4500
Dao 3,5
4000 Dao 10
Dao 0
3500 Outflow
Presin Psi
3000
2500
2000
1500
1000
500
0
0,0 5,0 10,0 15,0 20,0 25,0 30,0 35,0 40,0 45,0 50,0
Caudal de Gas MMPCD
Figura 1.5 (Sensibilidad al Dao de Formacin)
Sensibilidad de la Permeabilidad
9000
8000
permeabilidd 23 md
7000
Permeabilidad 50 md
6000
Presin Psi
outflow
permeabilidad 100 md
5000
4000
3000
2000
1000
0
0,0 20,0 40,0 60,0 80,0 100,0 120,0 140,0
Caudal de Gas MMPCD
Figura 1.6 (Sensibilidad a la Permeabilidad del Reservorio)

1.4.5 Optimizacin de la tubera de produccin

Uno de los componentes ms importantes en un sistema de produccin, es la sarta
de produccin. Debido a que cerca del 50 % de la prdida total de presin en un pozo de
gas puede ocurrir por la movilizacin de los fluidos desde el fondo del pozo hasta la
superficie. Un problema comn en los proyectos de completacin, es el seleccionar un
tamao de tubera de produccin basados en critrios totalmente irrelevantes, como por
ejemplo, el tamao que se tiene disponible en almacn. La seleccin del tamao del tubera
de produccin debe ser hecha en base a datos disponibles, ya sea pruebas de formacin o
datos de reservorio, lo cual no es posibles hacerlos en pozos exploratorios por falta de
informacin confiable.
A medida que el rea de flujo vertical se incrementa, las velocidades de flujo

disminuyen pudiendo llegar a generar que las condiciones de flujo sean inestables e
ineficientes, esto ocasiona que se forme un escurrimiento de lquido, formndose la
acumulacin de lquido en el fondo del pozo, que podra ocasionar el ahogo o muerte del
pozo. Una situacin similar se presenta en pozos de baja productividad y dimetros
excesivo de tubera, (Figura 1.7). Por el contrario, en las tuberas de produccin muy
pequeos el caudal de produccin es restringido a causa de la prdida excesiva de friccin.
Un problema comn que ocurre en la completacin de pozos de alto potencial, es el

de instalar tuberas de produccin con dimetros excesivos para mantener la seguridad.
Esto con frecuencia es contraproducente, ya que disminuye la vida til de los pozos, a
medida que la presin del reservorio decrece, los lquidos comienzan a escurrirse por falta
de velocidad del gas para arrastrar los lquidos en fondo.
La respuesta de la capacidad de produccin con la variacin del rea de flujo, es

muy importante para poder definir el dimetro de tubera que se deba bajar a un pozo, ya
que para dos dimetros distintos de tubera obtendremos distintos caudales. Por ejemplo, si
tenemos un dimetro d 2 mayor a d1 , el caudal q 2 aumenta un porcentaje con respecto al
caudal q1 ; quiere decir, que estamos frente a un pozo restringido por el comportamiento
de flujo de salda (outflow). La severidad de la restriccin, depender del porcentaje del
incremento del caudal con un cambio del tamao de la sarta. Por el contrario, para un
d 2 d1 el caudal q 2 es aproximadamente igual al caudal q1 , no se justificarn el costo
de una inversin para un cambio de tamao de tubera ver ( Figura 1.7.)
Optimizacin de Tubera Vertical y Lnea Horizontal
5000
4500
4000
Outflo
Presin (psia)
3500
Diametro: 2,445
3000
2500 Diametro: 2
2000 Diametro: 3
1500
Diametro: 4
1000
Inflow Diametro: 4
500
0 Diametro: 3
0 5 10 15 20 25 30 35 40 45 50 55 60
rocio 65 70
Caudal (MMPCD)
Figura 1.7 (Sensibilidad a los dimetros tubera y lnea de produccin)
1.4.6- Efecto de Agotamiento del Reservorio
Al aislar los componentes de las tuberias tanto vertical como horizontal, podemos
observar el efecto de Agotamiento del reservorio, con su disminucin de su capacidad
productiva, conforme transcurre el tiempo. Teniendo en cuenta los cambios de la relacin
gas-condensado RGC y el corte de agua.
Las intersecciones de las curvas aguas arriba y aguas abajo para las misma
condicioes de la presin de reservorio, da como resultado las capacidades de produccin
para esta relacion ( ver figura 1.8 ).
Mantener la produccin en un caudal constante, implicaria una disminucin de la

presin de fondo fluyente a medida que la presin del reservorio declina. Existen dos
formas para lograr esto:
la primera, es instalando un compresor para reducir la presin del separador.

la segunda, es instalando una lnea de flujo y tuberas de mayor dimetro para
disminuir la cada de presin en el sistema de tuberas.
Agotamiento de Reservorio
4000
Presin: 3694
Presin: 2000
3000 Presin: 1500
Presin (psia)
Presin: 1300
Rel. Agua:
2000
203 004
1000
0
0 2 4 6 8 10 12 14
Caudal (MMPCD)
Figura 1.8 ( Agotamiento de Reservorio)
1.5.- Anlisis del Sistema Nodo en Cabeza de Pozo
Con la ubicacin del nodo de solucin en la cabeza del pozo ( nodo 5 ), la lnea de
flujo horizontal esta aislada facilitando el anlisis de cambio de dimetro de la misma y de
la cada de presin en la lnea o conducto.
Nuevamente el sistema total es dividido en dos componentes, constituyendo el

separador y la lnea de flujo horizontal como un componente, y el reservorio ms la sarta de
tubera vertical como un segundo componente; ver la figura 1.10 muestra,.el primer
componente. La lnea de flujo empieza con la presin de separacin incrementandose, la
presin en la lnea de acuerdo a la prdida de presin debido a los efectos de friccin y
aceleracin, determinandose la presin final en cabeza de pozo para mover el caudal
asumido. La Figura 1.11, muestra el segundo componente del sistema; la linea de flujo
empieza con la presion de reservorio, la cual va disminuyendo de acuerdo a las
restrinciones encontradas, primeramente, se debe descontar la prdida de presin obtenida
en las perforaciones en el caso que el pozo este completado, luego se descuenta la prdida
de presin por elevacin, friccin y aceleracin obtenida en la tubera vertical encontrando
la presin en cabeza para cada caudal asumido. Las tablas de clculos de las prdidas de
presin y los procedimientos de clculos los mostramos, ms explcitamente, en el captulo
5 de este libro.
La presin del nodo para este caso esta dada por:
Entrada (Inflow) al nodo:
Pwh = PR Pres Ptub EC.(1.5)
Salida (Outflow) de nodo
Pwh = Psep + Plf EC.(1.6)
Procedimiento de clculo:
Asumir varios valores de q sc , y determine el correspondiente Pwf de los mtodos

de inflow performance.
Determine la presin de cabeza del pozo, Pwh correspondiente para cada q sc y
Pwf determinada en el paso 1.
Trazar un grfico Pwh vs q sc .
Utilizando una presin fija de separador y las ecuaciones en las tuberas de flujo,
calcular Pwh para varios caudales de flujo asumidos.
Trazar un grfico Pwh vs q sc en el mismo grfico que en el paso 3. La
interseccin da solamente el valor de Pwh y q sc para un dimetro de lnea que ira
a satisfacer ambos subsistemas.
1.6.- Anlisis del Sistema en el Separador
Con la ubicacin del nodo en el separador se puede dividir el sistema en dos

componentes, para optimizar la presin de separacin, con los distintos dimetros de
choques en el caso de que existan. El primer componente del sistema es el separador. El
segundo componente del sistema muestra, el reservorio, tubera y lneas de flujo. La Figura
1.12 nos muestra el efecto de la presin de separacin para los distintos choques y el
mximo caudal que podramos obtener. La solucin es obtenida haciendo el grfico
Psep q sc , como Psep calculado para la relacin:
vs
Psep = PR Pres Ptub Plf EC.(1.7)
Procedimiento de clculo:
Comenzar con la presin de reservorio para calcular la presin de fondo fluyente

correspondiente para cada caudal asumido.
Determinar la presin de cabeza para cada Pwf y q sc del paso 1, haciendo uso de
una correlacin de flujo vertical.
Con la presin de cabeza del paso 2, establecer la presin del separador Psep
,respectiva y permisible para cada caudal.
Trazar un grfico Psep vs q sc y determinar q sc para varios valores de Psep .
Al igual que en el punto anterior 1.5, estos procedimiento y clculos, los mostramos
ms explcitamente en capitulo 5.
Figura 1.10 Componentes de separador y lnea de

flujo horizontal
Figura 1.11 Componentes del reservorio y tubera vertical

El incremento o reduccin de presin del separador, esta ligado al comportamiento

del sistema de tubera y en particular a la lnea de flujo. Al disminuir la presin del
separador se logra un incremento en el caudal del pozos y para los pozos de de alta
productividad se ve reflejado mucho mejor. Muchas veces existe el criterio errneo de
producir un pozo bajo condiciones de flujo subcrtico, siendo mejor producir bajo
condiciones crticas as eliminando el efecto de contrapresion del separador al reservorio,
dejando baches de lquido en el fondo.
En pozos con baja productividad, el componente restrictivo puede ser el mismo

reservorio y un cambio de presin del separador tendr un efecto insignificante sobre el
caudal, porque adicionales cadas de presin ofrecen pequeos incrementos en la
produccin.
Sensibilidad del choque en el Separador
3500
3000
2500
Presin (psia)
2000
1500
1000
500
0
0 1 2 3 4 5 6 7 8 9 10 11 12
Caudal (MMPCD)
Nro. Choque 10 Nro. Choque 30 Nro. Choque 36 Nro. Choque 40
Figura 1.12 (Efecto de presin en el separador para distintos choques)
1.7.- Seleccin del Compresor

La seleccin y el dimensionamiento de un compresor para aumentar la capacidad
productiva de un sistema de pozos, requiere conocer la presin de succin y descarga
requerida, adems del volumen de gas para la venta, y la distancia donde se debe entregar el
gas, que es usualmente fijada. En base a todos estos datos requeridos, determinamos la
descarga y succin del compresor que esta en funcin al caudal de gas.
La presin del separador controla la presin de succin del compresor y est

directamente relacionada con la potencia del compresor (HP) estimada de la siguiente
manera:
HP = 23(R )1 n (n )(q ) EC.(2.8)
Donde:
Pr esin de descarg aabsoluta

R = Re lacin de compresin =
Pr esin de succinabsoluta
n = nmero de etapas
Para R > 4.5 use 2 etapas
Para R > 20 use 3 etapas
q = MMpcsd
El siguiente procedimiento, es usado para determinar los parmetros de diseo

necesarios y la potencia requerida para entregar una cantidad de gas a una presin fijada en
la lnea de venta
Comenzar con la PR , determine Psep para varios valores de q sc usando el
procedimiento para determinar el efecto de la presin del separador.
Trazar un grfico Psep vs q sc
Comenzar con la presin de lnea de venta, determine la presin de descarga
requerida en el compresor, Pdis para varios valores de caudal de flujo.
Trazar un grfico Pdis vs q sc en el mismo grfico tal como fue usado en el paso
2. La interseccin de estas curvas da la capacidad de flujo o de entrega.
Seleccione valores de q sc y determine los valores de Pdis , Psep y
P = Pdis Psep para cada q sc .

Determinar la relacin de compresin requerida, R , y la potencia requerida por el
compresor, HP .
1.8.- Anlisis del Sistema para pozos con Restricciones Superficiales
La Figura 1.13 muestra una descripcin fsica del pozo con un choque de superficie
instalado. Puesto que el choque de cabeza esta usualmente representado por el nodo 2 de
acuerdo a la ubicacin de los nodos, mostrada en la Figura 1.3, es seleccionado para
resolver el problema y determinar los caudales posibles para diferentes dimetros de
choque.
La solucin es dividida en dos partes:
1. La primera parte, sigue exactamente el procedimiento descrito en la seccin 1.5

(anlisis del sistema nodo en cabeza de pozo). En este caso, el desempeo de la
curva vertical del IPR representar la presin aguas arriba del nodo 5, Pwh (presin
de cabeza que controla el caudal) y el desempeo de la curva del segmento
horizontal, la presin aguas abajo del nodo 5, PD (presin necesaria para mover el
fluido al separador). As mismo, hemos considerado que no existe cada de presin
en el nodo, y que el caudal que se predice es donde la presin aguas arriba es igual
a la presin aguas abajo (Pwh = PD ) , ver Figura 1.14. Sin embargo, sabemos que
el choque crear una cada de presin en el nodo funcional 5 para cada caudal.
Figura 1.13 (Solucin de Nodo en el choque de superficie)

2.- La segunda parte se aboca a encontrar esta cada de presin, P para luego hacer
un grfico P vs q elaborado sobre la base de los clculos del desempeo del
choque.
La cada de presin para diferentes choques y caudales, se obtiene de la figura 1.14

y se hace un grfico P vs q .
Para diferentes dimetros de choques, calcular la presin de cabeza asumiendo
varios caudales.
Tabular estos datos en una tabla, incluyendo adems los valores de cada de presin
entre presin de cabeza, requerida para mover el caudal asumido a travs del choque
y la presin downstream necesaria para mover el fluido al separador.
Los P tabulados son plasmados en coordenadas cartesianas, para mostrar el
comportamiento del choque sealaremos en la Figura 1.14 diferentes dimetros de
choques.
Las curvas del comportamiento del sistema nos indican el P requerido para varios
caudales, tomando en cuenta el sistema completo desde la salida al separador. Las curvas
de desempeo del choque revelan un P creado para un conjunto de caudales
considerando diferentes tamaos de choques. Los puntos de interseccin de las P creadas
y requeridas representan las soluciones posibles. Por ejemplo, el caudal obtenido a travs
de la configuracin de un pozo sin restricciones, caer en un cierto porcentaje con la
instalacin de un choque en cabeza de un dimetro particular. Al igual que en inciso 1.5,
detallados se muestran en el capitulo 6.
Determinacin de Caudales por diferentes Choques

Presin de Separacin 1200 psi
9000
8000
7000
Presion Psi
6000
5000
4000
3000
2000
1000
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65
Caudal Gas MMPCD
Delta de Presion ck 20 Achond ck 24 Achond

ck 28 Achond ck 32 Achond ck 36 Achond
Presion fondo fluyente Linea Horizontal Presion Cabeza
Figura 1.14 ( Comportamiento total del sistema incluyendo optimizacin de choque)

Referencias Bibliogrficas
Gas Production Operations H. Dale Beggs
The Technology of Artificial Lift Methods Kermit E. Brown, Volume 4

Production Optimization of Oil and Gas Wells by Nodal Systems Analysis
Production Optimization, Using Nodal Analysis H. Dale Beggs
Manual de Anlisis de Pozos Dowell-Schlumberger
Subsurface Engineering Exxon Company, U.S.A. Technical Services

Training.
Nodal System Analysis of Oil and Gas Wells, By Kermit E. Brown, and James F.
Lea, SPE 14014
A Nodal Approach for Applying Systems Analysis to the Flowing and

Artificial Lift Oil or Gas Well, By Joe Mach, Eduardo Proao, Kermit E. Brown,
SPE 8025
Nodal System Analysis of Oil and gas Wells, By Kermit E. Brown, and James F.
Lea, SPE 14014
1
The Technology of Artificial Lift Methods Kermit E. Brown, Volume 4 Production Optimization of Oil
and Gas Wells by Nodal Systems Analysis, pag. 87
2
A Nodal Approach for Applying Systems Analysis to the Flowing and Artificial Lift Oil or Gas Well, By
Joe Mach, Eduardo Proao, Kermit E. Brown, SPE 8025
3
Nodal System Analysis of Oil and gas Wells, By Kermit E. Brown, and James F. Lea, SPE 14014
i
The Technology of Artificial Lift Methods Kermit E. Brown, Volume 4 Production Optimization of Oil
and Gas Wells by Nodal Systems Analysis, pag. 87
ii
A Nodal Approach for Applying Systems Analysis to the Flowing and Artificial Lift Oil or Gas Well, By
Joe Mach, Eduardo Proao, Kermit E. Brown, SPE 8025
iii
Nodal System Analysis of Oil and gas Wells, By Kermit E. Brown, and James F. Lea, SPE 14014
iv
Application of Nodal Analysis in Appalachian Gas Wells, By R.M Frear Jr., Stonewall Gas Co., and J.P.
Yu and J.R. Blair, West Virginia U., SPE 17061
v
Production Optimization, Using Nodal Analysis, By H. Dale Beggs, 1991, pag. 7
Clasificacin de los Fluidos en el Reservorio----------------------------------------------------29
2 Clasificacin de los fluidos en el Reservorio
2.1 Introducin
Las acumulaciones de gas y de petrleo ocurren en trampas subterrneas formadas
por caractersticas estructurales, estratigrficas o ambas. Por fortuna, estas acumulaciones
se presentan en las partes ms porosas y permeables de los estratos, siendo estos
principalmente areniscas, calizas y dolomitas, con las aberturas nter granulares o con
espacios porosos debido a diaclasas, fracturas y efectos de soluciones. Por lo que un
yacimiento est definido, como una trampa donde se encuentra contenido el petrleo, el
gas, o ambas como mezclas complejas de compuestos, como un solo sistema hidrulico
conectado cuyas caracterstica no solo depende de la composicin si no tambin de la
presin y temperatura a la que se encuentra. Muchos de los yacimientos de hidrocarburos se
hallan conectados hidrulicamente a rocas llenas de agua, denominadas acuferos, como
tambin muchos de estos yacimientos se hallan localizados en grandes cuencas
sedimentarias y comparten un acufero comn.
La temperatura de un reservorio es determinada por la profundidad, y el

comportamiento del fluido en un reservorio es determinado por su composicin relacin
PVT. En un reservorio se tiene diferentes clases de fluido, las cuales mostramos en tabla
2.1. Las temperaturas crticas de los hidrocarburos ms pesados son ms elevadas que los
componentes livianos. De all la temperatura crtica de la mezcla de un hidrocarburo
predominantemente compuesto por componentes pesado, es ms alta que el rango normal
de temperatura en el reservorio.
Tabla 2.1 Caractersticas y composicin de los diferentes tipos de

Fluido en el reservorio
Componente Petrleo Petrleo Gas y Gas seco
Voltil Condensado
C1 45.62 64.17 86.82 92.26

C2 3.17 8.03 4.07 3.67
C3 2.10 5.19 2.32 2.18
C4 1.50 3.86 1.67 1.15
C5 1.08 2.35 0.81 0.39
C6 1.45 1.21 0.57 0.14
C7+ 45.08 15.19 3.74 0.21
PM C7+ 231.0 178.00 110.00 145.00
Dens. Relativa 0.862 0.765 0.735 0.757
Color del Negro Anaranjado Caf Ligero Acuoso

Verdoso Oscuro
Lquido
29
Cuando la presin de reservorio cae por debajo del punto de saturacin, el diagrama
de fase del fluido original no es representativo, ya que el gas y lquido son producidos a
razones diferentes a la combinacin original, resultando un cambio en la composicin del
fluido. La segregacin gravitacional de las dos fases con diferentes densidades tambin
podra inhibir el contacto entre las dos fases previendo el equilibrio en el reservorio.
Los reservorios de hidrocarburos son clasificados de acuerdo a:
La composicin de la mezcla de hidrocarburos en el reservorio.

La presin y temperatura inicial del reservorio.
La presin y temperatura de produccin en superficie.
El comportamiento termodinmico de una mezcla natural de hidrocarburos, puede

ser utilizado para propsitos de clasificacin, tomando como base del diagrama el
comportamiento de las fases.
2.2.- Diagrama de Fases (Presin- Temperatura)
Un tpico diagrama de Temperatura y Presin es mostrado en la Figura 2.1. Estos

diagramas son esencialmente utilizados para:
Clasificar los reservorios.

Clasificar naturalmente el sistema de hidrocarburos.
Describir el comportamiento de fases del fluido
La Figura 2.1 presenta los siguientes elementos: La curva llamada envolvente de

fases, que resulta de unir las curvas de punto de burbuja y punto de roco que muestra la
mezcla para diferentes temperaturas; curvas que se unen en el punto denominado crtico. La
envolvente de fases divide el diagrama en tres regiones, la primera llamada regin de
lquidos, est situada fuera de la fase envolvente y a la izquierda de la isoterma crtica. La
segunda llamada regin de gases, se encuentra fuera de la fase envolvente y esta a la
derecha de la isoterma crtica; La tercera y ltima, encerrada por la fase envolvente, se
conoce como regin de dos fases, en esta regin, se encuentran todas las combinaciones de
temperatura y presin en que la mezcla de hidrocarburo puede permanecer en dos fases en
equilibrio, existiendo dentro de ella, las llamadas curvas de igualdad, que indican un
porcentaje de total de hidrocarburo que se encuentra en estado lquido y gaseoso. Todas
estas curvas inciden en un punto crtico. Se distinguen, adems, en el mismo diagrama, la
cricondetrmica y la cricondenbrica, las cuales son la temperatura y la presin mximas,
respectivamente, que en la mezcla de hidrocarburos pueden permanecer en dos fases en
equilibrioi.
Para un mejor entendimiento de la Figura 2.1 se darn todas las definiciones y

algunos conceptos bsicos asociados con el diagrama de fase.
30
Figura 2.1 (Diagrama de fase (Presin Temperatura))
2.2.1- Propiedades intensivas.- Denominados a aquellos que son independientes de la

cantidad de materia considerada como ser: la viscosidad, densidad, temperatura, etc.
funcin principal de las propiedades fsicas de los lquidos.
2.2.2- Punto Crtico.- Es el estado a condicin de presin y temperatura para el cual las
propiedades intensivas de las fases lquidas y gaseosas son idnticas, donde cuya
correspondencia es la presin y temperatura crtica.
2.2.3- Curva de Burbujeo (ebullicin) .- Es el lugar geomtrico de los puntos, presin

temperatura, para los cuales se forma la primera burbuja de gas, al pasar de la fase
lquida a la regin de dos fases, siendo este estado el equilibrio de un sistema
compuesto de petrleo crudo y gas, en la cual el petrleo ocupa prcticamente todo
el sistema excepto en una cantidad infinitesimal de gas.
El yacimiento de punto de burbujeo se considera cuando la temperatura normal est

debajo de la temperatura crtica, ocurriendo tambin que a la bajada de la presin
alcanzar el punto de burbujeo.
2.2.4- Curva de roco (condensacin) .- Es el lugar geomtrico de los puntos, presin

temperatura, en los cuales se forma la primera gota de lquido, al pasar de la regin
de vapor a la regin de las dos fases. El punto de roco es anlogo al punto de
burbuja, siendo el estado en equilibrio de un sistema que est compuesto de petrleo
y gas, lugar en la cual el gas ocupa prcticamente todo el sistema dando excepcin a
cantidades infinitesimales de petrleo.
31
2.2.5- Regin de dos fases.- Es la regin comprendida entre las curvas de burbujeo y roco
(cricondenbara y cricondenterma). En esta regin coexisten en equilibrio, las fases
lquida y gaseosa.
2.2.6- Cricondenbar .- Es la mxima presin a la cual pueden coexistir en equilibrio un

lquido y su vapor.
2.2.7- Cricondenterma .- Es la mxima temperatura a la cual pueden coexistir en

equilibrio un lquido y su vapor.
2.2.8- Zona de Condensacin Retrgrada .- Es aquella cuya zona est comprendida

entre los puntos de las curvas cricondenbar y cricondenterma (punto crtico y punto
de roco), y que a la reduccin de presin, a temperatura constante, ocurre una
condensacin.
2.2.9- Petrleo Saturado .- Es un lquido que se encuentra en equilibrio con su vapor

(gas) a determinada presin y temperatura. La cantidad de lquido y vapor puede ser
cualesquiera. En este sentido la presin de saturacin es la presin a la cual lquido
y vapor estn en equilibrio. En algunos casos la presin de burbujeo o presin de
roco puede usarse sinnimamente como presin de saturacin.
2.2.10- Petrleo Bajo Saturado .- Es el fluido capaz de recibir cantidades adicionales de

gas o vapor a distintas condiciones de presin y temperatura, en un fluido no
saturado, la disminucin de la presin no causa liberacin de gas existentes en
solucin en el fluido.
2.2.11- Petrleo Supersaturado .- Es aquel fluido que a condiciones de presin y

temperatura que se encuentra, tiene una mayor cantidad de gas disuelto que el que le
correspondera en condiciones de equilibrio.
2.2.12- Saturacin crtica de un Fluido .- Es la saturacin mnima necesaria para que

exista escurrimiento de dicho fluido en el yacimiento.
Inicialmente toda acumulacin de hidrocarburos tiene su propio diagrama de fases

que depende slo de la composicin de la mezcla. De acuerdo a esto, los yacimientos de
hidrocarburos se encuentran inicialmente, ya sea en estado monofsico (A, B, y C) o en
estado bifsico (D), de acuerdo con la composicin relativa de sus presiones y temperaturas
en los diagramas de fases.
Cuando la presin y la temperatura iniciales de un yacimiento caen fuera de la

regin de dos fases pueden comportarse:
1.- Como yacimientos normales de gas (A), donde la temperatura del yacimiento excede
el cricondentrmico.
32
2.- Como yacimiento de condensado retrgrado (de punto de roco) (B), donde la
temperatura del yacimiento se encuentra entre la temperatura crtica del punto
cricondentrmico.
3.- Como yacimientos de petrleo bajo-saturado (de punto burbujeo) donde, la

temperatura del yacimiento est debajo de la temperatura crtica.
Cuando la presin y la temperatura iniciales del yacimiento caen dentro de la regin

de dos fases pueden comportarse:
1.- Como yacimientos de petrleo saturado, donde, existe una zona de petrleo con un
casquete de gas.
2.- Como yacimiento de petrleo saturado sin estar asociados a un casquete de gas, esto
es, cuando la presin inicial es igual a la presin de saturacin o de burbujeo. La presin y
temperatura para este tipo de yacimientos se localizan exactamente sobre la lnea de
burbujeo (E).
2.3.- Clasificacin de los reservorios

Se aclara que el estado fsico de un fluido de yacimiento generalmente vara con la
presin, pues la temperatura es esencialmente constante. Es prctica comn clasificar a los
yacimientos de acuerdo a las caractersticas de los hidrocarburos producidos y a las
condiciones bajo las cuales se presenta su acumulacin en el subsuelo. As, tomando en
cuenta las caractersticas de los fluidos producidos, se tienen reservorios de:
Reservorio de Petrleo
Reservorio de Gas
2.3.1.- Reservorio de Petrleo
Si la temperatura del reservorio T es menor que la temperatura crtica Tc del fluido

del reservorio, el reservorio es clasificado como reservorio de petrleo. Dependiendo de la
presin inicial del reservorio P1 , los reservorios de petrleo pueden ser subclasificados en
las siguientes categoras:
2.3.1.1. Reservorio de Petrleo Subsaturado
Si la presin inicial del reservorio Pi, es igual est representada en la Figura 2.2 por
el punto 1, y mayor que la presin del punto de burbuja, Pb, y la temperatura esta por bajo
de la temperatura critica del fluido del reservorio.
33
2.3.1.2.- Reservorio de Petrleo Saturado

Cuando la presin inicial del reservorio esta en el punto de burbuja del fluido del
reservorio, como mostramos en la Figura 2.2, punto 2, el reservorio es llamado reservorio
saturado de petrleo.
2.3.1.3.- Reservorio con Capa de Gas

Si la presin inicial del reservorio es menor que la presin en el punto de burbuja
del fluido del reservorio, como indica en el punto 3 de Figura 2.2, el reservorio es
predominado por una capa de gas en la zona de dos fases, la cual contiene una zona de
lquido o de petrleo con una zona o capa de gas en la parte superior.
Figura 2.2 (Diagrama de Fase (Presin y Temperatura))
En general el petrleo es comnmente clasificado en los siguientes tipos:
Petrleo negro
Petrleo de bajo rendimiento
Petrleo de alto rendimiento (voltil)
Petrleo cerca al punto crtico
34
2.3.2.- Petrleo Negro

El diagrama de fase nos muestra el comportamiento del petrleo negro en la Figura
2.3, en la cual se debe notar qu lneas de cualidad son aproximadamente equidistantes
caracterizando este diagrama de fase de petrleo negro. Siguiendo la trayectoria de la
reduccin de presin indicada por la lnea vertical EF, la curva de rendimiento de lquido
mostrado en Figura 2.4, es el porcentaje de volumen lquido en funcin de la presin. La
curva de rendimiento de lquido se aproxima a la lnea recta, excepto las presiones muy
bajas. Cuando el petrleo negro es producido normalmente se tiene una relacin gas
petrleo entre 200 1500 PCS/STB y la gravedad del petrleo esta entre 15 40 API. En
el tanque de almacenamiento el petrleo normalmente es de color marrn a verde oscuro.
Figura 2.3 (Diagrama de Fase petrleo negro (Presin y Temperatura))
Figura 2.4 ( Curva del rendimiento liquido para petrleo negro)
35
2.3.3.- Petrleo Negro de bajo rendimiento
El diagrama de fase para un petrleo de bajo rendimiento es mostrado en la Figura

2.5. El diagrama es caracterizado por las lneas de calidad que estn espaciadas
estrechamente cerca de la curva de roci. En la curva de rendimiento de lquido (Figura
2.6) se muestra las caractersticas de rendimiento de esta categora de petrleo. Las otras
propiedades de este tipo de petrleo son:
Factor volumtrico de la formacin de petrleo menor que 1,2 bbl/STB

Relacin Gas Petrleo menor que 200 pcs/STB
Gravedad del petrleo menor que 35 API
Coloracin negro
Recuperacin substancial de lquido a condiciones de separacin como es
indicado por el punto G sobre o 85% de lnea de cualidad de la Figura 2.5
Figura 2.5 (Diagrama de fase para petrleo de bajo Rendimiento)
Figura 2.6 (Curva de Rendimiento para bajo rendimiento de Petrleo)
36
2.3.4.- Petrleo Voltil
El diagrama de fase para un petrleo voltil (alto rendimiento) es dado en la Figura

2.7. Observndose que las lneas de calidad estn juntas y estrechas cerca del punto de
burbuja y estn ms ampliamente espaciadas a bajas presiones. Este tipo de petrleo es
comnmente caracterizado por un alto rendimiento de lquido inmediatamente por debajo
del punto de burbuja como es mostrado en la Figura 2.8. Las otras propiedades
caractersticas de este petrleo comprende:
Factor volumtrico de la formacin menor que 2 bbl/STB

Relacin Gas Petrleo entre 2000 3200 PCS/STB
Gravedad del petrleo entre 4,5 55 API
Baja recuperacin de lquido a las condiciones de separador como es indicado
en el punto G en Figura 2.7.
Color verdoso para naranja
Figura 2.7 (Diagrama de fase para petrleo voltil de alto rendimiento)
Figura 2.8 (Curva de rendimiento de liquido para petrleo voltil)
37
2.3.5.- Petrleo Cerca al punto crtico

Si la temperatura de reservorio Tr esta cerca de la temperatura Tc del sistema de
hidrocarburo mostrado en la Figura 2.9, la mezcla de hidrocarburos es identificada como
petrleo cerca al punto crtico. Porque todas las lneas de calidad convergen al punto
crtico, una cada de presin isotrmica (como se muestra en la lnea vertical EF, Figura
2.9), puede llevar del 100% de petrleo del volumen poral de hidrocarburo a condiciones
iniciales al 55 % de petrleo al punto de burbuja si decae la presin en un valor de 10 a 50
psi por debajo del punto de burbuja, el comportamiento caracterstico de encogimiento de
petrleo cerca al punto crtico es mostrado en la Figura 2.10. Este petrleo es caracterizado
por un alto GOR ms de 3000 PCS/STB con un factor volumtrico mayor a 2.0 bbl/STB.
Las composiciones de este tipo de petrleo son normalmente caracterizado por 12,5 a 20
%mol de heptano plus, 35% o ms de etano a travs de hexano y el resto en metano.
Figura 2.9 (Diagrama de fase para petrleo cerca al punto crtico)
Figura 2.10 (Curva de rendimiento de lquido para petrleo cerca al punto crtico)
38
2.4.1.- Reservorio de Gas
Con el advenimiento de las perforaciones profundas han sido descubierto

yacimientos de gas a alta presin con propiedades materialmente diferentes de aquellos
yacimientos de gas seco anteriormente encontrados. El fluido del yacimiento esta
compuesto predominantemente por metano, pero se encuentra cantidades considerables de
hidrocarburos pesados.
Si la temperatura de reservorio es mayor que la temperatura crtica del fluido de

hidrocarburo, el reservorio es considerado un reservorio de gas. Los reservorios que
producen gas natural pueden ser clasificados, esencialmente, en cuatro categoras y estas
son:
2.4.2.- Reservorio de Condensacin Retrograda de Gas
Si la temperatura del reservorio Tr esta entre la temperatura crtica Tc y la

cricondetrmica Tct del fluido el reservorio, es clasificado como reservorio de condensacin
retrgrada.
El fluido existe como un gas a las condiciones iniciales del reservorio, cuando la
presin de reservorio declina a una temperatura constante, la lnea del punto de roco es
cruzada y se forma el lquido en el reservorio. Este lquido tambin se forma en el sistema
de tubera en el separador debido al cambio de presin y temperatura. ii.
Considrese que las condiciones iniciales de un reservorio de condensacin

retrgrada de gas es presentado por el punto 1 del diagrama de fases (presin
temperatura) de la Figura 2.11, la presin del reservorio esta por encima de la presin del
punto de roco, el sistema de hidrocarburo, el reservorio muestra una fase simple (fase
vapor). Cuando la presin de reservorio declina isotrmicamente durante la produccin, la
presin inicial (punto 1) cae al (punto 2) que es la presin declinada y esta por encima del
punto de roco; existe una la atraccin entre molculas de los componentes livianos y
pesados, ocasionando su movimiento por separado, esto origina que la atraccin entre los
componentes ms pesados sean ms efectivos de esta manera el lquido comienza a
condensarse.
Este proceso de condensacin retrgrada, continua con la precisin decreciente

antes de que llegue a su mximo condensacin de lquido econmico en el punto 3. La
reduccin en la presin permite a las molculas pesadas comenzar el proceso de
vaporizacin normal. Este es un proceso para lo cual pocas molculas de gas golpean la
superficie lquida y causan que ms molculas entren a la fase lquida. El proceso de
vaporizacin continua cuando la presin de reservorio esta por debajo de la presin de
roci.
39
2.4.3.- Reservorio de Gas-Condensado cerca al punto crtico

Si la temperatura de reservorio esta cerca de la temperatura crtica, como es
mostrado en la Figura 2.12, la mezcla de hidrocarburo es clasificado como reservorio de
gas condensado cerca del punto crtico. El comportamiento volumtrico de esta categora
de gas natural es descrita a travs de la declinacin isotrmica de presin como se muestra
en la lnea vertical 1 3 en la Figura 2.12. Todas las lneas de calidad convergen en el
punto crtico, un aumento rpido de lquido ocurrir inmediatamente por debajo del punto
de roco como la presin es reducida en el punto 2, este comportamiento puede ser
justificado por el hecho de que varias lneas de calidad son cruzadas rpidamente por la
reduccin isotermal de presin.
Figura 2.11 (Diagrama de fase para reservorio de gas con condensacin retrograda)
Figura 2.12 (Diagrama de fase para reservorio de gas condensado cerca del punto crtico)
40
2.4.4.- Reservorio de Gas-Hmedo
El diagrama de fase correspondiente a un reservorio de gas hmedo, se presenta en

la Figura 2.13, en ella se puede observar que la temperatura del reservorio es mayor que la
cricondetrmica de la mezcla, por tal razn nunca se integran las dos fases en el reservorio,
nicamente existe la fase gaseosa en el reservorio, si el reservorio es agotado
isotrmicamente a lo largo de la lnea vertical A B.
El gas producido fluye hacia la superficie, y por ende, la presin y la temperatura de

gas declinar..El gas entra en la regin de dos fases, en la tubera de produccin debido a
los cambios de presin y temperatura y a la separacin en la superficie. Esto es causado por
una disminucin suficiente en la energa cintica de molculas pesadas con la cada de
temperatura y su cambio subsiguiente para lquido a travs de fuerzas atractivas entre
molculas.
Cuando estos fluidos llevados a superficie entran en la regin de dos fases,

generando relaciones gas petrleo entre 50000 y 120000 PCS/ BBLS, l liquido
recuperable tiende a ser transparente, con densidades menores de 0.75 gr/cm3 .iii y los
contenidos de licuables en el gas son generalmente por debajo de los 30 Bbls/MMPC.
Estos yacimientos se encuentran en estado gaseoso cuya composicin predomina un

alto porcentaje de metano que se encuentra entre 75-90 % aunque las cantidades relativas
de los componentes ms pesados son mayores que en el caso del gas seco.
Figura 2.13 (Diagrama de fase para reservorio de gas hmedo)
2.4.5.- Reservorio de Gas-Seco

Este ltimo tipo de reservorio es lo que se conoce como reservorio de gas seco, cuyo
diagrama se presenta en la Figura 2.14. Estos reservorios contienen principalmente metano,
con pequeas cantidades de etano, propano, y ms pesados, el fluido de este reservorio
entran en la regin de dos fases a condiciones de superficie, durante la explotacin del
41
reservorio. Tericamente los reservorios de gas seco no producen lquido en la superficie,

por ende, la diferencia entre un gas seco y un gas hmedo es arbitraria y generalmente en
sistemas de hidrocarburos que produzcan con relaciones gas petrleo mayores de 120000
PCS/ Bbls se considera gas seco. 1
Figura 2.14 (Diagrama de fase para reservorio de gas Seco)
2.5 Correlaciones para determinar el punto de Roco
En un desarrollo o explotacin de un campo gasfero es muy importante conocer la

presin de roco para evitar los problemas de condensacin retrgrada, ya que el mismo
sobre lleva una mala explotacin del reservorio y por ende una baja recuperacin de
condensado con incidencias econmicas no recomendable. Por lo tanto, para explotar un
reservorio gasifero la presin de reservorio no debera caer por debajo de la presin de roco
debido a la condensacin del gas en el reservorio. Si la presin de reservorio es igual a la
presin de roco se debera realizar una inyeccin de gas seco para bajar el punto de roco.
Para la determinacion del punto de roco existen dos correlaciones existente en la

industria petrolera una correlacin esta hecha en base a la composicion de fluido y a las
propiedades del c7+ La segunda correlacin basada en los datos de produccin de
reservorio usualmente disponible. Pero ninguna de esta correlaciones remplazara al estudio
PVT de los fluidos si se dispone de ellas, las misma que deberan si analizadas para ver el
grado de representatividad del fluido.
2.5.1.- Determinacin del punto de roco con la composicin del gas
La preedicin de la presin de roco no es ampliamente practicado debido a la

complejidad del comportamiento de la fase retrgrada, es necesario la determinacin
experimental de la condicin del punto de rocoiv.Sage y Olds, y Et al presentaron distintas
correlaciones para determinar la presin de roci para varios sistemas de condensado.
42
La presin de punto de roco es estimada utilizando la correlacin generada por

Nemeth y Kennedy, que utiliza la composicin y temperaturav. Esta se describe como esa
presin en la cual los fluidos condensados iniciaran la cada de la primera gota de lquido
fuera de la fase gaseosa.
A[0,2 * %N2 + CO2 + %H2S + 0,4 * %Meth + %Eth + 2(% prop + %IBut + %N

[
pd = exp %But + %IPen + %NPen + %NHex ] + B * DenC7 + C * %Meth %C7+ + 0,2 ) + ]

D * T + E * L + F * L2 + G * L3 + H * M + I * M 2 + J * M 3 + K
Donde:
A = 2,0623054 x10 2
B = 6,6259728
C = 4,4670559 x103
D = 1,0448346 x10 4
E = 3,2673714 x10 2
F = 3,6453277 x103
G = 7,4299951x105
H = -0,11381195
I = 6,2476497 x10 4
J = 1,0716866 x106
K = 10,746622
L = ( )(
C 7+ MWC 7+ )
M = MWC7+ (DenC +
7 )
+ 0,0001
0,6882 * % NHep + 0,7068 * % NOct +
DenC7+ = %C7+
0,7217 * % NNon + 0,7342 * % NDec
100,2 * % NHep + 114,2 * % NOct + 128,3 * % NNon +
MWC7+ = %C7+
142,3 * % NDec
%C7+ = % NHep + % NOct + % NNon + % NDec
C7+ = %C7+ 100
La correlacin de Nemeth y Kennedy, es muy sensible a la concentracin de los

compuestos de gas ms pesados. Muchos anlisis de gas normalmente agrupan los
componentes ms pesados en un solo valor. El usuario conseguir un clculo mucho mejor
de la presin del punto de roco utilizando una suposicin adecuada para propagar
componentes ms pesados y repetir ms estrechamente el verdadero anlisis de gas.
43
El rango de propiedades usada para desarrollar esta correlacin incluyen presiones de

roci que van de 1000 a 10000 psi y temperatura que van de 40 a 320 o F y un amplio
rango de composicin de reservorio. La correlacin nos pueden predecir la presin de roci
en un rango de seguridad de +/- 10% para condensado que no contienen gran cantidad de
no hidrocarburo.
Ejemplo Prctico No1. Se tiene la composicin del gas y se desea conocer la presin de
roco. Se tiene una muestra recombinada cuya composicin presentamos en la tabla 2.1 la
presin inicial de reservorio 3916 psi gravedad API en el tanque es 58.
Tabla 2.1
Componente Fraccin
Molar
Metano 87,54 %
Etano 4,54 %
Propano 2,51 %
Iso-Butano 0,35 %
Butano Normal 0,97 %
Iso-Pentano 0,31 %
Pentano Normal 0,39 %
Hexano 0,56 %
Heptano 0,6 %
Octano 0,28 %
Nonanos 0,23 %
Decanos 0,6 %
Nitrgeno 1,03 %
Dixido de Carbono 0,09 %
Gas Sulfhdrico 0,00 %
Total 100
+ 0 , 6882 * % NHep + 0 , 7068 * % NOct +

Den C 7 = % C 7+ = 0 . 71189
0 , 7217 * % NNon + 0 , 7342 * % NDec
2.5.2.- Determinacin del punto de roco basados en datos de produccin de campo.
Esta correlacin esta basada a un paper presentado en Calgary Canada (SPE 75686)
denominada Correlacin para determinar la presin de roco y C7+ para reservorio de Gas
100,2 * % NHep + 114,2 * % NOct + 128,3 * % NNon +
MWC7+ = %C7+ = 121.04
142,3 * % NDec
[
A 0,2 * %N2 + CO2 + %H2S + 0,4 * %Meth+ %Eth + 2(% prop+ %IBut + %N

[
pd = exp %But + %IPen+ %NPen+ %NHex ] + B * DenC7 + C * %Meth %C7+ + 0,2 ) + = 1866 ]

D *T + E * L + F * L + G * L + H * M + I * M + J * M + K
2 3 2 3

44
Condensado en base a pruebas de produccin. y parametros que usualmente se dispone.
Este mtodo primeramente se basa en calcular el %C7+ en funcin a la relacin de

Gas/Condensado en la teoria el autor presenta dos correlaciones las cuales son:
Primera Correlacin %C7+ =f(GCR)

%C7+ =(GCR/70680)-0.8207
Segunda Correlacin %C7+ =f(GCR, SGg)
%C7+ =10260*(GCR*SGg)-0.8207
Correlacin del punto de Roco Pd = f(GCR, %C7+,API, Tr)
Pd = K1* GCRK2 / C7+K3 K8*API(K4*TrK5-K6* C7+K7) )
Los valores de las constantes son las siguientes:
K1= 346,77647
K2= 0,09741
K3= -0,29478
K4= -0,04783
K5= 0,28126
K6= 0,00068
K7= 1,90633
K8= 8,41763
Nomeclatura
%C7+ Porcentaje de heptano superior

Pd Presin de roco ( psi )
GCR Relacin Gas Condensado (pc/bbl)
SGg Gravedad especifica del gas del separador aire=1
Tr Temperatura de Reservorio (oF )
Ki Coeficiente de regrecin
Ejercicio No2 determinar la presin de roco con los siguientes datos de produccin
Tr =183 F Relacin Gas/Condensado 42711pc/bbls, API 58.8, SGg=0.65
%C7+ = (GCR/70680)^-0.8207 1.51194
%C7+ = 10260*(GCR*SGg)^ -0.8499 1.59012
Pd = K1*(GCR^K2/C7^K3 * K8 * API ^ (K4*Tr^K5 - K6*C7^K7))= 4052 psi
45
2.6 Pruebas PVT
Los fluidos encontrados en yacimientos petrolferos son esencialmente mezclas

complejas de compuestos hidrocarburos, que contienen con frecuencia impurezas como
nitrgeno, dixido de carbono y sulfuro de hidrgeno. La tabla 2.1 presenta la composicin
en porcentajes molar de varios lquidos tpicos encontrados en yacimientos, junto con la
gravedad del petrleo fiscal, la razn gas petrleo de la mezcla de yacimientos y otras
caractersticas de tales fluidos. La composicin del petrleo fiscal es completamente
diferente a su composicin a condiciones del yacimiento, debido principalmente a la
liberacin de la mayor parte del metano y etano en solucin y a la vaporizacin de
fracciones de propanos, butanos y pentanos a medida que la presin disminuye al pasar de
condiciones del yacimiento a condiciones atmosfricas normales.
Existen dos mtodos de obtener muestras de fluidos del yacimiento:
Muestreo de Fondo
Se baja un equipo especial de muestreo dentro y hasta el fondo del pozo, sujetado
por un cable con el muestrador, a pozo cerrado, luego se deja fluir el pozo a bajos
caudales para muestrear a condiciones de reservorio.
Muestreo de Superficie
Tomando muestras de gas y petrleo en la superficie y mezclndolas en las debidas

proporciones de acuerdo con la razn gas petrleo medida a tiempo de muestreo. Las
muestras deben obtenerse al comienzo de las operaciones de produccin del yacimiento,
preferiblemente en el primer pozo, para que en esta forma la muestra se representativa
del fluido original que se encuentra en el yacimiento. La composicin del fluido
obtenido en el saca muestras depende de la historia del pozo, anterior de la operacin de
muestreo. Si el pozo no ha sido acondicionado adecuadamente antes de obtener la
muestra, ser imposible obtener muestras respectivas de fluidos del yacimiento.
Kennerly y Reudelhumber, recomiendan un procedimiento para acondicionar
debidamente el pozo. La informacin obtenida del anlisis de una muestra de fluido
incluye generalmente los siguientes datos:
a. Razones Gas en solucin Petrleo y Gas liberado Petrleo y los

volmenes de las fases lquidas.
b. Factores volumtricos, gravedad del petrleo fiscal y razones Gas Petrleo

del separador a condiciones fiscales, para diferentes presiones del separador.
c. Presin del punto de burbujeo de los fluidos del yacimiento.
d. Compresibilidad del petrleo saturado a condiciones del yacimiento.
46
e. Viscosidad el petrleo a condiciones del yacimiento como funcin de

presin.
f. Anlisis fraccional de una muestra de gas obtenida de la cabeza del pozo y

del fluido saturado a condiciones de yacimiento.
Para un anlisis preliminar de un yacimiento, y si no se disponen de datos de

laboratorio, generalmente puede hacerse estimaciones razonables a partir de correlaciones
empricas basadas en datos fciles de obtener. Estos datos incluyen gravedad del petrleo
fiscal, gravedad especfica del gas producido, razn gas petrleo al comienzo de la
produccin, viscosidad del petrleo fiscal, temperatura del yacimiento y posicin inicial del
mismo.
Las variaciones en las propiedades de un fluido del yacimiento, de varias muestras

obtenidas en diferentes partes del yacimiento, son pequeas y no exceden a las variaciones
inherentes a las tcnicas de muestreo y anlisis, esto sucede en la mayora de los
yacimientos. Por otra parte, en algunos yacimientos, particularmente en aquellos con
grandes volmenes de arena, las variaciones en las propiedades de fluidos son
considerables.
2.6.1 TIPOS DE PRUEBAS PVT
Las pruebas de PVT pueden realizarse de 3 maneras:
1. Proceso a composicin constante (masa constante).
2. Proceso a volumen constante.
3. Proceso de liberacin diferencial (petrleo negro).
2.6.1.1.- Proceso a composicin constante: La composicin global no cambia, se

carga a la celda una cantidad de fluido, se expande la celda o el mercurio, se agita para
alcanzar equilibrio y al aumentar el volumen el gas se va liberando. Luego se miden las
variaciones de lquido y volmenes de gas.
GAS
OIL Proceso
GAS
OIL GAS GAS
OIL
Hg Proceso Proceso Proceso Proceso
OIL
Hg OIL
Hg
Hg
Hg
Pb=Pr P2 P3 P4
47
2.6.1.2 Proceso a Volumen constante: (Procedimiento para gas y petrleo voltil). Se

carga cada celda con un volumen suficiente de fluido, primero aumentamos el tamao de la
celda,(sacamos un volumen de mercurio) a ese gas de expansin se lo retira y se mide su
masa su composicin.
Removemos Gas
GAS
GAS Proceso GAS
GAS
OIL
OIL
OIL
OIL Proceso Proceso Proceso
OIL
Hg Hg Hg Hg
Hg
Pb PL--Pc
2.6.1.3. Proceso de Liberacin diferencial:

(Para petrleo negro). En este tipo de prueba se baja la presin, de cada celdas se extrae
todo el gas que se expanda. Para que la prueba tengas valores de la ecuacin de estado hay
que calibrar con la ecuacin de estado.
Removemos todo el Gas
GAS Proceso GAS
OIL OIL OIL

OIL OIL
Proceso Proceso Proceso
Hg Hg Hg
Hg Hg
Pb Presin constante
48
Reservoir Engineering - Tarek Ahmed, 1946

Manual de Explotacin de Yacimientos de Gas y Condensado y de Petrleo
voltil - SPE filial Bolivia, 2000
Gas Production Operations - H. Dale Beggs, 1984
Ingeniera Aplicada de Yacimientos Petrolferos - B.C. Craft y M. F. Hawkins,
1997
Gas Production Engineering - Sunjay Kumar, 1987
Phase Behavior Monograph Volume 20 SPE, Curtis H. Whitson and Michael R.
Brule.
Petroleum Engineering Tool Kit , Programs for Spreadshee,Software, Doug
Boone & Joe Clegg.
i
Manual de Explotacin de Yacimientos de Gas y Condensado y de Petrleo voltil, SPE filial Bolivia,
2000, Pg. 22
ii
Gas Production Operations, H. Dale Beggs, 1984,
iii
Manual de Explotacin de Yacimientos de Gas y Condensado y de Petrleo voltil, SPE filial Bolivia,
2000, pag. 24
iv
Phase Behavior Monograph Volume 20 SPE, Curtis H. Whitson and Michael R. Brule,
v
Petroleum Engineering Tool Kit , Programs for Spreadsheet Software,Doug Boone & Joe Clegg,
49
Propiedades del Gas Natural____________________________________________________________ 50
3 Propiedades del Gas Natural
3.1.-INTRODUCCIN
Las leyes que nos describen el comportamiento de los gases en trmino de presin,
volumen y temperatura, que han sido definidas algunos aos atrs, son relativamente simples
para un fluido hipottico conocido como un gas ideal. Este captulo revisa la ley de los gases
perfectos y como pueden ser modificados para describir el comportamiento real de los gases, los
cuales podran desviarse significativamente.
Un gas es definido como un fluido homogneo de baja densidad y viscosidad, la cual no

tiene un volumen definido, pero se expande completamente llenando el recipiente donde esta
contenido. El conocimiento de la presin-volumentemperatura (PVT) esta relacionado con otras
propiedades fsicas y qumicas de los gases en condiciones de reservorio. Es esencial conocer las
propiedades fsicas del gas natural obtenida en laboratorio a condiciones de fondo, para predecir
y conocer el comportamiento de los gases referida a las mezcla.
En el anlisis del comportamiento total de un sistema de produccin es mucho ms

importante el conocimiento de la composicin de las propiedades fsico qumicas del fluido que
ser producido para aplicar el mtodo adecuado en la optimizacin del pozo.
3.2.- Gas Natural
Es una mezcla de hidrocarburos gaseoso, presentes en forma natural en estructuras

subterrneas o trampas. El gas natural esta compuesto principalmente de metano (80%) y
proporciones significativas de etano, propano, butano, pentano y pequeas cantidades de hexano,
heptano y fracciones ms pesadas. Habr siempre alguna cantidad de condensado y/o petrleo
asociado como el gas.
El trmino gas natural tambin es usado para designar el gas tratado que abastece la
industria y a los usuarios comerciales y domsticos, y tienen una cualidad especfica. Esta mezcla
de hidrocarburos gaseosos presentan algunas impurezas, principalmente de: nitrgeno ( N 2 ) ,
dixido de carbono (CO 2 ) y gas sulfhdrico (H 2 S ) .
3.2.1.- Composicin del Gas Natural
Composicin tpica:
Metano (CH 4 ) usualmente > 80%
Etano (C 2 H 6 ) 2 a 10 %
otros hidrocarburos:
Propano (C 3 H 8 ) , Isobutano (i C 4 H 10 ) , Butano normal ( n C 4 H 10 ) , Isopentano
(i C 5 H 12 ) , Pentano normal (n C 5 H 12 ) , Hexano (C 6 H 14 ) , fracciones ms pesadas
(C7 H16+ );
hidrocarburos cclicos y aromticos: ocasionalmente pueden ocurrir en pequeas
proporciones.
Impurezas comunes: nitrgeno (N 2 ) , dixido de carbono (CO2 ) y gas sulfhdrico
(H 2 S )
3.2.2. Comportamiento de los gases ideales
La teora cintica de los gases postula que el gas ideal esta compuesto de una cantidad de
partculas llamadas molculas, cuyo volumen es insignificante comparado con el volumen total
ocupado por el gas. Tambin se asume que estas molculas no tiene una fuerza de atraccin o
repulsin entre ellas y as se asume que todas las colisiones de las molculas son perfectamente
elsticas.
3.2.3 . Ley de los Gases Ideales
Un gas ideal es un fluido en que:

El volumen ocupado por las molculas es pequeo con respecto al volumen ocupado por
el fluido total.
Las colisiones intermoleculares son enteramente elsticas, sin prdida de energa en la

colisin.
No tienen fuerzas atractivas o repulsivas entre las molculas.

La base para describir el comportamiento ideal de los gases viene dada de la combinacin
de las llamadas leyes de los gases ideales.
3.2.3.1. Ley de Boyle
Boylei observ experimentalmente que el volumen de un gas ideal es inversamente

proporcional a la presin absoluta, si la temperatura del gas es constante,.
1
V o pV = cons tan te
p
P1 V2
= o P1V1 = P2V2
P2 V1
Reagrupando:
p1
V2 = V1 Ec..(3.1)
p2
3.2.3.2. Ley de Charles
La ley de Charles expresa en sus dos partes:

1. A presin constante, el volumen variar directamente con la temperatura absoluta,
expresado en la ecuacin,
T
VT o = cons tan te
V
V1 T1 T1 T2
= o =
V2 T2 V1 V2
Reagrupando:
T2
V2 = V1 Ec..(3.2)
T1
2. A volumen constante, la presin absoluta varia con la temperatura expresada en la
ecuacin:
T
pT o = cons tan te
p
P1 T1 T1 T2
= o =
P2 T2 P1 P2
Reagrupando:
T2
p 2 = p1 Ec.(3.3)
T1
3.2.3.3. Ley de Charles y Boyle
Las relaciones de las leyes de Charles y Boyleii pueden ser combinadas para dar la siguiente
relacin:
p1V1 p 2V2
= = cons tan te Ec.(3.4)
T1 T2
3.2.3.4 . Ley de Avogadro
A volmenes iguales todos los gases en las mismas condiciones de presin y temperatura
tienen el mismo nmero de molculas,iii cuyo valor es 2.733 x10 6 moles en una lb mol . La ley
de Avogadro , menciona que el peso de un volumen de gas, es una funcin del peso de las
molculas. El volumen y el peso del gas en libras es igual al valor numrico del peso molecular y
es conocido como el volumen molar. Una libra mol de un gas ideal ocupa 378 pi3 a 60F y
14.73psia, estas condiciones de presin y temperatura son comnmente referidas a las
condiciones normales.
3.2.4 Ecuacin para los gases ideales
La ecuacin de estado para un gas ideal se puede reducir de una combinacin de las leyes
de Boyle, Charles / Gay Lussac y Avogadro.
pV = nRT Ec.(3.5)
Donde
p = Presin absoluta, psia
V = Volumen, ft3
T = Temperatura absoluta, R
n = Nmero de libras-mol, donde 1 lb-mol es el peso molecular del gas
(lb)
R = Constante universal de los gases, para las unidades decimales
tiene el valor de:
(14.7 psia )(379.4 ft 3 )

= 10.73 psia ft 3 lb mol R
(1lb mol )(520 R )
Los valores de la constante de los gases R en diferentes unidades, mostramos en la Tabla 3.1
El nmero de lb-mol de un gas es igual a la masa de gas dividido por el peso molecular del
gas, la ley ideal de gas puede ser expresada como:
m
pV = RT Ec.(3.6)
M
Donde
m = masa de gas, lb
M = peso molecular de gas, lbm lb mol
Tabla 3.1
Valores de la Constante de Gas, R
Unidades R
atm, cc g mole, K 82.06
BTU lb mole, R 1.987
psia, cu ft lb mole, R 10.73
lb sq ft abs, cu ft lb mole, R 1544
atm, cu ft lb mole, R 0.730
mm Hg , liters g mole, K 62.37
in. Hg , cu ft lb mole, R 21.85
cal g mole, K 1.987
kPa, m 3 kg mole, K 8.314
J kg mole, K 8314
3.2.5.- Mezclas de Gases Ideales
El comportamiento de la mezcla de los gases ideales es usualmente interesante, al igual que

los componentes puros. Debido a que la mezcla de componentes de hidrocarburos, y las otras
propiedades fsicas y qumicas pueden ser determinadas en forma individual.
El tratamiento anterior del comportamiento de los gases no slo se aplica a gases de un

nico componente. El ingeniero de gas rara vez trabaja con gases puros, y el comportamiento de
una mezcla de varios componentes se debe tratar como una mezcla. Para esto se requiere la
introduccin de dos leyes adicionales.
3.2.5.1.- Ley de Dalton
En una mezcla gaseosa cada gas ejerce una presin igual que aquella que ejercera si este
ocupa el mismo volumen como una mezcla total. Esta presin es llamada presin parcial. La
presin total es la suma del las presiones parciales. Esta ley es valida solamente cuando la mezcla
y cada componente de la mezcla obedece a la ley de los gases ideales.
La presin parcial ejercida por cada componente de la mezcla de gas puede ser calculado
utilizando la ley de los gases ideales. Considrese una mezcla conteniendo n A moles de
componente A, n B moles de componente B y nC moles de componente C. La presin ejercida
por cada componente de la mezcla gaseosa puede ser determinada con la ecuacin del gas ideal.
RT RT RT
pA = nA , pB = nB , p C = nC
V V V
De acuerdo con la ley de Dalton, la presin es la suma de las presiones parciales.
p = p A + p B + pC
RT RT RT
p = nA + nB + nC
V V V
RT RT
p=
V
ni =
V
n
La relacin de presin parcial del componente i, pi, a la presin total de la mezcla p es:
pi n n
= i = i = yi Ec.(3.7)
p ni n
Donde:
yi = Fraccin molar del componente i
ni = Nmero de moles del componente i, siendo i = 1,2,......,N
ni = Nmero total de moles en la mezcla.

3.2.5.2 Ley de Amagat

En una mezcla gaseosa el volumen total es la suma de los volmenes de cada componente
que ocupara en una presin y temperatura dada. Los volmenes ocupados por los componentes
individuales son conocidos como volmenes parciales. Esta ley es correcta slo si la mezcla y
cada componente obedecen a la ley de los gases ideales.
El volumen parcial ocupado por cada componente de una mezcla de gas de n A moles de
componente A, n B moles de componente B e nC moles de componente C, y as sucesivamente
puede ser calculado utilizando la ley de los gases ideales.
RT RT RT
VA = nA , VB = n B , VC = n C .......
p p p
De acuerdo con la Ley de Amagat, el volumen total es:
V = V A + V B + VC + ...
RT RT RT
V = nA + nB + nc + ...
p p p
RT RT
V=
p
ni =
p
n
La relacin de volumen parcial del componente i al volumen total de la mezcla es:

RT
ni
Vi p n
= = i = yi Ec.(3.8)
V RT n
n
p
Donde
Vi = Volumen ocupado por el componente i a las condiciones normales.
V = Volumen total de la mezcla medido a las condiciones normales,
3.2.5.3.- Fraccin de Volumtrica
La fraccin volumtrica de un componente especifico en una mezcla esta definido

como el volumen de un componente dividido el volumen total de la mezcla.
Vi = vi / Vt
Vi = fraccin volumtrica del componente i en el gas

vi = volumen ocupado por el componente i
Vt = volumen total de la mezcla
Es conveniente en cualquier clculo de ingeniera convertir de fraccin molar a fraccin

de peso o viceversa. El procedimiento de conversin de la composicin de la fase de gas de
la fraccin molar o fraccin de peso esta mostrada en el siguiente procedimiento.
1. Asuma que el nmero total de moles de la fase gasfera es uno n = 1

2. de ecuacin de la fraccin molar se tiene que ni = yi
3. el nmero de moles n de un componente es igual a el peso de los componentes m
dividido por el peso molecular de el componente MW, tambin puede ser
expresado por la ecuacin 3.9, y el peso de los componentes puede ser expresado
como ecuacin. 3.10
n = m/MW Ec.(3.9)
mi = yi* M Wi Ec.(3.10)
4. por lo tanto wi = mi/mt
El procedimiento es ilustrado a travs del siguiente ejemplo Tabla 3.2
Tabla 3.2
Determinacin de la composicin en peso
Componente yi Mwi Mi=yi*Mwi Wi= mi/mt
C1 0,65 16,04 10,426 0,3824
C2 0,10 30,07 3,007 0,1103
C3 0,10 44,10 4,410 0,1618
C4 0,10 58,12 5,812 0,2132
C5 0,05 72,15 3,608 0,1323
Sumatoria 1.00 27,263
3.2.5.4.- Peso Molecular Aparente
Una mezcla gaseosa se comporta como si fuera un gas puro con un peso molecular
definidoiv. Este peso molecular es conocido como un peso molecular aparente y es definido
como:
M a = yi M i Ec.(3.11)
Donde:
Ma = Peso molecular aparente de la mezcla.
yi = Fraccin molar del componente i.
Mi = Peso molecular del componente i.
EL peso molecular de cada componente i puede ser encontrado en la Tabla 3.4
3.2.5.5.-Fraccin Molar
La fraccin molar de un componente en particular es definido como el nmero de

moles de un componente dividido, el nmero de moles totales de todos los componente de
la mezcla.
yi = ni /nt
yi = fraccin molar del componente i en la mezcla
ni = nmero de moles del componente i
nt = nmero de moles total de la mezcla
En la Tabla 3.3 mostraremos el ejemplo 3.2 de la determinacin de la composicin en

fraccin molar
Tabla 3.3
Determinacin de la composicin en fraccin molar
Componente wi Mwi ni=wi/Mwi yi= ni/nt
C1 0,40 16,04 0,025 0,6628
C2 0,10 30,07 0,003 0,0884
C3 0,20 44,10 0,005 0,1205
C4 0,20 58,12 0,003 0,0915
C5 0,10 72,15 0,001 0,0368
Sumatoria 1.00 0,038
3.2.5.-Densidad relativa del Gas
La densidad del gas por definicin es la relacin entre las masas especficas del gas y
del aire, ambas medidas en las mismas condiciones de presin y temperatura, esto es:
g
g = Ec.(3.12)
ar
Admitindose comportamiento de gas ideal, en la ecuacin 3.5, el nmero de moles n
es la relacin entre la masa de gas m y su masa molecular M, ecuacin 3.6. La masa
especfica es definida conociendo la relacin entre la masa y el volumen, o sea:
m pM
g = = Ec.(3.13)
V RT
Y la masa especfica del aire es:

p * 28.97
aire = Ec.(3.14)
RT
Por tanto, la densidad de un gas es:
pM
gas RT M
g = = = Ec..(3.15)
aire p * 28.97 28.97
RT
Donde
g = Densidad del gas.
M = Peso Molecular, lbm lb mol

28.97 = Peso Molecular del aire.
Tabla 3.4 Propiedades de los Hidrocarburos
Freezing Point, F 14.696 psia

Boiling Point F, 14.696 psia
Vapor Pressure, 100F, psia
cu ft/lb
Critical Temperatura F
Critical Pressure,psia
Molecular Weights
Compound
N Formula
Critical Volume
1 Methane CH 4 16.043 -258.73 (5000) -296.44 666.4 -116.67 0.0988
2 Ethane C2 H 6 30.070 -127.49 (800) -297.04 706.5 89.92 0.0783
3 Propane C3 H 8 44.097 -43.75 188.64 -305.73 616.0 206.06 0.0727

4 Isobutane C4 H10 58.123 10.78 72.581 -255.28 527.9 274.46 0.0714
5 n-Butane C4 H10 58.123 31.08 51.706 -217.05 550.6 305.62 0.0703
6 Isopentane C5 H12 72.150 82.12 20.445 -255.82 490.4 369.10 0.0679
7 n-Pentane C5 H12 72.150 96.92 15.574 -201.51 488.6 385.8 0.0675
8 Neopentane C5 H12 72.150 49.10 36.69 2.17 464.0 321.13 0.0673
9 n-Hexane C 6 H 14 86.177 155.72 4.9597 -139.58 436.9 453.6 0.0688
10 2-Methylpentane C 6 H 14 86.177 140.47 6.769 -244.62 436.6 435.83 0.0682
11 3-Methylpentane C 6 H 14 86.177 145.89 6.103 ---------- 453.1 448.4 0.0682
12 Neohexane C 6 H 14 86.177 121.52 9.859 -147.72 446.8 420.13 0.0667
13 2,3-Dimethylbutane C 6 H 14 86.177 136.36 7.406 -199.38 453.5 440.29 0.0665
14 n-Heptane C7 H16 100.204 209.16 1.620 -131.05 396.8 512.7 0.0691

15 2-Methylhexane C7 H16 100.204 194.09 2.272 -180.89 396.5 495.00 0.0673
16 3-Methylhexane C7 H16 100.204 197.33 2.131 ---------- 408.1 503.80 0.0646
17 3-Ethylpentane C7 H16 100.204 200.25 2.013 -181.48 419.3 513.39 0.0665
18 2,2-Dimethylpentane C7 H16 100.204 174.54 3.494 -190.86 402.2 477.23 0.0665

cu ft/lb
N
Molecular Weights
Compound
Formula
Critical Volume
21 Triptane C7 H16 100.204 177.58 3.375 -12.81 428.4 496.44 0.0636
22 n-Octane C8 H18 114.231 258.21 0.5369 -70.18 360.7 564.22 0.0690
23 Diisobutyl C8 H18 114.231 228.39 1.102 -132.11 360.6 530.44 0.0676
24 Isooctane C8 H18 114.231 210.63 1.709 -161.27 372.4 519.46 0.0656
25 n-Nonane C9 H 20 128.258 303.47 0.1795 -64.28 331.8 610.68 0.0684
26 n-Decane C10 H 22 142.285 345.48 0.0608 -21.36 305.2 652.0 0.0679
27 Cyclopentane C5 H10 70.134 120.65 9.915 -136.91 653.8 461.2 0.0594

28 Methylcyclopentane C6 H12 84.161 161.25 4.503 -224.40 548.9 499.35 0.0607
29 Cyclohexane C6 H12 84.161 177.29 3.266 43.77 590.8 536.6 0.0586
30 Methylcyclohexane C7 H14 98.188 213.68 1.609 -195.87 503.5 570.27 0.0600
31 Ethene (Ethylene) C2 H 4 28.054 - -------- -272.47 731.0 48.54 0.0746
154.73
32 Propene (Propylene) C3 H 6 42.081 -53.84 227.7 -301.45 668.6 197.17 0.0689
33 1-Butene (Butylene) C4 H 8 56.108 20.79 62.10 -301.63 583.5 295.48 0.0685
34 Cis-2-Butene C4 H 8 56.108 38.69 45.95 -218.06 612.1 324.37 0.0668
35 Trans-2-Butene C4 H 8 56.108 33.58 49.89 -157.96 587.4 311.86 0.0679
36 Isobutene C4 H 8 56.108 19.59 63.02 -220.65 580.2 292.55 0.0682
37 1-Pentene C 5 H 10 70.134 85.93 19.12 -265.39 511.8 376.93 0.0676
38 1,2-Butadiene C4 H 6 54.092 51.53 36.53 -213.16 (653) (340) (0.065

39 1,3-Butadiene C4 H 6 54.092 24.06 59.46 -164.02 627.5 305 0.0654
40 Isoprene C5 H 8 68.119 93.31 16.68 -230.73 (558) (412) (0.065
41 Acetylene C2 H 2 26.038 - -------- -114.5 890.4 95.34 0.0695
120.49
42 Benzene C6 H 6 78.114 176.18 3.225 41.95 710.4 552.22 0.0531
43 Toluene C7 H 8 92.141 231.13 1.033 -139.00 595.5 605.57 0.0550
44 Ethylbenzene C8 H10 106.167 277.16 0.3716 -138.96 523.0 651.29 0.0565
45 o-Xylene C8 H10 106.167 291.97 0.2643 -13.59 541.6 674.92 0.0557

cu ft/lb
Molecular Weights
Compound
Formula
Critical Volume
N
46 m-Xylene C8 H10 106.167 282.41 0.3265 -54.18 512.9 651.02 0.0567

47 p-Xylene C8 H10 106.167 281.07 0.3424 55.83 509.2 649.54 0.0570
48 Styrene C8 H 8 104.152 293.25 0.2582 -23.10 587.8 703 0.0534
49 Isopropylbenzene C9 H12 120.194 306.34 0.1884 -140.81 465.4 676.3 0.0572
50 Methyl alcohol C H4 O 32.042 148.44 4.629 -143.79 1174 463.08 0.0590
51 Ethyl alcohol C 2 H 6O 46.069 172.90 2.312 -173.4 890.1 465.39 0.0581
52 Carbon monoxide CO 28.010 -312.68 -------- -337.00 507.5 -220.43 0.0532
53 Carbon dioxide CO 2 44.010 -109.257 -------- -69.83 1071 87.91 0.0344
54 Hydrogen sulfide H2S 34.08 -76.497 394.59 -121.88 1300 212.45 0.0461
55 Sulfur dioxide SO 2 64.06 14.11 85.46 -103.86 1143 315.8 0.0305
56 Ammonia NH 3 17.0305 -27.99 211.9 -107.88 1646 270.2 0.0681
57 Air N 2 + O2 28.9625 -317.8 -------- --------- 546.9 -221.31 0.0517

58 Hydrogen H2 2.0159 -422.955 -------- -435.26 188.1 -399.9 0.5165
59 Oxygen O2 31.9988 -297.332 -------- -361.82 731.4 -181.43 0.0367
60 Nitrogen N2 28.0134 -320.451 -------- -346.00 493.1 -232.51 0.0510
61 Chlorine Cl 2 70.906 -29.13 157.3 -149.73 1157 290.75 0.0280
62 Water 18.0153 212.000 0.9501 32.00 3198.8 705.16 0.0497

H 2O
63 Helium He 4.0026 -452.09 -------- -------- 32.99 -450.31 0.2300
64 Hydrogen H Cl 36.461 -121.27 906.71 -173.52 1205 124.77 0.0356

Chloride
3.3.1.- Gases Reales
Bsicamente la magnitud de desviacin de los gases reales con respecto a los gases
ideales incrementa cuando incrementamos la presin y temperatura, variando tambin con
la composicin del gas. El comportamiento de un gas real es diferente a un gas ideal, la
razn para esto es que la ley de los gases perfectos fue derivada bajo la asumisin que el
volumen de molculas es insignificante y no existe atraccin o repulsin entre las
molculas, lo cual no es el caso para gases reales.
En la prctica los gases no se comportan de acuerdo con la ley definida por la

ecuacin 3.5 para las presiones y temperaturas de trabajo. Para expresar de forma ms real
la relacin entre las variables P, V e T, un factor de correccin, denominado factor de
compresibilidad de gas Z, es introducido en la ecuacin 3.5:
pV = ZnRT Ec..(3.16)
Donde para un gas ideal, Z = 1.
El factor de compresibilidad vara con el cambio de presin, volumen, temperatura y

la composicin del gas. Esto debe determinarse experimentalmente. Los resultados de la
determinacin experimental del factor de compresibilidad son normalmente dados
grficamente y normalmente toman la forma en la Figura 3.1
3.3.2.- Mtodo de obtencin del factor de compresibilidad Z

Con la aparicin del teorema de los estados correspondiente, desarrollado por Van
der Waals (1873) posibilito la elaboracin de bacos universales para la obtencin del
factor de compresibilidad Z. Siguiendo este teorema, los gases exhiben el mismo
comportamiento cuando lo sometemos a las mismas condiciones de presin, temperatura y
volumen reducidos. El trmino reducido traduce la razn entre la variable y/o su valor
crtico:
p
pr = Ec.(3.17)
pc
T
Tr = Ec.(3.18)
Tc
V
Vr = Ec.(3.19)
Vc
Donde:
pr = Presin reducida.
Tr = Temperatura reducida.
Vr = Volumen reducido.
pc = Presin crtica.
Tc = Temperatura crtica.
Vc = Volumen crtico.
Figura 3.1 Grfico tpico del factor de compresibilidad como una funcin de la
presin y temperatura constante.
El teorema de los estados correspondientes no es perfecto, mas aun cuando es

aplicado a gases con estructuras qumicas similares (por ejemplo: hidrocarburos
parafnicos) ofrece un mtodo de correlacin con precisin satisfactoria para trabajos de
ingeniera.
3.3.3. Correlaciones de Standing y Katzv
La correlacin de Standing y Katz, favorece valores de Z en funcin de presiones y

temperaturas reducidas, y fue desarrollada con base en datos experimentales para gases
naturales sin impurezas. Su aplicacin para gases cidos requiere el uso de factores de
correccin para la presencia de CO 2 y H 2 S .
El procedimiento para la determinacin de Z sigue los siguientes pasos:

Paso 1 - Determine las propiedades pseudo crticas
a. Composicin conocida.
De Tabla 3.4 determinar el Peso Molecular, presin y temperatura

pseudo crticas para cada componente.
Obtener la Masa molecular aparente (M a ) de la mezcla, ecuacin 3.11
Obtener las coordenadas pseudo crticas. Estas pueden ser calculadas a

travs de media ponderadas de las coordenadas crticas de cada componente y su fraccin
molar en la mezcla:
nc
PPc = y i p ci Ec.(3.20)
i =1
nc
TPc = y i Tci Ec.(3.21)
i =1
Donde
Pci = Presin Pseudo crtica del componente i.

Tci = Temperatura Pseudo crtica del componente i.
yi = Fraccin molar del componente i.
nc = Nmero de componentes.
b. Composicin desconocida.
Con la densidad del gas conocida usar la Figura 3.2, donde la presin y
temperatura pseudo crticas son dadasvi. O a travs de las siguientes correlaciones
presentadas por Standing (1981):
Ppc = 677 15,0 g 37,5 g2 Ec.(3.22)
T pc = 168 + 325 g 12,5 g2 Ec.(3.23)

Estas ecuaciones estn limitadas por el contenido de impurezas presentadas en la
mezcla gaseosa, los mximos porcentajes son de 3% H 2 S y 5% N 2 , o un contenido total
de impurezas de 7%.
Para salvar este obstculo las propiedades obtenidas pueden ser corregidas. Una
alternativa es el uso de la Figura 3.2 para el clculo de las propiedades pseudo crticas de
una mezcla gaseosa de hidrocarburos cuando hay la presencia de contaminantes, y el uso
de la Figura 3.3.
Figura 3.2 Correlaciones para las Coordenadas Pseudo Crticas
(Tomado de Brown, Katz, Oberfell y Alden Natural Gasoline Association of America )

Figura 3.3 Propiedades seudo crticas del gas Natural
( Tomado de Standing y Katz AIME)

Para gas natural seco, es:

Ppc = 706 51,7 g 11,1 g2 Ec.(3.24)
T pc = 187 + 330 g 71,5 g2 Ec.(3.25)
Paso 2 - Correcciones de las propiedades pseudo crticas.
Debido a la presencia de gases que no son hidrocarburos, utilizamos el

factor de correccin de Wichert y Aziz dado por:
( ) (
= 120 A0,9 A1,6 + 15 B 0,5 B 4,0 ) Ec.(3.26)
Donde
= Factor de ajuste de las propiedades pseudo crticas
A = Suma de las fracciones molares de H 2 S + CO2 .
B = Fraccin molar de H 2 S .
Calcule la temperatura pseudo crtica ajustada por la siguiente expresin:
T Pc = TPc Ec.(3.27)
Y a la presin pseudo crtica ajustada por medio de:
PPcT Pc
P Pc = Ec.(3.28)
TPc + B(1 B )
El parmetro puede de ser tambin obtenido en la Figura 3.4
Paso 3 - Factor de compresibilidad Z
El factor de compresibilidad Z, de la mezcla gaseosa es obtenido de la carta de

Standing & Katz, Figura 3.5 en funcin de las coordenadas pseudo reducidas, o sea:
p
pPr = Ec.(3.29)
pPc
T
T pr = Ec.(3.30)
TP c
O, si hay presencia de contaminantes:
p
pPr = Ec.(3.31)
p P c
T
T pr = Ec.(3.32)
T P c
3.3.4. Correlaciones de Brill & Beggs
Una modificacin de las ecuaciones publicadas por Brill & Beggs1 (1974) favorece
valores del factor Z con precisin suficiente para la mayora de los clculos de
ingeniera:
1 A D
Z = A+ + C p Pr Ec.(3.33)
exp B
Donde:
A = 1.39(T pr 0.92) 0.5 0.36T pr 0.101 Ec.(3.34)
0.066
( )
B = 0.62 0.23T pr p pr +
T pr 0.86
0.037 p 2pr +

0.32
p6
9(T pr 1) pr
Ec.(3.35)
10
C = 0.132 0.32 log T pr Ec.(3.36)
( 2
D = anti log 0.3106 0.49T pr + 0.182 4T pr ) Ec..(3.37)
3.3.5. Correlaciones de Dranchuk, Purvis y Robinson2
Este mtodo se basa en la ecuacin de estado desarrollada por Benedicto, Webb

and Rubin, para representar el comportamiento de hidrocarburos leves. Utilizando esa
ecuacin, conteniendo ocho constantes caractersticas de cada sustancia, , Dranchuk et
a,. obtuvieron una ecuacin explcita de la siguiente forma:
A2 A3 A5 2 r5

Z = 1 + A1 + + + A4 + r + A5 A6 +
T 3 r T T
r Tr r r
Ec.(3.38)
(1 + A )EXP( A )
2
A7 r 2 2
3 8 r 8 r
Tr
Figura. 3.4 Factor de ajuste de temperatura pseudo crticas

Donde
Pr
r = 0,27 Ec..(3.39)
ZTr
y las constantes de la correlacin son dadas por:
A1 = 0,31506237 A5 =-0,61232032
A2 = -1,04670990 A6 =-0,10488813
A3 = -0,57832729 A7 = 0,68157001
A4 = 0,53530771 A8 = 0,68446549
Figura. 3.5 Factor de compresibilidad para gases naturales
( Tomado de Standing y Katz AIME)

3.3.6. Correlaciones de Hall-Yarboroughvii
La ecuacin de Hall-Yarborough, desarrollada usando la ecuacin de estado de

Starling-Carnahan:
2
0,06125 PPr t e 1,2(1t )
Z= Ec..(3.40)
Y
Donde
PPr = Presin pseudo reducida.
t = Recproco, da temperatura reducida (TPc T )
Y = Densidad reducida con la cual puede ser obtenida como la
solucin de la ecuacin.
Y +Y 2 +Y 3 Y 4
F K = 0,06125PPr t e 1,2(1t ) +
2
( )
14,76 t 9,76t 2 + 4,58t 3 Y 2
(1 Y )3 Ec. 3.41
( )
+ 90,7t 242,2t 2 + 42,4t 3 Y (2,18+ 2,82t ) = 0 eq.(3.41)
Esta ecuacin no lineal puede ser convenientemente y resuelta usando las tcnicas
simples de iteracin de Newton-Raphson. Los pasos son:
K
1. Tomar una estimacin inicial de Y , donde K es un contador de iteracin.
1
Y = 0,001
2. Substituya este valor en la ecuacin 3.38; a menos que el valor correcto de Y tenga
que ser inicialmente seleccionado, en la ecuacin 3.38 sera un valor muy pequeo, el valor
K
de F diferente de cero.
3. Utilizando la primera serie de expansin de Taylor , una mejor estimacin de Y

puede ser determinada por la ecuacin:
K +1 K FK
Y =Y Ec.(3.42)
dF K dY
Donde la expresin general para dF dK puede ser obtenido como la derivada de

la ecuacin 3.38.
dF 1 + 4Y + 4Y 2 4Y 3 + Y 4
dY
=
(1 Y )4
(
29,52t 19,52t 2 + 9,16t 3 Y
Ec.(3.43)
)
(
+ (2,18 + 2,82t ) 90,7t 242,2t 2 + 42,4t 3 Y (1,18 + 2,82t ) )
4. Iterando, usando la ecuacin 3.38 y 3.40 antes que haya convergencia dentro de
K
una aproximacin satisfactoria, F 0.
5. Substituya este valor correcto de Y en la ecuacin 3.40, para determinar el valor
de Z.
3.4.1. Factor Volumtrico del Gas Naturalviii
Relaciona el volumen de gas evaluado a condiciones de reservorio y el volumen del

mismo gas medido a condiciones superficiales de, p sc y Tsc . Generalmente, se expresa en
pies cbicos o barriles de volumen en el reservorio por pi cbico de gas las condiciones
normales, o sus recprocos, en pies cbicos a las condiciones normales por pi cbico o
barril de volumen en el reservorio.
V p ,T
Bg = Ec.(3.44)
V sc
De acuerdo con la ecuacin de estado de los gases reales, la ecuacin 3.16, el volumen
de una determinada masa de gas m, equivalente a un nmero de moles n es dada en
condiciones de reservorio por:
ZnRT
V= Ec.(3.45)
p
Donde T y P son la temperatura y la presin del reservorio. En las condiciones normales

esta misma masa de gas ocupa el volumen:
Z sc nRTsc
V= Ec.(3.46)
p sc
El factor volumtrico de la formacin en la ecuacin 3.44 es:
ZnRT
p
Bg =
Z sc nRTsc
p sc
ZTp sc vol
Bg = Ec.(3.47)
Z sc Tsc p std vol
Utilizando, Tsc = 520 R , p sc = 14.7 psia y Z sc = 1, la ecuacin 3.47 se
hace:
ZT (14.7) ZT 3
Bg = = 0.0283 ft scf Ec.(3.48)
1(520) p
ZT
B g = 0.00504 bbls scf Ec.(3.49)
p
p
B g = 35.35 scf ft 3 Ec.(3.50)
ZT
p
B g = 198.4 scf bbls Ec.(3.51)
ZT
3.5.1. Compresibilidad Isotrmica de Gas Naturalix
La compresibilidad isotrmica del gas natural, es definida como el cambio de

volumen por unidad de volumen para una unidad de cambio en presin y temperatura
constante,
1 V
C= Ec.(3.52)
V p T
La compresibilidad es requerida en muchas ecuaciones del reservorio de flujo de

gas y puede ser evaluado de la siguiente manera:
3.5.1.1. Compresibilidad para un gas ideal
nRT
V=
p
y
V nRT
=
p T p2
por tanto
p nRT
Cg =
nRT p 2
1
Cg = Ec.(3.53)
p
3.5.1.2. Compresibilidad para un gas real
nRTZ
V= y
p
V 1 Z Z
= nRT Por tanto
p T p p p 2

p 1 Z Z
Cg = nRT
2
o
nRTZ p p p
1 1 Z
Cg = Ec.(3.54)
p Z p
La evaluacin de C g para gases reales, requiere determinar cmo el factor Z varia con la
presin y temperatura. La mayora de las grficas y ecuaciones que determinan Z estn
expresadas como una funcin de presin y temperatura reducida. La compresibilidad
reducida viene a ser definida como: C r = C g p c . Este puede ser expresado como una
funcin de p r en un valor fijo de Tr por:
1 1 Z
Cr = Ec..(3.55)
pr Z pr Tr
Los valores de la derivada ( Z p r )T puede tener una pendiente constante de

r
Tr ver (Figura 3.5) para el factor Z. Los valores de C r Tr como una funcin de p r y
Tr vienen a ser presentadas grficamente por Mattar, et al. En la Figuras 3.6 y 3.7. El
cambio de Z con p puede tambin ser calculado utilizando una expresin analtica,
calculando el factor Z a las presiones ligeramente encima y abajo de la presin de
inters.
Z Z Z2
= 1 Ec..(3.56)
pr p
Tr r1 p r 2 Tr
3.6.1. Viscosidad del Gas Naturalx
Es una medida de resistencia al movimiento de fluido y esta definida como la relacin

entre el esfuerzo cortante por unidad de rea y el gradiente de velocidad en un punto
determinado.
F

=
A
Ec..(3.57)
dV

dL
La viscosidad normalmente es expresada en centipoises o poises, pero se puede convertir

para otras unidades:
2
1 poise = 100 centipoise = 6 . 72 x 10 lbm ft sec
3 2
= 2 . 09 x 10 lbf sec ft = 0 . 1 kg m sec
La viscosidad absoluta de un gas natural vara con la composicin y con las

condiciones de presin y temperatura:
= f ( p, T , Composici n)
3.6.1.1. Determinacin de la viscosidad: Mtodo de Carr, Kobayashi y Burrows
Es la relacin ms ampliamente utilizada en la industria del petrleo.
Paso 1 - Determine la viscosidad del gas a la presin atmosfrica utilizando la

Figura 3.8
1 = f (M , T )
El peso molecular M, lleva en consideracin la influencia de la composicin del

gas en la determinacin de la viscosidad 1 .
Paso 2 - Haga correlaciones para la presencia de N 2 , CO2 y H 2 S .
( 1 )cor = 1 + correlacin H 2 S + correlacinCO 2 + correlacin N 2
Ec.(3.58)
Paso 3 - Calcule la presin y temperatura pseudo reducidas: p r y Tr .
Paso 4 - Obtenga la relacin 1 en funcin de p r y Tr , utilizando los
grficos de las Figuras 3.9 y 3.10.
Paso 5 - Determine la viscosidad del gas por medio de:

= (1 )corr Ec..(3.59)
1
3.6.1.2. Determinacin de la viscosidad: Mtodo de Lee, Gonzlez y Eakin xi
La viscosidad del gas puede ser obtenida por medio de:
= 10 4 K EXP X Y ( ) eq.(3.60)
Donde
K=
(9,4 + 0,02 M )T 1,5
; X = 3,5 +
986
+ 0,01M ; Y = 2,4 0,2 X
209 + 19M + T T
En estas ecuaciones las unidades utilizadas son: T = R , g = cp ,

M = peso molecular , g = gr cm 3 .
Figura 3.6 Variacin de CrTr con la presin y temperatura reducida. Extrada de

Engineering Data Book Gas Processors Suppliers Association, 1987. (1,05 < = Tr
< = 1,4 ; 0,2 < = Pr < = 15,0)
Figura 3.7 Variacin de CrTr con la presin y temperatura reducida. Extrada de

Engineering Data Book Gas Processors Suppliers Association, 1987. (1,4 < = Tr <
= 3,0 ; 0,2 < = Pr < = 15,0)
Figura 3.8 Viscosidad del gas natural a la presin de 1 atm.

. (Tomado de Carr, Kobayashi y Burrows, AIME)
Figura 3.9 Razn de viscosidad del gas natural. (Tomado de Carr, Kobayashi y
Burrows, AIME)
. (Tomado de Carr, Kobayashi y Burrows, AIME)

Figura 3.10 Razn de viscosidad del gas natural
3.7.1.- Factor de Compresibilidad para un sistema bifsico para gases retrgrado

Las correlaciones de Standing y Katz para el factor de compresibilidad son validas
solamente para un sistema monofsico o gas seco. En los reservorios de Gas-Condensado
retrgrado existe una condensacin del fluido durante el agotamiento del reservorio por
debajo del punto de Roci, obtenindose un sistema bifsico con una fase lquida. En este
caso se deber utilizar el factor de compresibilidad para dos fases. El factor de
compresibilidad es usado en las ecuaciones de balance de materiales, para estimar el
volumen inicial In-Situ, y las reservas recuperables en los programas del libro.
El factor de compresibilidad del gas es normalmente usado cuando no se tiene

disponible los estudios de fluido del reservorio, esta practica es aceptable para un sistema
de gas condensado con condensacin retrgrada, sin embargo si el gas es rico la reserva
podra ser sobre estimada, si no se utiliza el factor de compresibilidad para las dos fases.
El comportamiento de la relacin del factor de compresibilidad de un gas condensado

frente a un gas seco como funcin de la presin, es uniforme hasta la presin de roci,
luego el factor de compresibilidad para el gas seco cae por debajo del factor de
compresibilidad de un reservorio de gas condensado.
En general, las correlaciones para la obtencin del factor de dos fases estn basadas
en las propiedades pseudo reducidas del gas en el reservorio. Reyes et Alii, presentaron una
correlacin emprica para determinar el factor de compresibilidad de dos fases cuando el
gas es rico (en porcentaje en mol del C7+ mayor e igual al 4 % y su densidad del gas
mayor a 0.911) relativamente puro (porcentaje de impureza H2S o CO2 menor o igual al 5
% en mol) .
r
Z 2 f = Ao + A1* Pr + A2 / Tr + A3 * Pr 2 + A4 / Tr 2 + A5 *
Tr Ec. 3.61
Vlida para los siguientes rangos (0.7< =Pr < =20) y (1.1 < =Tr < =2.1)
A0 = 2.24353 A3 = 0.000829231
A1= -0.0375281 A4 = 1.53428
A2 = -3.56539 A5 = 0.131987
La suma del coeficiente 0.308 en la ecuacin 3.61, puede ser interpretada como el
factor de compresibilidad crtica. Los coeficientes fueron obtenidos utilizando para un
ajuste de 67 muestras con fraccin de C7+. superior al 4 %, en porcentaje de impurezas
de H2S o CO2 menor al 5 % en mol. con un total de 478 puntos de anlisis que
presentan un error absoluto medio de 3.38 %. En este mismo porcentaje de impurezas,
mayor al 5 %, la correlacin muestra ser vlida. Por lo tanto para determinar el factor de
compresibilidad de dos fases debemos seguir el siguiente procedimiento:
1. Calcular las coordenadas seudo reducidas del fluido producido a partir de

su composicin o su densidad.
2. Para concentracin de C7+ superior al 4 % o su densidad del fluido

producido es mayor a 0.911 se utiliza la ecuacin para dos fases
3. Si su composicin y densidades son desconocidas, se sugiere que se utilice
la ultima composicin y densidad conocida
3.7.1.1. Coordenadas Seudo crticas del C7+
Las coordenadas crticas (Pc , Tc) de los componentes puros del gas natural son
fcilmente encontrados en la literatura para la fraccin del C7+, en tanto esas
propiedades deben ser estimada con las ecuaciones de Kessler and Lee, las cuales
pueden ser usadas para estimar las coordenadas seudo crticas de la fraccin C7+, la
expresin es la siguiente :
Ppc..(psia)..= exp.(Ao+ A1/ d (A2 + A3/ d + A4/ d2) *103 *Tb

+ (A5 + A6/ d + A7/ d2) *107 *Tb2 (A8 + A9/ d2) *1010 *Tb3) Ec. 3.62
Donde d es la densidad de la fraccin C7+, en relacin al agua Tb es una temperatura

normal de ebullicin (o R ) de la fraccin C7+ , los coeficientes de la ecuacin son:
A0= 8.3634 A1= 0.0566

A2 = 0.24244 A3 = 2.2898
A4 = 0.11857 A5 = 1.4685
A6 = 3.648 A7 = 0.47227
A8 = 0.42019 A9 = 1.6977
La temperatura seudo crtica es estimada de la siguiente ecuacin :
Tpc..(.R.)..= 341.7 + 811* d + (0.4244+ 0.1174* d) *Tb+ (0.4669 3.2623* d) *105 /Tb
Ec. 3.63
Las coordenadas seudo crticas estn calculadas a la temperatura normal de ebullicin.

Para determinar la temperatura normal de ebullicin se debe utilizar la ecuacin
propuesta por Whitson, donde M es la masa molecular del C7+
Tb..(.R.)..= (4.55579* M 0.15178* d 0.1542)3

Ec. 3.64
3.7.1.2. Coordenadas seudo crtica de la mezcla
Una manera simple de presentar la temperatura y la presin seudo crtica de la

mezcla, esta dada por las siguientes ecuaciones de Sutton :
Ppc (psi) = 756.8 131.0*SGg 3.6*SGg2 Ec. 3.65

Tpc ( R )= 169.2 + 349.5* SGg 74.0*SGg2 Ec. 3.66
Ejemplo No 3 Determine el factor de compresibilidad para un fluido bibsico
Pr= 4987 psi , Procio = 3155 psi , SGg = 0.67, API =55, Sgmezcla = 0.828 Tr=
Ppc (psi) = 756.8 131.0*SGg 3.6*SGg^2 = 645.51
Tpc ( R )= 169.2 + 349.5* SGg 74.0*SGg^2 = 408.60
Z 2fases = 0.6145
Gas Production Operations H. Dale Beggs, Second printing, Nov.1985.

Ingeniera Aplicada de Yacimientos Petrolferos B. C. Craff y M. F.
Hawkins, Jr, 1997
Natural Gas Production Engineering Chi U. Ikoku
Engineering Data Book Gas Processors Suppliers Association, Volume 2,
Tenth Edition, 1987
i
Natural Gas Production Engineering Chi U. Ikoku, pag,
ii
iii
iv
Gas Production Operations H. Dale Beggs, 1985, pag. 17
v
Engenharia do Gs Natural Oswaldo A. Pedrosa Jr., 2001 by Pontifcia Universidade Catlica, 2002
vi
Previso de Comportamento de Reservatrios de Petrleo, Adalberto Jos Rosa, Renato de Souza
Carvalho, pag. 296
vii
Engenharia do Gs Natural Oswaldo A. Pedrosa Jr., 2001 by Pontifcia Universidade Catlica, 2002
viii
Ingeniera Aplicada de Yacimientos Petrolferos B. C. Craff y M. F. Hawkins, Jr, 1997, pag.
ix
x
Gas Production Operations H. Dale Beggs, 1985,
xi
ANLISIS DEL RESERVORIO__________________________________________________________________87
4 Anlisis Del Reservorio
4.1 Introduccin
El estudio del comportamiento del reservorio es muy importante para optimizar la

capacidad de produccin. El anlisis de las caractersticas y los factores que afectan al flujo
de fluido a travs del reservorio, y el sistema de tubera, nos lleva a optimizar e incrementar
la capacidad de produccin, siendo esta la base para la seleccin de mtodos de prediccin
del comportamiento de flujo en todo el sistema. (Analizando como una sola unidad).
Los reservorios pueden ser petrolferos y gasferos, pero nos abocaremos a los que
son de inters para nuestro anlisis de acuerdo a su composicin y relacin gas-petrleo.
Sabemos que al viajar el fluido desde el reservorio hacia la caera de produccin existen
prdidas de presin, debido a la resistencia al flujo que ejercen la roca y las tuberas de
produccin. Estas prdidas de presin dependen principalmente del caudal de flujo,
propiedades del fluido, propiedades de la roca y los factores de friccin.
El ingeniero de optimizacin en la produccin de gas debe ser capaz de prever no

slo el caudal de un pozo o un campo productor, si no tambin debe tener muy definido el
concepto de reservorio, la reserva original In-Situ, reserva recuperable y el caudal
econmico de produccin, relacionando las reservas remanentes con la presin de
reservorio.
La Figura 4.1 nos muestra un esquema de caudal versus presin fluyente en el fondo
de pozo, llamada relacin del comportamiento de flujo de entrada (IPR inflow performance
relationship) la cual nos permite visualizar el caudal de produccin versus la presin de
flujo. La curva A nos muestra el comportamiento de un ndice de productividad constante,
debido a que la presin fluyente se encuentra por encima del punto de roco en un sistema
monofsico. En la curva B nos muestra un sistema combinado; primeramente, observamos
un sistema monofsico para luego tener un sistema bifsico con el ndice de productividad
variable, ya que la presin fluyente se encuentra por debajo de la presin de roco. La curva
C nos muestra un comportamiento de un sistema bifsico con un ndice de productividad
variable, debido a que la presin de reservorio se encuentra por debajo de la presin de
roco.
Para calcular la cada de presin que ocurre en un reservorio, es necesario tener una
ecuacin que represente este comportamiento y exprese las prdidas de energa o prdidas
de presin debido a la fuerzas de friccin que es una funcin de velocidad o rgimen de
flujo. La forma de la ecuacin puede ser bastante diferente para los varios tipos de fluido,
las ecuaciones bsicas en todas las formas estn basadas en la ley de Darcy.
Figura 4.1 Curvas IPR Tpicas

5.2 Ley de Darcy
4.2 Ley de Darcy
Esta es simplemente una relacin emprica que se derivo para el flujo de fluido a
travs del filtro de arena no consolidada. Darcy, propuso una ecuacin que relaciona la
velocidad aparente del fluido con el gradiente de presin dp/dx, la cual es valida para flujo
vertical, horizontal e inclinada y tambin demostr que la velocidad del fluido es
inversamente proporcional a la viscosidad, ( ).
Se bebe tomar en cuenta que los experimentos de Darcy, fueron hechos tomando el
agua como fluido base. El filtro de arena fue saturado completamente con agua.. Ya que los
filtros de arena de Darcy son de rea constante, la ecuacin no calcula los cambios de la
velocidad con respecto a la posicin, siendo escrita la Ley de Darcy en forma diferencial de
la siguiente manera:
1 p'
v = k Ec.(4.1)
x'
El signo negativo se agrega porque si x se mide en la direccin del flujo, la presin

p declina en la misma direccin (gradiente de presin negativo), de esto resulta que el
signo menos debe agregarse para hacer la velocidad v positiva. Si sustituimos la

velocidad aparente v la expresin Q= v` * A, tenemos:
kA p'
Q = Ec. (4.2)
x'
Donde
Q= el caudal en cc/seg.
A = rea en cm2.
p ' = Gradiente de presin en atmsfera por centmetro.
x '
= Viscosidad en centipoises
Las unidades de la constante resultante, k, son diferente dependiendo de las

unidades usadas. La ley es vlida para un sistema homogneo de flujo laminar a valores
bajos de nmero de Reynolds
4.2.1- Flujo Lineal
Para el flujo lineal, el rea de flujo es constante, debiendo integrar la ecuacin de

Darcy para obtener la cada de presin que ocurre en una longitud L dada:
p2
kdp qL

= dx
Ao
Ec.(4.3)
p1
Si se supone que k , , y q son independientes de la presin o que pueden ser

evaluados con una presin promedio del sistema, la ecuacin viene a ser:
p2
q L
dp = dx
kA 0
Ec.(4.4)
p1
Integrando la ecuacin da:

qu
p 2 p1 = L Ec.(4.5)
kA
o
CkA( p1 p 2 )
q= Ec.(4.6)
L
Donde C es un factor de conversin de unidades. El valor correcto para C es 1.0

3
para las unidades Darcy y 1.127 x10 para las unidades de campo.
TABLA 4.1
Unidades de ley de Darcy
Variable Smbolo Unidad de Unidad de

Darcy Campo
Caudal de flujo q cc seg bbl dia
Permeabilidad k darcys md
rea A cm 2 ft 2
Presin p atm psi
Viscosidad cp cp
Longitud L cm Pies
La geometra del sistema lineal es ilustrada en la Figura 4.2
Figura 4.2 Geometra para flujo lineal
Se puede observar la ecuacin 4.5, en un esquema de coordenadas cartesianas

de p vs L que producir una lnea recta de pendiente constante, q kA . Donde la
variacin de la presin con la distancia es lineal.
Si el flujo de fluido es compresible, el caudal de flujo de masa q debe ser

constante y es expresada en trminos de presin, temperatura y gravedad especfica de gas,
entonces la ecuacin ser:
8.93ZTL
p12 p 22 = q sc Ec.(4.7)
kA
Donde:
p = psia
o
T = R
= cp
L = ft
k = md
A = ft 2
q sc = scf dia
Para flujo de altas velocidades en la cual existe turbulencia la ley de Darcy, debe
modificarse para calcular la cada de presin causada por la turbulencia. Aplicando la
correccin de turbulencia en la ecuacin para flujo de gas, esta viene a ser:
8.93Z g LT 1.247 x10 10 ZTL g

p12 p 22 = q sc + 2
2
q sc Ec.(4.8)
kg A A
Donde:
Z = Factor de compresibilidad del gas, obtenido a partir T , p .

o
T = Temperatura de flujo, R.
g = Gravedad del gas.
q sc = Caudal de flujo de gas, a 14,7 psia, 60 F, scf dia .
g = Viscosidad de gas, a T , p, cp .
kg = Permeabilidad del gas, md .
2
A = rea de flujo, ft .
Se puede obtener una aproximacin al coeficiente de velocidad a travs de:
2.33 x1010
= Ec.(4.9)
k 1.2
Donde:
= ft 1
k = md
4.2.2. Flujo Radial
Aunque el flujo lineal raramente ocurre en un reservorio, nosotros usaremos estas

ecuaciones para calcular la cada de presin a travs de la formacin, siendo esta:
= Pwfs Pwf Ec. (4.10 )
Para flujo radial, tambin se puede usar la Ley de Darcy para calcular el flujo dentro
del pozo donde el fluido converge radialmente a un cilindro relativamente pequeo. En este
caso, el rea abierta al flujo no es constante, por tanto deber incluir en la integracin de la
ecuacin 4.2, la geometra de flujo de la Figura 4.3, en la que se puede ver que la seleccin
de rea abierta al flujo en cualquier radio es:
A = 2 r h Ec.(4.11)
Definiendo el cambio en la presin con la ubicacin como negativa con respecto a la

direccin de flujo, dp dx se vuelve dp dr . Haciendo estas substituciones en la
ecuacin 4.2 da:
Figura 4.3 Flujo Radial
k (2 r h) dp
q= Ec..(4.12)
dr
Donde:
r = Distancia radial.
h = Espesor del reservorio.
Para un flujo de gas, antes de la integracin de la ecuacin 4.12 ser combinada con
la ecuacin de estado y la ecuacin de la continuidad.
4.2.3.- Flujo de Gas
El flujo de gas para un flujo radial esta basado en la ley de Darcy, la cual considera que
el fluido es compresible y esta basado en la ecuacin de estado real de un gas, donde el gas
es medido bajo condiciones estndar de superficie. La ecuacin para un fluido monofsico
la definiremos de la siguiente forma:
La ecuacin de la continuidad es:
1 q1 = 2 q 2 =cons tan te Ec.(4.13)
La ecuacin de estado para un gas real es:
pM
= Ec.(4.14)
ZRT
El rgimen de flujo para un gas es normalmente dado en algunas condiciones

Standard de presin y temperatura, p sc y Tsc , usando estas condiciones en las
ecuaciones 4.13 y combinando en las ecuaciones 4.13 y 4.14.
q = sc q sc
o:
pM p sc M
q = q sc
ZRT Z sc RTsc
Resolviendo para q sc y expresando q con la ecuacin 4.12 muestra:
pTsc 2 rhk dp
q sc =
p sc ZT dr
Las variables en esta ecuacin son p e r . Separando las variables e integrando:

pR re
q p T Z dr
pdp = sc sc
pwf Tsc 2 kh rw r
(p R2 p wf2 ) = q sc p scTZ ln(r rw )

e
2 Tsc 2 kh
q sc =
(
khTsc p R2 p wf
2
) Ec.(4.15)
p sc T Z ln (re rw )
La ecuacin 4.15 es aplicable para cualquier grupo consistente de unidades. En

las unidades llamadas convencionales, de campo la ecuacin vendr a ser:
q sc =
(
703 x10 6 kh p R2 p wf
2
) Ec.(4.16)
T Z ln (re rw )
La ecuacin 4.16 incorpora los siguientes valores de presin y temperatura estndar,

p sc = 14.7 psia y Tsc = 520 R . Modificando esta ecuacin para flujo
estabilizado con presin media del reservorio:
q sc =
(
703 x10 6 kh p R2 p wf
2
) Ec.(4.17)
T Z [ln (re rw ) 0.75 + S ]
Donde:
q sc = Caudal de flujo de gas, Mscfd
k = Permeabilidad, md
h = Espesor del reservorio, ft
pR = Presin media del reservorio, psia
p wf = Presin fluyente en el fondo, psia
T = Temperatura del reservorio, R
= Viscosidad, cp
Z = Factor de compresibilidad del gas.
re = Radio de drenaje, ft
rw = Radio de pozo, ft
S = Factor de dao.
4.3. Rgimen de Flujo en estado estable
Rgimen de flujo en estado estable existe cuando no hay cambio de presin en el borde
externo en funcin al tiempo. Prcticamente, tambin esto significa que el gradiente de
presin se mantenga con el tiempo ver Figura 4.4. que nos muestra esquemticamente la
distribucin radial de presin en torno de un pozo productor, en rgimen permanente.
Las condiciones que proporcionan el rgimen permanente de presin en determinadas

reas del reservorio son usualmente atribuidas a:
Influjo natural de agua proveniente de un acufero capaz de mantener la presin

constante en la frontera externa del reservorio.
Inyeccin de agua en torno del pozo productor de modo de contrabalancear la salida

de los fluidos del reservorio.
La relacin desarrollada por la ley Darcy para flujo de estado estable para un pozo de
gas natural es la ecuacin 4.16, introducindose un factor de dao s en la regin prxima
del fondo de pozo, la forma de rescribir la ecuacin 4.17 es:
1424q Z T re
p e2 p wf
2
= ln + s Ec.(4.18)
kh rw
Esta ecuacin sugiere que el rgimen de produccin de un pozo de gas es

aproximadamente proporcional a la diferencia de las presiones al cuadrado. Las
propiedades de y Z son propiedades media entre p e y p wf .
4.3.1. Rgimen de flujo de estado semiestables, (Pseudo-Steady State)
El estado pseudo-estable significa que la presin en el borde externo no se mantiene,

y al momento que el rgimen de flujo llega a tocar las fronteras, genera el agotamiento lo
que significa que la presin en el borde externo cae en funcin del caudal que sale del
yacimiento y esa cada de presin se refleja en todo el gradiente de presin en la misma
manera, en otras palabras 5 psi que caen en un da en el borde externo son 5 psi que caen en
cualquier punto del reservorio, por eso vemos esas 5 psi en un da en el pozo. Esto hace que
el gradiente de presin vaya cayendo sistemticamente tal como muestran en las

figuras 4.4 y 4.5.
El rgimen semi estable o rgimen seudo permanente de presin, usualmente ocurre

en las siguientes situaciones:
Pozo produciendo a un caudal constante de un pequeo reservorio cerrado.
Reservorio drenado por muchos pozos, con capa pozo aislado

hidrulicamente .
Figura 4.4 Distribucin radial de presin en rgimen permanente
Figura 4.5 Distribucin radial de presin en rgimen pseudo permanente

Para un sistema de geometra radial representado en la Figura 4.3, la condicin de

rgimen pseudo permanente puede ser expresado por:
P
= cons tan te
t
La Figura 4.5 ilustra las distribuciones radiales de presin en diferentes tiempos en un

reservorio cilndrico cerrado con un pozo en el centro produciendo a un mismo caudal
volumtricamente constante.
Matemticamente el escurrimiento del gas en rgimen pseudo permanente o semi

estable es tratado con una secuencia de rgimen permanente.
1424q Z T re
p e2 p wf
2
= ln 0.75 Ec.(4.19)
kh rw
Introduciendo el factor de dao incorporando el trmino 0.75 dentro de la expresin

logartmica, tenemos:
1424qZT r
p e2 p wf
2
= ln 0.472 e + s Ec.(4.20)
kh rw
Las ecuaciones 4.19 y 4.20 no son solamente aproximaciones en trminos de

propiedades, si no porque ellas asumen flujo de Darcy en el reservorio. Para caudales de
flujo de gas bastante pequeos esta aproximacin es aceptable. Una forma de
presentacin de las ecuaciones 4.20 es:
(
q = C p R2 p wf
2
) Ec.(4.21)
Para caudales de flujo ms grandes donde el flujo en Darcy es evidente en el

reservorio,
(
q = C p R2 p wf
2
)n Ec.(4.22)
Donde: 0.5 < n < 1

4.3.2.- Ecuaciones para flujo Radial en funcin al Seudo potencial (psuedo

presin o potencial de gas real)
Otra forma de presentar las ecuaciones bsica del flujo de Darcy`s, esta expresada de la
siguiente manera en base al seudo potencial:
pR r
0.3964 khTbZb / qbPbTr ( p / uz )dp = dr / r
p wf rw Ec(4.23)
Remplazando las constante Tb = 520 o R, Pb = 14,7 , Zb = 1 la ecuacin se puede escribir
pR
(0.703kh / qbT ) * 2 ( p / uz )dp = ln r / rw

p wf Ec.(4.24)
El termino (p/uz)*dp en la ecuacin 4.24 puede ser escrita
pR pr po
2 ( p / uz ) dp = 2 ( p / uz ) dp 2 ( p / uz ) dp
p wf po pw Ec. (4.25)
El trmino anteriormente mencionado, es la expresin de Kirchhoff de la

transformacin integral, y el contexto es llamado potencial del gas real que esta
usualmente representado por el m(p), se la siguiente manera:
pr
m( p ) = ( p / uz )dp
pw
Ec. (4.26)
Usando la seudo presin de un gas real en la ecuacin 4.24 se tiene
(0.703kh / qbT ) * m( p)r m( p ) f = ln r / rw

Ec. (4.27)
si el caudal de produccin esta en Mpcd y la presin base en 14.7 psi y Tb en 60 o F se

tiene la siguiente ecuacin :
m(p)r-m(p)f = 1422 * Qg *(Tr / kh) *ln (re/rw) Ec. (4.28)

4.3.3. Capacidad de entrega de un pozo de gas con flujo no-darciano
Una relacin ms precisa para un flujo estable de gas fue desarrollada por Aronofsky
e Jenkins que da la solucin de la ecuacin diferencial para un flujo de gas a travs de
medios porosos, usando la ecuacin de flujo de Forchheimer. Esta solucin es:
(
kh p R2 p wf
2
)
q= Ec.(4.29)
rd
1424 Z T ln + s + Dq
rw
Donde D es el coeficiente no Darcy y rd es el radio de drene efectivo de Aronofsky y

Jenkins, rd = 0.472re. Por otro lado:
rd
= 1.5 t D Ec.(4.30)
rw
Donde
0.000264 k t
tD = Ec.(4.31)
C t rw2
tD = Tiempo requerido para estabilizar el flujo.
El trmino, llamado con frecuencia efecto de turbulencia se da en los pozos de altos

caudales o potencial los cuales pueden ser substanciales. El coeficiente de turbulencia de
Darcy D, est en el orden de 10-3 y para caudales de gas se lo interpreta en trminos de Dq,
prximo al valor del logaritmo natural de la relacin ln rd / rw. Los valores pequeos de
caudal q resultaran proporcionalmente valores pequeos de Dq. Rescribiendo la ecuacin
4.29 se tiene:
1424 Z T 0.472re 1424 Z T D 2

p R2 p wf
2
= ln + s q + q Ec.(4.32)
kh rw kh
El primer trmino, de lado derecho de la ecuacin 4.32, es idntico al desarrollado

por Darcy. El segundo trmino, nos muestra el efecto de fluido no Darciano. Todos los
multiplicadores de q y q2 pueden ser considerados constantes, por tanto la ecuacin 4.32
puede tomar la siguiente forma:
p R2 p wf
2
= Aq + B q 2 Ec.(4.33)
El coeficiente no Darcy puede ser obtenido:

6 x10 5 k s 0.1 h
D= Ec.(4.34)
rw h 2perf
Donde:
= Gravedad del gas
ks = Permeabilidad prxima al fondo de pozo, md
h = Espesor neto, ft
h perf = Espesor perforado, ft
= Viscosidad del gas , cp
4.3.4.- Flujo Transiente pozo de gas
Flujo transiente de gas en un reservorio puede ser aproximado por la ley de Darcy y
la ecuacin de la continuidad, en general,
k
= p Ec.(4.35)
t
En coordenadas radiales se reduce a:
1 k p
= r Ec.(4.36)
t r r r
De la ley de los gases reales,
m pMW
= = Ec.(4.37)
V ZRT
Y por lo tanto:
p 1 k p
= rp Ec.(4.38)
t Z r r Z r
Si la permeabilidad k es considerada constante, luego la ecuacin 4.38 puede ser

aproximada:
p 1 p p
= r Ec.(4.39)
k t Z r r Z r
Realizando la diferenciacin en el lado derecho de la ecuacin 4.39 y asumiendo que

Z y son constantes diferentes y ellos cambian uniformemente y lentamente con la
presin. Tenemos:
1 p p 2 p p
2
+p + = RHS Ec.(4.40)
Z r r r 2 r
Reorganizando y recordando que
2
1 2 p2 2 p p
=p + Ec.(4.41)
2 r 2 r 2 r
Se convierte en
1 1 p 2 2 p 2
+ = RHS Ec.(4.42)
2 Z r r r 2
Por tanto la ecuacin 4.42 puede ser escrita como:
p 2 2 p2 1 p 2
= + Ec. (4.43)
k p t r 2 r r
para un gas ideal, Cg = 1 / p es como un resultado, la ecuacin 5.36 conlleva,
2 p2 1 p2 c p2
+ = Ec.(4.44)
r 2 r r k t
Esta aproximacin en la forma de la ecuacin de la difusividad es su solucin, bajo

las suposiciones listadas en esta seccin pudiera tener la forma de las soluciones de la
ecuacin para petrleo, presumiendo que p2 es usado en lugar de p. Diferencia de presiones
al cuadrado pueden ser usadas como una aproximacin razonable, por ende, las
suposiciones usadas que derivan de la ecuacin 4.44 estn limitadas y de hecho, ellas
pueden llevar grandes errores en pozos de altos caudales con variaciones grandes en la
presin fluyente. Por tanto una mejor solucin puede ser desarrollada, usando la funcin
pseudo presin para gas real de Al-Hussainy y Ramey, definida como la funcin de m(p),
p
p
m( p ) = 2 dp Ec.(4.45)
po
Z
Donde po algunas veces es una presin de referencia arbitraria (puede ser cero). La
pseudo-presin diferencial m( p ) , definida como m(p) m (pw f ), entonces es la fuerza
impulsora en el reservorio. Para bajas presiones puede ser mostrado que:
pi
p p i2 p wf
2
2 Z
dp
Z
Ec.(4.46)
pwf
Considerando presiones altas ( pi y pw f superiores que 3000 psi ),
pi
Z dp 2 Z ( pi )
p p
2 p wf Ec.(4.47)
pwf
La pseudo-presin de un gas real puede ser usada en lugar de las diferencias de

presiones al cuadrado en cualquier relacin de capacidad de entrega de un pozo de gas,
La ecuacin 4.19 tendra la siguiente forma:
[ (
kh m( p ) m p wf )]
q= Ec.(4.48)
0.472 re
1424T ln + s + Dq
rw
Puede ser usada de forma ms apropiada la pseudo-presin de un gas real como un

factor de integracin para una solucin exacta de la ecuacin de la difusividad para un gas,
la ecuacin 4.39 puede ser la base para este anlisis.
La definicin de pseudo-presin de un gas real puede ser escrita,
m( p ) m( p ) p 2 p p
= = Ec.(4.49)
t p t Z t
Similarmente
m( p ) 2p p
= Ec.(4.50)
r 2 Z r
Por tanto, la ecuacin diferencial se convierte en:
2 m( p ) 1 m( p ) c t m( p )
+ = Ec.(4.51)
r 2 r r k t
La solucin de la ecuacin es exactamente similar a la solucin para la ecuacin de la

difusividad en trminos de presin. El tiempo y presin adimensionales estn definidos
como:
0.000264 k t
tD = Ec.(4.52)
( C t ) i rw2
[
kh m( p i ) m p wf ( )]
pD = Ec.(4.53)
1424qT
La ecuacin de difusividad generalizada para pozos de petrleo es:
q
p r .t = p i Ei (x ) Ec.(4.54)
4 k h
Donde E i ( x ) es la funcin integral exponencial de x y es dado por:
ct r 2
x= Ec.(4.55)
4k t
Para x < 0.01, la integral exponencial E i ( x ) puede ser aproximada por el ln ( x )
donde es la constante de Euler y es igual a 1.78.
(
Por tanto la ecuacin 5.47 se convierte p r .t p wf )
q 4k t
p wf = p i ln Ec.(4.56)
4 k h ct rw2
La aproximacin logartmica para la integral exponencial, conduce a una expresin

anloga para un gas natural en trminos de pseudo-presin.
[ (
k h m ( p i ) m p wf )]
log t + log
k
1
q= 3 .23 Ec.(4.57)
1638 T ( c t )i rw2
4.4.- Tipos de pruebas
La habilidad de analizar el comportamiento y los pronsticos de productividad de un

pozo de gas, nos dan las distintas pruebas de produccin que se pueden realizar en el pozo,
dndonos un mejor entendimiento del comportamiento del reservorio, con un grado de
seguridad que es de suma importancia en la industria del gas natural.
Un completo anlisis y entendimiento de los resultados de una prueba de pozo nos

determina el comportamiento del caudal para los distintos dimetros de tubera y el
comportamiento de flujo con la reduccin de la presin de reservorio.
4.4.1.-Prueba de flujo tras Flujo (Flow-After-Flow tests)
Llamada tambin pruebas convencionales de contrapresin (Conventional

Backpresure Test). En este tipo de prueba, el pozo se fluye a un determinado caudal
midiendo la presin fluyente de fondo la cual normalmente se mantiene en estado transiente
(no alcanzando el estado pseudo-estable). Luego el pozo cambia su flujo a un nuevo
rgimen, normalmente en estado transiente sin llegar estado pseudo estable. Ver figura 4.6.
La presin puede ser medida con un medidor de presin de fondo de pozo.
Este proceso es repetido para diferentes rgimenes de flujo estabilizados. Esto se

puede realizar para un nmero indeterminado de perodos de flujo, normalmente es
recomendable que sean cuatro, al final de la prueba de multi-flujo se efecta un cierre de
pozo cuyo tiempo depende del tiempo de estabilizacin.
Figura. 4.6 Prueba de flujo tras flujo

4.4.2.- Prueba Iscronal (tiempo de flujo tiempo de cierre)
En un reservorio de baja permeabilidad es muy frecuente que el cierre despus de un

flujo no llegue a la estabilizacin y es imprctico extender por mucho tiempo el cierre si las
condiciones de pozo no han llegado a un estado semi estable. El objetivo de la prueba
isocronal, es obtener datos representativos para establecer una curva de capacidad de
entrega estable produciendo el pozo a un flujo estable con el tiempo de cierre suficiente
para obtener datos estabilizados en cada prueba. El Radio de investigacin alcanzado en la
prueba a un determinado tiempo es independiente del caudal de flujo. Por tanto, si una serie
de pruebas de flujo son ejecutadas en un pozo, para cada uno por el mismo periodo de
tiempo, el radio de investigacin ser el mismo al fin de cada prueba.
La Figura 4.7 nos muestra el comportamiento de la prueba de flujo isocronal en

funcin al caudal y a la presin de fondo fluyente. Observndose que en el periodo de
cierre despus de cada periodo de flujo debe estar en funcin al tiempo de estabilizacin
determinado para la prueba, llegando a restituir a su presin de reservorio. Obsrvese que
tambin es necesario que haya un periodo de flujo estabilizado al fin de la prueba.
Figura. 4.7 Prueba de Flujo Isocronal
4.4.3.- Prueba Iscronal Modificada (tiempo de flujo = tiempo de cierre)
Este tipo de prueba esta diseada principalmente a reservorios de baja permeabilidad,

ya que el tiempo de estabilizacin del flujo radial es elevado tanto para los periodos de flujo
como para los periodos de prueba, y la variante que presenta frente a las pruebas isocronal
es que el periodo de flujo es igual al periodo de cierre y no se requiere alcanzar las
condiciones estabilizadas de presin entre cada etapa de flujo ( ver Figura 4.8) la cual nos
muestra un diagrama esquemtico del caudal y las presiones de fondo fluyente.
Figura 4.8 Prueba de Flujo Isocronal Modificada
4.4.4.- Pruebas de Produccin
Estas clases de prueba se realizan continuamente en el campo para determinar los

volmenes producidos por pozo y as poder controlar la produccin acumulada con el
tiempo. Para que estas pruebas sean vlidas para el anlisis nodal, es importante que los
pozos produzcan con velocidades mayores a las crticas para arrastrar el condensado y el
agua que se nos acumula en el fondo del pozo y las condiciones de produccin y presin
tienen que estar estabilizadas, la cual nos distorsiona las presiones fluyente, y por ende los
potenciales de pozo, las velocidades crticas para el arrastre de los fluidos en el fondo
mostraremos en captulos prximos.
Con los datos de produccin y las presiones fluyentes de fondo obtenido en una
prueba de produccin, podemos determinar los ndices de productividad del pozo y el
comportamiento de entrega como mostraremos en la siguiente ecuacin:
IP = Qg/ Pr^2-Pwf^2 Ec.(4.58)
4.5 Mtodo de Interpretacin de prueba
La habilidad de analizar el comportamiento y pronstico de los pozos de gas y el

ndice de productividad de los mismos, no los dan los mtodos de interpretacin de prueba,
ya que podemos realizar un anlisis completos de los resultados que arrojan las pruebas de
pozos tanto de produccin como as tambin de cierre. Previniendo los distintos problemas
que se pueden presentar con la declinacin de la presin de reservorio y los efectos de dao
del pozo. Los resultados de las pruebas de pozo son a menudos utilizados para
optimizar o maximizar la produccin, previniendo el desarrollo del campo, y las facilidades
en superficie como plantas de procesamientos.
Bsicamente existen dos tipos de datos para la determinacin de la capacidad de

entrega:
Datos de pruebas (Isocronales, Flujo tras Flujo, Prueba de Produccin)

Datos de Reservorio
En la bibliografa actual, existen varios mtodos de interpretacin de pruebas de los cuales

tomaremos los tres ms principales para nuestro propsito.
4.5.1. Mtodo Simplificado
En 1936, Rawlins y Schellhardt, presentaron la ecuacin 4.22 , como la ley de Darcy

para un fluido compresible, donde C contiene todos los trminos diferentes de la presin;
como la viscosidad del gas, permeabilidad al flujo de gas, la temperatura de la formacin,
etc. Rawlins y Schellhardt, describen que la ecuacin 4.22 Krg era responsable por la
turbulencia normalmente presente en pozos de gas, entonces modificaron la ecuacin con
un exponente n, ecuacin 4.23.
El exponente n puede variar de 1.0 para flujo completamente laminar y 0.5 para un
flujo completamente turbulento.
Si los valores para el coeficiente de flujo C y exponente n puede ser determinado por el
rgimen de flujo, para cualquier valor de p wf , puede ser calculado, el caudal y se puede
construir la curva del comportamiento de flujo de entrada. Un parmetro comnmente
usado para ver el potencial cuando la p wf = 0 , es llamado Potencial Absoluto de Flujo
Abierto (AOF), el cual es definido como el mximo caudal que un pozo de gas producira
sin contrapresin.
Considerando el mtodo clsico, se tiene dos constantes para determinar C y n.

La teora indica que C es una funcin de radio de investigacin que significa que si dos
periodos de flujo poseen un mismo radio de investigacin, ellas tendrn el mismoC. Las
razones de flujo poseen un mismo intervalo de tiempo, entonces tendr un mismo radio de
investigacin y por tanto un mismo C. Para perodos estables de flujo, el C ser el
C estabilizado, que es el que estamos tratando de determinar. Para una serie de periodos
de flujo iguales que no son largos o suficientes para alcanzar la estabilizacin, los Cs de
cada prueba sern los mismos, mas no sean los C estabilizados.
Si el pozo ha fluido a un caudal estabilizado, como se muestra en el esquema log-log

(Figura 4.9) podemos determinar un mximo potencial transiente de la prueba, si tenemos
un flujo extendido como se muestra en la (Figura 4.10 ) obtendremos un punto estabilizado
por la cual pasamos una lnea paralela a la lnea de los puntos transiente. De modo que el
valor de C es el valor estabilizado, y tambin el mximo potencial determinado. La

ecuacin 4.22 puede ser escrita de la siguiente manera:
(
log p R2 p wf
2
=
1
n
) 1
log q g log C
n
Ec.(4.59)
2 2
La grfica logartmica log-log de la diferencial de presin p R p wf ( ) versus q g , nos
muestra una lnea recta (Figura 4.9), el factor de turbulencia expresado por ( n ) es inversa a
la pendiente de esta lnea. La figura tambin nos muestra, una prueba de produccin con
cuatro caudales de flujos, que estaran sobre una misma lnea recta mostrando una
condicin de flujo estabilizado. El valor del exponente n relacionada al coeficiente de
turbulencia se puede determinar grficamente de la diferencia de caudal dividido en
relacin de la diferencial de presin en cualquier punto de la lnea recta. En forma analtica
mostraremos en la ecuacin 4.60
log q 2 log q1
n=
( ) ( )
Ec.(4.60)
log p R2 p wf
2 2 2
2 log p R p wf 1
Una vez determinado el valor del exponente n , el valor C se puede determinar usando la
siguiente ecuacin:
qg
C= Ec.(4.61)
( p R2 2 n
p wf )
El valor de la constante C en base a los datos de reservorio pude ser representado por la
siguiente ecuacin:
0.000703Kgh g MPCD
C= Ec. (4.62)
UgZTr (ln(0.472rd / rw) + s ) PSI 2
Figura. 4.9 Anlisis Simplificado de Prueba

Figura. 4.10 Anlisis Simplificado para un caudal estabilizado
4.5.2. Mtodo Jones Blount and Glaze
En 1976 Jones, Blount y Glaze, sugieren un procedimiento de anlisis que permite

determinar el efecto de turbulencia o no, que se presenta en la completacin de pozos
independiente del efecto de dao y flujo laminar. El procedimiento tambin avala el
coeficiente de flujo laminar A, y el efecto de dao si el producto k g h es conocido.
La ecuacin presentada por Jones, et al,. para flujo de estado estable (steady-state
flow) incluyendo el factor de turbulencia es:
1424T Z q re
p R2 p wf
2
= ln + s .......
kh rw
Ec.(4.63)
3.161x10 12 g Z q 2T 1 1
... +
h2 rw re
El primer trmino de lado derecho es la cada de presin de flujo laminar o flujo

Darciano, en cuanto al segundo trmino, es la cada de presin adicional debido a la
turbulencia. El coeficiente de velocidad , es obtenido en la ecuacin 4.9
Algunas veces es conveniente establecer una relacin entre dos parmetros que
indican el grado de turbulencia que ocurre en un reservorio de gas. Estos parmetros
son: el coeficiente de velocidad , y el coeficiente de turbulencia . La ecuacin 4.63
se puede describirse para un flujo de estado semi estable o pseudo estable como:
1424T Z 0.472re
p R2 p wf
2
= ln + s q.......
kh rw
Ec.(4.64)
3.161x10 12 g Z T
... + q2
2
rw h
Los trminos de la ecuacin 4.64 son agrupados en dos coeficientes de la siguiente

manera:
1424T Z 0.472re
Coeficiente Laminar A= ln + s Ec.(4.65)
kh rw
3.161x10 12 g Z T
Coeficiente Turbulencia B= Ec.(4.66)
rw h 2
Por tanto, la ecuacin 4.64 toma la forma de la ecuacin 4.33, dividiendo esta por q
toma la forma de la ecuacin general propuesta por Jones, Blount y Glaze.
p R2 p wf
2
= A + Bq Ec.(4.67)
q
Para determinar los dos coeficientes existen dos formas: La primera hace uso de las
pruebas convencionales con dos o ms valores de flujo estabilizado, por lo menos un flujo
estabilizado en pruebas de flujo isocronal. Los datos de caudal y presin obtenidos en la
conduccin de estas pruebas son producidos en coordenadas cartesianas como
(p R2 p wf2 )/ q , en el eje de las coordenadas y q , en el eje de las abscisas. En la figura
4.10 el diagrama resultante muestra una lnea cuya pendiente es el coeficiente B que indica
el grado de turbulencia. Prolongando la recta hasta el ejes de las coordenadas se tiene el
coeficiente laminar A, adoptando en este caso el valor de (p R2 p wf2 )/ q para un caudal
igual que cero, resultado que muestra la existencia o no, de dao a la formacin.
El segundo camino es de simples substituciones de los parmetros, previamente

determinados en las ecuaciones 4.65 y 4.66 .Una vez determinados los coeficientes A Y B
se procede a la construccin de la curva del comportamiento de IPR, asumiendo diferentes
valores de presin de fondo fluyente, p wf , determinando los caudales para estos mismo
valores. Tambin podemos asumir los caudales de produccin y determinar las

presiones fluyentes indiferentemente. Las ecuaciones presentadas son:
q=
A + (
A 2 + 4 B p R2 p wf
2
) Ec.(4.68)
2B
p wf = (
p R2 Aq + B q 2 ) Ec.(4.69)
Figura 4.11 Anlisis Mtodo Jones, Blount & Glaze
Sabemos que el valor de A = P^2/Qgcs Ec. (4.70)

B = (P/Qg) /Qg Ec. (4.71)
Para tener alguna medida cualitativa de importancia en la contribucin del efecto de

turbulencia en una estabilizacin, Jones et al. Sugiri la comparacin de el valor de A
calculado del Potencial AOF y el valor de A` estabilizada calculada con la siguiente
formula:
A = A + B (AOF) Ec. (4.72)
Jones et al. Sugiri que si la razn de A y A es mayor que 2 o 3, existe alguna

restriccin en la completacin del pozo. Ellos tambin sugieren que si el espesor de la
formacin h es usado en el calculo de B, podra ser remplazado por el espesor perforado,
teniendo una cada de presin por efecto de turbulencia. Por lo tanto, para optimizar la
produccin y evitar este efecto de turbulencia se puede calcular un nuevo valor de B2

cuya formula es la siguiente:
B2 = B1 (hp1/hp2) Ec. (4.73)
Donde
B2 = coeficiente de turbulencia despus de la completacin.

B1 = coeficiente de turbulencia antes de la completacin.
Hp1 = espesor perforado inicial
Hp2 = espesor perforado nuevo
Anlisis de la terminacin del pozo despus de una prueba:
1. Si el valor de A es bajo y menor a 0.05 no existe dao en la Formacin

2. Si el valor de A/A es bajo menor que 2 existe una pequea turbulencia en el
pozo.
3. Si el valor de A y A/A son bajos el pozo tiene una buena completacin
4. Si el valor de A es bajo y A/A es alto existe insuficiencia de rea de perforacin
y la estimulacin cida no es recomendable.
5. Si el valor de A es alto y A/A es bajo es recomendable una estimulacin de pozo
4.5.3.-Mtodo Brar y Aziz
La metodologa de clculo presentada difiere de los ya estudiados, no es necesario

contar con pruebas hasta alcanzar por lo menos un dato estabilizado, a partir de pruebas
transientes se puede estimar la capacidad de entrega de un pozo cuya formacin es poco
permeable, la cual no afecta a la construccin de la curva de comportamiento del IPR.
Como el periodo transiente esta ligado con las variaciones del tiempo, las constantes
de la ecuacin 4.67 estn determinadas para distintos periodos de pruebas por lo tanto, el
valor de A se convierte en At , valor que crece hasta un mximo, mantenindose constantes
en este punto. Se debe ignorar el cambio del coeficiente B y tomar en cuenta slo el que
corresponde a la ltima etapa de flujo.
Para condiciones de estado semi estables (pseudo-steady state),
(p R2 p wf2 ) = 2mlog 0.472

r
re
+
s
2.303
q sc ...
w Ec.(4.74)
... + 0.869 mDq sc 2
Donde
1637 T Z
m= Ec.(4.75)
kh
Para flujo transiente, (transient flow),
(p R2 p wf2 ) = mlog kt
2

3.23 + 0.869s q sc ...
C rw Ec.(4.76)
2
.... + 0.869mDq sc
Comparando las ecuaciones 4.74 y 4.75, los coeficientes A y At son representados por
las siguientes frmulas:
0.472re s
A = 2m log + Ec.(4.77)
rw 2.303
kt
At = m log 3.23 + 0.869s Ec.(4.78)
C rw2
y el coeficiente indicador del grado de turbulencia es igual para ambas expresiones:
B = 0.869mD Ec.(4.79)
El objetivo de analizar y determinar los valores de A y B para flujo estabilizado,

ecuacin 4.33, es el de determinar la curva de comportamiento de la capacidad de entrega
del pozo (IPR.). El factor de dao s, el coeficiente de turbulencia o no Darcy D se puede
determinar con la siguiente ecuacin:
La ecuacin 4.76 es posible escribirla como:
(p R2 p wf2 ) = A + Bq sc Ec.(4..80)
t
q sc
Donde At y B son definidos en las ecuaciones 4.78 y 4.79 respectivamente. El valor

de At se incrementar hasta que el tiempo de flujo estabilizado sea alcanzado. Un
( 2 2
)
esquema p R p wf / q sc vs q sc en coordenadas cartesianas resultar en una serie
de lneas rectas, lneas paralelas teniendo pendientes iguales para B y la interseccin At
( 2 2
)
igual a p R p wf / q sc para cada flujo de tiempo. La pendiente y la intercepcin
pueden ser tambin determinadas usando el anlisis de los mnimos cuadrados. La
ecuacin puede ser expresada como:
kt
At = m log 3.23 + 0.869s + m log t Ec.(4.81)
C rw2
La Figura 4.12, nos muestra el comportamiento de At versus t, en escala semi-log

resulta una lnea recta teniendo una pendiente de recta igual a m que tienen de unidades
Mpsi^2/MMpcd/ciclo y una interseccin en t = 1hr (log 1 = 0) igual a At1. El
procedimiento para analizar una prueba de flujo isocronal o isocronal modificado es:
1. Determinar At y B a partir de pruebas transientes para varios tiempos de

flujo utilizando la ecuacin 4.76 mnimos cuadrados.
2. Haga un grfico At versus t en escala semi-log para determinar m e At1 ,

similar a la figura 4.11.
3. Utilizando el valor de m , determinar el valor de la permeabilidad k haciendo

uso de la ecuacin 4.75.
4. Obtener el valor de s con la ecuacin 4.78 utilizando los valores de m , k , y
At1..
5. Determine un valor estabilizado para A utilizando la ecuacin 4.77 .
6. Utilizando el valor de B determinado en el paso 1, calcular D utilizando la

ecuacin 4.79.
7. La curva que muestra la relacin de comportamiento del influjo construida

siguiendo el mismo procedimiento descrito por Jones, Blount y Glaze al igual que el
potencial absoluto del pozo, AOF, usando los valores de A y B.
El mtodo de mnimos cuadrados puede ser usado para determinar A y B con N

periodos de flujo transientes.
p2
q
( ) q 2 (p 2 ) q
At = Ec.(4.82)
N q 2 q q
( )
N p2
( ) q
p2
q
B= Ec.(4.83)
2
N q q q
Figura 4.12 Mtodo Brar e Aziz
Valores para At y B sern obtenidos para cada tiempo en la cual p wf fue medido. El
valor de B ser constante, Brar y Aziz sugirieron utilizar el valor de B obtenido para pruebas
de flujo ms largas como el valor representativo.
4.5.1.1 Mtodo de Anlisis LIT (Pseudo-presiones)
La importancia de considerar las variaciones de viscosidad y el factor de

compresibilidad con la presin, en reservorios muy compactos donde el gradiente de
presin es pocas veces pequeo, debido a los problemas de las variaciones del factor de
compresibilidad y viscosidad del gas en determinados rangos de presin se utilizo la pseudo
presin para obtener ms correctos, la ecuacin 4.45, logrando as mismo un anlisis ms
riguroso de los fenmenos de flujo, conocido tambin como Anlisis Laminar Turbulento
(LIT). La ecuacin diferencial de Pseudo presiones para un rgimen semi-estable es
modificado para dar una expresin equivalente a la ecuacin 4.67 .
( )
m( p ) = m( p R ) m p wf = Aq sc + Bq sc
2
Ec.(4.84)
Los coeficientes A y B indican tambin el tipo de flujo: laminar y turbulento

respectivamente, estos coeficientes se obtienen mediante ponderacin utilizando el
concepto de mnimos cuadrados siguiendo el mismo procedimiento descrito en el mtodo
de Brar y Aziz, excepto que en las ecuaciones 4.82 y 4.83 se entregan valores de diferencial
de Pseudo-presin en lugar de diferencial de presin al cuadrado.
m( p )
q 2 m( p ) q
q
At = Ec.(4.85)
N q 2 q q
m( p )
N m( p ) q
q
B= Ec.(4.86)
N q 2 q q
Una ves encontrados los coeficientes de A y B se substituyen en la ecuacin 4.84

encontrando de esta manera la ecuacin general para este mtodo, visualizando el
comportamiento del influjo, construyndose siguiendo el mismo procedimiento descrito
por Jones, Blount y Glaze, excepto que se entregan valores de diferencial de Pseudo-
presin en lugar de diferencial de presin al cuadrado.
Ejemplo No 4.1 Se tiene una prueba isocronal y se desea conocer el potencial de

reservorio y el IPR para los mtodos: simplificado, Jones Blount Glaze y LIT
Datos
Pr = 10477 psi Prof. = 14300 pies
Tr = 270 oF API = 59
Tc = 95 oF SGg = 0.65
Prueba Isocronal
Perodo de flujo Choque Duracin Presin (Pwf) Presin (Pwh) Caudal Caudal Caudal
Periodo Fondo Cabeza Gas Condensado Agua
CK/64" hr psia psia MMscf/d BPD BPD
esttica inicial 24 10477 9044
flujo 1 12 12 7815 5950 4,817 72 6
flujo 2 16 12 6009 4380 6,296 88 11
flujo 3 20 12 4865 3350 7,337 98 17
flujo 4 24 12 3978 2460 8,08 109 17
Flujo extendido 16 38 6188 4500 6,281 89 10
Del
Pseudo M(p)-
Diferencial DP^2/Q Potencial Diferencial Del M(p)/Qg Qg2 bQg2 Caudal
Presin
(P^2) M(p) Del M(p) Gas
MMPSI^2 Mmpsia^2/MMscfd Mmpsi2/cp Mmpsi2/cp Mmpsi2/cp/Mmscfd (MMscf/d)^2 MMscf/d
109,77 3653
48,69 10,11 2373 1280 265,78 23,20 752 4,817
73,66 11,70 1626 2027 321,93 39,64 1124 6,296
86,10 11,73 1204 2449 333,75 53,83 1223 7,337
93,94 11,63 904 2749 340,17 65,29 1262 8,08
71,48 11,38 1696 1957 311,61 39,45 1059 6,281
Sumatoria 8504 1261,63 181,96 4360 26,53
Prueba Isocronal mtodo Simplificado
1000,00
Diferencial de Presin PSI^2
100,00
10,00
1,00 ISOCRONAL
0,10
1 10 100 1000
CAUDAL DE GAS MMpcd
log q1 log q 4
n= = 0.787

log p 2 p 2 log p 2 p 2
R wf 1 R wf 4
q g = 0.0000041 MMPcd/Psi^2
C =
(p R2 p 2
wf ) n
Qg = C (Pr^2-Pwf^2)^n = 9.13 MMpcd
Mtodo de Blount Glaze y = 0,4563x + 8,3157

R2 = 0,6753
13,91
11,91
Dp / Qg (psi^2/Mmpcd)
9,91
7,91
5,91
3,91
Metodo Jones Blunt and Glaze
1,91
Lineal (Metodo Jones Blunt and
-0,09 Glaze)
0 1 2 3 4 5 6 7 8 9
Caudal de Gas Mmpcd
Determinacin de coeficientes "A" e "B" e AOF

Mtodo de Blount and Glaze
A 0,465209185
B 8,3157 Mmpsia^2/Mmscfd
B' 12,43 Mmpsia^2/Mmscfd
B'/B 1,494211175
AOF Tras. 8,834 Mmscfd
Mtodo Anlisis Lit
10000,00
Del m(p)- b Qg^2
1000,00
100,00
10,00
1,00 Metodo Analisis Lit

1 10 100
Qg MMPCD
Coeficientes Laminar e Turbulento "A" e "B"
Coeficiente
Laminar A 164,3435929 MMpsia^2/cp/MMscfd
A' 168,56 MMpsia^2/cp/MMscfd
Coeficiente
Turbulento B 22,77641183 MMpsia^2/cp/(MMscfd)^2
Potencial
transiente AOF trans. 9,56 Mmscfd
Potencial
Estabilizado AOF Est. 9,49 Mmscfd
m( p R ) m( p wf ) = 164.344 * Qg + 22.76 * Qg ^ 2
Resumen del comportamiento del IPR con los tres mtodos expuestos anteriormente:
Mtodos Simplificado Blount Glaze Lit

Presin Pwf Caudal Caudal Caudal
Qsc
Asumida (Mmscfd) Qsc (Mmscfd) Qsc (Mmscfd)
10477 0,00 0,00 0,00

9429 2,47 2,23 2,41
8382 4,09 3,90 4,01
7334 5,38 5,21 5,25
6286 6,43 6,26 6,25
5239 7,28 7,09 7,10
4191 7,96 7,74 7,82
3143 8,48 8,23 8,44
2095 8,84 8,57 8,97
1048 9,06 8,77 9,42
0 9,13 8,83 9,56
Ejemplo No 4.2 Se desea conocer el mximo potencial del reservorio y el IPR con
datos de reservorio para los siguientes Mtodos Darcy`s Jones Blount Glaze y Pseudo
potencial.
Datos de pozo son:
Pr = 10477 psi SGg = 0.65 Rw = 0.35 Lp = 0.88 pies

Tr= 270 o F API = 59 Rd = 1000 pies Rp = 0.021 pies
Tc = 95 o F RGC = 71603 pc/Bbl Kg = 1.23 md 12 Balas/pie
Prof.= 14300 pies RGA = 475224 pc/Bbl S = 17 Bajo balanceado
Hn = 62 pies U= 0.67 Krg = 0.25 md Porosidad 0.18 %
Hp = 48 pies Sw = 0.60 Swi = 0.25 Rc= 0.063 pies
Mtodo Darcy Pseudo Presin

PRESION CAUDAL Caud./Perf. PRESION SEUDOPRESION CAUDAL Caud./Perf. PRESION
ASUMIDA CALCULADO Fluy.Fondo ASUMIDA (M(p)) CALCULADO Fluy.Fondo
PSIA MMPCD MPCD PSI PSIA/cp MMPCD MPCD PSI
10477 0,00 0 10477 3653 0,00 0 10477

9429 0,97 2 9417 3123 1,61 3 9407
8382 2,00 3 8353 2627 3,11 5 8333
7334 3,16 5 7283 2164 4,52 8 7255
6286 4,49 8 6204 1734 5,82 10 6172
5239 6,06 11 5109 1338 7,02 12 5080
4191 7,94 14 3981 974 8,13 14 3973
3143 10,14 18 2773 644 9,13 16 2830
2095 12,56 22 1272 347 10,03 17 1562
1048 14,82 26 83 10,83 19
0 16,31 28 11,08 19
METODO DE JONES BLOUNT & GLAZE
RESERVORIO PERFORACIONES (BLOUNT GLAZE)
Pr^2-Pwf^2 = Aq^2+Bq Pr^2-Pwf^2 = Aq^2+Bq

A=3.16*10^-12 Be SGg Tr Zr / hp^2 Rw A= 3.16*10^-12 Be SGg Tr Zr (1/Rp - 1/Rc) / Lp^2
B= 1.424* 10^3 u Tr Zr (ln (0.472 Re/RW) +S) / K h B= 1.424*10^3 u Tr Zr (ln Rc/ Rp) / Kp Lp
Be= 2.33*10^10/ k ^1.201 Be = 2.33*10^10/ Kp^1.201
Beta psi^-1 1,8171E+10 Beta psi^-1 56934663738
Pr^2-
PRESION Term. Flujo CAUDAL Trmino Trmino Caud./Perf. PRESION
Pwf^2
ASUMIDA Darcy CALCULADO Fluy.Fondo
PSIA B PSI^2 MMPCD A B MPCD PSI
10477 17590,892 0 0,00 5401 140619 0,00 10477

9429 16659,766 20855831 1,25 5231 133176 2,17 9413
8382 15295,586 39516310 2,56 4979 122271 4,45 8343
7334 13740,436 55981440 4,02 4690 109839 6,98 7266
6286 12123,515 70251219 5,68 4387 96914 9,86 6175
5239 10523,772 82325647 7,61 4089 84125 13,20 5061
4191 9004,412 92204724 9,84 3816 71980 17,09 3901
3143 7635,658 99888451 12,38 3594 61038 21,50 2628
2095 6502,630 105376828 15,05 3450 51981 26,12 824
1048 5684,756 108669854 17,37 3402 45443 30,16
0 5214,631 109767529 18,80 3432 41685 32,64
En realidad no existe gran variacin entre los mtodos, y la variacin de caudal disminuye
si consideramos las cadas de presin por las perforaciones
4.6.- PRODUCCIN POZOS HORIZONTALES
En los aos 1980s, los pozos horizontales empezaron a tener una mayor importancia
en el sector petrolero debido a muchos problemas encontrados con los pozos
convencionales verticales, principalmente en las perforaciones marinas por el alto costo de
produccin y recoleccin de los mismos. Debido a los problemas de conificacin de agua,
gas y arenamiento por su alta diferencial de produccin. Por lo tanto, se decidi optimizar
la produccin con pozos horizontales.
Un pozo horizontal de longitud L que penetra un reservorio con permeabilidad

horizontal Kh y la permeabilidad vertical Kv, crean un modelo de drenaje diferente a un
pozo vertical con una mayor rea de flujo y una menor presin diferencial. La Figura. 4.13
nos muestra el patrn de flujo junto con las ms importantes variables que afectan el
comportamiento del pozo. La forma del drenaje es elipsoidal, con el eje a la mitad de
drenaje del elipsoide, relativa a la longitud horizontal del pozo ver Fig. 4.13.
Figura 4.13 Patrn de Flujo Formado alrededor pozo Horizontal
Es muy importante considerar la anisotropa de la permeabilidad vertical en los pozos

horizontales. Mientras ms grande sea permeabilidad vertical, ms alto es el ndice de
productividad de un pozo horizontal. La permeabilidad vertical baja puede dar pozos
horizontales poco atractivos para cualquier inversin.
A menudo ignorado, pero importante, es el problema de la anisotropa de

permeabilidad horizontal. Un pozo normal perforado a lo largo de la permeabilidad
horizontal, sera mucho mejor que un pozo perforado en cualquier direccin arbitraria o
normal a la permeabilidad horizontal ms pequea.
Siempre se dan medidas o parmetros antes que un pozo horizontal sea perforado,
estas medidas son hechas en base a un pozo piloto vertical. Para la anisotropa de
permeabilidad horizontal, y la medida de los esfuerzos o tensin en el agujero del piloto se
puede identificar el mximo y mnimos de las tensiones horizontales. Normalmente, ellos
coinciden con el mximo y las direcciones de permeabilidad horizontales mnimas. Por
consiguiente, un pozo horizontal que no se piensa que es fracturado hidrulicamente debe
perforarse a lo largo de la direccin de tensin horizontal mnima.
Joshi, present una relacin de entrega de un pozo horizontal que fue aumentada por
et y Economides en el 1990. La relacin muestra (mezcla del pseudo estado en un plano
horizontal y pseudo estado estabilizado en un plano vertical) es:
k H hp
q= Ec.(4.87)
a + a 2 ( L / 2) 2 I h I ani h
141.2 ln + ani ln

L/2 L [rw ( I ani + 1)]
Donde I anillo, es una medida de la anisotropa de la permeabilidad vertical y

la permeabilidad horizontal que esta dada por:
kH Ec.(4.88)
I ani =
kV
En la ecuacin. 4.87 a es el medio eje grande del elipsoide de drenaje formado por un
pozo horizontal de longitud donde la expresin para este elipsoide es:
0.5
L
0.5
re H
4
L Ec.(4.89)
a = 0.5 + 0.25 + para < 0.9re H
2
L / 2 2

La relacin entre el ndice de productividad de un pozo horizontal y un pozo vertical

en un reservorio especfico podra ser muy grande. (Asumiendo que un candidato apropiado
es seleccionado, el pozo es perforado en una direccin ptima, y es estimulado muy
efectivamente). Esta razn del ndice de productividad puede ser manifestada por un
incremento en el caudal de produccin, o un decremento en la cada de presin o ambos.
Por consiguiente, los pozos horizontales pueden ser excelentemente manejados donde los
problemas de agua y cono de gas y arena estn presente.
4.6.1 Impacto del efecto de dao en el comportamiento de un pozo Horizontal
El efecto de dao en un pozo horizontal es adicionado al denominador Eq. (4.90) de

la manera siguiente:
k H hp
q=
a + a 2 ( L / 2) 2 I h I ani h Ec. (4.90)
141.2 ln + ani
ln + S '
L [rw ( I ani + 1)]
eq
L/2

Este efecto de dao, denominado como S'eq, es caracterstico de la forma de dao en

pozos horizontales, tomando en cuenta la anisotropa de la permeabilidad y probabilidad de
penetracin de dao ms profundo, o ms cercano a la seccin vertical.
El impacto de este efecto de dao en la reduccin de los caudales puede ser muy
grande. La primera expresin logartmica en el denominador de la Ec. (4.90) esta en el
rango entre 1.5 y 3, para la mayora de las aplicaciones. La segunda expresin logartmica
va entre 2.5 y 4.5, donde S'eq puede ser hasta 50, con valores comunes aproximadamente
20. Aun cuando es multiplicado por Iani *h / L, en un rango de 0.02 y 0.3, este efecto
en los caudales de produccin puede ser sustancial.
4.6.2.-Efectos de Produccin de Agua y Permeabilidades Relativas
En un reservorio gasfero, el agua est siempre presente por lo menos como

connata, o irresidual denotada como Swc. As, en todas las ecuaciones anteriores la
permeabilidad debera ser considerada como efectiva, y debera ser invariable (en ciertos
casos significativamente menor) que la obtenida de las muestras de ncleos o de otras
tcnicas de laboratorio usando un simple fluido.
Si ambos fluidos gas y agua estn fluyendo, la permeabilidad efectiva debera ser
usada, como la suma de estas permeabilidades que es invariable y menos que la
permeabilidad absoluta de la formacin (para cualquier fluido).
Esta permeabilidad efectiva esta relacionada a la permeabilidad relativa (tambin para

las distintas propiedades de la roca).
Kg=K*Krg Ec. (4.91)
Kw=K*Krw Ec.. (4.92)
Las permeabilidades relativas son determinadas en laboratorio y son caractersticas de

la roca reservorio que esta saturada con fluidos. No es una buena prctica usar las
permeabilidades relativas obtenidas en un reservorio para predecir el comportamiento de
otro reservorio similar.
Usualmente, las curvas de permeabilidades relativas son presentadas como funcin de

las saturaciones de agua Sw, como se muestra en la figura 4.14. Cuando la saturacin de
agua Sw es la saturacin de agua connata, Swc no habra un flujo de agua libre, y por
consiguiente, la permeabilidad efectiva Kw, debera ser igual a cero. Similarmente cuando
la saturacin de gas empieza con saturacin residual de gas Sgr,. no existe flujo y la
permeabilidad efectiva debera a ser igual a cero.
Figura 4.14 Efectos de las Permeabilidades Relativas en produccin de agua
Luego, en un reservorio de gas, en un sistema bifsico la ecuacin para el gas y agua es:
k rg kh(Pr Pwf ) gas

qg = Ec. 4.93
(
141.2 gg ln(re / rw) + S ' )
k rw kh(Pr Pwf ) w
qw = Ec. 4.94
(
141.2ww ln(re / rw) + S ' )
4.6.3 Relacin del ndice de productividad para un pozo horizontal de gas
La ecuacin 4.95, es anloga a la de darcys, para estado sostenido y pseudo estado

sostenido, la cual puede ser escrita para un pozo horizontal, permitiendo efectos de
turbulencia, el comportamiento de la entrada relaciona para un pozo horizontal en un
reservorio de gas son los siguientes:
Para estado sostenido:
k H h( p e2 p wf2 )
q= Ec.(4 .95)
a + a 2 (L / 2 ) 2 I h
I ani h
1424ZT ln + ani
ln + Dq
L/2 L rw (I ani + 1)

para pseudo-estado sostenido:

_ 2
k H h( p e p wf2 ) Ec.(4.96)
q=
a + a 2 (L / 2 ) 2 I h I ani h
3
1424 ZT ln + ani
ln + Dq
L/2 L rw (I ani + 1) 4

En ambos casos, el efecto de dao puede ser adicionado dentro del segundo corchete
del denominador de las ecuaciones Eqs. (4.95) y (4.96). As mismo de la aproximacin
_2
p p wf2 / Z la diferencia real de la pseudo presin puede usarse, en la ecuacin

anteriormente descripta. Los efectos de Turbulencia en un horizontal puede ser
despreciado, debido a que Dq no tiene influencia en la multiplicacin del Iani*h/L. As el
coeficiente de turbulencia puede ser desechado, ecepto cuando el pozo horizontal esta
parcialmente abierto o comunicado.
4.6.4. Efecto de dao en un pozo horizontal
En la Ecuacin. 4.90, el efecto de dao S'eq, fue adicionado a la ecuacin de caudal

del pozo horizontal. Frick y Economides, desarrollaron la ecuacin para el efecto de dao
que reflejan el dao cerca del pozo horizontal. La figura 4.14 describe la forma de dao a
lo largo y normal de un pozo horizontal.
La anisotropa de la permeabilidad debera generar una forma elptica normal al pozo.

La forma de dao depende de la anisotropa de la permeabilidad. I anillo dada por la Ec.
(4.88), la simulacin grfica de la respuesta est mostrada en la Fig 4.14, para tres
diferentes valores de I ani.
El tiempo de exposicin durante la perforacin y la completacin del pozo,

producira un cono elptico truncado con la base ms grande cerca de la seccin vertical del
pozo. Tambin, durante la produccin el perfil de presin en el pozo implicara una alta
gradiente normal en la trayectoria del pozo cerca de la seccin vertical. Por consiguiente, el
dao inducido por produccin sera elptico. La geometra de la forma de dao resulta de un
efecto de dao anlogo a Hawkins frmula para un pozo vertical:
:
k 1 4 aH2 max aH max

Seq' = 1 ln 2 + + 1 Ec.(4.97)
ks I ani + 1 3 rw rw
Figura 4.15 Distribucin de los daos y los efectos de las permeabilidades
Donde a Hmax es el eje horizontal ms grande (cerca de la seccin vertical) de dao

del cono (fig. 4.15). Ecuacin (4.97), no asume dao al final del pozo. El efecto de dao
puede ser adicionado a el denominador de un pozo horizontal de la ecuacin de caudal de
un pozo horizontal pero deber ser multiplicado por h* Iani / L, llamados relaciones de
anisotropa, como muestra en Eq. (4.90).
La estimulacin de la matriz de un pozo horizontal es ms amplia que el proceso de

un pozo vertical. Primero, por que los volmenes requeridos son ms grandes, como se
esperara. Segundo, la colocacin de fluidos para la estimulacin a lo largo del pozo es
mucho ms difcil. Por consiguiente, se han desarrollado tcnicas para manejar ambas
dificultades. stos incluyen mtodos mecnicos de estimulacin con la distribucin de
fluido (coild tubing, zona de aislacin levantando el collar de dao alrededor de la zona
estimulada) o no completada del pozo.
1 k a2 a 1 k k a2 a
Seq' = 1 ln s H 2max + s H max + 1 + ln i H 2max + i H max + 1
2 ks 2 k k r
rw rw i s w rw
k 3 Ec.(4.98)
1 ln (I ani + 1)
ki 4
4.7 Factores que afectan la curva del comportamiento del pozo con el tiempo
Como resultado del agotamiento del reservorio, con el transcurso del tiempo, existe
una cada en la produccin de los pozos debido al comportamiento o eficiencia del
reservorio que es afectado por algunos parmetros sujetos al cambio, cuyos parmetros son:
4.7.1Coeficiente C y exponente n
En reservorios de alta permeabilidad donde el flujo es de gas, estabiliza rpidamente,

el valor de C no tiene una variacin significativa con el tiempo, esto implica que el valor
del caudal mximo (AOF) obtenido de la curva de IPR permanece relativamente invariable
durante la vida del reservorio. De cualquier manera, es sabido que el caudal y la presin del
reservorio varia con el tiempo afectando el comportamiento del coeficiente, esto se acenta
ms en un reservorio de baja permeabilidad donde el caudal de produccin de gas durante
periodos de flujo relativamente cortos disminuye con el tiempo a una presin fija de
cabeza.
El coeficiente C contiene varios parmetros sujetos a cambios, a medida que la

presin del reservorio declina antes el agotamiento de la cada, por esta razn es necesario
recalcular el coeficiente C cada cierto perodo bajo nuevas condiciones de reservorio.
703 x 10 6 k h
C= Mpcd/psi^2 Ec.(4.99)
0.472 re
T Z ln + s
rw
Como ya se dijo anteriormente, n se encuentra generalmente entre 0.5 y 1, para bajas

permeabilidades en las pruebas de pozos el valor n aproximado a 1 y para altas
permeabilidades el valor de n tiende a 0.5, esto es evidente porque la velocidad del caudal
de flujo de gas es mayor en las formaciones con alta permeabilidad y la acumulacin de
lquidos en paredes del pozo durante la prueba de contrapresin dan como resultado
pendientes menores a 0.5 y exponente aparentemente mayores a 1.
4.7.2 Permeabilidad del gas
El nico factor que tiene un efecto apreciable es la permeabilidad del gas, kg ,a la

saturacin de lquido en el reservorio, como consecuencia de una condensacin retrgrada,
por la disminucin de presin de reservorio por debajo de la presin de roco formando una
fase lquida de condensado o por la presencia de agua en el reservorio, esto afecta
apreciablemente en la permeabilidad y saturacin del gas reduciendo cuantitativamente el
valor de ambos parmetros.
En reservorios de gas seco la permeabilidad permanece constante, pero, el grado

de turbulencia se incrementa debido a las altas velocidades requeridas para mantener un
caudal de flujo de masa constante.
4.7.3 Espesor de la formacin
Durante la vida de un reservorio, el espesor de la formacin es considerado invariable

en los casos que estemos frente a reservorios volumtricos sin empujes de agua, pero
cuando se tiene reservorios con un fuerte empuje del acufero el espesor del reservorio
disminuye dependiendo del avance del nivel acufero, tambin puede variar cuando los
intervalos de completacin son modificados por la necesidad de perforar intervalos mucho
ms grandes. En estas circunstancias el pozo debe ser probado para determinar las nuevas
condiciones.
4.7.4.- Viscosidad del gas y el factor de compresibilidad
Estos trminos son dependientes de la presin medida del reservorio. Los cambios de
viscosidad y del factor de compresibilidad del gas afectan los coeficientes de C y A de las
ecuaciones 4.21 y 4.67.
Las siguientes ecuaciones muestran la relacin que existe entre los valores actuales y
futuros de C y A mantenindose el radio de drene, el efecto dao y el espesor de la
formacin son constantes.
( Z )P
CF = CP Ec.(4.100)
( Z )F
( Z )F
AF = A P Ec.(4.101)
( Z )P
Donde P y F indican las condiciones al tiempo presente y tiempo futuro

respectivamente.
4.7.5 Radio de pozo y radio de drene
El radio del pozo es constante y slo el radio efectivo del pozo varia cuando existe
estimulacin, esto es notado en el resultado cuantitativo del factor dao.
El radio de drene depende del espaciamiento del pozo y puede ser considerado
constante una vez alcanzado un flujo estable.
4.7.6 Factor dao
Cuando el pozo es expuesto a una fractura o tratamiento cido, el factor dao o efecto
superficial de dao, varia incidiendo en los valores de C y n o los coeficientes de A y B
siguiendo el mtodo de anlisis que se utiliza debiendo recalcularlos, para lo cual se tiene
que realizar nuevas pruebas de pozo y evaluar los nuevos valores de los coeficientes
afectados.
4.8 Cada de presin a travs de las perforaciones
Hacemos una consideracin en la evaluacin prctica de las perforaciones hechas por

el Dr. Harry McCleod, en una de sus artculos publicados. Esto muestra que la
compactacin de la zona estudiada ocurre, alrededor de la perforacin en condiciones
normales. Puede ser que este problema difiera entre los valores de rea de los pozos cuando
estos estn comunicados con una formacin no consolidada, y son de inters en un rea
abierta al flujo. Para formaciones impermeables no solamente interesa un rea abierta de
flujo, sino tambin la longitud de la perforacin. Ambos tienen un efecto sobre el caudal
dentro del borde de la perforacin.
La Figura 4.15 muestra una tpica perforacin y la nomenclatura que se utiliza en este
anlisis. El orden para analizar el efecto de estas perforaciones es la capacidad de flujo,
utilizando diferentes trminos asumidos teniendo como base los trabajos de numerosos
autores, referencia en la Figura 4.16.
La Figura 4.16 muestra que para efectuar una perforacin de 90, esta debe ser
analizada como una seccin mucho ms pequea en las paredes del pozo. Adems, se ha
asumido en este anlisis que esta es una zona no daada. Varias otras suposiciones son
hechas tales como:
La permeabilidad de la zona compactada es:
a) 10% de permeabilidad de formacin si es perforado sobre balanceado.
b) 40% de permeabilidad de formacin si es perforado desbalanceado.
El espesor de la zona daada es pulgadas.
La pequea seccin de la pared del pozo puede ser analizada como un

reservorio infinito, esto es, si p wfs permanece constante en el borde de la zona
compactada, eliminando as mismo l - de la ley de Darcy para un lmite exterior
cerrado.
La ecuacin presentada anteriormente por Jones, Blount y Glaze puede

ser utilizada para evaluar las prdidas de presin a travs de las perforaciones.
Figura 4.15 Tubo perforado tpico
Figura 4.16 Tubo perforado 90

La ecuacin de Jones, Blount y Glaze, es modificada como sigue:
5.8.1 Cadas de presin a travs de las perforaciones
2 2
p wfs p wf = Aq + Bq Ec.(4.102)
3.16 x10 12 g TZ (1 r p 1 rc ) 2
p 2
p 2
= q .....
wfs wf
L 2
p
Ec.(4.103)
1.424 x10 3 TZ (ln rc r p )
+ q
k L
p p
Donde:
3.16 x10 12 g TZ (1 r p 1 rc )
B= Ec.(4.104)
Lp2

(
1.424 x10 3 TZ ln rc r p ) Ec.(4.105)
A=
k pLp

q = Caudal de flujo / peroracin (q/perforacin), Mcfd

= Factor de turbulencia.
2.33x1010
=
k 1p.201
g = Gravedad del gas (sin dimensin)
T = Temperatura del reservorio, R (F+460)
Z = Factor de compresibilidad (sin dimensin)
rc = Radio de la zona compactada, ft
rc = r p + 0.5 in
rp = Radio de perforacin., ft
Lp = Longitud de la perforacin, ft
= Viscosidad, cp
kp = Permeabilidad de la zona compactada, md
= 0.1 K (disparo sobre balanceado)
= 0.4 K (disparo desbalanceado)
4.8 Perforaciones de Pozos y Efectos de Dao
La perforacin moderna de pozos es realizada con caones que estn adheridos, ya

sea al cable de acero o a una tubera de serpentin. La Figura 4.17, muestra un dibujo
esquemtico de un sistema de can con las cargas de conformacin ordenadas en un
patrn helicoidal. Este patrn permite una buena densidad de perforacin con una pequea
puesta de fase de Angulo entre las perforaciones continuas.
Figura 4.17 Sistema Esquemtico del can

La sarta de perforacin contiene un cabezal de cable, un dispositivo de

correlacin, un dispositivo posicionador y los caones de perforacin. El cabezal de cable
conecta a la sarta de cable y al mismo tiempo proporciona un punto dbil en el cual se
conecta el cable si surgen problemas de aprisionamiento. El dispositivo de correlacin es
usado para identificar la posicin exacta con un registro de correlacin previamente
dispuesto localizando frecuentemente los collares de las tuberas de revestimiento. El
dispositivo posicionador orienta los disparos hacia la tubera de revestimiento para una
optima geometra de perforacin. Los caones de perforacin son cargados con cargas de
conformacin, que consisten en la caja, el explosivo, y el forro, tal como se muestra en la
Figura 4.17, la corriente elctrica inicia una onda explosiva, las fases de la detonacin son
mostradas en la Figura 4.17. las perforaciones con un dimetro entre 0.25 y 0.4 pulgadas
nos da una penetracin que esta entre los 6 a 12 pulgadas segn diseo de los fabricantes.
Figura 4.18 Proceso de detonacin de cargas

La perforacin es generalmente realizada bajo condiciones de equilibrio, luego se

debe reducir la presin de la columna para que exista una reaccin inmediata del reservorio
acarreando los deshecho de la perforacin efectuada, las dimensiones de las balas y el
enfasamiento de las perforaciones tienen un papel controlante en el comportamiento del
pozo.
4.8.1 Clculo del efecto de Dao de la Perforacin
Karakas y Tariq (1988), han presentado una solucin semianaltica para el clculo
del efecto de dao en las perforaciones , ellos dividen en tres componentes los cuales son:
el efecto de circulacin plano SH , el efecto de convergencia vertical Sv, y el efecto de
agujero Swb , el efecto del dao total de la perforacin es:
Sp = SH + SV + SWb Ec. (4.106)
La figura 4.19 da todas las variables pertinentes para el clculo de la capa externa de
perforacin. Estas incluyen el radio del pozo, rW, el radio de perforacin, rperfo, la longitud
de la perforacin, ngulo de enfasamiento de la perforacin y la distancia entre las
perforaciones, hperf, que es exactamente inversamente proporcional a la densidad de
perforacin. la Figura 4.19 (nos muestra las variables y fase de la perforacin) con cual
delineamos el mtodo para estimar los componentes individuales de la capa externa de la
perforacin.
Figura 4.19 Variables y Fases de las Perforaciones

Variable del pozo para el clculo de la capa externa de perforacin

Clculo de SH:
rW
S H = ln Ec. (4. 107)
rW (0 )
l perf
para = 0
rW (0 ) = 4
a 0 (rW + l perf ) Ec. (4.108)
para = 0
La constante a0 depende del enfasamiento de la perforacin y puede ser obtenida de

la tabla 4.3. Este efecto de dao es negativo (excepto para 0 = 0), pero su contribucin total
es generalmente pequea.
TABLA 4.3
Constante para el Clculo de efecto de dao de perforacin

Enfasamiento
de la a0 a1 a3 b0 b1 c0 c1
perforacin
00(3600) 0.250 -2.091 0.0453 5.1313 1.8672 1.6E-1 2.675
180 0.500 -.0250 0.0943 3.0373 1.7115 2.6E-2 4.532
120 0.648 -2.018 0.0643 1.6136 1.7770 6.6E-3 5.320
90 0.726 -1.905 0.1038 1.5674 1.6935 1.9E-4 6.155
60 0.813 -1.898 1.1023 1.3654 1.6490 3.0E-5 7.509
45 0.860 -1.788 1.2398 1.1915 1.6392 4.6E-6 8.791
a
De Karakas y Tariq, 1988
El clculo Sv esta en funcin de dos variables adimensinales que deben ser
calculadas de la siguiente manera:
h perf kH
hD = a + Ec. (4.109)

l perf kV

Donde KH y KV, son las permeabilidades horizontal y vertical, respectivamente
r perf k
rD = 1 + H Ec. (4.110)

2h perf kV

la pseudo-capa externa vertical es entonces:

S V = 10 a .hD
b 1
.rD
b Ec. (4.111)
con
a = a1 Logr D + a 2
y
b = b1 rD + b2
Las constantes a1, a2 y b2 son tambin funciones del enfasamiento de perforacin y

pueden ser obtenidas. En efecto, el dao vertical SV, es potencialmente el ms grande
contribuidor del dao total Sp; para pequeas densidades de perforacin. Si el hperf, es muy
grande, el SV puede ser muy grande. Para el clculo del SWb una cantidad adimensinal es
calculada primeramente.
rw
S wD = Ec.(4.112)
I perf + rw
Luego
S wb = c1* e c 2*rwD Ec.(4.113)
Las Constante C1 y C2 pueden ser obtenidas de la Tabla 4.3
4.8.2 Daos Cerca al Pozo y Perforaciones
Karakas y Tariq 1988, tambin han demostrado que los daos y las perforaciones
pueden ser caracterizadas por un efecto de dao compuesto.
Sd p = [ K
KS
] r k
1 * Ln s + S p + S p = (S d )o + 1 * S p Ec. (4.114)
rw ks
Si las perforaciones terminan dentro de la zona de los daos (Iperf.< rs). En la

ecuacin 4.113 (Sd)o, es el efecto de dao equivalente de agujero abierto dado por la
frmula de Hawkins. Si las perforaciones terminan fuera de la zona de dao, entonces:
Sd p = S p Ec. (4.115)
Donde Sp es evaluada en una longitud de perforacin modificada I perf. y un radio

modificado rw, estos son:
K
I perf = I perf 1 s * rs Ec.(4.115)
K
K
rw = rw 1 s * rs Ec.(4.116)
K
Estas variables son usadas en las ecuaciones 4.107 y 4.112, para los efectos de dao
que contribuyen al efecto de dao compuesto en la ecuacin 4.115.
The Technology of Artificial Lift Methods Kermit E. Brown, Volume 4

Production Optimization of Oil and Gas Wells by Nodal Systems Analysis.
Production Optimization, Using Nodal Analysis H. Dale Beggs, 1991.
Applied Reservoir Engineering Charles R. Smith, G. M. Tracy, R. Lance

Farrar, Volume1, 1992.
Petroleum Production System Michael J. Economides, A. Daniel Hill,

Christine Ehlig-Economides, 1994.
Petroleum Engineering Handbook Society Of Petroleum Engineers, third

printing, feb. 1992
Ingeniera Aplicada de Yacimientos Petrolferos B. C. Craff y M. F. Hawkins,

Jr, 1997
Anlisis de la prdida de Presin de la columna de Produccin ______________________________ 139
5 Prdidas de Presin en Lneas de Produccin
5.1 Introduccin
La mayora de los componentes de un sistema de distribucin (Figura. 5.1) son

atravesados por el fluido de produccin en condiciones de flujo multifsico, es decir, gas
condensado y agua, por esta razn resulta imprescindible adelantarse en el estudio de los
mecanismos y principios del fluido multifsico, para estimar adecuadamente las prdidas de
carga que se producen en la tubera o lnea de conduccin. El flujo en tuberas se define
como el movimiento de gas libre, mezcla de fluidos o una combinacin de algn modelo de
flujo en tuberas sobre diferentes condiciones de operacin. El gas proveniente del medio
poroso pasa a la etapa de transporte por tubera, ya sea con movimiento vertical o
direccional, hasta la superficie donde cambia a un sentido horizontal o inclinado hasta el
separador. (Figura 5.1).
Figura 5.1 Direcciones de flujo en la tubera

Si bien la bibliografa acostumbra dividir al flujo multifsico en flujo vertical y flujo

horizontal inclinado, los mecanismos y principios son los mismos, a tal punto que existen
numerosas correlaciones que son aplicables al flujo vertical, horizontal o inclinado.
Durante ese proceso, el fluido puede pasar por algunas restricciones (vlvulas, choques)
antes de ingresar a la lnea principal que conecta con el separador, los mismos que
interponen una resistencia en el canal de flujo provocando una perdida adicional de la
energa inicial disponible. Lo expuesto anteriormente, se puede visualizar en la Figura
6.1(Esquema del choque) del Capitulo 6 del Anlisis de Flujo de Gas A travs de los
Choques; que permite apreciar las posibles zonas de mayor cada de presin en un sistema
de produccin.
Todas las perdidas de carga son funcin del caudal de produccin y de las
propiedades o caracterstica de los componentes del sistema. En el caso de un flujo en fase
simple, sea este lquido o gas, la prdida de presin puede ser calculada fcilmente, sin
embargo durante la produccin de un pozo se tiene un fluido multifsico lo que complica el
clculo de dichas prdidas de presin, de los tres trminos que contribuyen a la prdida de
carga total del sistema (Figura. 5.2.). Son los Componentes de Friccin, Aceleracin y
Elevacin
5.2 Ecuacin de energa
La base terica para la mayora de las ecuaciones de flujo de fluidos es la ecuacin

general de energa, la cual es una expresin del balance o conservacin de energa entre dos
puntos en un sistema. Considerando un sistema de estado estable, el balance de energa se
puede escribir como:
mu 2 mgZ1 mu22 mgZ 2
U1 + p1V1 + 1 + + Q w = U 2 + p2V2 + + Ec. (5.1)
2 gc gc 2 gc gc
donde:
U1 = Energa interna
PV = Energa de expansin o compresin.
2
mu
= Energa cintica.
2 gc
mgZ
= Energa potencial.
gc
Q = Transferencia de calor.
w = Trabajo desenvuelto por el fluido.
Dividiendo la ecuacin 5.1 por m para obtener un balance de energa por unidad de
masa y escribiendo la ecuacin resultante en forma diferencial.
p u du g
dU + d + + dZ + dQ dw = 0 Ec. (5.2)
gc gc
La ecuacin 5.2 se puede convertir en un balance de energa mecnica usando la

siguiente relacin termodinmica:
dp
dh = Tds + Ec. (5.3)

y
p
dU = dh d Ec. (5.4)

o
dp p
dU = Tds + d Ec. (5.5)

Donde:
h = Entalpa.
s = Entropa.
T = Temperatura.
Substituyendo la ecuacin 5.5 en la ecuacin 5.2, obtenemos:

dp u du g
Tds + + + dZ + dQ dw = 0 Ec. (5.6)
gc gc
Para un proceso irreversible, los estados de desigualdad de Clausius son:
dQ
ds Ec. (5.7)
T
Tds=-dQ + d(lw) Ec. (5.8)
Donde
lw = Prdida de trabajo, debido a la irreversibilidad.
Substituyendo la ecuacin 5.8 en la ecuacin 5.6 tenemos:
dp u du g
+ + dZ + d (lw ) dw = 0 Ec. (5.9)
gc gc
Si ningn trabajo es hecho en o por el fluido, dw = 0, obtenemos:
dp u du g
+ + dZ + d (lw ) = 0 Ec. (5.10)
gc gc
Considere una tubera inclinada por donde escurra un fluido en rgimen permanente,
conforme se muestra en la Figura 5.2, ya que dZ 0 dL sen, la ecuacin de energa se
vuelve:
dp u du g
+ + dL sen + d (lw ) = 0 Ec. (5.11)
gc gc
Figura 5.2 Flujo de fluido en tubera inclinada
La ecuacin 5.11 se puede escribir en trminos de gradientes de presin,

multiplicando la ecuacin por / dL:
dp u du g d (l )
+ + sen + w = 0 Ec. (5.12)
dL g c dL g c dL
Considerando la cada de presin positiva en la direccin del flujo, la ecuacin 5.12 puede
ser escrita como:
dp g dp u du
= + sen + + Ec. (5.13)
dL gc dL f g c dL
Donde el gradiente de presin debido al esfuerzo viscoso o las prdidas por atraccin
son expresados como:
dp d (l )
= w Ec. (5.14)
dL f dL
Podemos definir un factor de atraccin:
w 2 w g c
f = = Ec. (5.15)
u / 2 gc
2
u2
Donde
w = Esfuerzo cortante.
u 2
= Energa cintica por unidad de volumen.
2 gc
La ecuacin 5.15; define la relacin entre el esfuerzo cortante sobre la pared, y la

energa cintica por unidad de volumen, refleja la importancia relativa del esfuerzo cortante
en la pared sobre la prdida total. El esfuerzo de corte sobre la pared puede ser vlido
fcilmente de un balance de fuerzas entre las fuerzas de presin y las fuerzas de viscosidad.
Figura 5.3.
dp d 2
p
1 p dL = w ( d )dL Ec. (5.16)
dL 4
1

Figura 5.3 Equilibrio de fuerzas
Re escribiendo la ecuacin 6.16,

d dp
w = Ec. (5.17)
4 dL f
Substituyendo la ecuacin 5.17 en la ecuacin 5.15 y resolviendo para el gradiente

de presin debido a la atraccin, obtenemos.
dp 2 f 'u2
= Ec. (5.18)
dL f gc d
Esta es la ecuacin de Fanning, y f es llamado el factor de friccin Fanning. En

trminos del factor de friccin Darcy Weisbach o Moody, f 0 4f, la ecuacin 5.18 toma
la siguiente forma:
dp f u 2
= Ec. (5.19)
dL f 2 g c d
Donde
u = Velocidad de flujo, ft / seg.
gc = Factor de conversin gravitacional 32.17 lbmft / lbf seg2
d = dimetro de la tubera, ft
f = Factor de friccin Moody.
5.3.- Gradiente de presin total
El gradiente de presin total puede ser considerado como la composicin de los tres
dp
componentes del sistema; El componente por el cambio de elevacin ,
dL elevf
dp
componente por la prdida por friccin , y el ltimo componente es por los
dL friccion
dp
cambios en la energa cintica o aceleracin . Rescribiendo la ecuacin en
dL acele
funcin al gradiente de presin total se tiene:
dp dp dp dp
= + + Ec. (5.20)
dL dL dL f dL
dp
, Es la prdida de presin debido al peso de la columna de fluido
dL elevf
hidrosttico esta es cero si se trata de un flujo horizontal. Es la prdida de
presin debido al cambio de la energa potencial entre las posiciones 1 y 2 de

la Figura 5.2 la expresin viene dada :
dp g . * sen
= Ec. (5.21)
dL elev gc
dp
, es la prdida de presin por friccin, que se origina por el
dL friccion
movimiento del fluido contra las paredes de la tubera de produccin que estn
en funcin del dimetro de la tubera y a las propiedades del fluido, como se
muestra en la ecuacin siguiente:
dp f u 2
= Ec. (5.22)
dL f 2 g c d
dp
, es la cada de presin que resulta de un cambio en la velocidad del
dL acele
fluido entre las posiciones 1 y 2 de la Figura 5.2 la ecuacin para la cada de
presin por aceleracin es la siguiente:
dp vdv
= Ec. (5.23)
dL accl g c dL
La ecuacin 5.13 se aplica perfectamente al flujo multifsico, debiendo modificarse

ligeramente algunos trminos:
La densidad, es la densidad de la mezcla Gas-Lquido (m)

La velocidad, es tambin la velocidad de la mezcla (vm)
El factor de friccin, es ahora factor de friccin bisfico (f m)
Por lo tanto, la ecuacin general para flujo multifsico se convierte en:
dp m um dum g fmmvm
= + m sen + Ec. (5.24)
dL g c dL gc 2 gcD
5.4 Nmero de Reynolds
Reynolds, aplic el anlisis adimensional al movimiento del fluido en relacin a la

viscosidad para dos flujos, definiendo los parmetros adimensionales. El nmero de
Reynolds, parmetro no dimensional, utilizado para distinguir entre un flujo laminar y
turbulento, se define como la relacin entre el momento de fuerzas del fluido y las fuerzas
viscosas o de corte.
du
N Re = Ec. (5.25)

Donde
d = Dimetro de la tubera.
= Densidad del fluido, lbm / ft3.
u = Velocidad del fluido, ft / seg.
= Viscosidad del fluido, lbm / ft seg.
La viscosidad dinmica del fluido , frecuentemente indicada en cp, puede ser

convertida en lbm / ft seg utilizando el factor de conversin de 1 cp = 6.7197x10-4 lbm / ft
seg. As mismo:
du
N Re = 1488 Ec. (5.26)

En trminos de caudal de gas o nmero de Reynolds para condiciones bases de
temperatura y presin se puede escribir como:
20 q g
N Re = Ec. (5.27)
d
Donde q esta en Mscfd, esta en cp, y d esta en pulgadas. El rgimen de flujo es

relacionado al nmero de Reynolds como sigue.
Tabla N. 5.1
TIPO DE FLUJO NRe
Laminar <2000
Crtico 2000 3000
Transicin 3000 4000
Turbulento >4000
5.5 Rugosidad Relativa
El factor de friccin a travs de una tubera es afectado por la rugosidad en la pared

de la tubera. En tanto, la rugosidad de la tubera no es fcilmente ni directamente medible.
La rugosidad es una funcin del material tubular, el mtodo de manufactura y del medio
ambiente al que es expuesto, por tanto, no es uniforme a lo largo de la longitud de la
tubera.
El efecto de rugosidad en tanto, no es debido a las dimensiones absolutas definida

como la medida de la altura y el dimetro existentes en una zona de mayor protuberancia y
una distribucin relativamente uniforme, donde los grados proporcionan el mismo
comportamiento del gradiente de presin. Ms as, las dimensiones relativas del dimetro
interno de la tubera, as mismo, la rugosidad relativa es la relacin de la rugosidad absoluta
y el dimetro interno de la tubera.
e
rugosidad relativa = Ec. (5.28)
d
Figura 5.4 Factores de friccin para cualquier tipo de tubera comercial.
Donde
e = Rugosidad absoluta, ft o in.
d = Dimetro interno, ft o in.
Como la rugosidad absoluta no es directamente medible, la seleccin de la

rugosidad de la tubera es difcil, especialmente si no se tienen los datos disponibles como
el gradiente de presin, factor de friccin y el nmero de Reynolds para hacer uso del
diagrama de Moody, (ver Figura 5.4), en este caso el valor de 0.0006 es recomendado para
tuberas y lneas. La Tabla 5.2 muestra valores tpicos de rugosidad absoluta utilizados en
problemas de flujo de gas natural y la Figura 5.5 presentado por Moody; permite obtener la
rugosidad relativa en funcin del dimetro, material de tubera y el factor de friccin para
un flujo completamente turbulento.
Figura 5.5 Valores de rugosidad relativa y factor de friccin.

TABLA N. 5.2
Valores de Rugosidad absoluta
MATERIAL DE TUBERA e (in)
Vidrio bronce 0.00006
Tubera de Pozo 0.0006
Aluminio 0.0002
Plstico 0.0002 0.0003
Acero comercial 0.0018
Fierro comn 0.0048
Fierro Galvanizado 0.006
Fierro Dulce 0.0102
Revestimiento de Cemento 0.012 0.12
5.6 Determinacin del factor de Friccin
Los valores de flujo normalmente se encuentran entre dos extremos: flujo laminar y
flujo turbulento, dentro de estos valores se distinguen cuatro zonas: laminar, crtico,
transicin y totalmente turbulento, para cada zona las ecuaciones varan en funcin del
nmero de Reynolds y la rugosidad (ver Figura 5.4).
5.6.1 Flujo laminar de fase simple
Existe flujo laminar cuando se presenta un movimiento estacionario permanente en

cada punto de la trayectoria del fluido, que dice, que las lneas de corriente se deslizan en
forma de capas con velocidades suficientemente bajas sin causar remolinos. El factor de
friccin para flujo laminar puede ser determinado analticamente. La ecuacin de Hagen-
Poiseuille para flujo laminar es:
dp 32 u
= Ec. (5.29)
dL f gc d 2
El factor de friccin de Moody para flujo laminar se puede determinar combinando

las ecuaciones 5.22 y 5.29 para luego combinar con la ecuacin 5.23, las expresiones son:
fu 2 32 u
= Ec. (5.30)
2 gc d gc d 2
y
64
f = Ec. (5.31)
N Re
En su forma equivalente para el factor de friccin de Fanning:

16
f '= Ec. (5.32)
N Re
5.6.2 Flujo turbulento de fase simple
Existe un movimiento turbulento cuando la velocidad media lineal excede la

velocidad crtica y las partculas siguen trayectoria errtica. Estudios experimentales de
flujo turbulento, han mostrado que el perfil de la velocidad y el gradiente de presin son
muy sensibles a las caractersticas de la pared de tubera de produccin y lneas de
surgencia, caractersticas que se pueden clasificar como: tuberas lisas y tuberas rugosas.
5.6.2.1 Tuberas lisas
Las ecuaciones que se presentan son vlidas para valores especficos de nmero de
Reynolds, como la propuesta por Drew, Koo y McAdams en 1930, utilizado para intervalos
de 3x103 < NRe < 3x106.
f = 0.0052 + 0.5NRe-0.32 Ec. (5.33)
Y la ecuacin presentada por Blasius, para nmero de Reynolds mayores a 105 en tuberas
lisas:
f = 0.31NRe-0.25 Ec. (5.34)
5.6.2.2 Tuberas rugosas
En flujo turbulento, la rugosidad tiene un efecto determinante en el factor de friccin

por consiguiente, la gradiente de presin debido a su dependencia con la rugosidad relativa
y nmero de Reynolds. Colebrook y While en 1939, propusieron una ecuacin aplicable a
tuberas lisas como tubera de flujo de transicin y totalmente rugosas en las zonas de flujo
turbulento, la dificultad radica en que la ecuacin no es lineal su solucin requiere de un
procedimiento iterativo, dicha ecuacin es:
1 2e 18.7
= 1.74 2 log + Ec. (5.35)
f d N f e
Re
La solucin iterativa consiste en estimar un valor inicial de fe, se aconseja utilizar el

determinado a travs de la ecuacin de Drew, Koo y McAdams como primer parmetro.
Con este dato se clcula el nuevo factor de friccin f , que substituira al primer valor, en el
caso de que numricamente f y fe estn muy lejos del valor de tolerancia considerado
normalmente 0.001.
Para valores considerados de rugosidad relativa entre 10-6 y 10-2 , nmero de

Reynolds entre 5x103 y 108, se puede utilizar la ecuacin de Jain, 1976.
1 e 21.25
= 1.44 2 log + 0.9 Ec. (5.36)
f d N Re
La ecuacin 5.35 da un error de +/- 1% cuando se compara como la ecuacin de Colebrook,

la ventaja sobre esta ltima es la obtencin directa del factor de friccin.
5.7 Flujo de fase simple
Ahora que ya se ha presentado ecuaciones y procedimientos para validar el factor de

friccin en flujo de simple fase, la ecuacin de gradiente de presin derivada anteriormente
se puede desenvolver ms extensamente. Combinando las ecuaciones 5.13 y 5.22, la
ecuacin da el gradiente de presin, y aplicable para cualquier fluido en cualquier ngulo
de inclinacin de tubera, esta llega a seriii.
dp g f u 2 u du
= sen + + Ec. (5.37)
dL g c 2 gc d g c dL
Donde el factor de friccin, f, es una funcin del nmero de Reynolds y la rugosidad
de la tubera. Esta relacin se muestra en el diagrama de Moody (Figura 5.4). La gradiente
de la presin total puede ser considerada para ser compuesta de tres componentes distintos:
dp dp dp dp
= + + Ec. (5.38)
dL dL el dL f dL ace
Donde:
(dp / dL)el = g sen / gc , es el componente debido a la energa potencial o a al cambio de

elevacin. Es tambin referido como el comportamiento hidrosttico, ya que es el nico
componente el cual se aplicar en condiciones estticas de flujo.
(dp / dL)f = f u2 / 2 gc d , es el componente debido a las prdidas por friccin.
(dp / dL)ace = u du / gc dL , es el componente debido al cambio de energa cintica.
La ecuacin 5.37 se aplica para cualquier fluido en estado estable, un flujo dimensional
para el cual f , , y u pueden ser definidos.
El cambio de elevacin al componente hidrosttico es cero, solamente para flujo

horizontal. Se aplica para fluido compresible o incomprensible, flujo pseudo estable y
transiente en las tuberas verticales o inclinadas. Para flujo descendente el seno del ngulo
es negativo, y la presin hidrosttica aumenta en direccin del flujo.
La prdida de friccin de los componentes se aplica para cualquier tipo de flujo en

cualquier ngulo de inclinacin de tubera. Esto siempre causa una cada de presin en
direccin de flujo. En flujo laminar las prdidas de friccin linealmente son proporcional a
la velocidad del fluido. En flujo turbulento, las prdidas de friccin son proporcionales a un
, donde 1.7 n 2.
El cambio de energa cintica o la aceleracin del componente es cero para rea

constante, flujo incomprensible. Para cualquier condicin de flujo en la cual ocurre un
cambio de velocidad flujo compresible, la cada de presin suceder en una cada de
presin en direccin que la velocidad aumenta.
5.8 Flujo de dos fases
Introduciendo una segunda fase dentro de una corriente de flujo complica el anlisis
de la ecuacin de gradiente de presin. El gradiente de presin es incrementado para la
misma masa de flujo, y el flujo podra desarrollar pulsaciones naturales. Los fluidos pueden
separarse debido a la diferencia de densidades y pueden fluir en velocidades diferentes en la
tubera. Una interfase de separacin puede existir entre la fase de lquido y gas.
Propiedades como la densidad, velocidad y la viscosidad, las cuales son relativamente

simples para fluidos individuales, se vuelven muy difciles para determinar. Las cadas de
presin para sistemas multifsico.iii Antes de modificar la ecuacin de gradiente de presin
para las condiciones de flujo de dos fases, se debe definir y evaluar ciertas variables nicas
para una mezcla de dos fases gas lquido.
5.8.1 Variables de flujo de dos fases
Para calcular el gradiente de presin se necesita los valores como densidad,

viscosidad, y en algunos casos la tensin superficial para las condiciones de flujo. Cuando
estas variables son calculadas para flujo de dos fases, se encuentran ciertas normas de
mezclas y definiciones nicamente para esta aplicacin.
Analizaremos algunas de las propiedades ms importantes, las cuales deben ser

estudiadas antes de adaptar las ecuaciones de gradiente de presin para las condiciones de
dos fases.
5.8.1.1 Escurrimiento de Lquido (Holdup), HL,
Es definido como la fraccin de un elemento de tubera que es ocupado por el

lquido en algn momento, esto es:
Volumen de lquido en un elemento de la tubera

HL =
Volumen de elemento de tubera
Evidentemente, si el volumen es muy pequeo, el escurrimiento del lquido ser cero

o uno. Siendo necesario determinar el escurrimiento del lquido, si se desea calcular tales
valores como densidad de la mezcla, velocidad a travs del gas y lquido, viscosidad
efectiva y transferencia de calor. En caso de fluctuaciones de flujos, tal como el flujo tapn,
el escurrimiento de lquido peridicamente cambia de punto a punto y es tomado como el
tiempo medio del valor.
El valor de escurrimiento del lquido vara desde cero, para flujo de gas de una sola
fase, a uno para flujo de lquido de una sola fase. El escurrimiento de lquido puede ser
medido experimentalmente por varios mtodos, tales como pruebas de resistividad o
capacidad, densitmetros nucleares, o por entrampar un segmento de corriente de flujo
entre vlvulas de echado rpido y midiendo el volumen de lquido apaado.
No se puede calcular analticamente un valor de escurrimiento de lquido. Se debe

determinar de correlaciones empricas en funcin de variables, tales como gas, propiedades
lquidas, flujo constante, y dimetro e inclinacin de tubera.
El Volumen relativo instantneo de lquido y gas algunas veces expresado en

trminos de fraccin de volumen ocupado por el gas, llamado altura del gas , Hg , o fraccin
al vaco. La altura de gas es expresada como:
Hg = 1 - H L
5.8.1.2 Suspensin de lquido, L
Llamado algunas veces lquido de entrada, es definido como la razn de volumen de

lquido en un elemento de tubera, el cual existir si el gas y el lquido viajan a la misma
velocidad (sin escurrimiento) dividido por el volumen de elemento de tubera. Esto se
puede calcular conociendo los caudales de flujo de gas y lquido in situ, utilizando:
qL
L = Ec. (5.39)
qL + q g
donde
qL = Suma de caudales de petrleo y agua in situ.
qg = Caudal de flujo de gas in situ.
La elevacin de gas no voltil o fraccin de gas es definido como:
qg
g = 1 L =
qL + q g
5.8.1.3 Densidad
Toda ecuacin de flujo de fluido requiere que un valor de densidad de fluido este
disponible. La densidad esta envuelta en la evaluacin de los cambios de energa debido a
la energa potencial y los cambios de energa cintica. Para calcular los cambios de
densidad con los cambios de presin y temperatura, se necesita tener una ecuacin de
estado para el lquido sobre consideraciones. Las ecuaciones de estado son fcilmente
disponibles para fluidos de fase simple.
Cuando dos lquidos inmiscibles como petrleo y el agua fluyen simultneamente, la

definicin de densidad se vuelve mas complicada. La densidad de una mezcla fluyente de
gas lquido es muy difcil de evaluar debido a la separacin gravitacional de las fases, y la
volatibilidad entre las mismas. La densidad de una mezcla de petrleo agua se puede
calcular de forma aproximada de la siguiente manera:
L = O fO + w f w Ec. (5.40)
donde
qO
fO = Ec. (5.41)
qO + qW
y
fw = 1 fO Ec. (5.42)
Para calcular la densidad de una mezcla de gas lquido, se necesita conocer el

escurrimiento de lquido utilizando tres ecuaciones para densidad de dos fases hecha por
varios autores de flujo de dos fases:
S = L H L + g H g Ec. (5.43)
n = L L + g g Ec. (5.44)
L 2L g 2g
K = + Ec. (5.45)
HL Hg
La ecuacin 5.43 es utilizada por la mayora de los autores para determinar el

gradiente de presin debido al cambio de presin. Algunas correlaciones estn basadas en
suposicin de no-volatilidad, y por lo tanto, los autores utilizan la ecuacin 5.44 para
densidad de dos fases. La ecuacin 5.45 es utilizada por algunos autores para obtener la
densidad de la mezcla utilizada para calcular los trminos de perdidas de friccin y nmero
de Reynolds.
5.8.1.4 Velocidad
Muchas correlaciones estn basadas en una variable llamada velocidad superficial. La

velocidad superficial de la fase de un lquido es definida como la velocidad en la cual esa
fase existira si este fluido pasa a travs de toda seccin transversal de tubera. La velocidad
de gas superficial es calculada por:
q
vsg = g Ec. (5.46)
A
El rea real por la cual los flujos de gas estn reducidos por la presencia de lquido al
A Hg.
Por tanto, la velocidad real es calculada por:
qg
vg = Ec. (5.47)
AH g
Donde A es el rea de tubera.
Las velocidades de lquido real y la superficial son calculadas de forma similar:

q
vSL = L Ec. (5.48)
A
qL
vL = Ec. (5.49)
AH L
Desde que Hg y HL son menores que uno, las velocidades reales son mayores que las
velocidades superficiales.
La velocidad de dos fases o mezcla es calculada en base a los caudales de flujo instantneo
total de la ecuacin:
q + qg
vm = L = vsL + vsg Ec. (5.50)
A
Las fases de gas y lquido viajan a velocidades diferentes en la tubera. Algunos

autores prefieren evaluar el grado de desprendimiento y as mismo el escurrimiento de
lquido para determinar la velocidad de deslizamiento Vs. La velocidad de deslizamiento es
definida como la diferencia entre las velocidades de gas real y liquido por:
vsg vsL
vs = v g vL = Ec. (5.51)
Hg HL
Utilizando las anteriores definiciones para varias velocidades, formas alternativas de

las ecuaciones para no-escurrimiento de lquido real son:
v
L = sL Ec. (5.52)
vm
HL = s m
[
v v + (vm vs ) + 4vs vsL
2
]
1/ 2
Ec. (5.53)
2vs
5.8.1.5 Viscosidad
La viscosidad de un fluido es utilizada en la determinacin del nmero de Reynolds,

as mismo como otros nmeros sin dimensiones utilizados como parmetros de correlacin.
El concepto de viscosidad de dos fases es discutida y tiene que ser definida en forma
diferente por varios autores. Las siguientes ecuaciones utilizadas por los autores para
calcular la viscosidad de dos fases de gas lquido son:
n = L L + g g Ec. (5.54)
Hg
s = LH g
L
Ec. (5.55)
n = L H L + g H g Ec. (5.56)
La viscosidad de una mezcla petrleo agua es generalmente calculada utilizando

las fracciones fluyentes de petrleo y agua en una mezcla con los factores de peso. La
ecuacin normalmente utilizada es:
n = L fO + w f w Ec. (5.57)
Esta ecuacin no es vlida si se formase una emulsin de agua petrleo. La

viscosidad de petrleo crudo, gas natural y agua puede ser estimada por correlaciones
empricas.
5.8.1.6. Tensin Superficial
La tensin interfacial depende de otras propiedades de fluido como la gravedad de

petrleo, gravedad del gas y gas disuelto. Entonces, si la fase lquida contiene petrleo y
agua, los mismos factores de peso son utilizados para calcular la densidad y la viscosidad.
Esto es:
L = O fO + w f w Ec. (5.58)
Donde
L = Tensin superficial de petrleo.
w = Tensin superficial de agua.
5.8.2. Modificacin de la ecuacin de gradiente de presin para flujo de dos fases
La ecuacin de gradiente de presin, se aplica para cualquier flujo de fluido en una lnea
inclinada dado un ngulo de horizontal, dado previamente como:
dp dp dp dp
= + + Ec. (5.59)
dL dL el dL f dL ace
5.8.3.- Modelo simplificado para predecir velocidad mnima del gas para remover lquido
del fondo y la velocidad erosional
Turner et al hicieron un anlisis de la velocidad mnima necesaria para evitar el

resbalamiento del lquido en el pozo y para mover la pelcula de lquido en la pared del
tubera. Este estudio gener un criterio para determinar la velocidad crtica usando el
modelo de agua y comparando con de las velocidades crticas del gas producido en el pozo
a temperatura en cabeza y presin en cabeza.
Las velocidades crticas encontradas por Turner para el modelo de cada resulta una
ecuacin simple de la velocidad hacia arriba, ascendente para evitar que las gotas ms
grandes caigan, lo cual es considerada como un esfuerzo de corte. Las ecuaciones de Turner
fueron desarrolladas de la siguiente manera:
Bajada (peso) Subida (velocidad)
g CdgApVt
( l g )Vol = Ec. (5.60)
gc 2 gc
Donde
g = aceleracin local de la gravedad 32.2 pie/seg^2

gc = constante gravitacional 32.2 lb-pie/lb-seg^2
Cd = Coeficiente de cada adimensional
g = densidad del gas lbs/pie^3
l = densidad del lquido lbs/pie^3
Vf = velocidad final pie/seg
Vol = volumen del lquido cado pie^3
Ap = rea de la tubera pie^2
Resolviendo
Vl = 6.55 * ( L g ) D Ec.(5.61)
Cdg
La relacin entre la cada, el dimetro, velocidad y tensin superficial es:
Vt ^ 2 gD
Nwe = Ec.(5.62)
.gc
Donde:
Nwe = nmero de weber
= tensin superficial lb/pie
Usando la relacin anterior, Tarner recomend que el nmero de Weber sea igual a 30 para
cadas grandes, modificando la ecuacin 5.61.
1.59 ^1 / 4( l g )^1 / 4
Vl = Ec. (5.63)
g ^1 / 2
Propiedades Agua Condensado
Tensin Superficial (Dyna/cm) 60 20

Temperatura Superf. ( oR ) 580 580
Densidad (lb/pc) 67 45
Gravedad del Gas 0,65
Remplazando los datos de las tablas tenemos:
5.3( l 0.00279 P)^1 / 4

Vg..agua = Ec.(5.64)
(0.00279. * P)^1 / 2
4.03( l 0.00279 P)^1 / 4

Vg..cond = Ec. (5.65)
(0.00279. * P)^1 / 2
El caudal mnimo para prevenir el resbalamiento de lquido es:
3.06 * Vg * A * P
Qg.. min = Ec .(5.66)
TZ
Qg min = caudal mnimo del gas para remover el lquido del fondo MMPCD
Vg =velocidad del gas pie/seg
A =rea conducto o tubera pc
P = Presin psi
Z = Factor se compresibilidad
Ejemplo No 1. Se desea conocer si los caudales determinados en la prueba de produccin

estn por encima del caudal crtico para validar estos datos, la temperatura en Cabeza 95
oF, Gravedad especfica del gas 0.65, rea de la tubera 0.03260506 pie^2, API, 61.7 ,
SGmezcla= 0.711
Datos de Produccin
Presin
Prueba No Fluy. PSI Presin Cabeza psi Caudal MMPCD ck n/64
1 2133 1470 1,647 18
2 1948 1400 2,18 20
Verificacin caudales crticos
Caudales obtenido Velocidad Velocidad Caudal min. Caudal min.
en Pruebas Gas/ Agua Gas /Cond. Gas/ Agua Gas /Cond.
Gas MMpcd pie/ seg. pie/ seg. MMPCD MMPCD
1,65 7,23 5,01 2,45 1,70
2,18 7,42 5,14 2,37 1,64
Como podemos observar en el cuadro para la prueba 1 el caudal de gas no es lo

suficiente para arrastrar el condensado en el fondo ya que el caudal mnimo para arrastrar el
condensado en el fondo es 1.7 MMPCD pero en la prueba se tiene 1.65 MMPCD lo que
significa que no estamos produciendo en forma ptima.
5.8.4.- Cambio del Componente de elevacin
Para flujo de dos fases, el componente de cambio de elevacin se vuelve:
dp g
= s sen Ec. (5.67)
dL el g c
donde s , es la densidad de la mezcla gas lquido en la lnea. Considerando una porcin

de lnea, la cual contiene lquido y gas, la densidad de la mezcla puede ser calculada con la
ecuacin 5.43. Suponindose que no hay prdidas de fuerza de transmisin entre la fase de
lquido y gas, el trmino densidad es definido por la ecuacin 5.44. el uso de la ecuacin
5.43 implica la determinacin de un valor exacto de elevacin de lquido, en vista que la
densidad definida en la ecuacin 5.44 puede ser calculada desde los caudales de flujo
lquido y gas instantneo.
5.85. Cambio del Componente del factor de friccin
El componente del factor de friccin se vuelve:
dp (
f u2 ) f
= Ec. (5.68)
dL f 2 gc d
Donde f, y son definidos de maneras diferentes por diferentes autores. El

componente de friccin no puede ser analticamente pronosticable, excepto para el caso de
flujo laminar de fase simple, en tanto, el mtodo que recibe la mayor atencin es el factor
de friccin de dos fases, de manera que las definiciones ms comunes son las siguientes:
dp f u2
= L L sL Ec. (5.69)
dL f 2 gc d
dp f g g usg2
= Ec. (5.70)
dL f 2 gc d
dp f tp f um2
= Ec. (5.71)
dL f 2 gc d
En general, el mtodo del factor de friccin de dos fases difiere nicamente en la

determinacin de su grado de extensin sobre el flujo constante. Por ejemplo, en flujo
constante, la ecuacin 5.70, es basado en gas normalmente utilizado. En cuanto que en
rgimen de burbuja, la ecuacin 5.69, es basada en lquido utilizado frecuentemente. La
definicin f en la ecuacin 5.71 puede diferir dependiendo del autor.
Muchos autores tienen construidas correlaciones del factor de friccin con algunas
formas de nmero de Reynolds. Las variaciones de nmero de Reynolds utilizados para
validar el factor de friccin, las cuales son definidas, cuando las correlaciones de este factor
son analizadas para una correlacin. Una variacin, se manifiesta de forma diferente,
habiendo muchas correlaciones de numerador del nmero de Reynolds.
5.8.6.-Cambio del Componente de aceleracin
El componente de aceleracin para flujo de dos fases es representado por:
dp ( u du )K
= Ec. (5.72)
dL ace gc d L
El componente de aceleracin es completamente ignorado por algunos autores e

ignorado en algunos modelos de flujo por otros. Cuando los autores consideran esto, hacen
varias suposiciones acerca de la magnitud relativa de sus parmetros para llegar a algn
procedimiento para determinar la cada de presin debido al cambio de la energa cintica.
Del anlisis de los diversos componentes que contribuyen al gradiente de presin

total, y la consideracin principal para resolver las ecuaciones, desarrollando los mtodos
para la prediccin de elevacin de lquido y factor de friccin de dos fases. Es realizado
mediante la aproximacin.
5.9.- Modelo de flujo de dos fases
Cuando dos fluidos con diferentes propiedades fsicas simultneamente estn en una
tubera, se tiene una gama amplia de posibles modelos de flujo, que hacen referencia a la
distribucin de cada fase en la tubera. Muchos autores ponen su atencin en el pronstico
de modelo de flujo que podra existir para varias condiciones fijadas, a los cuales se tiene
diferentes nombres a cada uno de los modelos.
Si es conocido el modelo de flujo es ms confiable la determinacin de la correlacin

de las prdidas de presin. As mismo como resultado del incremento del nmero de lneas
de flujo de dos fases a partir de plataformas y facilidades de plataforma, se ha
incrementado la preocupacin con respecto a la prediccin de no solo el modelo de flujo, si
no del tamao de chorro de lquido y de su frecuencia.
La prediccin de modelos de flujo, para pozos horizontales es la ms problemtica

que para flujos de pozos verticales. Para flujo horizontal, las fases se tienden a separar por
efecto de la densidad provocando una especie de flujo estratificado que es comn.
Govier, present una serie de descripcin de los modelos de flujo horizontal aire
agua, y para flujo vertical aire agua. Estos son mostrados en las Figuras 5.6 y 5.7, las
cuales dependen de la extensin de algunas magnitudes relativas u sL y u sg ,cuando el flujo
ocurre en una tubera inclinada con un cierto ngulo el modelo de flujo toma otra forma.
Para flujos ascendentes inclinado en pozos, los modelos son casi siempre tapn o niebla, el
efecto de gravedad en lquido evita la estratificacin. Para flujo descendente inclinado el
modelo es normalmente estratificado, niebla o anular.
5.9.1 Clculo de la presin transversal
El clculo de la presin de transversal de dos fases envuelve el uso transversal de un

procedimiento iterativo de ensayo y error si existe cambio en la inclinacin y temperatura
de la tubera con la distancia. Calculando transversalmente en la lnea de flujo, el segmento
esta dividido en un nmero de incremento de presin o longitud y las propiedades de fluido
y la gradiente de presin son evaluados en condiciones medias de presin, temperatura e
incremento en la inclinacin de la tubera.
Figura 5.6 Modelo de flujo vertical para dos fases

Figura 5.7 Modelo de flujo horizontal para flujo de dos fases

5.9.2. Determinacin de la distribucin de temperatura
Cuando se requiere que un mtodo de determinacin de la distribucin de

temperatura sea ms preciso, es necesario considerar la transferencia calorfica que se
produce en los fluidos, el clculo de la transferencia calorfica para flujo de dos fases puede
ser muy importante cuando se determina el gradiente de presin en pozos geotrmicos,
pozos de inyeccin de vapor. En general, los clculos de transferencia de calor se realizan
asumiendo una distribucin de temperatura conocida, ms para sistemas multicomponentes
se requiere conocer tanto la temperatura de entrada como de salida del sistema, para los
clculos del comportamiento de transferencia de calor, es necesario primero transformar la
ecuacin de balance de energa en una ecuacin de balance de calor. La combinacin de las
ecuaciones 5.2 y 5.4 y asumiendo que no existe trabajo hecho sobre o por el fluido (dw =
0):
udu g
dh + + dZ + dq = 0 Ec. (5.73)
gc gc
Si la entalpa especfica y el factor adicionado son expresados como calor por

unidad de masa, entonces la energa mecnica equivalente de calor constante, J debe ser
introducida:
udu g
Jdh + + dZ + Jdq = 0 Ec. (5.74)
gc gc
Expresando la elevacin en trminos de ngulo y longitud de tuberas y resolviendo

para un gradiente tenemos:
dh dq g sen udu
= Ec. (5.75)
dL dL gc J g c JdL
dq
El calor adicionado al sistema por unidad de longitud es negativo, por tanto que se
dL
pierde calor cuando la temperatura del fluido es mayor que la temperatura que se encuentra
circundante. El gradiente de prdida de presin de calor puede ser expresado como:
dq U ( d )
= (T Tg ) Ec. (5.76)
dL wT
Donde:
T = Temperatura media del fluido.
Tg = Temperatura media circundante.
U = Coeficiente de transferencia de calor.
wT = Caudal de flujo msico.
5.10.- Flujo en pozos de Gas
Varios mtodos estn disponibles para calcular la cada de presin esttica y fluyente
en pozos de gas. El mtodo ms utilizado es el de Cullender y Smith, y presin y Temp.
Media para gas seco (sistema monofasico) y Grey para Gas condensado (sistema
multifasico). Todos los mtodos inician con la ecuacin 5.13, con modificaciones para la
geometra de flujo. La mayora de los casos, el gradiente de aceleracin es ignorada.
Frecuentemente, es necesario calcular la presin esttica de fondo en un pozo de gas para
cualquier clculo de pronstico y capacidad de entrega.
5.10.1.- Presin de fondo
Para clculos de ingeniera de reservorios y produccin, la presin esttica o cierre en

fondo de pozo, Pwss es frecuentemente exigida. En muchos casos puede ser difcil o caro
obtener la Pws con medidores de presin convencionales estacionando la herramienta
enfrente de la formacin productora. De esta manera, fueron desarrollando tcnicas para
calcular SBHP al considerarse como un caso especial de ecuacin general de flujo vertical
(ecuacin 5.37); las tcnicas de determinacin parten de las siguientes asunciones:
ngulo de inclinacin de 90, entonces sen = 1
Velocidad de flujo es cero.
Tomando estas consideraciones y combinando la ecuacin 4.32 de los gases perfectos

obtenemos:
dp g g
= Ec. (5.77)
dh gc
Donde
pM
g = Ec. (5.78)
ZRT
Combinando esta con la ecuacin 5.77,

dp g M dh
= Ec. (5.79)
p gc Z R T
5.10.1.1. Mtodo de presin y temperatura media (esttica)
Si Z es evaluado a la presin y temperatura media, la integracin de la ecuacin

5.79, es:
Pws H
dp gM

Pwh p
=
g c R Z T 0
dH
La cual
gM
Pws = Pwh EXP Ec. (5.80)
g
c R Z T
Esta ecuacin sostiene cualquier combinacin consistente de unidades. Para las

unidades convencionales de campo,
[
Pws = Pwh EXP (0.01875 g H )/ (T Z ) ] Ec. (5.81)
Donde
Pws = Presin esttica de fondo pozo, psia.
Pwh = Presin en cabeza de pozo, psia.
g = Gravedad del gas (aire = 1)
H = Profundidad del pozo, ft.
T = Temperatura media en tubera, R..
Z = Factor de compresibilidad, evaluado a las
Condiciones medias de presin y temperatura
P = ( Pws + Pwh) / 2 .
La estimacin de Z hace el clculo iterativo, para este mtodo propone la siguiente

solucin cuyos pasos a seguir son:
1) Estimar el valor inicial de la presin esttica de fondo Pws . La siguiente igualdad

proporciona un dato aproximado real:
(
Pws = Pwh 1 + 2.5 105 H ) Ec. (5.82)
2) Con Pws y Pwh obtener la presin y temperatura media.

3) Determinar la presin y temperatura pseudo crticas.
4) Calcular la presin y temperatura pseudo reducidas, para luego encontrar el valor
del factor de compresibilidad por medio de correlaciones o grficos.
5) Haciendo uso de la ecuacin 5.82 valorar la presin de fondo esttica, Pws, repetir
el procedimiento desde el paso 2, donde la nueva presin Pws ser la ltima Pws
obtenida. La iteracin continuar hasta que la diferencia absoluta entre Pws y Pws
sea de 0.001 como margen de error.
5.10.1.2 Mtodo de Cullender y Smith. (Esttica)
Este mtodo considera las variaciones del factor de compresibilidad con los cambios de
presin y temperatura, as mismo como las alteraciones de temperatura con la profundidad.
De la ecuacin 5.79
H
TZ M
p dp = R 0 0.01875 g dH Ec. (5.83)
ZT
I=
p
Integrando el lado derecho de la ecuacin se tiene:
Pws
Idp = 0.01875
Pwh
g H Ec. (5.84)
El lado izquierdo de ambas ecuaciones se resolvi mediante mtodos numricos como

expansin por series, evaluado as mismo (P T / Z) para algn nmero de incremento entre
las presiones estticas de fondo y superficie. Cullender y Smith, proponen una resolucin en
dos pasos dividiendo la profundidad total del pozo en 0, H / 2 y H , posibilitando la
aplicacin tanto en pozos someros como profundos con o sin presencia de gases sulfurosos.
La ecuacin resultante de la integracin es:
Pws
2 Idp = (Pm Pwh )(I ms + I wh ) + (Pws Pms )(I ws + I ms ) Ec. (5.85)
Pwh
Donde:
Pwh = Presin a la profundidad cero, H = 0, psia.
Pms = Presin a la profundidad media, H / 2, psia.
Pws = Presin a profundidad total, H, psia.
Iwh = Integral evaluada a Pwh y Tsup.
Ims = Integral evaluada a Pms y Tmed.
Iws = Profundidad a la cual se quiere evaluar, Th, ft.
Separando en dos expresiones la ecuacin 5.85 e igualando cada una al valor de

0.01875 g H para finalmente despejar la presin a la profundidad media y la presin a la
profundidad total respectivamente, se encuentra que:
0.01875 g H
Pws = Pwh + Ec. (5.86)
I ms + I ts
0.01875 g H
Pws = Pms + Ec. (5.87)
I ms + I ws
Conociendo la temperatura de superficie y de fondo de pozo se puede determinar la
temperatura a cualquier profundidad utilizando la relacin:
T f Ts
Th = Ts + h Ec. (5.88)
H
Donde
Th = Temperatura a cualquier profundidad, R.
Tf = Temperatura final, R.
Ts = Temperatura en superficie, R.
H = Profundidad final, pie.
h = Profundidad a la cual quiere evaluarse Th, pies.
Los siguientes pasos muestran la aplicacin de estas ecuaciones en el clculo de

presin esttica de fondo de pozo.
Paso A
1. Encontrar el factor de compresibilidad para condiciones de boca de pozo. Con ese
valor resolver la integral I.
2. Estimar un valor inicial de Pms para la profundidad H / 2haciendo uso de la ecuacin
5.82, teniendo cuidado de que H sea substituido por H / 2.
3. determinar la temperatura media (ecuacin 5.88).
4. calcular los parmetros pseudo reducidos a las condiciones de penetracin media
y establecer el factor de compresibilidad, Z.
5. con los datos de puntos 2, 3 y 4 calcular la integral Ims usando I = TZ/ P donde P =
Pms .
6. determinar Pms aplicando la ecuacin 5.86. si Pms Pms ir al paso B caso contrario
repetir desde el punto 4 considerando Pms como el nuevo valor inicial para el
reclculo; este procedimiento continua hasta conseguir una diferencia de 10-3 entre
presiones.
Paso B
1. Determinar el valor inicial de presin esttica, Pms, aplicando la ecuacin 5.82.
Previamente sustituir Pwh por Pms y H por H / 2 en lugar de H se debe que el
anlisis del paso B parte de la profundidad final del paso A, que dice, dar la
profundidad media del pozo.
2. Obtener Z a Pms y Tf.
3. Encontrar Iws.
4. Reclcular Pws, como los resultados de los puntos 2 y 3, si Pws Pms el
procedimiento finaliza obteniendo de esta manera la presin esttica de fondo
procurada, si no repetir hasta encontrar dos Pws prximas.
5.10.2.-Presin dinmica de fondo pozo
El pozo es puesto en produccin, con una presin de fondo fluyente que impulsa el
fluido hasta la superficie. Si la energa es suficiente como para vencer la resistencia
encontrada durante la trayectoria de flujo el pozo fluye por surgencia natural; en caso de ser
insuficiente el pozo es sometido a mtodos de recuperacin o levantamiento artificial hasta
obtener la presin adecuada para elevar el fluido.
La produccin depende de:

El rea de formacin
Perforaciones.
Sarta de produccin
Restricciones en boca de pozo (choques) y tuberas de produccin (vlvulas de
fondo).
Facilidades de separacin.
La presin dinmica de fondo de pozo de gas, es la suma de presin en cabeza, la

presin ejercida por el peso de la columna de gas, despreciando la variacin de energa
cintica (velocidad 0), la ecuacin 5.37 de conservacin de energa mecnica puede ser
rescrita de la siguiente forma:
dp g f u2
= cos + Ec. (5.89)
dL g c 2 gc d
5.10.2.1.- Mtodo de presin y temperatura media (dinmico)
Substituyendo la expresin para la densidad del gas en trminos de P, T y Z en la

ecuacin 5.89 resulta:
dp pM f u2
= (cos + Ec. (5.90)
dL ZRT 2 gc d
La integracin de la ecuacin 5.90 asumiendo una temperatura media en la sarta de
flujo y evaluando Z a las condiciones medias de presin y temperatura,
25 g q 2 T Z f ( MD)( EXP( S ) 1)
P = P EXP( S ) +
2
wf
2
wh Ec. (5.91)
Sd 5
Donde
Pwf = Presin dinmica de fondo de pozo, psia
Pwh = Presin de cabeza, psia
S = 0.0375 g (TVD) / T Z
MD = Profundidad media, ft
TVD = Profundidad vertical verdadera, ft
T = Temperatura media, R [(Twf + Twh /2)]
Z = Factor de compresibilidad media
f = Factor de Friccin, f(NRe, e/d) (Jain o Colebrook)
q = Caudal de gas MMscfd
d = Dimetro de la tubera, pulgadas
Para una mejor comprensin en la aplicacin del mtodo, se muestra una aplicacin
prctica.
Ejemplo No 2. Mediante el mtodo de la presin media y temperatura media, se desea

determinar la presin fluyente en cabeza de pozo para los siguientes datos:
Qg = 2.41 MMPCD Prof.=14331 pies SGg = 0.65 RGC = 71603 PC/BBL

Pwf = 9429 psi Diam. = 1.995 pulg. API = 59 Tr= 270 o F
Tc = 70 oF Rugosidad = 0.003 SGm = 0.687 Pr= 10477 psi
1.- Primeramente se debe determinar una presin en cabeza estimada de la siguiente manera
(
pwh = pwf / 1 + 2 , 5 x 10 5
H ) ( )
pwh = 9429 / 1 + 2 , 5 x 10 514331 = 6942
2.- Determinar la constante S
S = 0,0375 g (TVD ) T Z
Con los datos de le gravedad especfica de la mezcla SGM = 0.687
Tr + Tc
Temperatura Media Tmedia = = 170 oF , 630 oR
2
Pr + Pc
Presin Media Pmedia = = (9429+6942)/2 =8185 psi
2
Con estos datos procedemos a calcular el factor de compresibilidad medio Zmed= 1.254,
Ug= 0.037 cp
S = 0,0375* 0.687 * 14331 630 * 1.254 = 0.4672 E^s =1.5956
3.- Determinamos el nmero de Reynolds
20 q g 20 * 2410 * 0.687
N Re = N Re = = 449249
d 0.037 * 1.995
4.- Factor de friccin
1 e 21 . 25
= 1 . 14 2 log + 0 .9

Se despeja y se obtiene f=0.02655
f d N Re
5.- Se procedi a dividir esta ecuacin en dos partes para facilitar su clculo
25 g q 2 T Z f (MD ) ( EXP (S ) 1)
grupo1 =
S d5
25* 0.687 * 2.412 * 630 * 1.254 * 0.02655 * 14331(1.5956 1)

grupo1 = = 1213368
0.4672 *1.9955
6.- Determinamos la presin en cabeza
Pwf ^ 2 Grupo1 9429^ 2 1213368

Pwh2 = Pwh2 = = 7414
Exp( S ) 1.5956
7.- Se debe tomar esta nueva presin en cabeza para determinar un nuevo Z promedio, y
volver a repetir todos los pasos a partir del punto 2 hasta que la nueva presin de cabeza
determinada sea igual a la anterior, la ltima presin determinada 7444 psi .
5.10.2.2.- Mtodo de Cullender y Smith (dinmico)
La derivacin del mtodo de Cullender y Smith, para pozos fluyentes, comienza con
la ecuacin 5.90. Las siguientes substituciones son hechas para la velocidad:
q
u=
A
Psc T Z
q = qsc
Tsc P Zsc
La cual da:
dp pM cos MTZ Psc2 f qsc2
= +
dL ZRT R p Tsc2 2 g c d A2
o
p dp M p
2
= cos + C
ZT dh R ZT
Donde
8 psc2 qsc2 f
C= 2
Tsc g c 2 d 5
Que es la constante para un caudal de flujo en un dimetro particular de tubera.

Separando las variables da:
p
Pwf dp MD
ZT M
Ptf p 2 =
R 0
dL Ec. (5.92)
cos + C
ZT
Es aplicable para cualquier conjunto consistente de unidades. Substituyendo en

unidades de campo e integrando el lado derecho de la ecuacin 5.92 da:
p
Pwf dp
ZT
Ptf 2
p TVD
= 18,75 g MD Ec. (5.93)
0.001 + F2
ZT MD
donde
0.667 f qsc2
F2 = Ec. (5.94)
d5
y
TVD
= cos Ec. (5.95)
MD
Para abreviar la escrita ecuacin 5.93 y dividiendo el comportamiento total del pozo
en dos secciones, H / 2 .
Para encima de H / 2 la ecuacin se simplifica a la siguiente expresin:
18,75 g ( MD) = (Pmf Ptf )(Im f + Itf ) Ec. (5.96)
Para debajo de H /2 antes de la profundidad total del pozo es:
18,75 g ( MD) = (Pwf Ptf )(Iwf + Itf ) Ec. (5.97)

Donde
p
I= ZT Ec. (5.98)
2
p TVD
0,001 + F2
ZT MD
La ecuacin 5.94, se puede simplificar utilizando la ecuacin del factor de friccin

de Nikuradse, para flujo turbulento totalmente desarrollado, esto es:
0,10796 q
F= d < 4,277 in. Ec. (5.99)
d 2,612
0,10337 q
F= d > 4,277 in. Ec. (5.100)
d 2,582
Aplicando la regla de Simpson para obtener un resultado de presin ms exacto, esta es:
18,75 g ( MD) 2 =
(Pwf Ptf )
(Itf + 4 Im f + Iwf ) Ec. (5.101)
3
El siguiente procedimiento es recomendado:
1. Determinar el valor del lado izquierdo de la ecuacin 6.82 para encima de H /2.
2. Determinar F2 de la ecuacin 5.99 o 5.100.
3. determinar Itf de la ecuacin 5.98 y las condiciones de cabeza del pozo.
4. Asumir Imf = Itf para condiciones medias de profundidad del pozo o el punto medio del
tuno de produccin.
5. Determinar Pmf de la ecuacin 5.96.
6. Utilizando el valor de Pwf determinado en el paso 5 y la temperatura media aritmtica
Tmf determinar el valor de Imf de la ecuacin 5.98.
7. Recalcular Pmf de la ecuacin 5.96. Si este valor recalculado no es menos de 1 psi de
Pmf calculado en el paso 5, repita los pasos 6 y 7 hasta que el criterio anterior este
satisfecho.
8. Asuma Iwf = Imf para las condiciones de fondo del tubo de produccin.
9. Repita los pasos 5 hasta el 7, utilizando la ecuacin 5.97 para debajo de H /2 del tubo
de produccin, y obtener un valor de la presin de fondo de pozo, Pwf.
10. Aplicar la regla de Simpson como expresa la ecuacin 5.101, obtenemos un valor de
presin de fondo de pozo ms exacto.
Para un mejor entendimiento mostraremos una aplicacin prctica con los mismos
datos del ejemplo 2, para ver la diferencia entre mtodos.
Ejemplo No 3. Mediante el mtodo Cullender y Smith, se desea determinar la presin

fluyente en cabeza de pozo para los siguientes datos:
Qg = 2.41 MMPCD Prof.=14331 pies SGg = 0.65 RGC = 71603 PC/BBL

Pwf = 9429 psi Diam. = 1.995 pulg. API = 59 Tr= 270 o F
Tc = 70 oF Rugosidad = 0.003 SGm = 0.687 Pr= 10477 psi
1.- Se debe determinar el resultado del lado izquierdo de la ecuacin
18,75 g (MD ) = ( pmf ptf )(I mf + I tf ) 18,75 * 0.687 * 14331 = 184562

2.- Determinar el Valor de F para dimetro menor 4,277
0,10796 q 0,10796 * 2.41

F= d < 4,277 in. F =( )^ 2 = 0.00184
d 2,612 1.9952, 612
3.-Con los datos de le gravedad especfica de la mezcla SGM = 0.687
Tr + Tc
Temperatura Media Tmedia = = 170 oF , 630 oR
2
4.- Clculo para la seccin inferior de flujo
p
I= ZT TVD
2 = cos
p TVD MD
0,001 + F2
TZ MD
9429
1.366 * 630 Iwf=Imf
I= 2
= 103.620
9429 14331
0,001 + 0.001842
1 .366 * 630 14331
18.75 * SGm * MD 18.75 * 0.687 *14331

Pmf = Pwf Pmf = 9429 = 8539
Iwf + Im f 103.62 + 103.620
5.- Con esta presin calculada, se determina el nuevo factor de compresibilidad y el valor
de Imf.
Z= 1.287
8539
I= 1.287 * 630 = 93.40
2 Imf= 93.40
8539 14331
0,001 + 0.00184 2
1.287 * 630 14331
18.75 * 0.687 * 14331

Pmf = 9429 = 8493
103.62 + 93.40
6.- Este clculo se repite hasta que la presin Pwf coincida con Pwf anterior .
7.- Clculo para la seccin Superior de flujo, con la presin Pwf determinada en la parte
inferior se calcula el factor de compresibilidad de la parte superior del flujo.
8493
I= 1.287 * 630 = 93.53 Itf = Imf
2
8493 14331
0,001 + 0.001842
1. 287 * 630 14331
18.75 * SGm * MD 184562

Pmf = Ptf Ptf = 8493 = 7506
Iwf + Im f 93.53 + 93.53
8.- Con esta presin calculada, se determina el nuevo factor de compresibilidad y el valor
de Itf.
Z= 1.233
7506
I= 1.233 * 630 = 57.07
2
7506 14331
0,001 + 0.00184 2
1.233 * 630 14331
9.- Este clculo se repite hasta que la presin Ptf coincida con Ptf anterior
5.10.2.3.- Mtodo de Grey
La ecuacin de Grey (1978) fue especficamente desarrollada para pozos de gas

condensado. Esta correlacin asume que el factor de friccin es dependiente de la rugosidad
efectiva e independiente del nmero de Reynolds. Esta sumisin es razonable solamente
para el flujo de gas a una alta velocidad.
De la ecuacin 5.89, tenemos la ecuacin bsica de energa para una lnea de produccin
vertical
dp g f u2
= cos + Ec. (5.89)
dL g c 2 gc d
Ptotal = Phidrostatico + P friccin
Pasos a seguir para la determinacin de las prdidas de presin por friccin y hidrosttica
qg qsc Bg
Velocidad superficial del gas (pie/seg.) Vsg u sg = = Ec.(5.102)
A A
Velocidad superficial del Cond.(pie/seg.) Vsc 5.6144 qo Bo Ec.(5.103)

uso =
A
Velocidad superficial del Agua (pie/seg.) Vsa qw Bw Ec. (5.104)

usw =
A
Velocidad superficial de la Mezcla (pie/seg.) Vm Vm = Vsg+Vsc+Vsa
Densidad del Gas g (lb./PC) 0.09318 PPMg Ec. (5.105)

g =
ZT
Fraccin de lquido L qc + qw Ec. (5.106)

L =
(qc + qw) + qg
Fraccin de gas g g = 1 Ec. (5.107)
qc
Fraccin de Condensado c c = Ec. (5.108)
(qc + qw)
Fraccin de Agua w w = 1 c Ec.(5.109)
Densidad del Lquido l (lb./PC) l = 62.4SGcc + ww Ec. (5.110)
Densidad de la mezcla noslip (lb/PC) ns = gg + ll Ec. (5.111)
Tensin Interfacial Gas/Agua (Dinas/cm.) w( 280) = 53 0.01048P ^0.637 Ec. (5.112)
ns ^ 2Vm^ 4 Ec. (5.113)

N1 =
g ( l g )
gD ^ 2( l g )
N2 = Ec. (5.114)

1 + 730 Rv
N 3 = 0.0814[1 0.0554 ln[ ] ] Ec. (5.115)
Rv + 1
Vsl
Relacin de Velocidades Rv RV = Ec.(5.116)
Vsg
205
Fraccin de Holup Hl H l = 2.314[N1[1 + ] ]^ N 3 Ec. (5.117)
N2
Gm = l Vsl + g Vsg Ec. (5.118)
Viscosidad de la mezcla Um m = ll + gg Ec. (5.119)
Nmero de Reynolds de la mezcla NRE = m Vm D/m Ec. (5.120)
La rugosidad efectiva Kc es asumida dependiendo del valor de la relacin de velocidades

Rv
Si Rv >0.007 entonce, use la ecuacin Ke = Ko

Si Rv 0.007 entonce, use la ecuacin 5.121
Ke = k + Rv( ( kô-k)/0.0007) Ec. (5.121)
Kô = 28.5 / nsVm^2 Ec. (5.122)
En base a las ecuaciones presentadas las prdidas de presin por friccin y hidrosttica son:
mgdL
Phh = Ec. (5.123)
144 gc
2 ftpGm^ 2dL
Pfric = Ec. (5.124)
144 gcDns
Ejemplo No 4. Se desea determinar la presin fluyente para un sistema multifsico con los
siguientes datos:
Qg = 14.04 MMPCD Prof.=8181 pies SGg = 0.65 RGC = 48169PC/BBL

Qc = 292 BPD Pwf = 3267 psi SGm = 0.687 RGA = 66666 PC/BBL
Qw = 211 BPD Pr= 3709 psi API = 59 Diam. = 2.445 pulg.
Tr= 270 o F Tc = 70 oF Rugosidad = 0.0012
1.- Determinar las velocidades de los fluidos

qg 14.04 *1000000 * 0.0041
Vg (pie/seg.) usg = = = 20.57
A 0.03260 * 86400
Vc.(pie/seg.) 5.614 * 292o

uso = = 0.58
0.0326 * 86400
Vw(pie/seg.) 211 * 1.01

usw = = 0.08
0.0326 * 86400
Vm (pie/seg.) Vm = 20.57+0.58+0.08 = 21.23
2.- Determinamos las distintas propiedades del Gas
g (lb./PC) 0.09318 * 3267 * 20.375

g = = 13.01
0.81 * 588
292 + 211
Fracc.de lquido L L = = 0.0002
(292 + 211) + (14.04 *178126)
Fracc. de gas g g = 1 0.0002 = 0.9998
Fracc. de Cond. c 292

c = = 0.58
(292 + 211)
Fracc. de Agua w w = 1 0.58 = 0.42
l (lb./PC) l = 62.4 * 0.731 * 0.58 + 62.4 * 0.42 = 52.66
m noslip (lb/PC) ns = 13.01 * 0.9998 + 52.66 * 0.002 = 13.02
w( 280 ) = 53 0.01048P ^0.637 = 51.18
13 . 2 ^ 2 * 21 . 23 ^ 4
N1 = = 362672
32 . 17 * 0 . 0074 * ( 52 . 66 13 . 01 )
32 . 17 * 2 . 445 ^ 2 ( 52 . 66 13 . 01 )
N2 = = 102463
0 . 0074
(0.58 + 0.08)
RV = = 0.032
20.57
1 + 730 * 0.032
N3 = 0.0814[1 0.0554 ln[ ] ] = 0.0671
0.032 + 1
Fraccin del Holdup 205

H l = 2 . 314 [362672 [1 + ] ]^ 0 . 0671 = 0 . 004424
102463
Gm =52.66*0.58 +13.01*20.57 = 298.306
m = 0.55 * 0.0002 + 0.023 * 0.9998 = 0.0231

NRE = ((13.02*21.23*2.445)/0.0231)=29250
Rv=0.0320 > 0.007 entonces, use la ecuacin Ke = Ko
Rugosidad efectiva Kô = 28.5 / nsVm^2 = 0.00036
Factor de Friccin 1 e 21 . 25
= 1 . 14 2 log + 0 .9
= 0 . 0246
f d N Re
2 ftpGm^ 2dL
Pfric = =242 psi
144 gcDns
mgdL
Phh = =749 psi
144 gc
P = 242 + 749 = 991 psi Presin en Cabeza = 3267-991 = 2276 psi
5.11.-Flujo de gas en lneas de surgencia
La lnea que conecta la boca del pozo con el separador adquiere mayor importancia
cuando esta tiene una longitud considerable, donde la cada de presin es relevante para
determinar la capacidad productiva del pozo.
Los factores bsicos envueltos en un flujo horizontal son los mismos que se aplicarn
en flujo vertical. En ambos sistemas la cada de presin es la suma de las prdidas por
friccin y aceleracin. La principal diferencia radica en las consideraciones del balance de
energa debido a la posicin de la tubera.
Para propsitos prcticos, en pozos puramente gasferos, el hecho de considerar una

lnea de flujo completamente horizontal permite suprimir los componentes de elevacin y
aceleracin de la ecuacin de energa reduciendo a la siguiente forma:
dp dp
=
dL dL f
Substituyendo la ecuacin 5.19 en esta expresin resulta:
dp f u
2
p M f u2
= = Ec. (5.125)
dL 2 gC d Z R T 2 gC d
En La mayora de los casos, la temperatura T, es considerada como constante y el

factor de compresibilidad, Z se debe evaluar a una presin y temperatura media de la lnea.
Esto requiere una solucin iterativa si una de las presiones es desconocida. La integracin
de la ecuacin 5.125 para una distancia L, entre la presin aguas arriba P1 y la presin
aguas abajo P2 , se tiene:
25 g q 2 T Z f L
P1 P2 =
2 2
Ec. (5.126)
d5
donde
p = presin, psia
g = Gravedad del gas
T = Temperatura, R
q = Caudal de gas, MMscfd (14.7 psia, 60 F)
Z = Factor de compresibilidad a las condiciones de
P yT
L = Longitud de la lnea, ft
d = Dimetro interno, pulg.
f = Factor de friccin
El factor de friccin se puede determinar por la ecuacin de Jain 5.36 o por el

diagrama de Moody (figura 5.4). La ecuacin 5.126, fue derivada usando las condiciones
bases o estndar 14.7 psia y 60 F. Se puede por una forma ms general dejando las
condiciones estndar en la ecuacin como variables. Es frecuentemente ventajoso expresar
la ecuacin 4.55 en trminos de caudal de flujo, esto es:
0.5
CT P12 P22
q= b d
2.5
Ec. (5.127)
Pb g f T Z L
el valor de C depende de la unidades utilizadas en la ecuacin. En la tabla 5.3, estn los

valores de C para varias combinaciones de unidades.
Un factor llamado eficiencia es usado algunas veces en la ecuacin 5.127 para justificar de
hecho, que las tuberas frecuentemente entregan menos gas que el calculado. El factor de
eficiencia normalmente vara entre 0.7 y 0.92 y es normalmente obtenido de la experiencia.
TABLA N.5.3
Valores de C para Diversas unidades
P T D L q C
psia R in. Mi scfd 77,54
psia R in. Ft scfd 5634
psia R in ft MMscfd 5.634x10-3
psia R in M M3/d 1,149x106
Si la cada de presin o caudal de flujo es desconocida, la solucin es iterativa. Si el

dimetro es desconocido, la solucin tambin ser iterativa ya que el dimetro es necesario
para evaluar el factor de friccin. Este factor, incit a los autores para sustituir una ecuacin
especfica para f en la ecuacin general de flujo para hacer la solucin para cualquier q o d
no iterativo. El factor de friccin especificado es dependiente del dimetro o nmero de
Reynolds y no as de la rugosidad de la tubera. Las expresiones incorporadas en la
ecuacin 5.127 para f en varias de las ms populares ecuaciones en lneas de surgencia son
mencionadas a continuacin:
Ecuacin f
0,085
Panhandle A 0 ,147
N Re
0,015
Pamhandle B 0 ,183
N Re
0,187
IGT 0, 2
N Re
0,032
Weymouth 1
d 3
Utilizando estas relaciones para el factor de friccin en la ecuacin 5.127, la forma general
de la ecuacin de flujo de las lneas de surgencia es:
a a a4
T P 2 P22 1 a5
2 3
q = a1E b 1 d Ec. (5.128)

Pb T Z L g
Donde E es el factor de eficiencia, y los valores de las constantes ai usadas en las varias
ecuaciones son tabuladas, en la tabla 5.4.
TABLA N.5.4
Valores de las constantes ai
Ecuacin a1 a2 a3 a4 a5
Panhandle A 435,87 1,0788 0.5394 0,4604 2,618
Panhandle B 737,00 1,0200 0,5100 0,4900 2,530
IGT 337,90 1,1110 0,5560 0,4000 2,667
Weymouth 433,50 1,0000 0,5000 0,5000 2,667
Las unidades usadas en la ecuacin 6.91 son:
q = scfd medido a Tb , Pb
T = R
P = psia
L = miles
d = pulg.
Gas Production Operations H. Dale Beggs, Second printing Nov. 1985
Gas Production Engineering Sanjay Kumar, volume 4, 1987
Natural Gas Production Engineering Chi U. Ikoku
Engineering Data Book Gas Processors Supliers Association, Volume 2, Tenth

Edition, 1987
Applied Reservoir Engineering Charles R. Smith, G. M. Tracy, R. Lance Farrar,

Volume1 1992, OGCI publications.
Gas well testing theory and practice Luis Mattar
i
Production Obtimization Usin Nodal Anlisis, H Dale Beggs, 1991, pag 62-64
ii
Gas Production Operations, H. Dale Beggs, 1984, pag 100
iii
Production Optimization Using Nodal Anlisis, H. Dale Beggs, 1991, pag 64-68
Anlisis de Flujo a Travs de los Choques__________________________________________________ 184
6 Anlisis de Flujo de Gas a Travs De Los Choques
6.1.-Introduccin
El caudal de flujo de casi todos los pozos fluyentes es controlado con un choque en la
cabeza del pozo para controlar el caudal de produccin y asegurar la estabilidad del mismo.
El choque, es un instrumento de restriccin ms comnmente usado para efectuar una
variacin de presin o reduccin de caudal, este dispositivo normalmente se encuentra a la
salida del rbol de surgencia y la lnea de descarga. La Figura 6.1, nos muestra un esquema
grfico del choque y la variacin del dimetro de entrada y salida. La Figura 6.2., nos
muestra el arbolito de produccin y la ubicacin del choque con respecto al mismo,
comnmente el choque esta instalado corriente arriba del cabezal, para evitar dao.
Los choques, consisten en una pieza de metal en forma cilndrica y alargada con un
pequeo orificio para permitir el paso del fluido, tambin son los dispositivos de
restricciones ms comnmente usados para causar una cada de presin o reducir el caudal
de flujo. Son capaces de causar grandes cadas de presin: Un gas que entra en un choque a
5000 psia y salga a 2000 psia o menos. Los choques, entonces tienen, varias aplicaciones
como dispositivos de control en la industria del petrleo y gas. Algunas veces estas
aplicaciones pueden ser utilizados para:
Mantener un caudal de flujo permisible en la cabeza del pozo.
Controlar el caudal de produccin.
Proteger los equipos de superficie.
Controlar y prevenir los problemas de arenamiento al proporcionar suficiente

contra presin en la formacin productora.
Permite obtener informacin para calcular el ndice de productividad en

cualquier tapa de la vida productiva de un pozo.
Prevenir una conificacin de gas y agua.
p1 p2
q d1
d2
Figura 6.1 Esquema del choque

6.2.-Clasificacin de los Choques
Los choques pueden ser clasificados en dos grupos:
1. Choques superficiales.
2. Choques de fondo.
6.2.1 Choques superficiales
Existen diferentes tipos de choques superficiales los cuales se encuentran colocados

en la cabeza del pozo en las lneas de flujo.
La funcin principal del choque superficial es la de estrangular el flujo para

proporcionar estabilidad en las instalaciones superficiales. Los choques superficiales han
sido diseados para su fcil acceso.
Los choques superficiales pueden ser clasificados conforme a su desempeo:
Tipo positivo
Tipo ajustables.
6.2.1.1 Tipo positivo
Esta compuesto de un cuerpo o caja en cuyo interior se puede instalar o cambiar

manualmente diferentes dimetros de orificios, este tipo de choque esta representado en las
Figuras 6.3 y 6.4
En la mayora de los casos se utiliza el de orificio positivo, debido a su simplicidad y

bajo costo.
6.2.1.2.-Tipo ajustable:
Es similar al choque positivo, con la excepcin de que para ajustar el dimetro de la

apertura del orificio de flujo tiene una varilla fina con graduaciones visibles que indican el
dimetro efectivo del orificio. La Figura 6.5, muestra las partes del choque de tipo
ajustables.
6.2.2.-Choques de fondo
Los choques de fondo se encuentran en el interior de la tubera de produccin, mas

propiamente abajo de la vlvula de seguridad. Los choques de fondo son muy tiles para
lograr mayores velocidades en el extremo inferior de la columna de flujo. Utilizando
choques de fondo se pueden sostener contrapresiones ms bajas contra las formaciones
productoras, estimulando as mismo los regmenes de produccin.
Los choques de fondo son frecuentemente utilizados para reducir la presin fluyente
en la cabeza del pozo y prevenir la formacin de hidrato en las lneas superficiales y su
control, la experiencia nos indica que para una velocidad de flujo dada al tamao de la
apertura del choque de fondo es considerablemente menor que el de un choque superficial.
Figura 6.2 Ubicacin del choque

Los choques de fondo se pueden clasificar en tres grupos:
1. Tipo fijo
2. Tipo ajustables
3. Tipo Removibles
6.2.2.1.-Tipo fijo
El choque de fondo tipo fijo se coloca rgidamente adherido a la tubera de

produccin esta se puede ajustar o reposar solo tirando la tubera. Esto es molesto y de alto
costo, en caso de pozos de alta presin se torna ms complicado y difcil. Esta dificultad
con el choque de fondo fijo obliga a combatir su utilizacin, especialmente en pozos que
producen hidrocarburos con arena, que ocasionalmente tienen un desgaste rpido y
requieren reposicin frecuente del choque. Por que se pueden ajustar y compensar el
desgaste, por lo general girando la tubera de produccin o manipulando las varillas en
superficie.
Figura 6.3 Porta orificios tipo positivo
6.2.2.2.-Tipo Removible
Es el ms prctico en control de pozo en el fondo del mismo, el cual se puede tirar a

la superficie para su inspeccin, ajuste o reposicin e instalar de nuevo a cualquier
profundidad deseada en la columna de flujo sin alterar la tubera de produccin.
El choque de fondo removible Otis, esta montado en un empacador equipado con

cuas de modo que se puede colocar en la tubera de flujo a cualquier profundidad deseada.
Se puede tirar con herramientas especiales que se bajan con Wire Line para medir desde la
superficie. Consiste de un mandril con, un juego de cuas, un elemento de empaque y el

dimetro del orificio que ajusta libremente en una superficie cnica y tiene barras de
extensin que se extienden para la cima por una borda en el mandril y dentro de los
casquillos en un collar se puede deslizar libremente en el extremo superior del mandril.
Figura 6.4 Choques positivos
6.3 Factores que influyen en el choque
Existen diversos parmetros que influyen en la seleccin del choque ptimo. La

seleccin del dimetro del choque, esta influenciado por la presin fluyente, el ndice de
productividad, razn gas petrleo, etc. No menos importante, es conocer el sistema de
flujo para el cual sern seleccionadas las ecuaciones para este propsito. Diversos factores
pueden ser considerados para determinar la produccin o caudal en el sistema de flujo,
entre estos tenemos:
Comportamiento de entrada del flujo.

Sistema superficial.
Sistema subsuperficiales.
Instalaciones superficiales.
Instalaciones subsuperficiales.
Todos estos factores mencionados, son para considerar el mejor manejo de flujo de
fase simple o de dos fases.
Figura 6.5 Vlvula de agolla regulable
6.4.-Modelos de Flujo
Cuando el flujo de gas o la mezcla gas-lquido fluyen a travs del choque el fluido
puede acelerarse hasta alcanzar la velocidad del sonido en la garganta del choque, cuando
esta condicin ocurre el fluido es llamado critico, y el cambio de presin aguas abajo del
choque no afectan al caudal de flujo, porque las perturbaciones de la presin no pueden
viajar agua arriba ms rpido que la velocidad sonica. As para la prediccin de la cada de
presin para un caudal de flujo relacionada para los fluidos compresibles fluyendo a travs
del choque, deberamos determinar el flujo crtico, o subcritico, para aplicar su
correspondiente correlacin. La Figura 6.6, muestra la dependencia del caudal de flujo a
travs del choque o la razn de las presiones aguas arriba a aguas abajo para un fluido
compresible. El caudal es independiente de la relacin de presin aguas arriba/aguas abajo,
cuando el fluido es crtico.
Dependencia del Flujo a traves del choque
Qg Caudal de Gas MMPCD 0,3
0,2
Flujo Critico Flujo Subcritico

0,1
Flujo a traves del

choque
0
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 1,1
Relacion de Presion P2/P1
Figura 6.6 Dependencia del flujo a travs de los choques
El flujo a travs de Choques, en general, puede ser de dos tipos: Subcrticos y Crtico..
6.4.1.-Flujo Subcrtico
El flujo es llamado Subcrtico cuando la velocidad del gas a travs de las restricciones
es menor a la velocidad del sonido del gas, y el caudal depende tanto de la presin de
entrada como de la presin salida. Los choques subsuperficiales son normalmente
proyectados para permitir el flujo subcrtico.
6.4.2.-Flujo Crtico
El flujo es llamado crtico cuando la velocidad del gas a travs de las restricciones es
igual a la velocidad del sonido (1100 ft/sec para el aire) en el gas. La velocidad mxima en
la cual un efecto de presin o una perturbacin se pueden propagar a travs de un gas no
puede exceder la velocidad del sonido del gas. As mismo, una vez que la velocidad del
sonido sea alcanzada, un aumento mayor en la diferencial de presin no aumentar la
presin en la garganta del choque. Por lo tanto, el caudal de flujo no puede exceder al
caudal de flujo crtico conseguido cuando la razn de presin salida P2 y entrada P1 llega a
un valor crtico, indica que esta presin sea decrecida. Al contrario del flujo subcrtico, el
caudal de flujo, en flujo crtico, depende solamente de la presin entrada, porque las
perturbaciones de presin que trafican en la velocidad del sonido implican que una
perturbacin de presin en la extremidad salida no tendr efecto alguno sobre la presin
entrada o sobre el caudal de flujo. Los choques en superficies son normalmente desechados
para favorecer un flujo crtico.
6.4.2.1.- Pruebas de flujo crtico
Otro tipo de medicin de gas, es las pruebas de flujo crtico en la cual el flujo de gas
es medido en la descarga a bajas presiones a travs de un orificio en la atmsfera. Esta
medida no es bastante segura debido a que se tiene una variacin del 1 al 6 % del caudal de
gas pero es muy apropiada para lugares remotos, donde no existen lneas y el acceso es
muy difcil. Debido a estos problemas presentados es necesario realizar las mediciones de
caudales a condiciones de flujo crtico para poder determinar la potenciabilidad del
reservorio y ver la factibilidad de la explotacin.
Cuando el flujo de gas pasa a travs de un orificio, la velocidad aguas arriba es

controlada por la diferencial de presin a travs del orificio. Como la diferencial de presin
incrementa por un incremento en la presin agua arriba, una disminucin en la presin
aguas abajo o ambos la velocidad del gas aumenta y se acerca a un valor conocido como la
velocidad crtica. Despus de que esta velocidad crtica alcanza la presin diferencial a
travs del orificio no afecta a la velocidad. La velocidad crtica se alcanza cuando la
relacin de la presin aguas arriba / aguas abajo es aproximadamente 2 o ms. Sobre la
relacin de presin crtica no importa cunta presin est aplicada en el lado de las aguas
arriba.
El Bureau of Mines, nos da la siguiente tabla de los factores de orificio y las

formula para los clculos de caudal de gas.
Qg = C P / (SGg * T)^0.5 Ec. 6.1

Donde :
Qg = Caudal de Gas MPCD
C = Coeficiente de orificio dado en tablas
P = Presin Superficie psi
SGg = Gravedad Especfica del gas
T = Temperatura Fluyente ( o F +460 ) o R
Choques Coeficiente Mxima Choques Coeficiente Mxima

n/64 Promedio Variacin n/64 Promedio Variacin
(%) (%)
4 1,524 3,61 32 101,8 2,29
6 3,355 1,14 40 154 1,56
10 6,301 2,25 48 224,9 1,03
12 14,47 3,88 56 309,3 2,31
14 19,97 3,82 64 406,7 2,09
16 25,86 1,88 72 520,8 1,26
20 39,77 2,13 80 657,5 3,61
24 56,58 2,74 88 807,8 2,05
28 81,09 2,33 96 1002 6,32
Determinacion del Coeficiente Promedio
1200
4 3 2
y = 6E-06x - 0,0007x + 0,1249x - 0,2966x + 0,4642
2
1000 R = 0,9999
Coeficiente Promedio
800
600
400 Coeficiente Promedio
200 Polinmica (Coeficiente

Promedio )
0
0 20 40 60 80 100 120
Choque n/64
Ejemplo se desea determinar el caudal de gas para los siguientes datos superficiales
P.surgencia = 2500 psia, Superficie = 85 o F, SGg = 0.67, ck = 24
Qg = 56.58*2500 / (0.67 * 545)^0.5

Qg = 7.402 MMPCD
6.5.-Flujo de Gas
La ecuacin general para flujo a travs de restricciones puede ser obtenida

combinando la ecuacin de Bernoulli con la ecuacin de estado. Las prdidas irreversibles
o prdidas de fatiga son explicadas por un coeficiente de descarga, el cual depende del tipo
de restricciones.
6.5.1.- Flujo de gas simple fase
Cuando un fluido compresible pasa a travs de una restriccin, la expansin del fluido
es un factor muy importante. Para fluidos isotropitos de un gas ideal a travs de un choque ,
el caudal esta relacionado a la relacin de presin P1/P2. La ecuacin general es valida
solamente en el rgimen subcrtico, antes del flujo crtico, cuando la mxima velocidad de
flujo (igual a la velocidad del sonido) es obtenida.
0, 5
(k +1) k
k p 2
2k
2 1 p2
q sc = 974,61C d p1 d ch Ec.(6.2)
g T1 k 1 p1 p1

Donde
q sc = Caudal de flujo de gas, Mscfd ( a 14,7 psia y 520R).

d ch = Dimetro del Choque, in.
p1 = Presin entrada del Choque, psia.
p2 = Presin salida del Choque, psia
T1 = Temperatura de entrada (o upstream), R.
Cd = Coeficiente de descarga que depende de la relacin de
dimetros de tubera y de las restricciones, el nmero de Reynold y la geometra del
dispositivo, generalmente se considera 0,865.
k = ( )
Relacin de calores especficos c p c v , sin dimensin
g = Gravedad del gas (para el aire=1)
Si la razn de presiones en la cual el flujo crtico ocurre, es representado por

( p 2 p1 )c , entonces a travs del clculo elemental se puede demostrar que:
k (k 1)
p2 2
= Ec.(6.3)
p
1 c k + 1
El flujo es subcrtico para ( p2 p1 )> ( p 2 p1 )c y crtico para

( p2 p1 ) ( p 2 p1 )c , Si la razn de presin de la operacin es menor que la razn de
presin crtica, la razn de presin ser de ecuacin 6.1 debe ser substituida por la razn de
presin crtica, porque la razn de flujo mxima a travs del choque es aquella que
corresponde al flujo crtico.
El valor de k normalmente es asumido, como una constante, k 1.293 para una

gravedad de 0.63. Generalmente, el valor de k utilizado entre 1.25 y 1.31, implicando a
partir de la relacin representado por la ecuacin 6.2 que el flujo es crtico cuando la razn
de presin esta en la regin de 0.5549 (para k = 1.25) para 0.5439 (para k = 1.31). El flujo
es considerado crtico cuando la razn de presin es menor o igual a 0.55, el cual indica que
k sea a 1.275 aproximadamente. Substituyndose ( p 2 p1 )c = 0.55 , y k=1.275 en la
ecuacin 6.1, resulta en:
2
456.71C d p1 d ch
q sc = Ec.(6.4)
( g T1 ) 0.5
Donde:
q sc = Caudal de flujo de gas, Mscfd ( a 14,7 psia y 520R).

d ch = Dimetro del Choque, in.
p1 = Presin entrada del Choque, psia.
T1 = Temperatura de entrada (o upstream), R.
Cd = Coeficiente de descarga, 0.865
g = Gravedad del gas (para el aire=1)
El mximo caudal de flujo de gas ocurre cuando el flujo es crtico. Para cualquier valor
de la razn de presiones por debajo del valor crticoi.
6.5.2.-Flujo de dos fases
La siguiente ecuacin puede ser utilizada para determinar la relacin entre la pwh , ql
y d para flujo de gas lquido en rgimen crticoii.
b ql (GLR )c
pwh = Ec.( 6.6)
da
Donde,
pwh = Presin de cabeza del pozo, psia

ql = Caudal de flujo lquido, STB day
GLR = Razn gas - lquido, scf STB
d = Dimetro del choque, pu lg adas
Valores de a , b y c son propuestos por diferentes autores, que estn presentados en la

Tabla 6.1.
Tabla 6.1
Constantes del Choque
Autores A b c
Ros 2,00 4,25 x103 0,500
Gilber 1,89 3,86 x103 0,546
Baxendell 1,93 3,12 x103 0,546
Achong 1,88 1,54 x10 3 0,650
Ejemplo No 1. Se desea determinar la presin despus del choque para un fluido

monofsico
Qg= 2.41 MMPCD CK 16/64 Tc=70 oF

Pcab. Antes ck= 2501 psi SGg=0.687
1.- Determinamos la relacin de la presin crtica (Pch/Pwh)crit

K (cp/cv)=1.25
k ( k 1) 1.25 (1.25 1)
p2 2 p2 2
= = = 0.5549
p1 c k + 1 p1 c 1.25 + 1
2.- Determinamos el caudal crtico
Cd = 0.865
0,5
( k +1) k
k p2
2k
2 1 p2
qsc = 974,61Cd p1 d ch
g T1 k 1 p1 p1

0,5

q sc = 974,61* 0.865 * 2501* (16 / 64) 21 1.25
0.687 * 530 1.25 1
ch [
(0.5549)2 1.25 (0.5549)(1.25+1) 1.25 ]

= 3.214
g
3.-Definimos el rgimen de flujo Crtico o Subcrtico
Caudal de prueba 2.41< 3.214 caudal critico (El flujo es Subcritico)
4.- Debido a que la ecuacin es muy grande se subdividi la ecuacin
1 q sc2 k 1
grupo 2 = g * T1 = 0.0244
(0.97461* Cd * d ) 2 2
ch
p12 k
5.-Clculo de la relacin de presin
k k 1
p2 p2
2 k
= 1 grupo 2 *
p1 p1 c
p2
(
= 1 0.0244 * (0.5549 )c 2 1.25 )
1.25 1.25 1
= 0.7234
p1
6.- Con esta nueva relacin de presin determinada, tenemos que remplazar la relacin
crtica determinando una nueva relacin de presin crtica hasta igualar esta nueva relacin
con la anterior. Con esta relacin se procede a calcular la presin despus del choque.
Gas Production Engineering Sanjay Kumar, volume 4, 1987
Petroleum Produciton Systems, Michael J. Economides, A. Daniel Hill,
Christine Ehlig-Economides, 1994
Production Optimization, Using Nodal Anlisis, H. Dale Beggs, 1991
i
Petroleum Produciton Systems, Michael J. Economides, A. Daniel Hill, Christine Ehlig-Economides, pag
229
ii
Production Optimization, Using Nodal Anlisis, H. Dale Beggs, 1991, pag 124
Balance De Materiales Para Reservorios de Gas______________________________________________197
7 Balance De Mat eria e n Reser vorios de Gas
7.1.-INTRODUCCIN
El balance de materia en reservorios de gas, es un balance de moles de gas existente

en el reservorio. Matemticamente el balance de materia esta representada a travs de una
ecuacin denominada ecuacin de balance de materiales (EBM). Los moles existente en un
reservorio en un determinado instante son la diferencia entre los moles original en el
reservorio y los moles producidos. Como los volmenes de los fluidos producidos son
generalmente medidos a una determinada condicin estndar de presin y temperatura. La
ecuacin de balance de materiales es comnmente escrita de manera que para cualquier
instante, el volumen de los fluidos existente en el reservorio sea la diferencia entre el
volumen inicial y el volumen producido, ambos medidos a esa condicin de presin.
La ecuacin de balance de materiales para reservorio de gas esta sujeto a los dos
principales mecanismos de produccin, tales como la expansin de fluido y la entrada de
agua. Como principales utilizaciones prcticas de la ecuacin de balance de materiales se
pueden citar:
Determinacin de volumen original de gas.

Determinacin de volumen original de condensado.
Determinacin de la entrada de agua proveniente de acuferos.
Prevencin del comportamiento de reservorios.
Inicialmente la EBM ser utilizada para la determinacin de los volmenes originales

de gas existente en un reservorio, en cuanto en captulos posteriores este dato ser la base
para los estudios de previsin y comportamiento futuro utilizando mtodos analticos de
entrega de gas.
Conocindose el volumen de reservorio, la porosidad de roca y la saturacin de agua

irresidual, pueden ser calculados los volmenes originales de gas a travs del mtodo
volumtrico. El volumen de gas, medido en condiciones estndar, esta dado por:
V r Sgi
G= Ec.(7.1 )
Bgi
V r (1 Swi ) Ec.(7.2 )
G=
Bgi
Donde Vr es el volumen total del reservorio, es la porosidad de la roca. Sg la
saturacin de gas, Bg el factor volumtrico de gas, Sw la saturacin de agua y el ndice i se
refiere a las condiciones iniciales. El factor volumtrico de gas (Bg) , as mismo como otras
propiedades de los fluidos de la roca reservorio pueden ser determinados en laboratorio o
Balance De Materiales Para Reservorios de Gas______________________________________________198
estimados a travs de correlaciones empricas. Como se demuestra en el captulo 3 de

propiedades de los fluidos.
En muchos casos la porosidad, la saturacin de agua irresidual y el volumen del

reservorio no son conocidos con la precisin necesaria para el mtodo volumtrico, por
tanto no puede ser aplicado, en esas situaciones la EBM puede ser empleada para calcular
el Volumen In-Situ y la entrada de agua, en el caso que el reservorio tenga el empuje del
acufero.
La utilizacin del mtodo de balance de materiales exige un histrico de produccin

del reservorio, la calidad de los resultados a ser obtenidos dependen mucho de la calidad de
los datos registrados en el histrico de produccin. Eso significa que las cantidades de
agua, y gas producidos en un campo gasifero dependen de las presiones de reservorio y sus
medidas, ya que son instrumentos valiosos para el estudio de los reservorios portadores de
fluidos.
Los reservorios de gas pueden ser clasificados de acuerdo a su energa con la cual
producen, los cuales pueden ser:
Reservorios Volumtricos
Reservorios con empuje de agua.
7.2.-ECUACIN DE BALANCE DE MATERIAL
La ecuacin de balance de materal puede ser obtenida a partir del principio de

conservacin de masa en el interior del reservorio, representado por la siguiente expresin:
Masa producida = masa inicial masa actual.
Si la composicin del gas producido es constante, los volmenes producidos y

volmenes remanente en el reservorio son directamente proporcional a la masa, ya que para
la ecuacin de estado de los gases se tiene:
Masa ZRT
V= Ec.(7.3 )
Masa molecular p
Donde Z es el factor de compresibilidad del gas, R la constante universal de los gases,

T la temperatura y P la presin. En el capitulo 3 nos presenta los mtodos, y correlaciones
para la determinacin de Z. Dada la masa molecular del gas (obtenida a partir de su
composicin), el volumen de gas (medido bajo una cierta condicin de presin y
temperatura). La ecuacin de conservacin de masa puede ser expresada a travs de un
balance de volmenes, bajo una condicin de P y T cualquiera de referencia:
Volumen producido = volumen inicial volumen actual Ec.(7.4 )

Explotación Del Gas y Optimización de La Producción

Cargado por

Copyright:

Formatos disponibles

Explotación Del Gas y Optimización de La Producción

Cargado por

Información del documento

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

Explotación Del Gas y Optimización de La Producción

Cargado por

Copyright:

Formatos disponibles

Explotacin del Gas y

Balance de Materiales con entrada de agua

Jos Luis Rivero S.

Software y Mtodos Analticos

Por: Jos Luis Rivero Sandoval

Debido al crecimiento de la demanda de la energa, en los distintos campo de la

La publicacin de este libro es muy importante en el mercado del gas para

Los tres programas presentados WELLGAS, BALANGAS, PRONOSGAS ,

Ing. Jos Luis Rivero

Santa Cruz de la Sierra Bolivia

Primeramente quiero agradecer a la Universidad Gabriel Rene Moreno y YPFB

Al Dr. Hugo Araujo por la revisin criteriosa del libro, comentarios y

Al Ing. Luis Kin por las contribuciones recibidas durante el proceso de

Finalmente quiero agradecer a mi esposa Amaly por la comprensin, apoyo e

Captulo 1: Anlisis Nodal Pozos de Gas

1.1.- Introduccin al Anlisis Nodal10

Captulo 2: Clasificacin de los Fluidos en el Reservorio

2.2.12 Saturacin Crtica del Fluido32

Captulo 3 : Propiedades del Gas Natural

3.1 Introduccin a las Propiedades del Gas Natural.50

3.3.6 Correlacin de Hall Yarborough71

Captulo 4 : Anlisis de Reservorio

4.9 Perforaciones de Pozo y Efecto de Dao..........................................................132

Captulo 5: Perdida de Presin en Lneas de Produccin

Captulo 6 : Anlisis de Flujo de Gas a travs de los Choques

Captulo 7: Balance de Materiales Reservorio de Gas

Captulo 8 : Intrusin de Agua

8.6.2 Modelo Pseudo Permanente Modificado (MPSS).236

Captulo 9: Desarrollo de un Campo de Gas

Captulo 10: Manual y Aplicacin de los software

El mayor crecimiento registrado fue dado en el campo de la generacin de

En los pases industrializados, dadas estas ventajas econmicas ambientales, la

En base a la experiencia obtenida con el transcursos de los aos en

El clculo de reserva es muy importante ya que nos da una idea general de la

Otro factor importante que se debe tomar en cuenta en el balance de materiales

En el captulo 9, mostramos el proceso de un desarrollo ptimo de campo de gas

En el captulo 1, sealamos los conceptos bsicos del anlisis nodal para la

En el captulo 2-3, clasificamos los fluidos en el reservorio en funcin al

En los captulos 4-5-6, se explica el anlisis de reservorio, los distintas tipos de

En el capitulo 7-8, sugerimos un nuevo mtodo de Balance de materiales para

En general, el libro Explotacin del Gas y Optimizacin de la Produccin nos da

El desarrollo de un campo de gas no solo depende del reservorio y la

El anlisis nodal es presentado para evaluar efectivamente un sistema completo de

El objetivo principal del anlisis nodal, es el de diagnosticar el comportamiento de

Para que ocurra el flujo de fluidos en un sistema de produccin, es necesario que la

1. Flujo a travs del medio poroso.

p1 = pr pwfs = Prdidas de presin en medios porosos.

p7 = pwf pwh = Prdida de presin total en la tubera de produccin

Figura 1.1(Sistema de Produccin)

La prdida total de presin en un sistema de produccin es el punto inicial Pr

La presin en cada componente es dependiente del caudal de produccin, el caudal

extremadamente pequeas. Otro ejemplo de error en el diseo de terminacin es

La Interrelacin entre caudal y presin es aprovechada por el Anlisis Nodal para

1.2.-Anlisis del Sistema de Produccin

La optimizacin de la produccin de pozos de gas para un Sistema de Produccin,

Se darn algunas definiciones sobre el Anlisis del Sistema de Produccin de

Desde que el Anlisis Nodal TM fue propuesto la esencia del contenido se

Primero, la colocacin de los nodos, que se pueden encontrar en diferentes partes

Optimizacin del Dimetro de Tubera

Figura 1.2 (Determinacin de la capacidad de flujo)