130363

METODOLOGA PARA OBTENER UN MODELO DE POZO QUE PERMITA
PREDECIR EL COMPORTAMIENTO DE LA PRODUCCIN DE UN POZO

DE GAS: CASO ESTUDIO CAMPO CUPIAGUA.

NEYLA YOHANA ARIAS HOYOS

UNIVERSIDAD INDUSTRIAL DE SANTANDER
FACULTAD DE INGENIERAS FISICO-QUIMICAS
ESCUELA DE INGENIERA DE PETRLEOS
BUCARAMANGA
2009
METODOLOGA PARA OBTENER UN MODELO DE POZO QUE PERMITA

PREDECIR EL COMPORTAMIENTO DE LA PRODUCCIN DE UN POZO DE
GAS: CASO ESTUDIO CAMPO CUPIAGUA.

NEYLA YOHANA ARIAS HOYOS

Trabajo de Grado para optar por el titulo de Ing. de Petroleos

DIRECTOR
NICOLAS SANTOS SANTOS
Director Escuela de Ingenieria de Petroleos

UNIVERSIDAD INDUSTRIAL DE SANTANDER
FACULTAD DE INGENIERIAS FISICO-QUIMICAS
ESCUELA DE INGENIERIA DE PETROLEOS
BUCARAMANGA
2009

TABLA DE CONTENIDO

INTRODUCCIN 1
1. GENERALIDADES 3
1.1 CLASIFICACIN DE LOS YACIMIENTOS DE GAS 3
1.1.1 Yacimientos de gas seco 3
1.1.2 Yacimientos de gas hmedo 4
1.1.3 Yacimientos de gas condensado ...5
1.1.4 Problemas asociados a los yacimientos de gas condensado ..7
1.1.4.1 Acumulacin de liquido en el pozo .8
1.1.4.2 Acumulacin de lquido en la formacin ..9
2. BANCO DE CONDENSADO ..11
2.1 MODELO TERICO .12
2.1.1 Anlisis del modelo .14
3. ENTREGABILIDAD DE POZOS DE GAS CONDENSADO ..18
3.1 ECUACIN PROPUESTA ...23
3.1.1 Regmenes de flujo .25
3.1.1.1 Regin cercana al pozo .25
3.1.1.2 Regin de acumulacin de condensado .26
3.1.1.3 Regin de solo gas .26
3.1.2 Coexistencia de las regiones de flujo ..27
3.1.3 Calculo de pseudopresion .27
3.1.3.1 Regin 1 ...29
3.1.3.2 Regin 2 ...30
3.1.3.3 Regin 3 ...30
3.2 MTODO SEMIANALITICO PARA PREDECIR LA

ENTREGABILIDAD32
4. CLCULOS DE LA PRODUCTIVIDAD .39
4.1 MODELO BALANCE DE MATERIA ...40
4.1.1 Calculo de la razn de aceite fluyendo ( ) en la regin 1 41
4.1.2 Distribucin de presin alrededor del pozo 44
4.1.2.1 Pozo vertical sin fractura ...45
4.1.2.2 Pozo vertical fracturado .46
4.2 ESTIMACIN DE LOS PARMETROS DE CONTROL DE LA
PRODUCTIVIDAD EN POZOS DE GAS CONDENSADO .48
4.2.1 Estudio de simulacin 49
4.2.1.1 Simulacin establecida ..49
4.2.2 Influencia de la permeabilidad relativa en la productividad de un
pozo ..51
4.2.2.1 Permeabilidad relativa al gas a la saturacin de liquido
irreducible .52
4.2.2.2 Permeabilidad efectiva al gas a la saturacin de aceite
cerca a la cara del pozo .52
4.2.2.3 Exponente Corey del gas ..53
4.2.3 Estimacin de la permeabilidad relativa .53
4.2.3.1 Formulas de pseudopresion para una y dos fases ...54
4.2.3.1.1 Pseudopresion una sola fase .54
4.2.3.1.2 Pseudopresion dos fases 56
4.2.4 Generacin de una curva de permeabilidad estimada .58
4.3 PRODUCTIVIDAD DE LAS FRACTURAS EN POZOS DE GAS
CONDENSADO .60
4.3.1 Modelo representativo del yacimiento .62
4.3.2 Presentacin de los resultados en trminos del ndice de

productividad ...63
4.3.2.1 Modelo de fracturas hidrulicas ..64
5. PERMEABILIDAD RELATIVA 66
5.1 MODELO DE PERMEABILIDAD RELATIVA 67
5.1.2 Simulacin numrica..70
5.2 SATURACIN CRITICA DE CONDENSADO Y PERMEABILIDAD
RELATIVA AL GAS ...73
5.2.1 Modelo de red para condensacin en medio poroso 75
5.2.2 Longitud critica de liquido y saturacin critica de condensado en
un tubo capilar .76
5.2.3 Calculo de la saturacin critica de condensado ..77
5.2.4 Calculo de la conductancia del fluido ..78
5.2.5 Renormalizacion del mtodo y permeabilidades relativas
....81
6. FLUJ O COMPOSICIONAL .86
6.1 MODELO COMPOSICIONAL DEL RETROGRADO86
6.1.1 Ecuaciones que gobiernan los simuladores composicionales de
flujo en dominios multimecanicos ...87
6.2 EFECTOS DE LOS CAMBIOS COMPOSICIONALES EN EL
TRASIENTE DE FLUJ O 94
6.2.1 Desarrollo del modelo .96
6.3 INTERACCIONES NO-LINEALES ENTRE LAS RESTRICCIONES DE
FLUJ O CERCA AL POZO 100
6.3.1 Restricciones de flujo cerca al pozo .102
6.3.1.1 Representacin del pozo .103

6.3.2 Anlisis del modelo 106
6.4 FLUJ O EN LAS REGIONES PERFORADAS 107
6.4.1 Simulador de pozos perforados 3D ...108
6.4.1.1 Ecuaciones que gobiernan, una fase ..109
6.4.1.2 Mtodo de solucin matemtico, una fase ...112
6.4.1.3 Ecuaciones que gobiernan, dos fases ..111
6.4.1.4 Propiedades del fluido, dos fases ..113
6.4.1.5 Permeabilidad relativa, dos fases ..114
6.4.1.6 Mtodo matemtico de solucin, dos fases .115
6.4.2 Simulador de pozo a hueco abierto 1D 117
6.4.3 Anlisis del modelo ..117
6.5 FLUJ O EN POZOS FRACTURADOS HIDRULICAMENTE ..118
6.5.1 Simulador In-House .120
6.5.1.1 Ecuaciones que gobiernan 120
6.5.2 Sistema 2D para pozos fracturados hidrulicamente .121
6.5.3 Sistema 1D para hueco abierto ..123
6.5.4 Conductividad efectiva de la fractura 124
6.5.5 Sistemas de flujo no-Darcy en una sola fase ..125
6.5.6 Flujo en sistemas de gas condensado ..128
6.5.6.1 Relacin de movilidad en la fractura 128
6.5.6.2 Relacin de movilidad en la matriz 131
6.5.7 Factor dao de fractura y radio efectivo del pozo ...133
6.5.7.1 Mtodo iterativo para estimar el radio efectivo del
pozo.134
6.5.8 Anlisis del modelo 136
6.6 MODELO SIMPLIFICADO DE FLUJ O CERCA AL POZO ..137
6.6.1 Enfoque computacional y resultados.138
6.6.2 Ejemplo de aplicacin ..138
6.6.2.1 Datos Occidente de frica .138
6.6.2.2 Datos Govier-Fogarasi .140

6.6.2.3 Datos comisin de Railroad de Texas ..140
6.6.3 Incertidumbre de los datos de entrada para calcular la presin de
fondo fluyendo ..140
6.6.3.1 Precisin absoluta de dos medidores ..140
6.6.3.2 Tasa de liquido o CGR tendencia/eficiencia del
separador141
6.6.3.3 Rangos de la gravedad del condensado ..142
6.6.3.4 Pared rugosa .142
6.6.4 J ustificacin del modelo homogneo 142
6.6.4.1 Comparacin con la correlacin de Gray 143
6.6.4.2 Comparacin con el modelo de Ansari et al. .144
6.6.4.3 Comparacin con el modelo de Aziz et al. ...144
6.6.5 Anlisis del modelo ..145
7. METODOLOGA PARA OBTENER UN MODELO DE POZO 146
7.1 METODOLOGA APLICADA AL CAMPO170
7.1.1 Generalidades del campo Cupiagua .170
7.1.2 Aplicacin ..170
8. CONCLUSIONES ..180
RECOMENDACIONES .182
BIBLIOGRAFA ...183

LISTA DE FIGURAS

Figura 1: Diagrama presin-temperatura para yacimientos de gas condensado
....6
Figura 2: Diagrama presin-temperatura para un yacimiento de gas condensado
pobre .....7
Figura 3: Diagrama presin-temperatura para un yacimientos de gas condensado
rico .8
Figura 4: Formacin del bloque de condensado .12
Figura 5: Perfil de presin de un yacimiento homogneo, que se encuentra
fluyendo por encima de la presin del punto de roco.15
fluyendo por debajo de la presin del punto de roco..16
Figura 7: Diagrama esquemtico del mecanismo de condensacin cerca al pozo
..33
Figura 8: Diagrama de fases para la composicin del fluido original y la
composicin del banco de condensado 34
Figura 9: Perfil tpico del banco de condensado con numero capilar junto con la
aproximacin por el mtodo semianalitico....37
Figura 10: Saturacin de la fase aceite, alrededor de un pozo vertical sin
fracturas..42
Figura 11: Relacin gas-aceite fluyendo en un pozo vertical sin fracturas

..43
Figura 12: Tasa de flujo de aceite a condiciones estndar en un pozo vertical
..44
Figura 13: Distribucin de la presin alrededor de un pozo sin fractura con dao a
la formacin ...45
Figura 14: Grafico log-log de derivadas de pseudopresion normalizada para 10
das de agotamiento con diferentes tiempos de produccin .55
Figura 15: Grafico log-log de datos de simulacin con y sin efectos del numero
capilar .57
Figura 16: Saturacin del condensado despus de 60 das. (Incluyendo los efectos
del numero de captura) ...72
Figura 17: Saturacin del condensado despus de 60 das (sin los efectos del
numero de captura) ..73
Figura 18: Normalizacin del ndice de productividad 74
Figura 19: Configuracin del liquido en tubo circular ..79
Figura 20: Esquema de renormalizacion para diferente numero de celdas
..82
Figura 21: Esquema del proceso de renormalizacion; ..83
Figura 22: Recorrido de la isoterma de agotamiento para un gas retrogrado
..88
Figura 23: Permeabilidad vs Radio para 1D ..105
Figura 24: Caso de un yacimiento multicapas ...106
Figura 25: Esquema de un pozo fracturado hidrulicamente .122

Figura 26: Variacin de la permeabilidad relativa al gas y de la permeabilidad
relativa al gas base dentro de la fractura ...130
Figura 27: Variacin de la permeabilidad relativa al gas y de la permeabilidad
relativa al gas base dentro de la matriz ..133
Figura 28: Porcentaje de error para los datos del Occidente de
frica..139
Figura 29: Porcentaje de error para los datos Govier-Fogarasi .141
Figura 30: Localizacin del rea del Campo Cupiagua 171

LISTA DE TABLAS

Tabla 1: Composiciones tpicas de los yacimientos de gas .4
Tabla 2: Parmetros de simulacin ...50
Tabla 3: Composicin del fluido .52
Tabla 4: Flujo fraccional de gas ( ) y el correspondiente gas condensado ( )
y la tasa de flujo de gas condensaco ( ) en Vol/Vol a todas las condiciones de
flujo ...114
Tabla 5: Propiedades del Campo Cupiagua ..172

RESUMEN

TITULO: METODOLOGA PARA OBTENER UN MODELO DE POZO QUE
PERMITA PREDECIR EL COMPORTAMIENTO DE LA PRODUCCIN DE UN
POZO DE GAS: CASO ESTUDIO CAMPO CUPIAGUA*

AUTOR: NEYLA YOHANA ARIAS HOYOS**

PALABRAS CLAVES: Yacimiento de Gas Condensado, Banco de Condensado,
Productividad, Saturacin Critica de Condensado, Permeabilidad Relativa, Modelo
de pozo.

DESCRIPCION:

Los yacimientos de gas condensado, cada da ganan ms importancia en el pas, mas ahora con el
retorno del Campo Cupiagua, se vuelcan las miradas hacia este tipo de yacimiento, que presentan
un comportamiento complejo debido a la existencia de un sistema fluido de dos fases cuando la
presin del yacimiento cae por debajo de la presin del punto de roco. El objetivo principal de este
estudio es desarrollar una metodologa que permita obtener un modelo de pozo que tenga en
cuenta los efectos del fluido composicional. Para tal fin, se realiza una revisin de algunos de los
modelo existentes para yacimientos de gas condensado, que permitirn entender y comprender
mejor el comportamiento de estos yacimientos.

En este documento se pueden encontrar conceptos tericos de los yacimientos de gas
condensado, como se forma el banco de gas condensado, las regiones de flujo que se presentan
en un yacimiento de gas condensado que se encuentra en deplecin, una recopilacin de trabajos
existentes relacionados con los yacimientos de gas condensado y de cmo los modelados,
finalmente se presenta una metodologa que puede servir de basa para realizar un modelo de pozo
que permita predecir el comportamiento de la produccin a nivel de pozo en un yacimiento de gas
condensado.

________________________________
* Proyecto de Grado.
** Facultad de Ingenieras Fisicoqumicas. Escuela de Ingeniera de Petrleos.
Tutor: Nicolas Santos Santos.
ABSTRACT

TITLE: METHODOLOGY TO OBTAIN A SINGLE WELL MODEL PREDICTING
THE BEHAVIOR OF THE PRODUCTION OF A GAS WELL: CASE STUDY
CUPIAGUA FIELD*

AUTHOR: NEYLA YOHANA ARIAS HOYOS **

KEY WORDS: gas field condensate, condensate Bank, Productivity, Critical
Condensate Saturation, Relative Permeability, Single Well Model.

DESCRIPTION:

Gas condensate deposits, daily gain more importance in the country, but now with
the return of Cupiagua Field, eyes are turning towards this type of deposit which
have a complex behavior due to the existence of a system of two fluid phases
when the reservoir drops below the dewpoint pressure. The main objective of this
study is to develop a methodology to obtain a model that takes into account well
the effects of compositional fluid. To this end, a review of some of the existing
model for condensate gas fields, which will understand and better understand the
behavior of these sites.

This document can be found theoretical concepts of deposits of gas condensate,
as the bank is formed of gas condensate, the regions of flow that are presented in
a condensate gas field is in depletion, a compilation of existing work related
deposits of gas condensate and the models finally presents a methodology that
can serve as a basis for a model to predict well the behavior of the production wells
at a gas field condensate.

________________________________
* Draft Grade.
** Faculty of Engineering Physical Chemistry. School of Engineering Petroleum
Tutor: Nicolas Santos.

INTRODUCCIN

Los yacimientos de gas condensado, cada da ganan ms importancia en el pas,
mas ahora que debido al bajo precio del petrleo, los yacimientos maduros y/o de
crudo pesado que resulten econmicamente no viables, posiblemente no se
puedan producir, se vuelcan las miradas hacia estos yacimientos de
comportamientos complejos, que an no son comprendidos en su totalidad, debido
a la presencia de un sistema fluido de dos fases cerca al pozo, cuando la presin
en esta zona cae por debajo de la presin del punto de roco.

Se considera que la presin en los yacimientos de gas condensado en el momento
de su descubrimiento, se encuentra por encima o cercana a la presin del punto
de roco, por lo que solo existe flujo de gas. Sin embargo, a medida que se lleva a
cabo la produccin, ocurre una cada de presin, que origina tericamente tres
problemas principales: (1) Reduccin de la productividad del pozo, (2) Disminucin
de la cantidad de gas para la venta y (3) Acumulacin de condensados que
bloquean la produccin de gas. Cuando la presin de fondo fluyendo cae por
debajo de la presin del punto de roco.

Este tipo de yacimiento tan particular, en los que las propiedades son muy
susceptibles a cambios durante la produccin. Tratar de modelarlos es un poco
complejo y una simulacin de todo el campo requerira un buen refinamiento en la
zona cercana a la cara del pozo, para describir la mejor forma del gradiente de
presin y los cambios composicionales, implicando esto un mayor tiempo de
cmputo y la implementacin de buenas tcnicas de refinamiento.
El objetivo principal de este estudio es desarrollar una metodologa que

permita obtener un modelo de pozo que tenga en cuenta los efectos del
fluido composicional. Para tal fin, se realiza una revisin de algunos de los
modelo existentes para yacimientos de gas condensado, que permitirn
entender y comprender mejor el comportamiento de estos yacimientos.

1. GENERALIDADES

1.1 CLASIFICACIN DE LOS YACIMIENTOS DE GAS

Los yacimientos de gas se dividen de acuerdo a su composicin en
yacimientos de gas seco, gas hmedo y gas condensado y sus
composiciones tpicas se observan en la Tabla 1.

1.1.1 Yacimientos de gas seco.

Como se observa en la Tabla 1, los yacimientos de gas seco contienen
principalmente metano ( ) con pequeas cantidades de pentano y
componentes pesados ( ). Debido al alto contenido de componentes
voltiles del gas seco, la condensacin del lquido solo se alcanza a
temperaturas bajo .

La temperatura de los yacimientos de gas seco es mayor que la temperatura
cricondentrmica y ni las condiciones de yacimiento ni las de superficie
entran en la regin de dos fases durante la cada de presin del yacimiento,
por lo que la mezcla de hidrocarburos siempre es gaseosa. Tericamente los
yacimientos de gas seco no producen lquidos en superficie, sin embargo, del
gas seco se pueden extraer cierta cantidad de lquidos por medio de
procesos criognicos (enfriamiento).
1.1.2 Yacimientos de gas hmedo

Como se observa en la Tabla 1, los gases hmedos se caracterizan por un
mayor contenido de componentes intermedios y pesados que los gases
secos. A condiciones de separacin en superficie la mezcla cae en la regin
de dos fases generando relaciones gas-liquido mayores que .
El lquido del tanque tiende hacer incoloro, con gravedad y un
contenido de liquido menor de .
Tabla 1: Composiciones tpicas de los yacimientos de gas.
COMPONENTE GAS
SECO
GAS
HMEDO
GAS
CONDENSADO

96.0 90.0 75.0

2.0 3.0 7.0

1.0 2.0 4.5

0.5 2.0 3.0

0.5 1.0 2.0

- 0.5 2.5

- 1.5 6.0

- 115 125
RGL, PCN/BN - 26000 7000
Liquido API
de
Tanque Color
- 60 55
- Incoloro
Amar. claro
Armar. Claro
Amarillo

Fuente: Rojas, Gonzalo Ingeniera de Yacimientos de Gas Condensado

Los yacimientos de gas hmedo, difieren de los gases condensados en que:
(1) no ocurre condensacin retrograda durante la cada de presin, (2) tiene
menos cantidad de componentes pesados y (3) la cantidad de lquido

condensado en el separador es menor.

1.1.3 Yacimientos de gas condensado

El gas condensado, es un fluido monofsico a condiciones de yacimiento
originalmente. Est compuesto principalmente por metano (>60% Tabla 1),
etano y otros hidrocarburos pesados. Bajo ciertas condiciones de presin y
temperatura, se convierte en un fluido bifsico, este comportamiento se
conoce como condensacin retrograda. La temperatura del yacimiento se
encuentra entre la temperatura crtica y la cricondentrmica de la mezcla
como se puede observar en la Fig.1. La condensacin retrograda se produce
durante la disminucin de la presin isotrmicamente en el yacimiento,
obtenindose en superficie la produccin de gas y liquido. El condensado
puede ser incoloro, amarillo, anaranjado, marrn, verdoso e incluso negro,
con una densidad entre y y una relacin gas condensado ( )
mayor que . El lquido condensado, en el yacimiento
generalmente no fluye o no se puede producir.

Durante la produccin de un yacimiento de gas condensado, las mayores
cadas de presin se presentan en las cercanas del pozo productor, por lo
que la formacin de lquido condensado se origina en primera instancia en
esa zona. El volumen de la fase liquida, no solo depende de la presin y la
temperatura sino de los componentes que conforman el gas, es decir, un gas
sin suficientes componentes pesados, no desprende liquido al caer la
presin, por el contrario un gas con suficientes componentes pesados si
generara liquido al caer la presin por debajo de la presin del punto de

roco.

Figura 1: Diagrama presin-temperatura para yacimiento de gas
condensado.

Fuente: yacimientos-de-gas-condensado.blogspot.com/2009/01/revison-de-los-yacimientos-de-gas.html

Con base a los componentes del gas condensado, este se puede clasificar
en dos:
Gas condensado pobre: genera un pequeo volumen de liquido,
aproximadamente menor a de liquido por de gas.

Gas condensado rico: genera un volumen de liquido superior a
por .
Figura 2: Diagrama presin-temperatura para un yacimiento de gas

condensado pobre.


1.1.4 Problemas asociados a los yacimientos de gas condensado.

Debido a las caractersticas especiales de los yacimientos de gas
condensado, en comparacin con los yacimientos de aceite negro, estos
yacimientos presentan ciertos problemas como:

Figura 3: Diagrama presin-temperatura para un yacimiento de gas

condensado rico.


1.1.4.1 Acumulacin de lquido en el pozo: la presencia de lquido en los
pozos de gas afecta negativamente sus caractersticas de flujo. Estos
lquidos pueden provenir de la condensacin de hidrocarburos (condensado)
o de agua producida conjuntamente con el gas. En ambos casos, la fase
liquida de alta densidad debe ser llevada a superficie por el gas, sino estos
se acumulan en el fondo del pozo, produciendo una contrapresin adicional
sobre la formacin que afecta la capacidad productiva del pozo. En pozos
con baja presin de fondo, el lquido puede matarlo.

Inicialmente, los pozos tienen presin y tasa de flujo elevadas, que impiden
la acumulacin de liquido en el fondo, pero a medida que la presin de fondo
disminuye y aumenta la produccin de liquido (agua y/o condensado), se
deben buscar mtodos artificiales para levantar ese fluido, los mas
recomendados son: levantamiento artificial con gas o gas lift, bombeo
mecnico cuando la presin de fondo fluyendo lo permite y hacer fluir el pozo
hacia la atmosfera para alcanzar altas tasas de flujo que levanten el liquido,
en la aplicacin de este ultimo mtodo se pueden producir incendios y
contaminacin ambiental.

1.1.4.2 Acumulacin de liquido en la formacin: una de las caractersticas
bsicas de los yacimientos de gas condensado es la formacin de
condensado retrogrado cuando la presin cae por debajo de la presin del
punto de roco, durante el agotamiento isotrmico de la presin del
yacimiento. Esta acumulacin de condensado retrogrado puede ocurrir en:

a. Zonas cercanas al pozo productor: se presenta cuando la presin
de fondo fluyendo es menor que la presin del punto de roco
( ) y la presin promedio del yacimiento es mayor o igual a la
presin de roco ( ). Un banco o anillo de condensado
retrogrado crece alrededor de un pozo de gas condensado, este
banco crece a medida que la presin del yacimiento cae,
disminuyendo la productividad del pozo y generando perdida de
componentes pesados en superficie. El comportamiento es rpido en
yacimientos de gas condensado rico, cercanos al punto de roco y
baja permeabilidad. En el peor de los caso, la saturacin de
condensado alrededor del pozo ( ) puede alcanzar valores de 50 a
10
60%, generando reduccin en la productividad de los pozos de 2 a 10

veces mayores que las existentes previamente.

b. Todo el yacimiento: se presentan cuando la presin promedio del
yacimiento cae por debajo de la presin del punto de roco retrograda
ocurre precipitacin de condensado en forma de pequeas gotas que
se acumulan en los canales porosos ms pequeos de la roca. Para
que el condensado retrogrado alcance cierta movilidad en el
yacimiento su saturacin debe ser mayor que su saturacin critica (30-
40%), lo que es difcil de alcanzar, generalmente el volumen del
condensado retrogrado no excede un valor de 20% vol.

11
12
2. BANCO DE CONDENSADO

Durante el proceso de produccin de un yacimiento de gas condensado, la
temperatura de formacin normalmente no cambia, pero la presin se
reduce. Las mayores cadas de presin tienen lugar cerca de los pozos
productores, cuando la presin se reduce hasta cierto punto, denominado
presin de saturacin o presin del punto de roco, una fase liquida rica en
componentes pesados se separa. La reduccin continua de la presin
incrementa la fase liquida hasta que alcanza un volumen mximo, luego el
volumen de liquido se reduce. Este comportamiento se ve claramente en un
diagrama presin temperatura, como se ilustra en la Fig. 4.

El volumen de la fase liquida presente depende de la presin, la temperatura
y la composicin del fluido. Una vez que los fluidos del yacimiento ingresan al
pozo, tanto las condiciones de temperatura como de presin pueden variar.
El lquido condensado puede producirse dentro del pozo; sin embargo,
tambin pueden acumularse en el fondo como resultado de los cambios
producidos en el pozo. Si el gas no tiene suficiente energa para transportar
el lquido a la superficie, se produce el retorno del lquido originndose un
aumento del nivel de lquido en el pozo que finalmente restringir la
produccin.

13
Figura 4: Formacin del bloque de condensado.

Fuente: Fan Li, Harris Billy W., J amaluddin A., Kamath J airam, Mott Robert, Pope Gary A., Shandrygin Alexander
and Whitson Curtis Hays. Revision de los yacimientos de gas condensado

2.1. MODELO TERICO

Wheaton y Zhang
1
desarrollaron un modelo terico para describir la
formacin del banco de condensado. Este modelo tiene como prioridad tratar
la dinmica del banco de condensado. Ellos definieron que la saturacin de
condensado alrededor del pozo productor, no solo depende de la presin,
1
WheatonR.J.andZhang,H.R.ChangesandCondensateAccumulationAroundProductionWells.
PaperSPE62930,Copyright2000.
14
sino que tambin se debe tener en cuenta la concentracin total de los

componentes pesados y como varan estos con el tiempo cerca del pozo.
Para entender la dinmica de la acumulacin de condensado y el potencial
de re-vaporizacin del condensado, es esencial entender el
comportamiento del pozo y del campo en general.

Para un volumen de porosidad media con lmites abiertos, la siguiente
ecuacin de conservacin composicional puede ser escrita para el i-simo
componente si la dispersin es despreciable:
(1)
Donde, es la porosidad, es la densidad molar del flujo total, es la
densidad molar del componente i-esimo, representa la fase y el vector
velocidad de Darcy en cada fase.

Usando la Ley de Darcy en trminos de coordenadas cilndricas y despus
de algunos procedimientos matemticos se puede llegar a:
(2)
Donde, es la fraccin molar total del componente i-esimo, es la tasa de
movilidad del componente i-esimo a travs de todas las fase y est dada por
la siguiente ecuacin.
(3)
Donde, es la viscosidad, es la permeabilidad y es la permeabilidad
relativa.
15
La movilidad est definida en la Ec. (2), en trminos de la funcin Chi se

puede igualar, a travs de la siguiente ecuacin.
(4)
donde, es la fraccin molar del componente i es la fase .

Para facilitar el anlisis de la Ec. (2), se puede generalizar as:
(5)
donde,
(6)
(7)

2.1.1 Anlisis del modelo

La Ec. 5 muestra el cambio en la fraccin molar total del componente con el
tiempo y a una distancia r del pozo durante la condensacin o reevaporacin
de condensado en un yacimiento. Esta ecuacin ofrece una nueva forma de
modelar el flujo de fluidos para un componente de gas condensado y
describe la misma fsica que las otras descripciones matemticas
convencionales de las componentes de conservacin de masa. Sin embargo,
se hace fcilmente un anlisis terico general para los cambios
composicionales durante el proceso de formacin del banco de condensados
16
el cual no necesita desarrollar una compleja simulacin composicional

completa.

fluyendo por encima de la presin del punto de roco.

Fuente: R.J . Wheaton and H.R. Zhang, BG International, Condensate Banking Dynamics in Gas Condensate Fields:
Compositional Changes and Condensate Accumulation Around Production Wells. Paper SPE 62930, (2000)

La saturacin y acumulacin de condensados cerca de la cara del pozo
depende de la composicin de la mezcla de componentes y de la presin en
el yacimiento. El cambio composicional depende de la variacin de las
propiedades del fluido (estas incluyen la viscosidad) as como de los
cambios en presin. Si el fluido del pozo de gas condensado tiene una
presin de fondo fluyendo constante considerando que la presin de roco
puede estar arriba o abajo (vea Fig. 5 y 6), es probable que sea muy
grande cerca del pozo (en especial en sistemas de baja permeabilidad) y
insignificante. Ms all del pozo, sin embargo, puede
convertirse en un trmino considerable y en insignificante.
17

fluyendo por debajo de la presin del punto de roco.

Fuente: R.J . Wheaton and H.R. Zhang, BG International, Condensate Banking Dynamics in Gas Condensate Fields:
Compositional Changes and Condensate Accumulation Around Production Wells. Paper SPE 62930, (2000)

En el primer termino de la Ec. (5), predominan los efectos de
almacenamiento (banking) cerca al pozo mientras que all los efectos del
segundo trmino son insignificantes. El segundo trmino es ms importante
en zonas lejanas del campo en donde se da el efecto de lquidos perdidos.
Esta ecuacin ofrece claramente una ventaja sobre cualquier otra ecuacin
composicional convencional por la separacin de los efectos la variacin de
la presin en tiempo y espacio en los cambios composicionales durante la
deplecin.

En la ecuacin, y son funcin de la movilidad del componente y
varan con la presin y la composicin. El trmino de movilidad relativa
aparece como una funcin de las permeabilidades relativas del gas y del
condensado, de la viscosidad de las dos fases y de la composicin molar de
18
cada fase. Por consiguiente es claro que es un parmetro crtico en el

comportamiento del banco de condensados. Adems, para las ecuaciones 5-
7, se ve ms afectado por las propiedades PVT de la mezcla de gas
condensado y est ms dominado por la movilidad total del fluido.

Los resultados mostrados en la Fig. 5 y 6 implican que, en un yacimiento
homogneo, el cambio de presin durante la deplecin se puede dar de la
siguiente manera:
y (8)

Wheaton y Zhang
2
concluyeron para un sistema de dos componentes
solamente que el incremento de la concentracin de los componentes
pesados cerca al pozo y el proceso de inhibicin era continuo. Una vez se
finaliza el proceso de incremento de pesados en esta regin debido al flujo
de fluidos en el pozo cuando la presin de fondo fluyendo est por debajo de
la presin de roco, no es posible reversar el proceso o intentar desplazar el
banco de condensados que se ha formado, acudiendo a mtodos como un
simple cierre del pozo para incrementar la presin por encima del punto
inicial de roco de la mezcla en el yacimiento.

Para yacimientos con baja permeabilidad, el gradiente de presin cerca al
pozo es generalmente grande, la tasa de crecimiento y de expansin del
banco de condensados ser relativamente alta, lo cual puede agravarse con
tasas ms altas de produccin.
2
Wheaton R. J . and Zhang H. R. Op.cit, 12
19
3. ENTREGABILIDAD DE POZOS DE GAS CONDENSADO

3

La entregabilidad del pozo (Well deliverability) se define como la relacin de
la tasa de produccin de un pozo en funcin de algunas limitaciones de
presin. Esta relacin de tasa-presin seria definida en la cabeza del pozo, y
de esta manera se define la relacin de entregabilidad de pozos, la cual
incluye todas las fuentes para las prdidas de presin desde yacimiento
hasta los separadores de superficie.

La contribucin de las diferentes fuentes de las prdidas de presin sobre la
entregabilidad del pozo depende principalmente de la permeabilidad del
yacimiento. Para yacimientos con una permeabilidad mayor de , las
prdidas de presin en el yacimiento son despreciables, comparadas con la
prdida total de presin. Para yacimientos con baja permeabilidad (menor a
) la entregabilidad del pozo se reduce debido a que el bloque de
condensado puede tener un efecto significativo en esta.

Los efectos del bloque de condensado depende de: (1) permeabilidades
relativas, (2) propiedades y (3) la forma en que el pozo es producido. Sin
embargo la reduccin de la entregabilidad del gas, debido al bloque de
condensado, solo es importante cuando la presin de fondo fluyendo ( )
alcanza un mnimo y el pozo se encuentra en declinacin. En el momento de
3
Fevang Oivind Gas Condensate Flow Behavior and Sampling. A Dissertation for the Partial
FulfillmentofRequirementsfortheDegreeofDoktorIngenior.DivisionofPetroleumEngineeringand
AppliedGeophysicsTheNorwegianInstituteofTechnologyUniversityofTrondheim.October(1995)
20
iniciar el desarrollo de esta investigacin, Fevang se baso en investigaciones

previas como las mencionadas a continuacin.

Muska
4
trato el problema del bloque de condensado en sus discusiones del
ciclaje de gas, donde el introduce un mtodo simple para estimar el radio del
bloque de condensado en funcin del tiempo, la tasa de gas y las
propiedades de la roca y el fluido. Fetkovich
5
uso los resultados de Muskat
para derivar la dependencia de la tasa y el tiempo sobre el dao por el
bloque con el fin de utilizarlo en la ecuacin estndar de la tasa de gas.

Kniazeff y Naville
6
y Eilerts
7
,
8
et al, fueron los primeros en modelar
numricamente la entregabilidad de pozos de gas condensado en
coordenadas radiales. Estos estudios muestran perfiles de presin y
saturacin radiales en funcin del tiempo y de otras variables operacionales,
confirmando que el bloque de condensado reduce el la entregabilidad del
pozo. Kniazeff and Naville tambin estudiaron el efecto del flujo no-Darcy (en
la fase gas) sobre la entregabilidad del pozo.

Gondouin et al
9
hicieron una significativa contribucin para entender el
fundamento de la entregabilidad del gas condensado. A travs de las
4
Muskat,M.PhysicalPrinciplesofOilProduction,McGrawHillBookCompany,Inc.(1949).
5
Fetkovich, M.J. The Isochranal Testing of Oil Wells. Paper SPE 4529 presented at the (1973) SPE
AnnualTechnicalConferenceandExhibition,LasVegas,Sep.30Oct3.
6
Kniazeff,V.J.andNaville,S.A.TwoPhaseFlowofVolatileHydrocarbonsSPEJ(1965)3744;Trans,
AIME,234.
7
Eilerts,C.K.,Sumner,E.F.andPotts,N.L.IntegrationofPartialDifferentialEquationforTransient
RadialFlowofGasCondensateFluidsinPorousStructuresSPEJ(1965)141152.
8
Eilerts, C. K., and Sumner, E. F. Integration of Partial Differential Equation for Multicomponent,
TwoPhaseTransientRadialFlow,SPEJ(1967)125135
9
Gondouin, M., Iffly, R. and Husson, J. An Attempt to Predict the Time Dependence of Well
DeliverabilityinGasCondensateFields,SPEJ(1967)112124;Trans,AIME,240.
21
simulaciones de aceite negro, ellos extendieron el trabajo de Kniazeff y

Naville, mostrando la importancia del bloque de condensado y los efectos del
flujo no-Darcy sobre el desempeo de la contra presin. Ellos tambin dan
procedimientos experimentales y medidas que cuantifican los efectos de la
permeabilidad relativa y el flujo no-Darcy multifasico.

ODell y Miller
10
presentaron la primera ecuacin de velocidad del gas
usando la funcin de psedopresin que describe el efecto del bloque de
condensado. La ecuacin es vlida cuando:
1. El wellstream producido corresponde al gas original del yacimiento, y
2. El radio del bloque de condensado es relativamente pequeo (la
presin en el yacimiento esta significativamente por encima del punto
de roco).

A partir de sus resultados, es claro que la entregabilidad del pozo puede ser
reducida para pequeas regiones de bloque de condensado.

Fussell
11
presenta simulaciones composicionales de las ecuaciones de
estado ( ) en coordenadas radiales, de pozos de gas condensado,
produciendo por deplecin de presin por debajo del punto de roco. El
muestra que el efecto de la acumulacin del condensado sobre la
productividad del pozo y que la ecuacin de ODell-Miller (con una pequea
correccin para contabilizar el gas disuelto en la fase de aceite fluyendo),
sobre predice las prdidas de entregabilidad debido al bloque de
10
ODell,H.G.andMiller,R.N.SuccessfullyCyclingaLowPermeability,HighYieldGasCondensate
Reservoir,JPT(1967)4147;Trans.,AIME240.
11
Fussell,D.D.SingleWellPerformancePredictionsforGasCondensateReservoir,JPT(1973)258
268,860870;Trans.,AIME255
22
condensado, comparado con resultados de simulaciones. El tambin evala

el efecto del comportamiento de fase y las caractersticas de la permeabilidad
relativa sobre el desempeo de la produccin.

J ones y Raghavan
12
,
13
estudiaron principalmente el comportamiento del
trasiente de presin (drawndown y buildup), de pozos radiales. Ellos utilizan
simuladores composicionales de las ecuaciones de estado, con mezclas
simples de gas condensado para tres componente ( ), la clave de
la observacin hecha sobre la entregabilidad del pozo, es que la funcin de
la pseudopresin presentada por Fussell es exacta durante todos los tiempos
de la deplecin. La integral debe ser evaluada usando presiones y
saturaciones conocidas como funcin del radio a un tiempo de deplecin
dado. Sin embargo esto no es muy til debido, a que ellos tienes que hacer
simulaciones composicionales para conocer las saturaciones y presiones a
un tiempo de deplecin dado.

Posteriormente a esa investigacin literaria, Fevang realizo una revisin de
los conceptos y ecuaciones para los yacimientos de gas seco y gas
condensado, a partir de esto tomo como punto de partida la ecuacin general
de tasa de gas en estado pseudoestable para cualquier geometra (radial,
pozos con fracturas vertical u horizontal) es:

12
Jones,J.R.andRaghavan,R.InterpretationofFlowingWellResponseinGasCondensateWells.
PaperSPE14204presentedatthe1985SPEAnnualTechnicalConferenceandExhibition,LasVegas,
Sept2225.
13
Jones, J. R. and Vo, D. T. and Raghavan, R. Interpretation of Pressure Buildup Response in Gas
Condensate Wells. Paper SPE 15535 presented at the 1986 SPE Annual Technical Conference and
Exhibition,NewOrleans,Oct58.
23
o
(9)
Donde

(10)

y para unidades de campo y para unidades SI.

La integral de presin incluye solo propiedades del fluido: viscosidad del gas
, factor volumtrico de formacin , (o equivalente factor de
comprensibilidad del gas ). La constante incluye propiedades bsicas del
yacimiento tales como la permeabilidad , espesor , radio de drenaje y
radio del pozo . El dao total efectivo es la suma del dao constante , y
una tasa dependiente del dao , es el coeficiente de la tasa
dependiente del dao.

La ecuacin de tasa de gas (Ec. 9) cambia de una simple relacin lineal
tasa-presin de la Ley de Darcy al modelo cuadrtico de Forchheimer
14
,
15
.
(11)
14
Golan,M.andWhitson,C.H.WellPerformance,2
nd
ed.PrenticeHall(1986)
15
Forchheimer,P.WasserbewegungdurchBodenZeitzVer.DeutschIng.(1901)17811788.
24
donde
(12)

La integral de presin en la ecuacin de tasa de gas recibe un nombre
especial, pseudopresin (Al-Hussainy et al.
16
), denominado utilizando el
smbolo actual SPE y definida por:
(13)

Sustituyendo , en la ecuacin de tasa de gas obtenemos,
(14)

3.1 ECUACIN PROPUESTA

Teniendo en cuenta todos los conceptos anteriores, el autor propone la
siguiente ecuacin general para pozos de gas condensado de cualquier
geometra (radial, pozos con fracturas vertical o horizontal) es, una formula
composicional,
16
AlHussainy,R.,Ramey,H.J.,Jr.,andCrawford,P.B.TheFlowofRealGasesThroughPorous
Media,JPT(1966)624.
25
(15)
La constante es la misma para el caso de la ecuacin tradicional de tasa
de gas, convierte el wellstream de base molar a una base en
volumen de gas estndar. determina la fraccin molar del wellstream que
es gas despus de la separacin en superficie. La primera parte de la
ecuacin representa la contribucin de la pseudopresin de la fase de gas
fluyendo en el yacimiento a la pseudopresin de la corriente de produccin
en el pozo. La segunda parte de la integral representa la contribucin de la
pseudopresin de la fase de aceite fluyendo en el yacimiento a la
pseudopresin de la corriente de produccin en el pozo.

En trminos de los parmetros de aceite negro, la ecuacin de tasa de
gas para gas condensado es
(16)
donde,
(17)

y son lo mismo como para la ecuacin general de tasa de gas, las
permeabilidades relativas

y son definidas relativas a la permeabilidad
absoluta y no relativa a la permeabilidad en saturacin de agua irreducible.

Las prdidas de entregabilidad debido al bloque de condensado son
particularmente importantes en yacimientos donde la parte principal de la
cada total de presin del pozo (yacimiento a separador) ocurre en el
26
yacimiento (usualmente en yacimientos de baja permeabilidad). La tasa de

produccin de los pozos petrolferos de yacimientos de baja permeabilidad
es tan baja que la cada de presin causada por los efectos del flujo no Darcy
pueden ser despreciables. Los pozos que producen de yacimientos de gas
condensado de baja permeabilidad estn usualmente fracturados, lo cual
elimina los efectos del flujo no Darcy.

3.1.1 Regmenes de flujo

El flujo de fluidos hacia un pozo productor de gas condensado en un
yacimiento que est en deplecin puede ser dividido en tres regiones
principales de flujo, extendiendo desde la cara del pozo hacia afuera.

Regin 1: una regin cercana a la cara del pozo, saturada con aceite y
gas, en la cual ambos estn fluyendo simultneamente.
Regin 2: una regin de acumulacin de condensado, donde el liquido
condensado esta inmvil y solo el gas esta fluyendo.
Regin 3: una regin que contiene solo una fase de gas.

Para unas condiciones dadas de produccin, una, dos o todas las tres
regiones podran existir. El desarrollo de las tres regiones de flujo podra
alcanzar una clase de estado pseudoestable. Las condiciones de flujo
pseudoestable, despus que se han establecido, cambian solo gradualmente
con el tiempo.

27
3.1.1.1 Regin cercana al pozo: El flujo composicional (GOR) dentro de la

Regin 1 se mantiene totalmente constante. Esto significa que el gas que
entra en la Regin 1 tiene la misma composicin que la mezcla del
wellstream producido. Por lo tanto si conocemos el wellstream producido,
entonces conoceremos el flujo composicional dentro de la regin. La Regin
1 es la principal fuente de prdida de entregabilidad en los pozos de gas
condensado, estas prdidas dependen principalmente de la permeabilidad
relativa al gas en esa zona y del tamao de la misma. La permeabilidad
relativa al gas en la Regin 1 es principalmente funcin de la distribucin de
la saturacin de lquido. La saturacin de liquido que se produce depender
de la riqueza del wellstream y de las propiedades PVT del yacimiento de gas
en la Regin 1. La distribucin de saturacin en la Regin 1 es determinada
especficamente para asegurar que todo el liquido condensado que entra,
tenga la suficiente movilidad para fluir a travs y fuera de ella sin ninguna
acumulacin. El tamao de la Regin 1 se incrementa gradualmente con el
tiempo. El crecimiento depende de las propiedades PVT del gas y la tasa de
produccin.

3.1.1.2 Regin de acumulacin de condensado: La Regin 2 define una
regin de acumulacin neta de condensado. El condensado formado en esta
regin tiene cero o muy baja movilidad, solo fluye gas. La acumulacin del
condensado es causada por dos mecanismos: (1) el condensado producido
por el yacimiento de gas, debido a la declinacin de la matriz del yacimiento
y (2) el condensado producido debido al gradiente de presin que se genera
en la corriente del yacimiento de gas en la Regin 2.

28
3.1.1.3 Regin de solo gas: La Regin 3 siempre existe en un yacimiento

de gas condensado que se encuentra subsaturado. La composicin en la
Regin 3 es constante e igual a la original del yacimiento de gas.

3.1.2 Coexistencia de las regiones de flujo

La Regin 1 siempre existe cuando la presin de fondo fluyendo es inferior a
la presin del punto de roco. Un periodo corto de transicin es requerido
para acumular saturaciones en estado estable en la Regin 1. El periodo de
acumulacin es ms largo para un pozo fracturado verticalmente, debido a
que el flujo lineal esta cercano al pozo, comparado con pozos radiales u
horizontales.

La Regin 2 siempre existe junto con la Regin 1 despus que la presin del
yacimiento cae por debajo de la presin del punto de roco, en este caso, la
Regin 3 puede desaparecer. No es posible que existan las Regiones 2 y 3
en ausencia de la Regin 1. Las tres regiones existen en yacimientos
ligeramente subsaturados y con presin de fondo fluyendo menor que la
presin del punto de roco. La Regin 2 puede desaparecer o tener un efecto
insignificante para yacimientos altamente subsaturados.

Para gas condensado muy rico, cercanos al punto crtico, la Regin 1 puede
existir en toda el rea de drenaje (en ausencia de las Regiones 2 y 3),
despus que la presin del yacimiento cae por debajo de la del punto de
roco.

29
3.1.3 Calculo de pseudopresin

Basado en las observaciones de las tres regiones para muchos sistemas de
gas condensado, se desarrollo un mtodo simple para calcular la integral de
pseudopresin en las ecuaciones 15 y 16. El enfoque es una extensin del
mtodo de pseudopresiones propuesto originalmente por Evinger y Muskat
17

para soluciones en pozos de gas de seco.

La integral de pseudopresin se dividi en tres partes, correspondiente a
cada regin de flujo.
Total
Regin 1
Regin 2 (18)
Region 3

Dada la produccin del y , se sabe inmediatamente , porque es
igual a la presin del punto de roco del wellstream producido. Usando las
propiedades de aceite negro, con

definido como la tasa de solucin de
aceite en el gas, se busca la presin en las tablas donde y se
define esta presin como . En un tratamiento composicional, la presin del
17
Evinger,H.H.andMuskat,MCalculationofTheoreticalProductivityFactor,Trans.,AIME(1942)
146,126139.
30
punto de roco del wellstream composicional producido, se define como . Si

, la integral de la Regin 1 debe ir solo desde hasta , en este
caso, las regiones 2 y 3 no existen.

3.1.3.1 Regin 1: en esta regin los lmites de presin varan desde la a
. En esta regin, las dos fases, gas y lquido condensado estn fluyendo.
En esta regin la integral de pseudopresin se resuelve modificando el
enfoque de Evinger-Muskat. A presiones las propiedades
, , , y se encuentran directamente de los datos para aceite
negro. La siguiente ecuacin define el de produccin para un sistema
con .
(19)

Esta expresin da el radio en funcin de la presin, con todas las
variables PVT conocidas en funcin de la presin.
(20)

La ecuacin 20 se puede expresar en trminos de volumen de aceite relativo
al wellstream producido durante la expansin a composicin constante,

(21)

31
El volumen de aceite relativo , puede ser expresado en trminos de las

propiedades PVT del aceite negro, combinando las ecuaciones 19 y 20.
(22)

Como se mostro por Evinger-Muskat, las permeabilidades relativas y
puede expresarse como una funcin de el radio cuando ambas fases
con mviles.

3.1.3.2 Regin 2: en esta regin los lmites se encuentran entre y ,
que es la presin del punto de roco inicial del yacimiento. En esta regin
estn presentes la fase liquida y gaseosa, pero solo fluye la fase gaseosa. La
existe cuando , la integral de la Regin 2 es evaluada usando
donde la saturacin de aceite es estimada como funcin de la presin a
partir de experimentos CVD. El volumen relativo de aceite
, la relacin presin-saturacin
. Si el valor de no se conoce de los datos
PVT para aceite negro, se calcula usando la siguiente ecuacin.
(23)
Donde,
(24)
(25)
32
Donde el subndice representa la presin en el estado actual y

representa la presin en el estado anterior.

3.1.3.3 Regin 3: los lmites de esta regin son desde a . No hay
lquidos condensados en esta zona, la fase de gas coexiste con el agua
connata. La permeabilidad de la fase gaseosa es constante e igual a
. Solo encontrando las propiedades se resuelve la integral en la
Regin 3. La tradicional funcin de pseudopresin para una sola fase puede
ser utilizada.


La ecuacin fue desarrollada, con base en el entendimiento fsico de las
condiciones de flujo en las tres regiones. Estas tres regiones de flujo se
presentan cuando los pozos de gas condensado producen con una presin
de fondo fluyendo por debajo de la presin del punto de roco. La Regin 1
tiene flujo composicional constante ( ), donde el gas y el condensado
fluyen simultneamente. En esta regin la perdida de entregabilidad es
causa por una reduccin de la permeabilidad relativa. Esta entregabilidad del
pozo puede ser determinada con exactitud mediante las Ecs.15 o 16.

La funcin de pseudopresin multifasica se puede calcular fcilmente a partir
del de produccin y las propiedades . El efecto de la reduccin de la
permeabilidad relativa al gas (Bloque de condensado) es incorporado en la
33
funcin de pseudopresin. La geometra del pozo, el factor dao, entre otros,

son tenidos en cuenta en la constante de productividad . Esta funcin es
calculada en tres partes, con base a las tres regiones de flujo.

3.2 MTODO SEMIANALTICO PARA PREDECIR LA ENTREGABILIDAD

La Fig. 7 muestra en detalles un diagrama esquemtico del mecanismo de
acumulacin de condensado cerca al pozo. La presin en la Regin III es
superior a la presin del punto de roco del yacimiento de gas, por lo tanto,
solo hay una fase de gas en esta regin. La composicin total es la misma
composicin del gas en esta regin. Una vez este gas entra a la Regin II, el
condensado sale y alimenta el crecimiento del banco. Puede que este gas
tenga una fase gaseosa y una fraccin de condensado asociado
dependiendo de la presin. Estas fracciones dependen del perfil de presin y
se puede obtener por el desarrollo de un clculo flash para la composicin
del alimento a diferentes presiones en la regin del banco. Los valores de
saturacin que se encuentran dentro del banco dependen de las
permeabilidades relativas as como tambin del flujo fraccional en el gas de
alimento. A cualquier presin dentro del banco de condensados, la fraccin
de volumen de gas y de aceite depende de y y sus composiciones y
pueden ser obtenidas por desarrollo de un clculo flash a la presin dada.

Para el gas de alimento dentro del banco de condensados, el balance molar
puede ser escrito como:
34
(25)
donde,
(26)

Figura 7: Diagrama esquemtico del mecanismo de condensacin cerca al
pozo.

Fuente: Chowdhury Nitin, Sharma Ravi, Pope Gary A., and Seperhrnoori Kamy, A Semi-Analytical Method to
Predict Well Deliverability in Gas-Condensate Reservoir SPE 90320. Copyright 2004.

A partir del detalle de las simulaciones de malla fina se ha encontrado que
las fracciones molares de los componentes en el gas y las fases
condensadas permanecen casi constantes incluso dentro del banco de
condensados. As, para el banco de condensados, un balance molar puede
ser escrito como:
35
(27)

Aunque las fracciones molares remanentes, y son casi las mismas
dentro del banco de condensados, las fracciones molares totales de los
componentes, dentro del banco son diferentes a los valores del gas
original. La acumulacin de los componentes ms pesados en el banco de
condensados resulta en una alta saturacin de aceite en esta regin
comparado con la fraccin del volumen de aceite en el gas de alimento. Este
cambio en la composicin permite un cambio en el diagrama de fase del
fluido como es mostrado en la Fig. 8.

Se ha desarrollado un mtodo semianaltico para predecir las tasas de
produccin de gas y aceite (condensado) de pozos con bancos de
condensados. El flujo cercano de estos pozos rpidamente alcanza un
estado estable o casi estable, cuando se puede asumir que tanto la presin
de fondo fluyendo ( ) como la presin promedio del rea de drenaje
pueden ser casi constante y se encuentre dominada por el comportamiento
de las permeabilidades relativas de las tres fase gas, aceite y agua cerca al
pozo.

Figura 8: Diagramas de fases para la composicin del fluido original y la
composicin del banco de condensado.
36

Las permeabilidades relativas de las fase de gas y aceite en el banco de
condensado, dependen de las saturaciones en el banco. Sin embargo, las
saturaciones son desconocidas a menos que se realice un grid de
simulacin. La solucin a este problema es usar la teora de flujo fraccional.

Para flujo radial, despreciando los efectos de capilaridad y gravedad y
usando la Ley de Darcy, se puede escribir el flujo de aceite como:
(28)
Similarmente, se puede escribir el flujo de gas como:
(29)
El flujo fraccional para cualquier fase es definido como:
(30)
37
Usando la definicin anterior y la Ley de Darcy se llega a:

(31)
Entonces, la relacin de permeabilidades relativas puede ser calculada a
partir de la siguiente ecuacin:
(32)
El flujo fraccional de la fase de gas y aceite puede ser calculado a partir de
los volmenes de gas y aceite del banco de condensado:
(33)
(34)

La Ley de Darcy se escribi en trminos de tasa de flujo molar, considerando
la tasa de flujo molar constante en estado estable incluso con densidad
variable, mientras que el volumen de flujo no.
(35)
Donde, en la tasa molar de flujo de gas y es la densidad
molar del gas .

La Ec. 35 puede ser integrado en un rango dado de presiones o a una
distancia radial para calcular la tasa molar de flujo de gas. La Fig. 9 muestra
un esquema del perfil del banco de condensado, con la aproximacin por el
38
mtodo semianaltico. La regin entre la cara del pozo y el radio equivalente

del pozo es dividido en nmeros pequeos de intervalos radiales, en la cual
es integrada la Ec. 35. Para cada uno de esos pequeos incrementos se
pueden aproximar como constantes las propiedades fsicas a una presin
conocida y la tasa de flujo molar es constante en estado estable. Este
procedimiento es anlogo al que usa un grid fino cerca al pozo en la regin
de inters y un grid con intervalos ms amplios en las regiones lejanas del
pozo. Integrando la Ec. 35 obtenemos:
(36)
(37)

Figura 9: Perfil tpico del banco de condensado con numero capilar junto con
la aproximacin por el mtodo semi-analtico.

39


Un grid fino de simulacin composicional fue construido utilizando el
simulador composicional de la Universidad de Texas UTCOMP. El modelo
semianaltico fue codificado dentro del simulador y comparado con un grid
fino de simulacin para un gas seco con un mximo abandono de liquido del
1%, con presin inicial del yacimiento de y de temperatura,
la presin del punto de rocio es de y un gas rico en condensado
con un mximo abandono de liquido del 20%, la presin inicial del yacimiento
es de y de temperatura, la presin del punto de rocio es de
, ambos fluidos se modelaron utilizando un grid fino con un angulo
de 36 y un grid basto similar a un grid cartesiano.
Los resultados se calcularon para el caso de un solo estrato y para mltiples
estratos, obteniendo los siguientes resultados:
Para un yacimiento con un solo estrato:
El grid fino y el grid semianaltico lograron captar el banco de
condensado cerca al pozo, mientras que el grid basto no muestra
ninguna acumulacin de condensado cerca al pozo.

La variacin de la viscosidad y la densidad molar del gas condensado,
que se presenta a medida que se alejan del pozo es muy pequea y
podra considerarse insignificante, esto es importante porque el
enfoque de pseudopresin, parte de la idea de que la variacin de la
40
densidad y la viscosidad es significativa y necesita ser tenidas en

cuenta para la exactitud de la tasa de flujo de gas.

La tasa de produccin obtenidas utilizando el mtodo semianaltico
son cercana a las del grid fino en estado estable y despus de
alrededor de 60 das. Para el caso del gas rico, el mtodo
semianaltico y el grid basto presentaron tasa de produccin cercanas
al grid fino en estado estable y despus de alrededor de 40 das.
Para mltiples estratos:
Los resultados del mtodo semi-analtico estn bien, de acuerdo con
los resultados del grid fino. Los resultados muestran que el mtodo
puede ser utilizado para predecir con exactitud las tasas de
produccin para cada estrato.

4. CLCULOS DE LA PRODUCTIVIDAD

La productividad en muchos pozos de gas condensado se reduce por la
formacin de un banco de condensados en la regin cercana al pozo.
Pronosticar la entregabilidad del pozo de gas condensado usualmente
requiere un grid fino de simulacin numrica para modelar la formacin del
banco de condensado y para darse cuenta de altas velocidades en
fenmenos tales como el flujo no-Darcy y los cambios relativos a la
permeabilidad de nmeros capilares grandes.

41
La productividad del pozo es un punto importante en el desarrollo de

yacimientos de gas condensado con bajas y medias permeabilidades. Sin
embargo obtener pronsticos exactos de la productividad puede ser difcil
porque necesitan comprender y contabilizar los procesos complejos que
ocurren en la regin cerca al pozo. Cuando la presin cae por debajo de la
presin del punto de roco, una regin altamente saturada de liquido se
acumula alrededor del pozo, alterando el flujo de gas y reduciendo la
productividad, por lo tanto, es de vital importancia tener en cuenta el efecto
del bloque de condensado, en los clculos de productividad. En la mayora
de los pozos de gas condensado el efecto neto del fenmeno de altas
velocidades es mejorar la productividad y reducir el deterioro debido a la
obstruccin de condensado, por lo anterior es importante incluir la alta
velocidad en la simulacin de pozos de gas condensado.

4.1 MODELO BALANCE DE MASA

Predecir el comportamiento de un yacimiento de gas condensado, se puede
predecir fcilmente, empleando un balance de masa combinado con una
relacin de inflow-performance relationship (IPR
18
). Las propiedades PVT del
fluido pueden ser calculadas a partir de un modelo de aceite negro ajustado,
donde propiedades tales como, la relacin gas-aceite y el factor volumtrico
de formacin son tabuladas contra la presin. Estudios de simulacin han
demostrado que no hay diferencias significativas entre los resultados de un
18
MOTT, Robert Engineering Calculations of Gas Condensate Well Production. Paper SPE 77551
42
simulador de aceite ajustado y uno composicional, para yacimientos de gas

condensado en deplecin
19
.

El IPR del pozo se puede basar en el mtodo de la integral de
pseudopresin, que puede tener en cuenta los efectos del bloque de
condensado y de las altas velocidades. El IPR del pozo es empleado para
calcular la mxima tasa de produccin de gas, para un valor particular de
presin promedio del yacimiento.

La integral de pseudopresin, para la produccin de pozos de gas est dada
por:
(38)
Donde, es el ndice de productividad, calculado de acuerdo al tipo de pozo,
y es la pseudopresin de las dos fases, definida por
20
:
(39)
El clculo de requiere las permeabilidades relativas y junto con
las propiedades PVT. El clculo de la integral de pseudopresin se divide de
acuerdo a las tres regiones, las cuales fueron resueltas en el Cap. 3.

19
ElBanbi, A.H. et al Producing Rich GasCondensate Reservoir Case History and Comparison
betweenCompositionalandModifiedBlackOilApproaches.PaperSPE58955,presentedattheSPE
InternationalPetroleumConferenceandExhibition,Villahermosa,Mexico,Feb13,2000.
20
Fevang, O. and Whitson, C.H.Modeling Gas Condensate Well Deliverability, SPERE, November
1996,p221230
43
4.1.1 Calculo de la razn de aceite fluyendo ( ) en la regin 1.

Un nuevo mtodo para calcular es basado en la estimacin del volumen
de la Regin 1. La Fig. 10, 11 y 12 muestran la tasa de flujo de aceite a
condiciones de superficie en funcin de la distancia al pozo. Estas muestran
que la tasa de flujo de aceite alcanza un valor mximo a una distancia
mayor de 2500 ft del pozo. La tasa de flujo de aceite es el flujo de aceite a
condiciones estndar, la cual se encuentra en dos fases en el yacimiento.
Para este estudio, el valor mximo ocurre cerca del lado externo de la
Regin 2, pero despus de la cada de presin por debajo del punto de roco,
la Regin 3 no es muy extensa y la tasa de flujo de aceite mxima se halla
en la mitad de la Regin 2.

Para calcular la expansin de la Regin 1, durante el periodo , se
considera que la presin del yacimiento declina de a
Figura 10: Saturacin de la fase aceite, alrededor de un pozo vertical sin
fracturas.
44

Fuente: MOTT, Robert Engineering Calculations of Gas Condensate Well Production. Paper SPE 77551

La cantidad de aceite necesario para expandir la Regin 1 en un delta de
volumen poroso es:
(44)
es la saturacin de aceite dentro de la Regin 1, estimada a partir de Ec.
20. es la saturacin de aceite en la Regin 2, que se asume es la misma
que la saturacin de un experimento de deplecin a volumen constante con
el fluido de yacimiento (corregido para permitir la presencia de agua
connata). El trmino dentro del corchete en la Ec. 44 representa la derivada
de la cantidad de aceite a condiciones de superficie con respecto al volumen
de la fase aceite en el yacimiento.
Figura 11: Relacin gas-aceite fluyendo en un pozo vertical sin fracturas.
45


Combinando las Ecs. 43 y 44, obtenemos una sola ecuacin con dos
incgnitas, el incremento en el volumen poroso y la presin en la Regin 1.
Esta puede ser solucionada si hay una manera de calcular el volumen poroso
dentro de una curva de contorno de presin. Esto se hace en trminos de la
relacin , definida como la fraccin total de la cada de presin que ocurre
dentro de la presin de contorno.
(45)

Figura 12: Tasa de flujo de aceite a condiciones estndar en un pozo
vertical.
46


4.1.2 Distribucin de presin alrededor del pozo

El clculo de la tasa de flujo de aceite requiere de un mtodo para estimar el
volumen poroso dentro de una presin de contorno dada. Adems, el rea de
una presin de contorno es necesaria para estimar la velocidad del gas a una
presin determinada. Este dato es obtenido de tablas de volumen poroso y
rea contra la fraccin de cada de presin . Estas tablas son calculadas
empleando algunas simplificaciones sobre la distribucin de presin
alrededor del pozo. Estos clculos asumen el flujo de una sola fase e ignoran
el cambio de propiedades con la presin.
4.1.2.1 Pozo vertical sin fractura: en este caso la presin de contorno es
cilndrica. Se asume que algn mecanismo de dao es debido a una zona
47
alrededor del pozo con permeabilidad diferente. Si el radio de esta zona es

, la relacin entre la presin y el logaritmo del radio, es modelado por dos
lneas rectas como se muestra en la Fig. 13. La pendiente de la lnea est
dada por en la zona externa y en la zona interna, donde:
(46)
(47)

Figura 13: Distribucin de la presin alrededor de un pozo sin fractura con
dao a la formacin.


48
4.1.2.2 Pozo vertical fracturado: Prats

21
da un esquema grfico de los
contornos del potencial de presin alrededor de un pozo fracturado con
conductividad infinita para una distancia adimensional ( , donde L es la
longitud promedio de la fractura) de 2.5. Los contornos de presin se asumen
elpticos. Fuera de esta regin, se asume que la presin vara linealmente
con el logaritmo del radio. Para cada valor de potencial, el eje mayor y menor
de la elipse (como las distancias adimensionales y ) son ledos en
la tabla derivada de la Fig. 4 del artculo de Prats. Entonces pueden ser
empleados para calcular el volumen dentro de la curva de contorno de
presin y el rea de dicho contorno.


El balance de materia y los modelos de productividad, pueden ser
combinados en una hoja de clculo de Excel para predecir el
comportamiento de un yacimiento de gas condensado. La hoja de clculo
puede modelar tres tipos de geometra, un pozo vertical sin fractura, un pozo
vertical fracturado y un pozo horizontal sin fractura. Los resultados obtenidos
se compararon con los resultados obtenidos de la simulacin composicional
utilizando un grid fino, consiguindose:

En el pozo vertical en un yacimiento de gas rico en condensado, el
modelo de simulacin de grid fino y la hoja de clculo del balance de
21
Prats,M.EffectofVerticalonReservoirBehaviorIncompressibleFluidCasa.SPEJJune(1961)105
118
49
materia, ignoraran los efectos de alta velocidad. Se realizaron dos

hojas de clculos una calculando la integral de pseudopresin por el
mtodo de CVDMB y otra por el nuevo mtodo de estimando del flujo
OGR. El mtodo de pseudopresin CVDMB sobrestima la
productividad del pozo, especialmente despus que la presin de
fondo fluyendo cae por debajo de la presin de roco. El nuevo
mtodo, el cual estima el flujo de OGR calculando el tamao de la
Regin 1, muestra mejores resultados de acuerdo con los resultados
de simulacin. Al momento de incluir los cambios de la permeabilidad
relativa con el nmero capilar se mejora la productividad, y los
resultados obtenindose con el clculo de la pseudopresin estn ms
cercanos a los resultados de la simulacin.

En un pozo vertical con fractura hidrulica, en un yacimiento de gas
seco, en este caso se utilizo un simulador 2D. La hoja de clculo sobre
estima ligeramente la productividad del pozo durante los tiempos
tempranos de produccin, pero la discrepancia es pequea
comparada otras incertidumbres en la prediccin del comportamiento
de un pozo fracturado.

En un pozo horizontal con baja permeabilidad en un yacimiento de gas
rico en condensado, se utilizo un simulador 3D. La hoja de clculo
muestra resultados muy cercanos a los obtenidos por la simulacin,
sin embargo la productividad es ligeramente sobre estimada en pozos
cortos.

El comportamiento de un pozo de gas condensado, puede ser
predicho con el modelo del balance de materia por deplecin del
yacimiento y por el mtodo de la integral de pseudopresin para el
influjo del pozo.
50

4.2 ESTIMACIN DE LOS PARMETROS DE CONTROL DE LA
PRODUCTIVIDAD EN POZOS DE GAS CONDENSADO

La capacidad para predecir la entregabilidad de un pozo es la clave para el
desarrollo de los yacimientos de gas condensado. La entregabilidad de un
pozo depende principalmente de la permeabilidad relativa en el punto final y
cerca a la cara del pozo, as como la permeabilidad del yacimiento. Esos
parmetros y el numero capilar se pueden obtener de los datos de presin
buildup utilizando la pseudopresin para una y dos fases simultneamente.
Adicionalmente esos parmetros pueden ser utilizados para estimar la curva
de permeabilidad relativa al gas
22
.

En los yacimientos de gas condensado, se forma un banco de condensado
cuando la presin de fondo fluyendo cae por debajo de la presin del punto
de roco. Esto crea tres regiones de diferente saturacin alrededor del pozo.
Cerca del pozo, una regin de alta saturacin de condensado que reduce la
permeabilidad efectiva al gas, disminuyendo la productividad. Esta
disminucin reduce las altas tasa de gas y/o bajas fuerzas capilares, que
reducen la saturacin del condensado en las inmediaciones del pozo, lo que
resulta en un aumento de la permeabilidad relativa del gas. Esto se llama
efecto del nmero capilar, acoplamiento positivo, extraccin de viscosidad o
extraccin de velocidad. Por el contrario altas tasa de gas inducen la inercia
(flujo turbulento o efectos no Darcy), que reducen la productividad. Por lo
tanto, la productividad es un equilibrio entre el nmero capilar y los efectos
inerciales.
22
Bozorgzabeh, M. and Gringarten, A. C. Estimating Productivity-Controlling Parameters in Gas/Condensate Wells
From Transient Pressure Data. Paper SPE 94018
51

Para predecir la entregabilidad en pozos de gas condensado generalmente
se utilizan simuladores numricos composicionales, pero estos simuladores
requieren un grid fino para modelar la formacin del banco de condensado
con mucha precisin. Los efectos del numero capilar y flujo no Darcy se
calculan a travs de correlaciones empricas, que requieren un numero
capilar base, y las curvas de permeabilidad absoluta y permeabilidad relativa.

4.2.1 Estudio de simulacin

La simulacin composicional se utiliza para identificar los parmetros claves
que determinar la entregabilidad del pozo. La simulacin genera un nuevo
pozo sinttico con datos de ensayo, que se interpretaran con el fin de
desarrollar una metodologa que permita obtener los parmetros utilizados en
la pseudopresin del gas seco y en los mtodos de pseudopresin para dos
fases en equilibrio. Se usa para que coincidan los datos de ensayo, primero
con los parmetros del yacimiento y a continuacin, con los parmetros
generados, para comparar los resultados y comprobar la metodologa.

4.2.1.1 Simulacin establecida. Se utilizo un simulador composicional
comercial (Eclipse 2003), con la ecuacin de Peng-Robinson (PR) de las
ecuaciones de estado (EOS), para modelar el comportamiento de un pozo
1D, con una sola capa y yacimiento de gas condensado radial. Se realizo un
anlisis de sensibilidad sobre el tamao de la celda y tiempo de paso
(timesteps) para evitar errores de truncamiento por discretizacion espacial y
52
temporal y as obtener el gradiente de presin lineal y el perfil de saturacin

en una escala semilog. Los tamaos de las celdas cerca a la cara del pozo,
se disminuyeron para aumentar la precisin en el radio exterior del banco de
condensado. El radio externo de aproximadamente 12000 pies es
significativamente ms grande que el rea tpica de drenaje en un pozo de
condensado, para garantizar que los efectos de frontera no se sientan en el
pozo. Los parmetros utilizados para la simulacin se muestran en la Tabla
2.

Las propiedades del fluido son del campo en Mar del Norte de gas
condensado, con un mximo de abandono de liquido del 2.1 vol% y una
presin de saturacin de . La Tabla 3 muestra la composicin del
fluido.
Tabla 2: Parmetros de la simulacin

PARMETROS VALOR
Presin inicial del yacimiento (psia) 3600
Espesor de la formacin (ft) 100
Porosidad promedio (%) 10
Permeabilidad absoluta (md) 10
Radio neto de la red 1
Saturacin de agua connata (%) 20
Radio del pozo (ft) 0.25
Compresibilidad del agua (1/psi) 3.0E-6
Compresibilidad de la roca (1/psi) 4.0E-6

53
Fuente: Bozorgzabeh, M. and Gringarten, A. C. Estimating Productivity-Controlling Parameters in Gas/Condensate

Wells From Transient Pressure Data. Paper SPE 94018
Los parmetros PR EOS fueron ajustados por regresin para obtener una
buena correlacin entre los valores reales y los predichos.

Las curvas de imbibicin de permeabilidad relativa usadas en la simulacin
se basan en las expresiones de Corey. (Liu et al. 2001):
(48)
(49)
El agua connata se asume inmvil y la saturacin de aceite residual se
asume cero. La saturacin critica de gas tambin se asume cero porque su
impacto en las curvas de permeabilidad relativa al gas y la saturacin de
aceite es insignificante comparada con el exponente de la permeabilidad
relativa de Corey para el gas Las carreras de simulacin fueron
diseadas para mostrar el impacto de los coeficiente de la formula de Corey
sobre la productividad de un pozo y as identificar los otros parmetros de
control.

4.2.2 Influencia de la permeabilidad relativa en la
productividad de un pozo.

A travs de un estudio de simulacin numrica se identificaron los efectos de
la permeabilidad en relacin a la saturacin de lquido irreducible, la
54
permeabilidad efectiva al gas cerca a la cara del pozo a la saturacin de

aceite y el exponente de Corey para el gas en la productividad.
Tabla 3: Composicin del fluido.
COMPONENTE Mol%

1.58

2.41

79.6

6.87

3.57

1.89

0.89

2.62

0.46

0.12

0.0026

Wells From Transient Pressure Data. Paper SPE 94018.

4.2.2.1 Permeabilidad relativa al gas a la saturacin de lquido
irreducible. Este estudio demuestra que la productividad es controlada por
la permeabilidad relativa al gas a la saturacin de liquido irreducible, ,
55
adems de la relacin entre y , segn Fevang y Whitson

23
y que
la reduccin del valor de , disminuye la produccin.

4.2.2.2 Permeabilidad efectiva al gas a la saturacin de aceite cerca a la
cara del pozo. Para analizar este efecto se, utilizaron dos modelos. Los
resultados mostraron que el modelo con menor a la , posee menor
productividad, por lo tanto, una baja permeabilidad relativa al gas a la
saturacin de aceite critica, disminuye la productividad.

4.2.2.3 Exponente Corey del gas. La productividad es mayor en un modelo
con menor exponente, debido a que un menor valor del exponente,
corresponde a un sistema ms miscible, con mayor permeabilidad relativa al
gas en la zona de saturacin inmvil, incrementando la productividad.

Se puede afirmar que y a la saturacin cerca a la cara del pozo y el
exponente de Corey son los parmetros claves para la funcin de
permeabilidad relativa y deben definirse correctamente.

4.2.3 Estimacin de la permeabilidad relativa.

23
Fevang,O.andWhitson,C.H.ModelingGasCondensateWellDeliverability.SPE30714.(1996)
56
La simulacin composicional con y sin efectos del nmero capilar se utilizo

para generar datos de pozos sintticos, que ilustran los cambios de la
funcin derivada en el perfil de saturacin cuando la pseudopresin de una
sola fase es utilizada. Estos datos sintticos fueron analizados con las
formulas de pseudopresin de una y dos fases para establecer una teora
justificada para la estimacin de los parmetros de permeabilidad relativa y el
numero capilar base. Los efectos no-Darcy no son incluidos en este anlisis,
ya que solo afectan el dao de la cara del pozo y no la estabilizacin de
lneas de flujo radial.
4.2.3.1 Formulas de pseudopresin para una y dos fases. Cuando la
presin de fondo fluyendo cae por debajo de la presin del punto de roco,
ocurre la condensacin retrgrada, y un banco de liquido se forma alrededor
del pozo. Este banco se expande a medida que aumenta el tiempo de
produccin. La interpretacin de los datos del trasiente de presin de una
prueba de pozo asume que la ecuacin de difusin lineal, es la verdadera
para una sola fase y un fluido ligeramente compresible. En los yacimientos
de gas seco y yacimientos de gas condensado, la pseudopresin de una y
dos fases deben ser utilizadas.

4.2.3.1.1 Pseudopresin una sola fase. La ecuacin de difusividad puede
ser linealizada por el gas seco para reemplazar la presin por la funcin de
pseudopresin del gas seco para una sola fase:
(50)
donde, es el factor compresibilidad del gas real.

57
El uso de la pseudopresin de una sola fase en yacimientos de gas

condensado equivale a considerar el gas como el fluido dominante y el
condensado como el fluido heterogneo. Pero no se usa la pseudopresin,
se utiliza una pseudopresin normalizada:
(51)

La Fig. 14 muestra tres lneas de estabilizacin; la estabilizacin inicial
corresponde a la mejora de la permeabilidad relativa al gas debido a los
efectos del numero capilar; la segunda estabilizacin, varia con el tiempo de
produccin y depende de la saturacin de aceite en el banco de condensado:
a mayor saturacin de aceite, menor permeabilidad relativa al gas y mayor
nivel de estabilizacin de flujo radial; la tercera estabilizacin en el mismo
nivel para todos los tiempos de produccin, porque corresponde a la
movilidad del gas en la regin del yacimiento cuando la presin est por
encima de la presin del punto de roco. El comienzo de la segunda lnea de
estabilizacin de una acumulacin de derivados, corresponde al final de los
efectos del numero capilar, por lo tanto, el tiempo correspondiente se puede
utilizar para calcular el nmero capilar base.

Figura 14: Grafica Log-log de derivadas de pseudopresin normalizada para
10 das de agotamiento con diferentes tiempos de produccin.
58

4.2.3.1.2 Pseudopresin dos fases. La ecuacin de difusividad para un
yacimiento de gas condensado, puede ser linealizada en una funcin de
pseudopresin de dos fases:
(52)
El clculo de la integral anterior, es difcil porque requiere el conocimiento de
la permeabilidad relativa al gas que es funcin de la presin, varios mtodos
han sido propuestos para tal fin.

Para evaluar con precisin la permeabilidad de la formacin se puede
calcular la pseudopresin de dos fases, utilizando los datos de simulacin de
pozos, primero sin efectos del numero capilar e investigando el impacto de
59
, , , y el dao en las derivadas de pseudopresin de dos fases de

los niveles de estabilizacin.

Las pseudopresiones de dos fase se calcularon solo para la primera regin
de movilidad, donde fluyen gas y condensado juntos en el pozo. En esta
regin ambas fases se encuentran en estado de equilibrio. Se utilizaron tanto
la relacin gas/aceite ( ) a la presin en el momento del cierre como la
sugerida por Raghavan (1976) y la presin del banco , para calcular
de la Ec. 20.

La Fig. 15 muestra la estabilizacin de la derivada de los datos que incluyen
los efectos del nmero capilar, , esta al mismo nivel que la estabilizacin
final de la derivada de la pseudopresin para una fase y corresponde a la
permeabilidad absoluta. Por otra parte dos estabilizaciones de la derivada
pueden verse en los datos que no incluyen los efectos del nmero capilar,
. La primera estabilizacin, corresponde a la velocidad de despojo
(stripping) de la zona, que origina una permeabilidad mayor que la
permeabilidad absoluta, y por lo tanto es incorrecta. La segunda
estabilizacin, es la regin del banco de condensado, donde el numero
capilar es inferior al nmero capilar base que corresponde a la permeabilidad
absoluta, porque el numero capilar no tiene ningn efecto sobre la
permeabilidad relativa al gas en esta regin.
Figura 15: Grafica Log-log de datos de simulacin con y sin efectos del
numero capilar.
60


4.2.4 Generacin de una curva de permeabilidad estimada.

El mtodo ms comn para la generacin de curvas de permeabilidad
relativa es a partir de datos de laboratorio, en estado estable. Sin embargo,
la permeabilidad relativa tambin puede determinarse del desempeo del
campo o de los datos del trasiente de presin.

61
Basados en los resultados anteriores de la simulacin, se propone el

siguiente procedimiento para generar una curva de permeabilidad estimada
como respuesta al incremento de presin, usando las pseudopresiones de
una y dos fases:
1. Determinar la movilidad de una y dos fases del anlisis
del incremento de presin, utilizando las pseudopresiones de una y
dos fases o la interpretacin de pruebas de pozo para un modelo
compuesto por tres regiones radiales
24
.
2. Determinar la relacin de permeabilidad relativa de una y
dos fases, proporcionando la movilidad calculada en el paso 1 y la
relacin de viscosidad.
(53)
3. Calcular la saturacin del gas a la presin en el momento
de cierre, usando la relacin presin-saturacin de Boe, que predice la
saturacin en el pozo durante el cierre:
(54)
donde, es la movilidad del gas y la movilidad total, es la
compresibilidad total, es el factor volumtrico de formacin de la
derivada parcial con respecto a de la constante .
4. Usando la funcin de Corey para la permeabilidad
relativa al gas, para calcular , utilizando la relacin de permeabilidad
relativa calculada en el paso 2 y la saturacin de aceite al momento
del cierre del paso anterior.
(55)
24
Bozargzadeh,M.andGringarten,A.C.CondensateBankCharacterizationFromWellTestData
andFluidPVTProperties.PaperSPE899004PA.(2006)
62
5. Seleccione un valor para .

6. Cree una tabla de - vs- saturacin, a partir de la Ec.
55.
7. A partir de la Ec. 54 cree una tabla de saturacin vs-
presin y genere una tabla -vs- presin.
8. Calcule para cada uno de los valores de presin,
a partir de la Ec. 20 y a la presin del banco, .
9. Calcule , dividiendo de para cada valor de
presin.
10. Calcule la permeabilidad absoluta:
(56)
Utilizando la permeabilidad efectiva a la condensacin inicial, obtenida
del anlisis de pseudopresin para una sola fase y el estimado
en el paso 5.
11. Se genera las pseudopresiones de dos fases a partir de
la Ec. 18 y y , calculadas entre el paso 9, analizar los datos
obtenidos para la permeabilidad absoluta, si no es igual a la obtenida
en el paso 10, repetir desde el paso 5 hasta el 1, hasta que convergen
las permeabilidades absolutas.

4.2.5 Anlisis del modelo.

Para verificar este modelo se utilizo un simulador numrico composicional y
los resultados obtenidos se compararon con una prueba Drillsteam (DST) a
un yacimiento de gas condensado pobre del Mar del Norte, concluyendo:
63

La curva de permeabilidad relativa utilizada en el
simulador composicional, est basada en una saturacin de agua
connata del 20%. El punto final de la curva permeabilidad relativa al
gas fue obtenido correctamente y es consistente de acuerdo al
obtenido de la prueba.

La permeabilidad estimada y la permeabilidad absoluta,
fueron sustituidas en el modelo de simulacin y se genero una nueva
historia de produccin. Esta nueva simulacin est bien de acuerdo
con la simulacin previa. Esto confirma la capacidad de proponer
tcnicas que generan la permeabilidad absoluta y la permeabilidad
relativa estimada, las cuales pueden ser utilizadas en la prediccin de
la entregabilidad del pozo en yacimientos de gas condensado.

4.3 PRODUCTIVIDAD DE LAS FRACTURAS EN POZOS DE GAS
CONDENSADO

La productividad de los pozos de gas condensado puede ser reducida
drsticamente como resultado de la acumulacin de lquido alrededor del
pozo durante el agotamiento. Este fenmeno ha sido descrito para pozos sin
fracturas. Para el anlisis de proyectos de ingeniera relacionados con la
presin de agotamiento, es necesario poder representar este fenmeno en
un modelo numrico. Uno de los mtodos para calcular la reduccin del
ndice de productividad (PI) en un simulador composicional fue descrito por
64
Flores
25
. Sin embargo, hasta hace poco no era posible estimar el efecto de la
fractura de los pozos debido a la complejidad del problema. Carlson y Myer
26

presentaron un estudio del efecto del condensado en sistemas de baja
permeabilidad y fracturados.

Para el modelado del Campo Smorbukk un gran campo en la parte norte de
Noruega y del Mar del norte. Hace parte del desarrollo Asgard, que
comprende tres campos: Midgard, un campo de gas pobre, Smorbukk Sor un
campo de gas condensado y Smorbukk que es el ms profundo y tiene
mayores reservas de liquido de los tres. Se han estimado las reservas de
Asgard en de gas y de hidrocarburos
lquidos. Se utiliz el software de fracturas Simtech junto con un modelo de
aceite negro extendido, permiten modelar un sistema con fracturas. Las
propiedades del gas condensado, se manejan en el simulador como un
aceite negro condensado opcin propuesta por Spivak y Dixon
27
en la que
se asume que la fase vapor contiene el gas y el vaporizado aceite,
mientras que la fase de lquidos hidrocarburos es el condensado y el disuelto
es el gas. El equilibrio de fase es entonces aproximado por la presin que
depende de las funciones y .
La simulacin con fracturas es inherentemente mas difcil, mientras el modelo
con fracturas podra ser la interfaz por lo menos de un simulador
composicional, este enfoque podra significar un aumento en el tiempo de las
25
Flores, J. Compositional Simulation to Develop an Optimum Gas Cycling Scheme Kaybob
BeaverhillLakeCpool,Alberta.PaperNo.782931presentedatthe(1978)AnnualMeetingofthe
PetroleumSocietyofCIM,Calgary,1316June.
26
Carlson, M. R. and Myer, J. W. G. The Effects of Retrograde Liquid Condensation On Single Well
Productivity Determined Via Direct (Compositional) Modeling of Hydraulic Fracture in a Low
PermeabilityReservoir.PaperSPE29561presentedatthe(1995)SPERockyMountainRegion/Low
PermeabilityReservoirSymposium,DenverColorado,1922March.
27
Spivak,A.andDixon,T.N.SimulationofGasCondensateReservoirs.PaperSPE4271presented
atthe(1973)SPENumericalSimulationSymposium,Houston,1112January.
65
corridas y problemas de estabilidad. La eficiencia computacional de las

propiedades PVT del aceite negro permite el estudio tiempo/costo efectivo en
la influencia de varios parmetros.

Debido a la complejidad y un gran nmero de importante parmetros
potenciales, el problema estudiado lo dividieron en varias fases. En la
primera fase, se desarrollo la aproximacin del aceite pseudonegro y su
exactitud se verifico por una comparacin detallada del simulador
composicional realizadas por Statoil. Despus de la viabilidad de los dos
componentes modelados se estableci, el efecto de la estimulacin sobre el
ndice de productividad (PI) y su estudio mediante fracturas geomtricas
idealizadas, yacimientos heterogneos con modelos basados en una
interpretacin petrofsica. Posteriormente a la confirmacin de la viabilidad
de las fracturas, simulaciones se llevaron a cabo con fracturas realistas
resultantes de simular el bombeo con diferentes tipos de tratamientos de
estimulacin.

4.3.1 Modelo representativo del yacimiento

Se utilizaron dos de las 18 capas del yacimiento para la construccin del
modelo. Debido a la falta de datos de permeabilidad vertical, asumieron
en todas partes. Debido a que se espera que la tensin
interfacial (IFT) tenga un gran efecto sobre el movimiento del condensado,
utilizaron dos escenarios detallados para el modelo.
Se defini la saturacin crtica de aceite , o alternativamente la saturacin
de aceite residual , es una funcin de la saturacin de agua inicial ,
66
que se supone que es la misma saturacin de agua connata (es decir

. La permeabilidad relativa al aceite a .
Todas las otras curvas de permeabilidad se comportan como se describe a
continuacin:
1. A baja tensin interfacial (< 0.5 min/m), o
y es una lnea recta desde hasta
.
2. Alta tensin interfacial (>0.5 min/m), o
y es una parbola entre hasta
.

El modelo representa un grid de uno a cuatro elementos de simetra en
coordenadas cartesianas, con un rea de 1107.7 x 1107.7m (equivalente a
un radio de drenaje de 1250m). Utilizaron un grid de 16 x 10. Todos los
modelos se iniciaron en equilibrio gravedad/presin capilar y con saturacin
de agua inicial.

4.3.2 Presentacin de los resultados en trminos del ndice de
productividad

Cada simulacin proporciona una gran cantidad de informacin, usualmente
se utiliza una simple medida para comparar diferentes casos. En el este caso
la medida de comparacin utilizada es el ndice de productividad (PI), cuya
definicin clsica es:
67
(57)
Esta ecuacin asume una sola fase de flujo, yacimientos de una sola capa y
en estado de equilibrio. Las dificultades que surgen al aplicar la ecuacin
anterior son flujo en varias capas y zona heterognea. En primer lugar, la
presin del radio externa (P
e
), puede variar considerablemente entre capa y
capa debido al agotamiento desigual. En segundo lugar el tiempo para
alcanzar el estado de equilibrio en las diferentes capas vara en orden de
magnitud en funcin de las propiedades de cada capa. Por lo tanto el
sistema siempre est en flujo transitorio. A partir de todo lo anterior se usa
una definicin nueva llamada ndice de productividad (PI) instantneo
(58)
donde es la presin promedio del yacimiento durante todo el tiempo t.

4.3.2.1 Modelo de fracturas hidrulicas: Dos tipos de fracturas se
utilizaron es el modelo. Una fractura idealizada representada por una fractura
rectangular de altura y anchura constante, que abarca todo el espesor bruto
del yacimiento. Una fractura realista representada en un modelo 3D de altura
y anchura variable, tal como lo predice el simulador de fractura. Ambas
integradas se integraron a un modelo de yacimiento con transmisibilidad
interna en el grid.


Los modelos homogneos muestran alrededor un valor de ndice de
productividad superior en un 50%. Esto demuestra que la heterogeneidad
68
juega un papel significativo y en el momento de simular, utilizar las capas fina

son necesarias para obtener respuestas precisas.

La necesidad de fracturas despus de la cada de presin puede resultar del
desarrollo de un programa o una estrategia deliberada. El diferencial de
deplecin crea cambios de stress sustanciales, que tienden a contener las
fracturas en las capas ms repletadas. Asumiendo que cada capa acta de
manera independiente, los cambios de stress en cada capa pueden
estimarse a partir de:
(59)

La heterogeneidad disminuye los beneficios relativos de la fractura. Las
tcnicas de anlisis y desarrollo de modelos numricos del ndice de
productividad en funcin del tipo de fluido, la heterogeneidad del yacimiento,
caractersticas de flujo multifase, conductividad de las fracturas, duracin y
contencin y turbulencia del gas, permiten tener un mejor conocimiento del
ndice de productividad y as utilizarlo de la manera ms adecuada.

69

5. PERMEABILIDAD RELATIVA

Muchos pozos de gas condensado muestran una disminucin significativa en
la productividad cuando la presin cae por debajo de la presin del punto de
roco. La causa ms aceptada para esta disminucin del ndice de
productividad ( ) es la disminucin de permeabilidad relativa al gas, debido
a la acumulacin de condensado en la regin cercana a la cara del pozo.
Las predicciones del rendimiento del flujo, requiere de modelos precisos de
permeabilidad relativa al gas como funcin de la tensin interfacial ( ),
entre las fases del gas y del condensado.

Tradicionalmente se ha asumido que la permeabilidad relativa es funcin
principalmente de: tamao y distribucin de los poros, mojabilidad de la roca,
saturacin, historia de saturacin y tensin interfacial ( ), la mayor
incertidumbre es el efecto del numero capilar ( ). Recientes
publicaciones indican que la permeabilidad relativa mejora a altos nmeros
capilares. Adems, la vaporizacin del agua cerca a la cara del pozo puede
tener un efecto pronunciado en la permeabilidad relativa gas/aceite.

Debido a que la permeabilidad relativa depende de la composicin del fluido,
la presin del pozo y la saturacin de agua y del condensado, un modelo es
necesario tanto para la interpretacin de datos de laboratorio y los predichos
70
en procesos de simulacin. La clave de un modelo de permeabilidad relativa

al gas depende de la saturacin crtica del condensado y del nmero capilar.

5.1 MODELO DE PERMEABILIDAD RELATIVA

El problema fundamental con la acumulacin del condensado en el
yacimiento es que las fuerzas capilares pueden mantener el condensado en
los poros a menos que las fuerzas de desplazamiento del condensado
exceden las fuerzas capilares. En la medida en que presin en la fuerza de
desplazamiento de la fase gaseosa y la fuerza de flotabilidad en el
condensado, pueden ser mayor que la fuerza capilar en el condensado, la
saturacin del condensado se reducir y la permeabilidad relativa al gas
aumentara.

Varias definiciones del nmero capilar se han hecho, los ms comunes se
han escrito en trminos de la velocidad de desplazamiento de los fluidos,
mediante el uso de la Ley de Darcy para reemplazar el gradiente de presin
por la velocidad. Sin embargo es la fuerza del fluido atrapado, la definicin
ms fundamental, y est definida como:
(60)
Donde el gradiente de el potencial de flujo esta dado por
es el desplazamiento de fase. En algunos casos, las
fuerzas de flotabilidad pueden contribuir de manera significativa en la fuerza
total de la fase atrapada. Para cuantificar este efecto, se introduce el nmero
de Bond definido como:
71
(61)
Para casos especiales como flujo vertical, los vectores fuerza son colineales
y solo puede sumarse los valores escalares de la viscosidad y la fuerza de
flotabilidad y correlacionar la saturacin residual de aceite con esta suma, o
en algunos casos esta fuerza es insignificante comparada con otras fuerzas y
tanto el numero capilar o el numero de Bond pueden ser utilizados por ellos
mismo. Sin embargo, en general, las fuerzas de la fase atrapada no son
colineales en el flujo del yacimiento y el vector suma debe ser utilizado. De
hecho, las fuerzas viscosas son frecuentemente casi horizontales y la fuerza
de flotabilidad es vertical, de este modo el escalar suma es un error. Una
generalizacin del nmero capilar y el nmero de Bond, obtenida por J in
28
,
llamada numero de captura, que para la fase desplazada por la fase se
define como:
(62)
Esta definicin explcitamente no tiene en cuenta los efectos de expansin y
mojabilidad en la captura de la fase residual. Por lo tanto, se debe utilizar la
correlacin de saturacin residual junto con la gran variedad de fases y
diferentes mojabilidades y la gran variedad de tipos de roca. La saturacin
residual modelada con base al nmero de captura est dada por:
(63)
Donde, es la saturacin residual de la fase y es la saturacin de la
fase . Los subndices alto y bajo se refieren a alto y bajo nmero de captura.
, es cero generalmente. y son los parmetros de captura, obtenidos
28
Jin,Minquan.AstudyofNonAqueousPhaseLiquidCharacterizationandSurfactantRemediation
Ph.D.Dissertation,U.ofTexas,Austin,1995.
72
del ajuste de los datos de saturacin residual de la fase . generalmente es

1.0. La transferencia de masa puede reducirse al valor de para valores
menores que .

Establecer la correlacin de saturacin residual con el numero de captura, es
el primer y fundamental caso paso para correlacionar datos de permeabilidad
relativa como funcin de la tensin interfacial. El siguiente paso es
correlacionar el punto final de la permeabilidad relativa de cada fase, que se
pueden correlacionar a partir de la siguiente ecuacin:
(64)
Donde es la saturacin residual de la fase conjugada, el condensado es
la fase conjugada del gas. El paso final es calcular la permeabilidad relativa a
cada fase como funcin de la saturacin. Una solucin a este problema es
asumir la simple funcin de la permeabilidad relativa de Corey
29
. Entonces se
requiere una correlacin de los exponentes de Corey con el nmero de
captura. La Ec. 64 est escrita en trminos del exponente ms que del punto
final y puede ser utilizada para este fin. Sin embargo, no todos los datos de
permeabilidad relativa pueden ajustarse con el modelo de Corey, por lo que
se generaliza este punto de vista usando las siguientes ecuaciones:
(65)
29
DelshadM.,M.Delshad,D.Bhuyan,G.A.PopeandL.W.Lake:EffectsofCapillaryNumberonthe
Residual Saturation of Three Phase Micellar Solution. Paper SPE 14911 Proceeding of the SPE/DOE
5
th
SymposiumonEnhancedOilRecovery,Tulsa,Oklohoma,Apr.2023,1986.
73
Donde, es la permeabilidad relativa y es la permeabilidad relativa en el

punto final dada por el nmero de captura y la saturacin. y
son la permeabilidad relativa y la permeabilidad relativa en el punto final a
bajo numero de captura en la fase. La saturacin normalizada ( ) dada en la
Ec. anterior se define como:
(66)
Donde es el nmero de fases presentes, es la saturacin y es la
saturacin residual de la fase .

5.1.2 Simulacin numrica

Se utilizaron las curvas de permeabilidad de las Fig. para investigar el efecto
del numero de captura sobre la productividad de un pozo en un yacimiento
de gas condensado. La ecuacin de estado (EOS) del simulador de
yacimientos composicional UTCOMP se utilizo para este anlisis
30
. La
descripcin del fluido y el comportamiento de fase es el mismo dado por Wu
et al
31
. Utilizando el punto de vista de Wang et al
32
para la precisin del
modelo se utilizaron seis componentes hidrocarburos para el comportamiento
del gas y del fluido condensado. Por simplicidad se utilizo la descripcin de
una capa para iniciar el anlisis.
30
Chang,YihBor.,G.A.PopeandK.Sepehrnoori.Ahigherorderfinitedifferencecompositional
simulator.JournalofPetroleumScienceandEngineering,5(1990)3550.
31
Wu,WeiJr,P.Wang,M.Delshad,C.Wang,G.A.PopeandM.Sharma.ModelingNonEquilibrium
MassTransferEffectsforaGasCondensateField.PaperSPE39746presentedattheAsiaPacific
ConferenceheldinKualaLumpur,Malaysia,March2324,1998.
32
Wang,P.,G.A.PopeandK.Sepehrnoori.DevelopmentofEquationsofStateforGasCondensate
forCompositionalPetroleumReservoirSimulation.SubmittedtoIndustrial&EngineeringChemistry
Research,Nov.1997.
74

La simulacin utiliza un dominio bidimensional vertical de seccin transversal
(x-z) con un ngulo de 36. El grid de simulacin tiene 8 capas con la
permeabilidad ms alta en el topo (90 md) y la permeabilidad ms baja en el
fondo (15 md). Se utilizaron 19 bloques en la direccin con un tamao de 1
a 500 pies y con bloques ms pequeos cerca a la cara del pozo. El pozo
produce a una tasa constante de y atraviesa todas las 8
capas. Posee presin constante a ( ) en la frontera que se aplica al
limite externo del yacimiento, la cual permite que el fluido fluya por todo el
grid de simulacin. La temperatura y presin del yacimiento son y
, inicialmente el lquido se encuentra en la regin de dos fases.


Se realizaron simulaciones con y sin efectos del nmero de captura en la
permeabilidad relativa, demostrando la importancia en la saturacin del
condensado y el ndice de productividad. Las Figs. (16 y 17) ilustran la
distribucin del condensado con y sin el nmero de captura. La Fig.16
muestra que la saturacin del condensado tiende a un mximo en estratos
de alta permeabilidad. Esto es debido a que la saturacin del condensado
aumenta cuando el nmero de captura disminuye con el incremento de la
distancia del pozo y finalmente disminuye con el aumento de la presin. La
Fig. 17 muestra que la saturacin del condensado en la regin cercana a la
cara del pozo para el caso sin nmero de captura, es alta cerca al pozo y
disminuye cuando se aleja del pozo en todos los estratos. La normalizacin
del ndice de productividad, con y sin nmero de captura y dos diferentes
75
puntos finales del gas, se muestran en la Fig. 18 y es funcin de la presin

promedio.
El ndice de productividad en este anlisis es calculado a partir de la
siguiente ecuacin:

(67)
Donde, es la tasa de produccin total en , es la presin
promedio del yacimiento en y es la presin de fondo fluyendo en
. La Fig. 18 muestra que el ndice de productividad disminuye
rpidamente como el condensado se acumula en el yacimiento, pero este
efecto se acenta cuando se reduce la saturacin de condensado a altos
valores de numero de captura. El efecto del nmero de captura es ms
significativo en el caso de punto final de permeabilidad relativa al gas bajo
que en el valor ms alto.

Figura 16: Saturacin de condensado despus de 60 das. (Incluyendo los
efectos del numero de captura).

76

Fuente: Pope G. A., Wu W., Narayanaswamy G., Delshad M., Sharma M. and Wang P. Modeling Relative
Permeability Effects in Gas-Condensate Reservoirs. Paper 49266, prepared for presentation at the 1996 Annual
Technical Conference and Exhibition held in New Orleans, Louisiana, 27-30 September 1996.

Este anlisis de simulacin muestra porque la permeabilidad relativa no debe
ser modelada basada solo en la tensin interfacial. La naturaleza
fundamental de los cambios de la permeabilidad relativa se ha demostrado
utilizando el concepto de nmero de captura. La ecuacin base de la tensin
interfacial indica claramente que no es el nico factor que afecta la
saturacin residual y la saturacin de condensado en particular.

Figura 17: Saturacin del condensado despus de 60 das. (sin los efectos
del numero de captura).

77
Permeability Effects in Gas-Condensate Reservoirs. Paper 49266. 1996.

5.2 SATURACIN CRITICA DE CONDENSADO Y PERMEABILIDAD
RELATIVA AL GAS

Cuando los yacimientos de gas condensado se desarrollan por agotamiento
de presin, la entregabilidad del gas se ve afectada por la cantidad y la
distribucin de la condensacin en torno al pozo. Comprender los parmetros
que afectan la distribucin y la cantidad de saturacin de condensado ( ) en
el pozo y el efecto de estos parmetros sobre el gas y el flujo de liquido, es
importante en el desarrollo de los mtodos para aumentar la entregabilidad
del gas. La fuerte reduccin de la entregabilidad del gas puede deberse a la
forma de permeabilidad relativa de la fase gas ( ) los mecanismos del
deterioro de la productividad del gas no estn claros. Los dos principales
parmetros que afectan el recobro del condensado y la entregabilidad de
pozos de gas son: la saturacin critica de condensado ( , que afecta el
recobro de liquido) y la permeabilidad relativa al gas ( , que afecta la
entregabilidad de gas del pozo).

Figura 18: Normalizacin del ndice de productividad.

78

Permeability Effects in Gas-Condensate Reservoirs. Paper 49266. 1996.

Varios autores han medido la saturacin crtica de condensado y han
estudiado la permeabilidad relativa al gas , llegando a diferentes
observaciones como las mencionadas a continuacin:

La permeabilidad relativa de ambas fases liquido y gas disminuye
cuando se incrementa la tensin interfacial y aumenta cuando se
incrementa la tasa de flujo.
La permeabilidad relativa al gas es ms sensible a la tensin
interfacial y la tasa de flujo que la permeabilidad relativa al lquido.
Una comprensin terica de los efectos de la tensin interfacial, la
gravedad, fuerzas viscosas y la notabilidad en la saturacin critica del
79
condensado, la permeabilidad relativa del gas y la permeabilidad

relativa del lquido.

5.2.1 Modelo de red para condensacin en medio poroso.

El medio poroso utilizado en este estudio est representado por un simple
modelo de red de dos dimensiones (2D), con tubos capilares de longitud L y
de intersecciones cubicas de ancho b, se asume que b es igual al mximo
dimetro del tubo capilar en la red, hay cuatro tubos conectados a cada
interseccin. La distribucin log-normal se asume con una desviacin
estndar de 5,74 y un radio promedio de 23 para el tubo capilar. El
tamao de distribucin de la red es de 20x20 intersecciones, el volumen de
poro de la red con una longitud de tubo de 6000 es de unos 8,93 , el
volumen total de las intersecciones con lado de 32,67 m de ancho en torno a
la red es de 0,03 y una permeabilidad estimada de la red de 1 Darcy. El
valor de umbral de radio utilizado en este estudio es 20 (para
condensacin de liquido, primero en forma de tubo con ).

5.2.2 Longitud critica de liquido y saturacin critica de condensado en
un tubo capilar

80
Fang et al.
33
Desarrollo un modelo fenomenolgico para la saturacin critica
de condensado en el medio poroso representado por una red. Su modelo se
baso en dos hiptesis: las fuerzas viscosas son insignificantes y la altura
critica de lquido ( ) formado in-situ en un tubo capilar vertical es la longitud
mxima de la columna de lquido antes que sea mvil. En este estudio, se
ampla el modelo de Fang et al. incluyendo las fuerzas viscosas y ofreciendo
la opcin . Durante el crecimiento del lquido en el tubo capilar, la
tensin interfacial, la gravedad, la histresis del ngulo de contacto, la fuerza
viscosa y el tamao del tubo capilar son los principales factores que afectan
la estabilidad de la nueva fase de lquido.

La altura crtica del lquido formado in-situ en un tubo vertical de radio r se
expresa como: (Ver la Fig.19)
(68)
donde es la presin en el tope y es la presin en el tope de la columna
de lquido; es el angulo de contacto de retroceso y es el angulo de
contancto de avance; es la densidad diferente entre el gas y la fase
liquido. En este estudio, para simplificar se asume que la densidad del gas
es igual a cero y la tensin interfacial ( ) no cambia con la presin. En un
tubo capilar horizontal, la longitud crtica del lquido es calculada
indirectamente a partir de:
(69)

33
Fang,F.etal.AphenomenologicalModelingofCriticalCondensateSaturation.PaperSPE36716
presented at the (1996) SPE Annual Technical Conference and Exhibition, Denver, Colorado, 6 9
October.
81
La ecuacin anterior no implica las fuerzas viscosas (cuando ), el

lquido no puede moverse, a menos que . Despus de determinar ,
la saturacin critica de condensado en un tubo vertical puede ser calculada a
partir de . En un tubo vertical la saturacin critica de condesado es igual
a cero cuando las fuerzas viscosas es mayor que las fuerzas
interfaciales . De lo contrario la saturacin critica de
condensado es igual a 1 o menos de acuerdo a la Ec. 69.

5.2.3 Clculos de la Saturacin critica de condensado S
cc
:

El primer paso para estimar la saturacin critica del condensado, es
determinar el diferencial de presin para cada tubo. La distribucin de
presin en toda la red es calculada a partir de la siguiente ecuacin:

(70)

donde, es la tasa de flujo de fluido en el tubo . La tasa de flujo en un tubo
capilar circular puede ser calculada:
(71)

Si el fluido es incomprensible. En la ecuacin anterior es la viscosidad de
fluido, es la densidad del fluido, y y son el diferencial de presin y el
radio del tubo , respectivamente. Una vez que la presin de entrada y salida
en los extremos de la red se conocen, las Ecs. 70 y 71 se pueden utilizar
para obtener la distribucin de presin en la red. El mtodo de iteraciones de
82
Gauss se puede utilizar para resolver el conjunto de ecuaciones lineales.

en cada tubo puede calcularse si se conoce la distribucin de la presin de
acuerdo a las Ecs. 68 y 69. Entonces el valor de en toda la red se calcula
fcilmente.

5.2.4 Calculo de la conductancia del fluido

Para este estudio se considero un fluido en una sola fase en un tubo capilar
de radio y de longitud . De acuerdo a la Ley de Darcy, la tasa de flujo a
travs del tubo es , donde es la permeabilidad
en el tubo . Una forma alternativa de la ecuacin de Darcy es
, donde es la conductancia en el tubo y puede ser
calculada a partir de:
(72)

Figura 19: Configuracin del liquido en tubo circular.

83

Tubo circular despus de la formacin de liquido

Tubo circular con liquido en el fondo.
Fuente: Li, Kewen and Firoozabadi, Abbas Phenomenological Modeling of Critical Condensate Saturation
and Relative Permeabilities in Gas/Condensate Systems. Paper SPE 56014

donde, y son la conductancia absoluta y la saturacin del fluido en un
tubo , respectivamente; es la viscosidad del fluido . La conductancia
efectiva del fluido en un tubo , se relaciona con la permeabilidad efectiva
de la fase , que puede calcularse usando la siguiente ecuacin cuando
:
(73)

Si se considera un tubo saturado con un poco de liquido (ver la Fig. 19)
que es mayor que . En este caso, el liquido se localiza en el fondo del tubo
84
y el liquido y el gas pueden fluir. La ecuacin para calcular la conductancia

efectiva del gas de un tubo en un flujo bifsico es:
(74)
donde, y son la viscosidad del gas y del liquido condensado,
respectivamente. es la saturacin del liquido condensado en el tubo .
Similarmente, la conductancia efectiva del liquido de un tubo en un flujo
bifsico es:
(75)
donde, es la conductancia de la fase liquido cuando . La
conductividad efectiva de la fase gas y liquido de toda la red se calcula con
base a la conductancia individual de los tubos. Una vez que la conductancia
efectiva de la fase gas y liquido en la red es calculada, la permeabilidad
relativa del gas y del condensado puede ser calculado fcilmente.

5.2.5 Renormalizacion del mtodo y permeabilidades relativas

La renormalizacion es un mtodo exacto y computacionalmente rentable en
comparacin con la simulacin numrica directa. Para una fase determinada
y conociendo la distribucin de saturacin, la conductancia efectiva de la fase
85
es calculada en pequeas regiones, en las llamadas celdas de

renormalizacion y despus a gran escala. En la Fig. 20 se observan con
unidades diferentes de la red el numero de celdas renormalizadas para
una red 2D, en la figura el numero de celdas varia de 2 a 5. Una celda con un
numero de celdas de renormalizacion , es una unidad compuesta de 2
nodos y dos tubos . La conductancia efectiva de una fase es el
conjunto de sistemas que puede obtener despus de la repetida aplicacin
del proceso de renormalizacion. La Fig. 21 muestra el proceso de
renormalizacion para . La red principal est integrada por cuatro
unidades celdas . Despus del primer paso de renormalizacion, la red se
reduce a una unidad celda . En el ltimo paso, la red primaria se convierte
en dos conductores, uno horizontal y otro vertical.


La saturacin critica de condensado se calcula mediante la modificacin del
modelo de en el que las fuerzas viscosas son consideradas y la altura
critica del liquido formado in-situ en los tubos capilares que puede ser
menor que la longitud del tubo. La saturacin critica de condensado es
funcin de la tensin interfacial, la tasa de flujo, la gravedad y el ngulo de
contacto de histresis.

Figura 20: Esquema de renormalizacion para diferente numero de celdas.
86

Fuente: Li, Kewen and Firoozabadi, Abbas Phenomenological Modeling of Critical Condensate Saturation and
Relative Permeabilities in Gas/Condensate Systems. Paper SPE 56014

La saturacin critica de condensado , aumenta a medida que
aumenta la tensin interfacial ( ) y decrece
cuando aumenta la gravedad . Cuando la tensin interfacial es lo
suficientemente grande, la saturacin critica de condensado llega a un
mismo valor mximo para diferentes valores de graveada , esto
significa que el efecto de la gravedad sobre es insignificante a
altos valores de . Cuando es pequeo, es decir, a bajas tasa de
flujo, es afectada por la gravedad, mientras que para valores altos
de , el efecto disminuye.

Figura 21: Esquema del proceso de renormalizacion;
87

Fuente: Li, Kewen and Firoozabadi, Abbas Phenomenological Modeling of Critical Condensate Saturation and
Relative Permeabilities in Gas/Condensate Systems. Paper SPE 56014.

La saturacin critica de condensado , disminuye con el aumento de
las fuerzas viscosas, es decir, aumenta la cada de presin.
En un sistema gas-liquido, se asume que la fase mojante es el liquido,
pero la mojabilidad de sistema roca-fluido en un yacimiento de gas
condensado cerca del pozo, podra ser modificado para que la fase
mojante sea el gas a travs de algunos qumicos, cambiando tambin
la saturacin critica de condensado y la permeabilidad relativa al
gas . La saturacin critica de condensado disminuye cuando el
angulo de contacto de retroceso aumenta de 0 a 85, a un valor de
88
tensin interfacial fijo. La mojabilidad del gas tiene un efecto ms

pronunciado en la saturacin critica de condensado. la saturacin
critica de condesado es sensible a la tensin interfacial si el sistema
roca-fluido, es mojado por la fase liquido. Sin embargo, para un
sistema donde la fase mojante es gas, la saturacin critica de
condensado es mnima y no aumenta mucho, incluso a muy altos
valores de tensin interfacial.
Despus de calculas valores de saturacin critica de condensado a
diferentes tamaos de red, se encontr que a una tensin interfacial
, con o sin fuerzas viscosas, la saturacin critica de
condensado es aproximadamente constante cuando el tamao de la
red es mayor que 15 x 15, demostrando que el tamao elegido para
el modelo de 20 x 20, es adecuado.

La permeabilidad relativa al gas es un parmetro muy importante ya que
cuando disminuye abruptamente, la entregabilidad del gas se reduce
drsticamente. La permeabilidad relativa al gas es afectada por la tensin
interfacial, la tasa de flujo, el ngulo de contacto de histresis y el bloqueo de
condensado.
La permeabilidad relativa tanto al gas como a la fase liquido,
disminuye cuando aumenta la tensin interfacial. La permeabilidad
relativa al gas es ms sensible a la tensin interfacial que la
permeabilidad relativa a la fase liquida.
La permeabilidad relativa al gas y la permeabilidad relativa a la fase
liquido, aumenta cuando se incrementa la gravedad , es decir, con
un aumento de la tasa de flujo. La permeabilidad relativa a la fase
liquido es menos sensible a las fuerzas viscosas que la permeabilidad
relativa al gas.
89
La permeabilidad relativa al gas y a la fase liquido aumenta cuando la

mojabilidad del liquido disminuye ( ) y
. Si la mojabilidad de la roca alrededor
de la cara del pozo puede ser alterada de mojada por la fase liquido a
mojada por el gas, la entregabilidad del pozo puede aumentar
considerablemente. La permeabilidad relativa al gas, disminuye muy
rpido para las condiciones de mojado por la fase liquido y se
aproxima a cero cuando la saturacin critica se incrementa un 22%.
En algunos yacimientos de gas, un pozo puede ser bloqueado por
completo y la produccin de gas puede detenerse altas presiones
34
.
Este fenmeno es el resultado de la formacin in-situ de liquido cerca
a la cara del pozo, cuando la presin del yacimiento cae por debajo de
la presin del punto de roco. El liquido recin formado en el medio
poroso, bloque los poros grandes y pequeos, reduciendo de manera
significativa la eficacia de la permeabilidad relativa al gas.
34
Gravier,J.F.etal.DeterminationofGasCondensateRelativePermeabilityonWholeCoresUnder
ReservoirConditionsSPEFEFebruary(1986)9;TransAIME,281.
90

6. FLUJO COMPOSICIONAL

6.1 MODELO COMPOSICIONAL DEL RETROGRADO
35

En un yacimiento tpico de gas natural, todos los hidrocarburos existen en la
fase gaseosa en el momento de su descubrimiento. Dependiendo de la
composicin inicial de la mezcla de hidrocarburos in-situ y su
comportamiento durante el agotamiento, se reconocen hasta tres tipos de
yacimientos de gas natural: yacimientos de gas seco, gas hmedo y gas
retrogrado o yacimientos de gas condensado. Este ltimo es el ms rico en
trminos de hidrocarburos pesados, por lo que es ms probable que se
desarrolle una segunda fase de hidrocarburos ms pesados (liquido
condensado) en una isoterma de agotamiento, como se ve en la Fig. 22. El
comportamiento del rendimiento de este tipo de yacimiento suelen utilizar
simuladores composicionales debido a que su rendimiento se ve muy
influenciado por los cambios en la composicin de fluidos.

Durante los aos 1970 y 1980, el aumento de la computacin, hizo posible la
aplicacin de simuladores de composicin a escala sobre todo el yacimiento.
Desde entonces una gran cantidad de esfuerzos se han dedicado a
desarrollar un sistema ptimo para la solucin de conjuntos de ecuaciones
35
Ayala Luis F., Ertekin Turgay and Adewumi Michael. Compositional Modeling of Retrograde Gas-Condensate
Reservoir in Multimechanistic Flow Domains. Paper SPE 94856
91
diferenciales de composicin no-lineal. Despus de muchos trabajos se ha

desarrollado un mtodo de solucin totalmente implcito y el mtodo IMPES
(presin implcita, saturacin y composicin explicita o IMPESC),
respectivamente. Young y Stephenson (1983) prefirieron clasificar los dos
mtodos para la solucin bsica del modelo de composicin como Newton-
Raphson y no Newton-Raphson.

En este anlisis, se implementa el flujo multimecanistico provisto de una
ecuacin diferencial en un balance de masa para la simulacin
composicional de yacimientos de gas-condensado. El mtodo de solucin
numrica Newton-Raphson IMPIS (presin y saturacin implcita) se utiliza.
En este caso, el flujo de fluido resulta de la accin combinada del transporte
tanto por presin como por la concentracin molecular del pozo. El objetivo
de este anlisis es explorar la aplicabilidad de este concepto para el estudio
multicomponente, condensacin de gases de composicin variable como en
el caso de gases retrgrados.

6.1.1 Ecuaciones que gobiernan los simuladores composicional de flujo
en dominios multimecnicos

Un balance de materiales composicionales tpico que los simuladores
numricos pueden resolver en un grid discreto para yacimientos
hidrocarburos se pueden escribir de la siguiente manera:

(76)
92

Figura 22: Recorrido de la isoterma de agotamiento para un gas retrogrado.

Fuente: Ayala Luis F., Ertekin Turgay and Adewumi Michael. Compositional Modeling of Retrograde Gas-
Condensate Reservoir in Multimechanistic Flow Domains. Paper SPE 94856

Esta ecuacin debe ser escrita para cada uno de los componentes que estn
en la mezcla de hidrocarburos. La ecuacin anterior se obtiene mediante la
sustitucin de un modelo adecuado de la velocidad en trminos de los
gradientes de los potenciales responsables del flujo, se puede escribir en
trminos de la ecuacin de continuidad de la siguiente manera:
(77)

93
El modelo que une el campo de velocidad a los gradientes del campo

potencial en la mayora de los yacimientos es la ley de Darcy. Sin embargo
en un flujo multimecnico, la velocidad de una especie ( , ) es el
resultado de las contribuciones de flujo causada por un componente Darcian
y un componente Fickian. Es decir, una especie que se mueve como
consecuencia de los efectos de difusin y otra como movimiento de la fase
Darcian. Por lo tanto, la velocidad total de una especie es el resultado de la
adiccin de la accin de un componente tipo Darcian ( , debido al flujo de
la fase macroscpica) y un componente tipo Fickian ( , a causa de los
efectos de difusividad de la especie dentro de la fase) como se muestra a
continuacin:

(78)

La mayor parte del movimiento de una fase en el medio poroso causada por
los gradientes de presin de gravedad se supone que es descrita por la ley
de Darcy. La mayor parte de la fase de velocidad descrita por la ley de
Darcy, est indicada en un valor promedio de toda la seccin del rea
transversal incluyendo granos y espacio poroso, dado por:

(79)
donde, y

94
La difusin o Velocidad de Fickian para cada una de las especies se hace

igual a la concentracin neta efectiva de la velocidad de la difusin en la
fase:
(80)
El fenmeno de difusin de la fase liquida se puede considerar menos
importante que los efectos de la difusin de la fase gas. Por lo tanto se
supone que el fenmeno multimecnico solo tiene lugar en la fase gas,
mientras que el flujo de la fase liquida es solo atribuible al gradiente de
presin. Entonces
(81)

En un medio poroso, la difusin de las especies dentro de una fase no se
produce a travs de toda la seccin transversal (roca y espacio poroso), sino
a travs de zonas abiertas a la difusin de esa fase. Para una seccin
rectangular de cualquier volumen, la superficie disponible para la difusin del
gas est dada por:
(82)
Si se utiliza el concepto de coeficiente de difusividad efectiva, el nmero total
de molculas de gas difundidas, en un lugar dado del yacimiento, por unidad
de rea esta dado por:
(83)
A partir de esto, la velocidad efectiva Fickian, por unidad de rea promedio
de la fase de gas, est dada por:
(84)
95
Porque todas las velocidades Fickian en las Ecs. 79 son valores promedios
de velocidad con respecto a las secciones de reas transversales de flujo
( ), transformndolo de la siguiente forma:
(85)
Por lo tanto,
(86)
Cuando se sustituyen todas las consideraciones anteriores en las Ecs. 79, se
obtiene las siguientes ecuaciones multimecnicas:

(87)

Sustituyendo estas ecuaciones en la ecuacin de continuidad Ec. 78, se
obtiene el balance de masa para composicin multimecnica, que el
simulador numrico debe resolver para cada uno de los componentes de
cada grid en el sistema:

(88)

Cuando hay agua presente en el yacimiento, se debe incluir una ecuacin
diferencial parcial. Si no se consideran los efectos multimecnicos que se
96
presentan en la fase agua, los movimientos de esta fase se rigen por la ley
de Darcy, como se ve en el siguiente balance de materiales:
(89)
La inclusin de la difusin mediante la ley de Fick con multiplicador de
saturacin y porosidad y la adveccin calculada por la ley de Darcy, ha
llamado el modelo advectivo-difusivo (ADM).


En el medio poroso, el flujo de fluidos ha sido relacionado con el fluido de los
gradientes de presin y de gravedad, a travs de una relacin conocida
como ley de Darcy. La permeabilidad, mide la facilidad con que un fluido se
mueve en el medio poroso debido a la creacin de gradientes de presin. En
sistemas con muy baja permeabilidad conocidos como sistemas apretados
(k<0,1md), el flujo de fluidos macroscpicos podra no darse, como est
previsto en la ley de Darcy, a pesar de la presencia de grandes gradientes de
presin. Ertekin et al. (1986) sugiri que en ese caso, el flujo molar causada
por los gradientes de concentracin (es decir difusin), podra asumir como el
principal mecanismo del fluido, en relacin con el flujo Darcian o flujo
macroscpico.

El modelo presentado en la Ec. 83, describe el funcionamiento eficaz de
Fickian o la difusividad del flujo de gas: en la Ec. 83 representa la
densidad molar o concentracin molecular del gas por unidad de volumen de
la fase y es el coeficiente eficaz de difusin del gas. El uso de un
97
coeficiente eficaz de difusin permite mirar el efecto combinado de los

diferentes mecanismos de difusin que podran tener lugar en el medio
poroso. La difusin de un gas a travs del medio poroso puede ser
impulsada por diferentes mecanismos: difusin libre molecular, difusin
Knudsen y difusin de superficie. La difusin libre implica la difusin de s
mismo a travs de una fase de desequilibrio debido a la aparicin de
molcula a molcula en las regiones de las colisiones moleculares de
diferentes concentraciones. La difusin de Knudsen tiene a dominar en el
flujo de los gases a travs de un estricto sistema de poros muy pequeos. En
este caso, el dimetro de poro es tan pequeo que las molculas del gas
tienden a chocar con las paredes de los poros y por lo tanto, no ve a la otra.
Entonces las molculas-pared dominan las colisiones y las interacciones
molecula-molecula son raras. La difusin de superficie es el transporte de las
molculas de gas adsorbidas a lo largo del gas y de la superficie o de las
interacciones gas/liquido de las interfaces.

Silvia y Belery (1989) justifica el uso de los coeficientes de difusin efectiva
en el modelo de difusin en yacimientos naturalmente fracturados y
presentaron un modelo para el clculo de los coeficientes de difusin efectiva
de un componentes en una mezcla multicomponente. Cussler explico que un
aumento al doble de la presin del gas reduce el coeficiente de difusin a la
mitad, mientras que el aumento al doble de la presin se debe al aumento de
la temperatura absoluta.

Se considera que el flujo multimecanico, desempea un papel importante en
la recuperacin en matrices muy apretadas de almacenamiento de lquido
condensado. Cuando la matriz est rodeada por una red de fracturas de alta
permeabilidad, el sistema es representado como un yacimiento de gas
98
condensado naturalmente fracturado. si los bloques de la matriz se

consideran rectangulares, el sistema de fracturas naturales puede parecer la
configuracin de cubo de azcar, estudiada por Warren y Root (1963).
6.2 EFECTOS DE LOS CAMBIOS COMPOSICIONALES EN EL
TRASIENTE DE FLUJO

Los yacimientos de gas condensado difieren en su comportamiento de los
yacimientos convencionales y de la optimizacin del recobro, necesidad de
un anlisis cuidadoso del yacimiento, planificacin y gestin. En el momento
de su descubrimiento, los yacimientos de gas condensado a menudo
contienen una sola fase de vapor de gas. A medida que son producidos, la
presin disminuye en el yacimiento y en las instalaciones de superficie, se da
la condensacin de liquido del gas. Esta condensacin isotrmica despus
que la presin del yacimiento cae por debajo de la presin de roco, se
conoce como condensacin retrograda.

La caracterizacin de los yacimientos de gas condensado es una difcil tarea,
ya que el flujo multifasico en el yacimiento y el cambio de composicin de la
mezcla de fluidos que fluyen hacia el pozo, complican la interpretacin de
las pruebas de pozo considerablemente. Los analistas de campos de gas
condensado, generalmente asumen una sola fase de flujo, fluido
comprensible, entre otras, estas suposiciones pueden dar una idea
equivocada del campo. El mtodo de pseudopresin del gas muestra buenos
resultados en el anlisis de campos de gas condensado, pero en campos
muy ricos los resultados obtenidos pueden ser inexactos. Esto se puede
asumir a diversos factores como por ejemplo: este mtodo asume una sola
fase de gas, el efecto del petrleo se descuida, la permeabilidad relativa al
99
gas no se tiene en cuenta, al igual que la transferencia de masa y los efectos

del cambio de composicin.

Olaberinjo
36
considero el comportamiento de flujo en un yacimiento de gas
condensado y desarrollo un modelo que permite calcular aproximadamente
la presin de agotamiento, de un yacimiento de gas condensado, a partir de
la composicin original del fluido del yacimiento, utilizando los clculos de
equilibrio liquido-vapor, las propiedades del fluido, propiedades
adimensionales y la acumulacin del liquido retrogrado como gua.

6.2.1 Desarrollo del modelo.

Toda interpretacin directa de la presin transitoria de una prueba de pozo,
se basa en la ecuacin de difusin lineal de Horne
37
.

(90)
La Ec. 90 se aplica a una sola fase ligeramente compresible, yacimiento de
aceite (es decir, con ).

En los yacimientos de gas condensado, lquidos de condensacin retrograda
y el gas pueden estar presentes juntos y puede ocurr flujo multifasico. El
flujo de estos fluidos en el medio poroso solo puede estar representado por
36
Olaberinjo, A. F. Modeling the Effects of Compositional Changes in Transient Flow of Gas
CondensateReservoirs.PaperSPE106339presentedatthe2006SPETechnicalSymposiumofSaudi
ArabiaSectionheldinDhahran,SaudiArabia,2123May(2006).
37
Horne,R.N.ModernWellTestAnalysis,PetrowayInc.PaloAltoCA(1990)
100
una ecuacin diferencial no lineal. Por lo tanto, la Ec. 90 no es vlida para

representar estos lquidos.

La ecuacin de continuidad para el componente es un sistema de gas
condensado es representado por:

(91)

La ecuacin de Darcy que describe el flujo de la fase de petrleo y gas
cuando cada una de las fases esta fluyendo a tasas constante , dadas
por:
(92)
(93)

Donde es el volumen de vapor de lquido que fluye o la movilidad de las
fases que fluyen, definida por:
(94)

Sustituyendo las Ecs. 92 y 93 en 91:
(95)
La compresibilidad se define:
(96)

101
Introduciendo el efecto de la transferencia de masa, en trminos del

concepto molar, la densidad molar est dada por:
(97)
Entonces,
(98)
La Ec. 95 se convierte:
(95)

Sin embargo a partir de la Ec. 94
(96)

Entonces
(97)

La introduccin de la ecuacin de estado para los gases y la relacin de
permeabilidad.
(98)

(99)

La Ec. 97 se convierte en:
(100)

102
La Ec. 100 es una firme ecuacin de difusividad del gas condensado con
componente en dos fases de flujo, en un yacimiento de geometra radial,
derivada de conformidad con la hiptesis de que la presin inicial del
yacimiento es igual a la presin del punto de roco, el pozo produce a una
tasa constante (es decir, la tasa de flujo de cada componente en su fase es
constante) y la formacin es penetrada completamente, los efectos
gravitacionales y los efectos capilares son ignorados y el yacimiento es
homogneo.
Un mtodo clsico para tratar la Ec. 100 es definir una variable llamada
pseudopresin y pseudo compresibilidad, :

(101)

Ahora definimos pseudo compresibilidad, :

(102)

103
Sustituyendo las Ecs. 101 y 102 en la Ec. 100, obtenemos la ecuacin que
gobierna la transmisin de presin en un yacimiento de gas condensado con
tendencia a cambios composicionales y efecto de transferencia de masa
como:
(103)


Se utilizaron tres caso de aplicacin, en el primer caso se analizo un gas
condensado con alto abandono de liquido, el efecto se traduce en una mayor
cada de presin con una tasa de declinacin de presin que varia
considerablemente de los otras casos. En el segundo caso, se analizo un
gas condensado con bajo abandono de liquido, y en el tercer caso, se
analizo un gas condensado con muy bajo abandono de liquido, para estos
gases, el cambio en la diferencia de presin con el tiempo no es significativa,
ya que el coeficiente de flujo de la fase gas a la fase aceite y la
permeabilidad relativa al gas son considerablemente altas. Esto se debe a
que en el caso donde el gas condensado es rico, si se afecta la
permeabilidad relativa al gas.

En conclusin, un nuevo mtodo para calcular el agotamiento de la presin
en un yacimiento de gas condensado es demostrado. Este mtodo considera
los cambios de composicin de los lquidos del yacimiento, los efectos de la
transferencia de masa durante la produccin, las impurezas que puedan
estar presentes en la mezcla, las dimensiones de las propiedades, el anlisis
104
de pseudopresiones y las propiedades de las fases vapor y liquido.

Adicionalmente, el modelo podra ser utilizado para interpretar los resultados
del flujo de las dos fases probado en yacimiento, sobre todo para un campo
de gas condensado o algunos sistemas con flujo de gas y aceite.

La composicin general de la mezcla que fluye en cualquier parte del
yacimiento cambia a medida que se avanza con la produccin, por lo tanto
es una funcin del tiempo y la presin. La composicin del vapor y del
lquido, junto con la fraccin molar del lquido y del vapor puede calcularse
directamente a partir de clculos de equilibrio liquido-vapor. El modelo
proporciona una base para la comprensin de los mecanismos y factores que
controlan la produccin del gas y del condensado durante la declinacin.

La distribucin de presin en un yacimiento de gas condensado expresada
en trminos de pseudopresin cumple la clsica ecuacin de flujo de lquido
es estado estable. La interpretacin de la pseudopresin conduce una
correcta estimacin de la capacidad de flujo, y el efecto dao. Tambin es
til para clculos de entregabilidad del pozo y anlisis de sensibilidad.
Adems, la correlacin entre dimensiones de pseudopresin y el tiempo
puede ser representada por la solucin de flujo de fluidos ligeramente
compresibles. El mtodo funciona mejor cuando las diferencias de presin
y son grandes. Tambin es aplicable a yacimientos de
aceite y gas con composicionales.

105
6.3 INTERACCIONES NO-LINEALES ENTRE LAS RESTRICCIONES DE

FLUJO CERCA AL POZO
38

La productividad de los pozos, se ve reducida por restricciones de flujo cerca
al pozo, estas restricciones de flujo incluyen: dao mecnico,
completamiento parcial, alteraciones debido a la formacin de lquido
residual o lquido de abandono y los efectos del flujo no-Darcy.
A pesar de todos los avances tecnolgicos que existen hoy en da, aun los
resultados obtenidos de un proceso de simulacin dependen en gran medida
del juicio subjetivo del Ingeniero de simulacin, esto se debe a varias
razones: (1) a menudo no se comprenden del todo los fenmenos fsicos
complejo que suceden en torno a un pozo; (2) no se puede realizar un
seguimiento de los caminos de flujo en la regin cercana al pozo, controlada
por la geometra del pozo, completamiento y prcticas de estimulacin y
operaciones de pozo; (3) no se tiene una descripcin exacta de la formacin
de las propiedades cerca al pozo; (4) los yacimientos de simulacin y las
herramientas utilizadas presentan limitaciones. Debido a todo lo anterior la
manera ms sencilla de calcular la entregabilidad de un pozo, es utilizar la
ecuacin de rendimiento del pozo. Considerando el estado de equilibrio
pseudoestable, la entregabilidad y el ndice de produccin pueden ser
representados en las ecuaciones 104 y 105 con potencial pseudo-gas
(104)
donde,
(105)
38
Hwang, M. K. Modeling Nonlinear Interaction Among Near-Well Flow Restrictions in Well-Deliverability Prediction
and Simulation. Paper SPE 65407.
106
donde es el espesor total de la formacin, permeabilidad del gas,

tasa de flujo de gas, presion promedio del yacimiento, presin de
fondo fluyendo, radio del pozo, radio del yacimiento,

es el dao total y
una constante de conversin.

Un simulador de yacimiento puede ser utilizado para calcular la
entregabilidad del pozo, la masa en solucin y el impulso de conservacin de
energa. En este caso, el clculo de la entregabilidad puede ser dividido en
dos partes: profundidad en el yacimiento de la regin cercana al pozo y
desde la regin cercana al pozo hasta el pozo. La correspondiente ecuacin
de flujo es dada a continuacin:
(106)
(107)
donde es el ndice del pozo y es el radio de distancia desde el centro
hasta el primer grid en un grid radial (o radio de Peaceman en un grid
rectangular), donde es definido.

6.3.1 Restricciones de flujo cerca al pozo

La precisin de la entregabilidad calculada, depende de la forma en que las
restricciones de flujo son representadas. Aquellas restricciones de flujo cerca
al pozo son comnmente llamados efectos de dao (skin). Existen varios
107
tipos de estos efectos y son fluidodinmicos, interactan los unos con los
otros. Algunos de estos efectos son estticos y otros dinmicos. Un efecto de
dao esttico, refleja una cada de presin adicional debido al dao en el
pozo, se llama dao mecnico ( ). tambin se puede presentar una cada
de presin adicional debido al cambio de la permeabilidad de la formacin en
la regin cercana al pozo, esto se presenta por el abandono de
condensados, se denomina dao alterado a la formacin ( ). Efectos de
dao dinmico, incluyen una cada de presin adicional, causados por la
convergencia de flujo, debido a un completamiento parcial y por el flujo no-
Darcy. Estos se denominan dao de completamiento parcial ( ) y dao no-
Darcy ( ), respectivamente.
Cuando un pozo es producido, a menudo se presentan mltiples efectos de
dao, generalmente estas interacciones son ignoradas y el dao total se
calcula sumando los distintos componentes:
(108)
Sin embargo, las interacciones de los daos son no lineales y a menudo
multiplicativas. Por lo tanto, el dao total es a menudo mucho mayor que el
calculado convencionalmente.

6.3.1.1 Representacin del pozo. En la simulacin de yacimientos, una
restriccin de flujo cerca al pozo, puede ser representada de dos maneras.
(1) introducir un factor de dao como un parmetro y (2) ajustar la red de
distribucin o la formacin de propiedades, de manera que la presin y la
distribucin de velocidad en la regin cercana al pozo sean adecuadamente
cambiadas. Otros problema que se debe considerar en la representacin de
un pozo, son el tamao de la red de distribucin y el efecto sobre el ndice de
productividad simulada en un modelo escala del yacimiento. Un sistema de
108
red con grandes grid a menudo origina una distribucin errnea de presin,
velocidad, saturacin de fluidos y composicin. Este efecto tiene mayor
incidencia en la regin cercana al pozo.

Para conservar la exactitud y la productividad y la productividad del modelo
de simulacin escala del campo, varios trabajos han propuesto ajustar el
grueso de la red del ndice de productividad (J) obtenido. Sobre la base de la
J correcta, conocida a priori, algunos parmetros en la ecuacin pueden ser
cambiados. Fevang y Whitson
39
proponen utilizar la integral de
pseudopresiones ( , equivale a ) en funcin de tres parmetros ( ,
,
), donde es la produccin acumulada de la relacin aceite/gas ( ).
La tabla de puede ser construida a partir de una serie de corridas en una
grid fina, usando la correcta
.
Similarmente, Malachowski et al.
40
, utilizaron
el pseudodao ( ) dado como funcin de dos parmetros ( , ),
donde es el caudal del componente clave y es la fraccin de molar
del wellstream. El ajuste ms amplio ha J ha sido incorporar los efectos el
numero capilar y el flujo no-Darcy en la regin cercana al pozo
41
,
42
. Estos
dos efectos son muy sensibles al tamao del grid.

39
Fevang,O.andWhitson,C.H.ModelingGasCondensateDeliverabilitySPERENovember(1996)
220
40
Malachowski,M.A.,Yanosik,J.L.,andBatten,A.H.SimulationofWellProductivityLossesDueto
NearNearWellCondensateAccumulationinFieldScaleSimulation.PaperSPE30715presentedat
the(1995)SPEAnnualTechnicalConferenceandExhibition,Dallas,2025October.
41
Lee,S.andChaverra,M.ModelingandInterpretationofCondensateBankingfortheNearCritical
CupiaguaField.PaperSPE49265preparedforpresentationatthe(1998)SPEAnnualTechnical
ConferenceandExhibition,NewOrleans,2730September.
42
GCSim:PredictionofGasCondensateReservoirRateTimeProductionPerformance,PARAA/S
InternalReport.September(1998)
109
Para un pozo vertical 1D radial, un yacimiento homogneo, la ecuacin del

dao total, se considera afectada por mltiples zonas como se ve en la Fig.
23 Cuando los mencionados cuatro daos estn presente, es necesario
calcular el correspondiente dao total efectivo para estimar su efecto en J. La
ecuacin total de dao puede ser obtenida de las interacciones de todos los
daos, como se menciono anteriormente, despus de unos artificios
matemticos se llega a la siguiente ecuacin algebraica que permite calcular
el dao total efectivo:
(109)
donde,

(110)

Figura 23: Permeabilidad vs. Radio. Para 1D.

Fuente: Hwang, M. K. Modeling Nonlinear Interaction Among Near-Well Flow Restrictions in Well-Deliverability
Prediction and Simulation. Paper SPE 65407.
110

Para el caso 2D, multicapas, consideramos un yacimiento como el que se
muestra en la Fig. 24. En la figura se observa que hay varias capas que se
encuentran en estado pseudoestable y se comportan como una nica capa
en el yacimiento. Para el caso 2D la ecuacin queda de la siguiente manera:
(111)
donde, i corre por todas las capas abiertas; es la permeabilidad al gas
en la celda una de la capa, y
(112)
(113)
donde, I es el nmero de zonas alteradas en la direccin radial del
yacimiento.
Figura 24: Caso de un yacimiento multicapas.

111

Fuente: Hwang, M. K. Modeling Nonlinear Interaction Among Near-Well Flow Restrictions in Well-Deliverability
Prediction and Simulation. Paper SPE 65407.


La Ec. 111 resume en su totalidad las interacciones entre todos los
componentes de dao y sus contribuciones al dao total. Para un yacimiento
homogneo la Ec. 111 queda reducida a la Ec. 109.

La Ec. 111 es muy til para la entregabilidad del pozo y el anlisis del
rendimiento del pozo. Proporciona una valiosa visin sobre cmo afectan
cada componente del dao a J, a travs de S
T
, especialmente con lquidos
de abandono o los efectos de fluidos residuales del completamiento.

En el anlisis de pruebas de pozo, la Ec. 111 reduce el nivel de
incertidumbre asociado con la estimacin de parmetros de cuantificacin y
cmo interactan entre ellos mismos. Las Ec. 109 y 111 se pueden usar para
calcular el efecto de las restricciones de flujo as como la estimulacin de la
productividad del pozo.

Un simulador de yacimiento puede ser utilizado para cuantificar directamente
los efectos del dao en el pozo, as como los efectos combinados de las
interacciones no lineales entre los distintos componentes del dao.
112

6.4 FLUJO EN LAS REGIONES PERFORADAS

La regin ms importante que afecta la productividad de un pozo es la regin
perforada. Muchos investigadores matemticos han centrado sus estudios en
este tema, pero lo han enfocado a la regin de una fase, prestndole menos
atencin a la regin de dos fases. Diversos investigadores como Danesh,
Henderson y Ali, demostraron que la permeabilidad relativa al gas y al
condensado, puede aumentar significativamente por un incremento de la
tasa de flujo. Este efecto, conocido como Acoplamiento Positivo complica el
flujo de gas y condensado en la cara del pozo.
El flujo de gas y condensado en la regin perforada se analizara mediante un
modelo de elemento finitos. Este modelo permite cambios en las propiedades
del fluido y contabiliza los acoplamientos positivos y efectos inerciales
negativos utilizando un flujo fraccional basado en correlaciones de
permeabilidad relativa.

En sistemas de gas condensado, los valores de tensin interfacial
relativamente bajos, los procesos de condensacin alrededor de la cara del
pozo cuando la presin cae por debajo de la presin del punto de roco,
originan una regin en la que fluyen las fases de gas y de condensado. Esto,
junto con la presencia de condensado en todos los poros, origina un
mecanismo de flujo diferente a los yacimientos convencionales de gas y de
113
aceite. Este comportamiento resulta en un mejoramiento en la permeabilidad

relativa, causada por un aumento en la velocidad.

6.4.1 Simulador de pozos perforados 3D

El modelo 3D utilizado en este anlisis consiste en un pozo con un radio
en una sola capa, yacimiento cilndrico con un radio externo , con un
medio poroso homogneo, con permeabilidad absoluta y espesor de la
formacin . Perforaciones cilndricas con un radio y una longitud , que
se distribuyen en espiral alrededor del pozo con una densidad de disparo de
(numero de perforaciones por pie). El espaciamiento entre perforaciones
es constante e inversamente proporcional a la densidad de perforacin. El
ngulo de fase, , indica el ngulo entre dos tneles de perforacin
sucesivos y es un valor constante.

6.4.1.1 Ecuaciones que gobiernan, una fase. En un simulador de flujo de
una sola fase, las ecuaciones que describen el flujo de un fluido compresible
se resuelven en condiciones de estado estable. La ecuacin de continuidad
para este caso se puede escribir como:
(114)
donde, . Es el operador de divergencia, es la densidad, y es la
velocidad.
114

Los diferenciales parciales en la ecuacin de flujo incluyen los efectos
inerciales (no-Darcy), propuesto por Forchheimer (1914). Comnmente se
aplica en un sistema de flujo 1D. La ecuacin de Forchheimer para un
sistema 3D puede escribirse de la siguiente forma:
(115)
donde, es el operador gradiente, es la presin , es la viscosidad, es la
permeabilidad absoluta, es el valor absoluto del vector velocidad, es el
coeficiente de Forchheimer, tambin se refiere al coeficiente no-Darcy o
inercial. El termino del lado derecho expresa la tasa o impulso de
transferencia en cada direccin. La Ec. 115 puede ser reorganizada de la
siguiente manera:
(116)
El termino de la velocidad absoluta, , puede estar relacionada con el
gradiente de presin y su solucin esta dada por una ecuacin polinomica de
segundo orden, la Ec. 115 para es:
(117)
La raz negativa del polinomio de segundo orden se descarta porque no
puede ser negativo.
Sustituyendo la Ec. 117 en la Ec. 116, obtenemos:
115
(118)
Para evitar la incoherencia dimensional, todos los parmetros en las
ecuaciones presentadas e implementadas en el simulador se encuentran en
unidades (SI).

6.4.1.2 Mtodo de solucin matemtico, una fase. La ecuacin de
diferenciales parciales no lineales ( ), Ec. 117 depende principalmente de
la variable de presin .
Las condiciones lmite para este sistema son:
La presin en el lmite exterior es conocida.
No hay entrada ni salida en la cara del pozo, excepto a travs de las
perforaciones.
La presin en el lmite interno es conocida.
Debido a la alta conductividad de las perforaciones comparada con el
medio poroso, la presin dentro de las perforaciones es constante e
igual a la presin en la cara del pozo.

La ecuacin de diferenciales parciales no lineales y las condiciones de
confinamiento se resuelven utilizando el software Comsul Multipysics (2002).
Este paquete matemtico utiliza el mtodo de elementos finitos.

116
En el sistema 3D, hay un gran contraste entre el tamao de las perforaciones

y el radio externo. Esto exige una estructura de malla que debe captar las
pequeas curvaturas de las perforaciones, reconociendo el gran agotamiento
de la zona. La principal dificultad en la simulacin de perforaciones es el
complejo flujo alrededor de las perforaciones asociados con los cambios
bruscos en los parmetros de flujo. Por lo tanto, el simulador numrico debe
adoptar una malla ms fina alrededor de las perforaciones con una adecuada
representacin geomtrica. Adicionalmente, el nmero de nodos o elementos
en la malla debe ser optimizado para disminuir el tiempo de CPU y requisitos
de memoria.

6.4.1.3 Ecuaciones que gobiernan, dos fases. En un sistema de flujo de
dos fases, hay transferencia de masa entre las dos fases de flujo que viajan
a travs del medio poroso, pero la masa total que fluye fuera del sistema es
la misma que fluyen en el medio poroso.

La ecuacin de continuidad para este sistema es:
(119)
donde, los subndices y denotan gas y condensado , respectivamente.

En un sistema de dos fases con efectos insignificante de presin capilar y
gravitacionales, se puede aplicar la misma ecuacin de flujo que para una
sola fase de flujo (Ec. 116) reemplazando la permeabilidad absoluta ( ) por la
permeabilidad efectiva ( ) y el factor inercial de flujo de una sola fase ( ) por
la de flujo bifsico para cada una de las fases de flujo. Obteniendo:
117
(120)
Donde

La ley de Darcy tambin puede extenderse a flujo de dos fases:
(121)
Las Ecs. 120 y 121 sugieren que ( ) en la Ec. 120 equivale a la
permeabilidad relativa ( ) para cada fase en la Ec. 121. Sustituyendo la Ec.
121 en la Ec. 119 y dividiendo por obtenemos:
(122)
Esta ecuacin es la combinacin de la ecuacin de continuidad y las
ecuaciones de flujo para sistemas de gas-condensado con presin capilar
insignificante.

Se ha asumido que podra haber transferencia de masa entre las dos fases
que fluyen, pero la composicin total del fluido ( ) permanece constante en
el fluido que viaja a travs del medio poroso. Es decir, no hay ninguna fuente
o sumidero en la masa total de cada componente. Por lo tanto,

(123)
118
donde, , y son las fracciones de masa del componente en la mezcla,

en el gas y el condensado, respectivamente; es la tasa de flujo
volumtrico; y es el flujo fraccional del gas. La tabla 4 contiene y el
correspondiente gas condensado ( ) y la tasa de flujo de gas condensado
( ) en vol/vol en todas las condiciones de flujo.

6.4.1.4 Propiedades del fluido, dos fases. Las propiedades del fluido en
equilibrio de fases dependen solo de la presin para una temperatura dada.
Las propiedades del fluido pueden calcularse utilizando un modelo
composicional o simplemente utilizando una relacin de correlaciones
empricas de presin. El subprograma matemtico del paquete Maltad (2003)
Permite cubico spines fue utilizado para obtener los datos de composicin,
densidad, viscosidad y tensin interfacial, sin oscilaciones. Cabe sealar que
una mezcla binaria se utilizo para simplificar los clculos y reducir el tiempo
de CPU.

Tabla 4: Flujo fraccional de gas ( ) y el correspondiente gas
condensado ( ) y la tasa de flujo de gas condensaco ( ) en Vol/Vol
a todas las condiciones de flujo.

1 Infinito 0
0.995 199 0.005
0.99 99 0.01
0.981 51.6 0.019
0.972 34.7 0.029
119
0.952 19.8 0.05

0.910 10.1 0.099
0.714 2.5 0.4
0 0 Infinito

Fuente: J amiolahmady, M., Danesh, A., Sohrabi, M. and Ataei, R. Gas-Condensate Flow in Perforated
Regions. Paper SPE 94072.

6.4.1.5 Permeabilidad relativa, dos fases. En la Ec. 122, la permeabilidad
relativa al gas ( ) es estimada utilizando el flujo fraccional, basada en
correlaciones desarrolladas recientemente (J amiolahmady et al. 2006c). en
esta formulacin, es la interpolacin entre la curva base y la curva de
fluidos miscibles utilizando una funcin de interpolacin , que tambin
representa el efecto de microporos. La curva base es medida a alto valor de
tensin interfacial y baja velocidad, ajustado por los efectos inerciales. La
curva de permeabilidad relativa del gas miscible, tambin es ajustada para
incluir el efecto inercial. La permeabilidad relativa al condensado esta
relacionada con la permeabilidad relativa al gas en la definicin de flujo
fraccional.

6.4.1.6 Mtodo matemtico de solucin, dos fases. Para dos fases existe
una ecuacin principal (Ec. 122) y una ecuacin auxiliar (Ec. 123).
En la Ec. 122:
Las propiedades del fluido, y , son solo funcin de la presin.
120
La permeabilidad relativa al condensado, , es obtenida a partir de la

siguiente ecuacin:
(124)
Este clculo implica , , y la permeabilidad relativa al gas .
La permeabilidad relativa al gas es interpolada a partir de la siguiente
ecuacin:
(125)
Este clculo requiere la permeabilidad de la curva base ( ), la
permeabilidad de la curva miscible ( ) y el parmetro de
interpolacin ( ) que se obtienen a partir de las siguientes
ecuaciones:
(126)
Depende de las propiedades de la roca, , y valores medidos
experimentalmente que es funcin de .
(127)
Depende de las propiedades de la roca, y .

En la ecuacin 123.
121
Las propiedades del fluido, y , y la composicin del fluido son solo

funcin de la presin. Por lo tanto, es funcin de la presin del
pozo.

En estas dos ecuaciones, la variable principal es dependiente de la presin,
pero para resolver la presin y el flujo fraccional de gas, se establecieron las
siguientes condiciones lmite:

1. La composicin total de la mezcla es conocida, es decir, el flujo
fraccional en el pozo se conoce. Para una presin dada, un cambio en
la composicin de la mezcla ocasiona un cambio en el flujo fraccional.
2. La presin en el lmite interno es conocido, es decir, se conoce la
presin en la cara del pozo.
3. No hay entradas ni salidas en la cara del pozo, excepto en las
perforaciones.
4. La presin en el lmite externo es conocida.
5. Debido a la alta conductividad de las perforaciones comparada con el
medio poroso, la presin dentro de las perforaciones es constante e
igual a la presin en la cara del pozo.

Similar que para el caso de flujo de una sola fase, para resolver el sistema se
utilizo el software Comsul Multipysics (2002)
6.4.2 Simulador de pozo a hueco abierto 1D

122
El modelo 1D de hueco abierto utilizado en este anlisis consiste en un pozo

con un radio en una sola capa, yacimiento cilndrico con un radio externo
. El simulador 1D utiliza una simetra existente en y un angulo de
orientacin de el flujo radial. Es decir, las ecuaciones se resuelven para una
sola variable independiente, , con una extra termino en el interior de los
parntesis del operador de divergencia, en la parte izquierda de las Ecs. 118
y 122, con el fin de incluir la variacin de este parmetro con la direccin.


Este estudio se centro en el flujo de gas condensado alrededor de las
perforaciones, teniendo en cuenta la transferencia de masa entre las fases
de gas y de lquido y la dependencia de la permeabilidad relativa del gas y
del condensado, con la velocidad y la tensin interfacial.

El xito del completamiento de las zonas perforada para flujo lineal de una
sola fase es afectado por los parmetros de perforacin (radio, longitud, fase
angular y numero de perforaciones) y el radio del pozo. el efecto de la
velocidad se tiene en cuenta mediante la inclusin del trmino inercial. En un
sistema de flujo de dos fases, las propiedades de los fluidos y el flujo
fraccional de las dos fases fluyendo, juegan un papel importante en la
eficiencia del completamiento de las zonas perforadas. El gran numero de
parmetros pertinentes requiere de una cuidadosa seleccin de la expresin
de la eficiencia de perforacin. Se selecciono la relacin de productividad
( ) para expresar la eficiencia de perforacin. La relacin de productividad
123
es igual a la relacin de la tasa de flujo volumtrico para el caso perforado

con la del hueco abierto sin perforar a las condiciones de la cara del pozo,
porque el flujo fraccional de los dos casos es igual:
(128)

Las Ec. 118 y 122, permiten analizar el efecto de las perforaciones en el flujo
del gas condensado. Este anlisis es realizado mediante un modelo de
elementos finitos. Estas ecuaciones tienen en cuenta el cambio de las
propiedades y contabiliza los efectos de acoplamiento positivo e inercial
mediante el flujo fraccional.

6.5 FLUJO EN POZOS FRACTURADOS HIDRULICAMENTE

El flujo de gas condensado alrededor de un pozo fracturado hidrulicamente,
es diferente a los sistemas convencionales de gas y aceite. Esto se debe
principalmente a los cambios de fase, acumulacin de condensado y
acoplamiento (es decir, incremento de la permeabilidad relativa, incremento
de la velocidad y/o disminucin de la tensin interfacial) y efectos inerciales
(es decir, disminucin de la permeabilidad relativa y de la velocidad).

La descripcin de los pozos fracturados hidrulicamente en yacimientos de
gas condensado utilizando los simuladores de yacimientos requiere de grid
muy finos para capturar los cambios significativos de flujo y de las
124
propiedades de la roca que ocurren alrededor de la fractura. Esta es una

tarea muy compleja, que consume mucho tiempo y poco practica. Por tal
razn, se desarrollo un simulador matemtico de dos dimensiones, basado
en el mtodo de diferenciales finitas, que combina los efectos de
acoplamiento e inerciales, utilizando correlaciones generalizadas
desarrolladas recientemente. Este modelo de pozo, que simula el flujo en
estado de equilibrio del gas y condensado alrededor del pozo fracturado
hidrulicamente, tambin tiene en cuenta el cambio de fase y los sistemas de
baja tensin interfacial.

El fracturamiento hidrulico de un pozo es la tcnica ms conocida para
mejorar la productividad y especialmente en yacimientos de gas condensado
apretados. La fractura reduce la resistencia al flujo alrededor del pozo
disminuyendo la presin drawdown y por lo tanto reduciendo el impacto
negativo del banco de condensado. Debido a la importancia y gran aplicacin
del fracturamiento hidrulico, los clculos de productividad en estos sistemas
ha sido de gran inters para muchos investigadores, cuyos estudios han
tenido por objeto determinar el mejoramiento de la productividad del pozo o
optimizar el diseo de la fractura, investigando para tal fin, el comportamiento
del flujo y la distribucin de presin alrededor de la fractura en estado estable
y en condiciones de estado pseudo estable. Los resultados obtenidos son
correlaciones para calcular la productividad del pozo, el factor dao y el radio
efectivo del pozo.

6.5.1 Simulador In-House

125
La Fig. 25 muestra un esquema de un pozo fracturado hidrulicamente,

cuadrado o cilndrico, con presin constante en el rea de drenaje. El
sistema puede dividirse en dos zonas, la zona fracturada y la zona sin
fractura. La zona fracturada es cilndrica y contiene la fractura y la matriz con
radio exterior igual a la longitud media de fractura. La zona sin fractura, solo
contiene la matriz.
En el desarrollo del simulador In-House se asume:
a. La fractura y la matriz son dos medios porosos uniformes diferentes.
b. El ancho de la fractura es constante.
c. La fractura ha penetrado en forma simtrica en ambas direcciones.
d. El flujo en la matriz del pozo es insignificante, comparada con el flujo
en la fractura del pozo.
e. El efecto del casing perforado en el rgimen de flujo es insignificante.
f. La fractura vertical ha penetrado a travs de toda la altura del
yacimiento, es decir, la fuerza de gravedad no se tiene en cuenta en
el sistema 2D.

6.5.1.1 Ecuaciones que gobiernan

Las ecuaciones que gobiernan son similares a las descritas para el flujo
alrededor de las regiones perforadas. (J amiolahmady et al.
43
).
La ecuacin de continuidad para sistemas de flujo de gas condensado en
estado estable es:
(128)
43
Jamiolahmady M., Danesh A. and Duncan D. B. Measurement and Modeling of Gas Condensate
Flow Around a Rock Perforation Transport in Porous Media Volume 63 Number 2, pp 323347,
(2006).
126
Las ecuaciones de flujo para la fase de gas y condensado son:

(129)
(130)
y son la peremabilidad relativa al gas y al condensado
respectivamente, afectada por la inercia y el acoplamiento. Se utilizo la
correlacin generalizada desarrollada por J amiolahamdy et al.
44

Combinando la Ec.128 con las anteriores se obtiene:
(131)

6.5.2 Sistema 2D para pozos fracturados hidrulicamente.

En los sistemas 2D con las suposiciones hechas anteriormente, para flujo en
dos dimensiones, la ecuacin anterior se puede escribir como:
(132)
Figura 25: Esquema de un pozo fracturado hidrulicamente

44
Jamiolahmady M., Danesh A., Tehnari D. and Sohrabi M. Variations of Gas/Condensate Relative
Permeability With Production Rate at NearWellbore: A General Correlation SPE 83960, SPE
ReservoirEvaluation&EngineeringJournal,Volume9,Number5,pp688697,(Dec2006).
127

Fuente: Mahdiyar H., J amiolahmady M. and Sohrabi M. A New Flow Skin Factor Formulation for Hydraulically
Fractured Wells in Gas Condensate Reservoir. Paper SPE 113595, 2008.

Las condiciones lmite que se requieren para solucionar la ecuacin
diferencial son las siguientes:

(133)
donde, es la longitud media del area de drenaje.

128
En la Ec. 132 la permeabilidad relativa y las propiedades de la roca y del

fluido, dependen del y de la velocidad del fluido y son funciones
complejas de la presin y de su gradiente. en condiciones de estado
estable, es decir, flujo total de la composicin constante, dependiente de la
presin local para una temperatura dada por la Ec. 123 (Cap. 6.4).

El mtodo de iteracin numrica de diferenciales infinitas, fue utilizado para
desarrollar la principal ecuacin diferencial no-lineales, Ec. 132 y la ecuacin
auxiliar, Ec. 123 (Cap. 6.4).
El modelo se dividi en 1000 bloques, 40 bloques en la direccin y 25
bloques en la direccin .
Los bloques del grid cerca al pozo y los que se encuentran cerca a la punta
de la fractura son ms finas que el resto de los bloques, para calcular con
precisin las variaciones bruscas de los parmetros de flujo en esta zona,
especialmente cerca a la cara del pozo, donde la velocidad del flujo es
mxima.

6.5.3 Sistema 1D, para hueco abierto

Como se menciono anteriormente, el objetivo principal de este anlisis es
desarrollar una frmula para calcular el radio efectivo del pozo en hueco
abierto para replicar el flujo en pozos fracturados hidrulicamente. La
aplicabilidad y precisin de la formula ha sido verificada y comparada por el
flujo del pozo en pozos fracturados hidrulicamente, predicho con el
129
simulador 2D para pozos fracturados hidrulicamente, con el sistema de

hueco abierto equivalente. Para simular el flujo de gas condensado en
sistemas hueco abierto, otro simulador in-house fue desarrollado. El sistema
1D de hueco abierto utilizado en este anlisis, consiste en yacimiento
cilndrico de una sola capa con radio externo .
La forma en una dimensin de la Ec. 131 es:
(134)
Presin interna y externa constante, son las condiciones lmite para resolver
esta ecuacin diferencial parcial no-lineal. El simulador 1D hueco abierto
posee un grid de 50 bloques, con finos bloques cerca del pozo.

6.5.4 Conductividad efectiva de fractura

En un pozo fracturado hidrulicamente, la conductividad de fractura, es
definida:
(135)
Donde, y son la permeabilidad absoluta de la fractura y de la matriz
respectivamente. es el espesor de la fractura y es la longitud media de
fractura. En sistemas de una sola fase y flujo Darcy el resultado de un pozo
fracturado hidrulicamente, es controlado por , asi como por el radio de
penetracin de fractura, . Sin embargo, para una sola fase, flujo
130
no-Darcy y sistemas de flujo de gas-condensado, la conductividad de fractura

absoluta es reemplazada por la conductividad de fractura efectiva.

6.5.5 Sistemas de flujo no-Darcy en una sola fase

Segn Forchheimer la ecuacin para la cada de presin es una funcin de
segundo orden de la velocidad superficial de la siguiente manera:
(136)
En esta ecuacin el vector es el factor inercial y es el valor absoluto del
vector de velocidad. La Ec. 136 puede ser ajustada de la siguiente manera:
(137)
Por lo tanto,
(138)
donde, es llamada funcin no-Darcy y se relaciona con el nmero de
Reynolds de la siguiente manera:
(139)
La comparando de la ecuacin 138 con la ecuacin de Darcy es:
(140)
131
La permeabilidad efectiva debido al flujo no-Darcy o efectos inerciales es

. En otras palabras, es la permeabilidad relativa de una sola fase
con rgimen de flujo no-Darcy.

En un pozo fracturado hidrulicamente, dentro de la matriz, debido a la baja
velocidad, el efecto inercial es insignificante. Sin embargo, dentro de la
fractura, el fluido (gas o aceite) comnmente fluye con altas velocidades, de
modo que el efecto inercial podra ser muy importante y puede disminuir la
conductividad efectiva de fractura significativamente.

Muchos investigadores han estudiado este tema y han presentado algunas
correlaciones para predecir el efecto inercial, como Guppy et al
45
que
desarrollo la siguiente correlacin:
(141)
donde, es al constante de flujo adimensional y se define como:
(142)
donde, es la tasa total de flujo de ambas extremos de la fractura.
La Ec. 141 puede reescribirse en trminos de de la siguiente manera:
(143)
Guppy et al.
46
, simularon la fractura de un pozo y presentaron la siguiente
correlacin:
45
Guppy K. H., CincoLey H., Ramey Jr. H. R. and Samaneigo V. F. NonDarcy Flow in Wells With
FiniteConductivityverticalFracturesSPE8281,SPEJournal,pp.681,(Oct.1982).
132
(144)
J . L. Gidely
47
propuso que el efecto inercial reduce la conductividad efectiva
de la fractura de la siguiente manera:
(145)
Donde, se estima basado en las propiedades de flujo en la cara del pozo.
A partir de la Ec.143 se indica que el promedio del efecto inercial en la
fractura (representado por el numero de Reynolds) es 62% de la entrada de
la fractura y la cara del pozo. Por lo tanto, en una sola fase y sistema de flujo
no-Darcy la conductividad efectiva de fractura, puede calcularse de la
siguiente manera:
(146)
donde, es el promedio de la funcin no-Darcy en la fractura, en otras
palabras, es la permeabilidad relativa promedio de flujo en una sola fase en
la fractura. por lo tanto, la Ec. 146 puede reescribirse:
(147)
donde, y son la permeabilidad efectiva de la fractura y la matriz
en un sistema que fluye en una sola fase, respectivamente.

En un sistema fluido multifase, la relacin de movilidad de las fases fluyendo
es muy importante. Si se reemplaza la movilidad absoluta por la movilidad
46
Ibid
47
Giddley J. L. A Method for Correcting Dimensionless Fracture Conductivity for NonDarcy Flow
EffectsSPE20710,SPEProductionEngineering,pp.391(Nov.1991).
133
efectiva, la conductividad absoluta de fractura cambia a la conductividad

efectiva de fractura.

La Ec. 147 puede ser reescrita de la siguiente manera:
(148)
En esta ecuacin, que se deriva del sistema de flujo en un sola fase, y
son las relaciones de movilidad en la fractura y en la matriz,
respectivamente.

6.5.6 Flujo en
sistemas de gas condensado

6.5.6.1 Relacin
de movilidad en la fractura

Dentro de la fractura la relacin de movilidad promedio se define como:
(149)
donde, los subndices y son gas y condensado, respectivamente. La
mayor dificultas para utilizar esta ecuacin es la estimacin de la
134
permeabilidad relativa promedia que es afectada por el acoplamiento y los

efectos inerciales.

Se evalu el impacto de la variacin de , el acoplamiento y la inercia a
travs de una serie de simulaciones en el flujo del gas condensado en el
interior de pozos fracturados hidrulicamente para diferentes diseos de
fractura y yacimientos. La Fig. 26 muestra la variacin de la permeabilidad
relativa al gas de la permeabilidad base relativa al gas, a lo largo de la
fractura. Esta variacin muestra dos puntos clave:

1. La
variacin de la permeabilidad relativa base dentro de la fractura es
insignificante, es decir, no hay cambio significativo de fase en el
interior de la fractura. en otras palabras, se puede asumir que es
constante a lo largo de la fractura y es igual a las condiciones del
pozo y la permeabilidad relativa base y las propiedades del fluido se
pueden estimar de las propiedades del pozo.
2. En el
interior de la fractura debido a las altas velocidades de flujo, el efecto
inercial es dominante y el efecto de acoplamiento es insignificante.

Con base en estos resultados, en lugar de la Ec. 149, para el clculo de la
relacin de movilidad promedio en la fractura se puede calcular a partir:
(150)
Donde,
135
(151)
Donde el subndice , se refiere a la base es la relacin de movilidad
base en la fractura, depende de la permeabilidad base relativa al gas, a las
condiciones del pozo.
Figura 26: Variacin de la permeabilidad relativa al gas y de la
permeabilidad relativa al gas base dentro de la fractura.

Fractured Wells in Gas Condensate Reservoir. Paper SPE 113595 prepared for presentation at the 2008 SPE
Europe/EAGE Annual Conference and Exhibition held in Rome, Italy, 9-12 J une 2008.

La permeabilidad base relativa al condensado, puede calcularse a partir:
(152)
136
Donde y las viscosidades son calculados a condiciones de pozo. En la

Ec. 150, el numero de Reynolds, para el flujo de gas condensado es
calculado as:
(153)
Donde, y , son la densidad y viscosidad promedio, son estimadas por
promedio volumtrico de las propiedades de dos fases.

(154)
En la Ec. 153, es el factor inercial en una sola fase y es la velocidad a la
entrada de la fractura en el pozo:
(155)
donde, es la tasa total de flujo en el pozo a condiciones de fondo del pozo.

6.5.6.2 Relacin
de movilidad en la matriz

La Fig. 27 muestra la variacin de la permeabilidad relativa al gas y la
permeabilidad base relativa al gas en la matriz. De esta figura se obtienen 4
puntos importantes:

137
1. En la
matriz de la zona fracturada, que contiene una fractura con una
conductividad efectiva aceptable, la variacin de la permeabilidad
base relativa, representa el efecto de cambio de fase en ese medio,
es insignificante. As la movilidad relativa en la matriz se puede
estimar basados en a condiciones de pozo.
2. Lejos del
pozo, la permeabilidad relativa es casi constante y paralela a la
permeabilidad relativa base. Este se debe al impacto de la tensin
interfacial, que debe contabilizarse.
3. Cerca del
pozo, el aumento de la velocidad mejora la permeabilidad relativa al
gas, pero su efecto en la permeabilidad promedio relativa al gas
dentro de la matriz es insignificante.
4. La forma
de la variacin de en la matriz de un sistema de pozos fracturados
hidrulicamente y en un sistema equivalente hueco abierto es el
mismo y la diferencia entre sus valores es insignificante.
De acuerdo a todo lo anterior, la relacin de movilidad en la matriz se puede
calcular a partir:
(156)
donde, y son la permeabilidad relativa al gas y al condensado
en la matriz respectivamente, a condiciones del pozo en el sistema
equivalente hueco abierto.

138
Por lo tanto, la conductividad efectiva de fractura en un sistema de gas

condensado, es estimada por la siguiente ecuacin:
(157)
donde, , y son calculados por las Ecs. M,L y ; respectivamente.

Figura 27: Variacin de la permeabilidad relativa al gas y de la
permeabilidad relativa al gas base dentro de la matriz

139
Fractured Wells in Gas Condensate Reservoir. Paper SPE 113595 prepared for presentation at the 2008 SPE
Europe/EAGE Annual Conference and Exhibition held in Rome, Italy, 9-12 J une 2008.

6.5.7 Factor
dao de fractura y radio efectivo del pozo

En sistemas de gas condensado, un pseudo factor dao de fractura ( ), no
tiene en cuenta el factor dao debido al dao, es definido a por la siguiente
ecuacin:
(158)
En esta ecuacin, es el diferencial de pseudopresiones y se define:
(159)
Donde, y son la permeabilidad relativa al gas y condensado en la
matriz. La Ec. 158 puede escribirse con base al radio del pozo de la siguiente
manera:
(160)
La relacin entre el pseudo factor de dao y el radio efectivo del pozo se
puede obtener de la siguiente manera:
(161)
140
Para sistemas de una sola fase y flujo Darcy, Mahdiyar et al.

48
, introdujo una
correlacin para estimar el factor dao de fractura en pozos fracturados
hidrulicamente con rea de drenaje cuadrada y en estado estable:
(162)
Donde
(163)
En esa ecuacin es el radio externo adimensional, que representa el
efecto del radio de penetracin en . es el pseudo factor dao,
representa la resistencia al flujo Darcy en la zona fracturada.
Combinando las Ecs. 161 y 162 obtenemos una ecuacin para el radio del
pozo:
(164)
Se propone ampliar el uso de las ecuaciones anteriores para sistemas gas-
condensado, reemplazando la conductividad absoluta de fractura por la
conductividad efectiva de fractura. Entonces la correspondiente ecuacin
para la estimacin del pseudo factor dao y el radio efectivo en condiciones
de estado estable son:
(165)
(166)
48
Mahdiyar H. Jamiolahmady M. and Danesh A. New Mechanical and Damage Skin Factor
Correlations for Hydraulically Fractured Wells SPE 107634, European Formation Damage
Conference,Scheviningen,Netherland(2007)
141

6.5.7.1 Mtodo iterativo para estimar el radio efectivo del pozo

1. La primera estimacin del radio efectivo del pozo, se calcula con base
al radio efectivo Darcy del pozo de la siguiente manera:

2. La distribucin de presin y el clculo del flujo del pozo se realiza con
base al radio efectivo del pozo.
3. La conductividad efectiva de fractura es calcula usado la siguiente
ecuacin:

Donde cada uno de los parmetros se puede calcular a partir de las
ecuaciones anteriormente descritas.
4. El nuevo
radio efectivo del pozo, es calculado utilizando la siguiente ecuacin:

5. Si la diferencia entre el sucesivo radio efectivo del pozo es un valor
significativo, los clculos se repiten desde el paso 2, de otra forma se
calcula la distribucin de presin y otros parmetros de flujo y se
reportan.


142
El mtodo propuesto es aplicable tanto para una fase como para sistemas de
gas condensado. Por ejemplo, para el gas seco, la permeabilidad relativa al
gas y al condensado son 1.0 y cero respectivamente y la formulacin de la
conductividad efectiva de fractura es simplificada de la ecuacin de Guppy.
Adems para sistemas de baja velocidad, el numero de Reynolds podra
aproximarse a cero y la formula se convierte en un sistema de flujo Darcy de
una sola fase.

Para verificacin de la precisin del mtodo desarrollado para sistemas de
gas condensado, se simulo la conductividad por cuatro diferentes sistemas
de fracturas hidrulicas, en el simulador in-House 2D para pozos fracturados
hidrulicamente y en el simulador 1D para hueco abierto. En estas
simulaciones se asume que el radio interno es el modelo de hueco abierto
1D es el radio efectivo del pozo en pozos fracturados hidrulicamente y es
calculado mediante el mtodo iterativo explicado anteriormente. Todos los
parmetros en el modelo hueco abierto y el pozo fracturado hidrulicamente
son los mismos.
Los resultados obtenidos de estas simulaciones mostraron una buena
exactitud del mtodo, que se puede atribuir a que este tiene en cuenta el
efecto de la velocidad dominante dentro de la factura como un efecto inercial,
mientras que el efecto de la velocidad en la matriz es insignificante. Muestra
adicionalmente que dentro de la matriz la principal fuerza es la tensin
interfacial.

6.6 MODELO SIMPLIFICADO DE FLUJO CERCA AL POZO
143

La prediccin del rendimiento de yacimientos a largo plazo con un modelo de
pozo realista est lleno de incertidumbres debido a la complejidad del flujo
bifsico. Esto se presenta a que incluso en modelos bifsicos bien
estructurados se aparta del desempeo esperado cuando ocurren cambios
en las condiciones de flujo. Estos cambios incluyen, tasa gas/liquido, presin
en la cabeza del pozo y presin en las lneas de flujo con el tiempo, entre
otros, el influjo de agua es el problema ms grave en la prediccin.

En este anlisis se explorara la posibilidad de utilizar mtodos simplificados
para calcular la presin de fondo fluyendo (BHP), la presin en la cabeza del
pozo (WHP), las tasas medidas, la gravedad de los fluidos producidos y las
dimensiones de la tubera.

6.6.1 Enfoque computacional y resultados

Los datos utilizados para este anlisis provienen de la literatura (Govier y
Fogarasi 1975; Preffer et al. 1988) y algunos campos de frica Occidental
con tasa de produccin de gas/condensado de media a alta. En total se
examinaron 167 pruebas independientes. Los mtodos utilizados fueron:
Gray (Aziz et al. 1972), homogneos (Ansari et al. 1994) estado estable
OLGA, y el enfoque para el gas hmedo desarrollado por Peffer et al (1988).
El termino aceleracin, la friccin y la cabeza hidrosttica estn implcitos en
todos los mtodos reportados para este anlisis.
144

6.2.1 Ejemplo de aplicacin

Se consideraron tres sets de datos para este anlisis. El primer set obtenido
del Occidente de frica, y los otros dos reportados por Govier and Fogarasi
(1975) y por Peffer et al (1988).

6.6.2.1 Datos Occidente de frica. Unos 25 conjuntos de datos estn
disponibles de pruebas drillstem. La tasa gas/condensado vario entre 50 y 37
, mientras que la presin oscilo entre 5100 y 7300 . Todos
los mtodos muestran resultados comparables con el mtodo de Aziz et al.
49

Se nota una gran similitud entre el modelo homogneo y el modelo de Ansari
et al.
50

Figura 28: Porcentaje de error para los datos del Occidente de frica

49
Aziz,K,Govier,G.W.,andFogarasi,M,PressureDropinWellsProducingOilandGas,Cdn.J.Pet.
Tech.(JulySeptember1972)38.
50
Ansari,A.M.etal.AComprehensiveMechanisticModelforUpwardTwoPhaseFlowin
WellboresSPEPF(May1994)143;Trans,AIME,207.
145

Fuente: Kabir C.S. and Hasan A. R. Simplified Wellbore-Flow Modeling Gas/Condensate Systems. Paper SPE
89754

La Fig. 28 muestra que el error aleatorio determinado por los diferentes BHPs
medidos y calculados por el mtodo del gas hmedo es aparente. El triple
cambio en el valor de la relacin gas/aceite ( ) no se espera que ocurran
en estos yacimientos, dada por la similitud en las propiedades presin
/volumen /temperatura (PVT). Los datos PVT medidos en el laboratorio y el
modelo de la ecuacin de estado muestran que ocurre un pequeo aumento
en el contenido de liquido al aumentar la profundidad. Se espera que el error
de BHP disminuya con el aumento del , porque se asume flujo constante
en una sola fase. Cualquier error en la medicin de la tasa de liquido puede
traducirse en grandes errores en el calculo de BHP. Las mediciones del
liquido estn sujetas a largos errores en comparacin con el gas, debido al
pequeo volumen medido en un sistema tpico gas/condensado.

6.6.2.2 Datos Govier-Fogarasi. Se utilizaron 96 pruebas de las 102
reportadas por Govier y Fogarasi. La Fig. 29 muestra la tendencia esperada,
es decir, los errores de la BHP tienen a disminuir con el aumento del .

146
6.6.2.3 Datos comisin de Railroad de Texas (TRRC). En este conjunto

de datos se realiza un estudio a fondo, debido a las grandes variaciones del
(5180 a 780000 ), la presin en la cabeza del pozo (3107 a
12313 ) y la presin de fondo fluyendo (4125 a 11250 ). En las
graficas del porcentaje de error de la presin de fondo fluyendo con , la
presin en la cabeza del pozo y la tasa, no se encontr la tendencia.
Posiblemente la medicin de la tasa es la responsable de que no haya
encontrado la tendencia.

6.6.3 Incertidumbre de los datos de entrada para calcular la presin de
fondo fluyendo

6.6.3.1 Precisin absoluta de dos medidores. Todas la presiones de
recorrido calculadas presuponen que los medidores en la cabeza del pozo
tienen las mismas caractersticas as como el lugar en el fondo del pozo,
especialmente en lo relacionado con la precisin absoluta. En la practicas
estas condiciones rara vez se cumple. Debido a lo anterior se examinaran.

Figura 29: Porcentaje de error para los datos Govier-Fogarasi

147

Fuente: Kabir C.S. and Hasan A. R. Simplified Wellbore-Flow Modeling Gas/Condensate Systems. Paper SPE
89754

6.6.3.2 Tasa de liquido o CGR tendencia/eficiencia del separador. La
produccin de liquido es relativamente pequea en un pozo tpico
gas/condensado. en consecuencia la separacin de este liquido de grandes
volmenes de gas producido, no es muy eficiente. La tendencia creciente de
la presin en cabeza de pozo se espera que se produzca a causa de un
almacenamiento trmico (Hasan et al. 2005) y no es una anomala. Al
parecerla separacin de altas tasas de liquido es inadecuada debido al
excesivo choque causando tendencias opuestas cuando se mide el flujo de
fluidos en el separador. Se garantizo que todas las tasa se encontraran por
encima de la tasa critica de descarga en los datos del Occidente de frica. El
uso del modelo de Li et al.(2002) mostro que la tasa critica vara entre 1.77 y
2.61 , esto prueba el continuo aumento del condensado.

6.6.3.3 Rangos de la gravedad del condensado. Moses (1986) reporto que
el rango de la gravedad del liquido debera ser de 40 a 60 . Otros
148
autores han sugerido como el limite superior. Sin embargo los datos
de los campos no pueden fijarse a este rango. Por ejemplo los 120 pruebas
reportadas de Govier y Fogarasi, 56 tenan mayor que 70 y 5 tenan
valores menores que 40. Se supone que la contaminacin del fluido es la
responsable de esta anomala, la presencia de gas disuelto es una condicin
que podra aligerar el liquido en el tanque. Para los clculos realizados se
establece que los valores lmite superior e inferior fueron violados.

6.6.3.4 Pared rugosa. Este tem es otra posible fuente de error. Para los
clculos se utilizo una tubera con una pared rugosa de 0.0018 . Al hacer
esto, se propone que todos tubos se realizaron de acero comercial, sin
ningn recubrimiento contra la corrosin o revestimiento. Esta suposicin
funciona muy bien en los pozos del Occidente de frica, pero los pozos de
los otros dos conjunto de datos han desconocido caracteres.

6.6.4 Justificacin del modelo homogneo

Los resultados presentados muestran claramente que las estimaciones de un
modelo simple homogneo son comparables con los datos obtenidos al
aplicar los modelos de Ansari et al
51
. , Gray, OLGA, and Aziz et al
52
. La
explicacin del rendimiento de un modelo homogneo, frente a los modelos
sofisticados, se muestra a continuacin:

51
Ansari,Op.cit,141
52
Aziz,Op.cit,141
149
6.6.4.1 Comparacin con la correlacin de Gray. La correlacin de Gray

es utilizada en enfoques pseudohomogeneos para el clculo de las perdidas
en cabeza por friccin, esttica y aceleracin, sin embargo Gray indico que
incluso para pozos de condensados, las gotas de liquido no se mueven tan
rpido como la fase gaseosa, es decir, hay algo de liquido estancado. Para
tener en cuenta este estancamiento del liquido, Gray propone la siguiente
expresin para la fraccin de volumen in-situ.
(167)
Donde
(168)
y y son velocidad adimensional y el numero dimetro.

En un modelo homogneo, . En otras palabras, la correlacion de
Gray es idntica a un modelo homogneo cuando el termino elevado al
exponente desaparece. En los clculos realizados el primer termino de la
Ec.167 vario entre 0.004 y 0.01. Por lo tanto, se puede concluir que la
correlacin de Gray colapse en el modelo homogneo para la mayora de los
pozos gas/condensado.
6.6.4.2 Comparacin con el modelo Ansari et al. En el flujo anular, la
nica diferencia entre el modelo Ansari et al. (1994) y el modelo homogneo
es la manera como se calcula el factor de friccin. El modelo Ansari et al., el
factor de friccin de la fase gaseosa, , que fluye solo en el pozo, es
calculado inicialmente usualmente a partir de la relacin numero de Reynold
( )/ factor de friccion., Posteriormente, Wallis (1969) sugiri, que
150
ese factor de friccion se multiplica por un factor que tiene en cuenta la

rugosidad, de la fase gaseosa que fluye en contra de una pelcula de liquido
a lo largo de la pared del tubo. Este factor multiplicador depende del
espesor adimensional de la pelcula de liquido, dada por ( ). Es decir,
(169)
Ansari et al., propuso un mtodo para estimar el espesor de la pelcula, pero
en los clculos realizados se utilizo una mtodo ms sencillo propuesto por
Lockhart-Martinelli, obtenindose un espesor adimensional de la pelcula de
liquido menor de 0.015 para los campos analizados. Sin embargo el clculo
de friccin por el mtodo de Ansari et al. fue ligeramente superior al obtenido
por el modelo homogneo. Por lo tanto, el clculo del factor de friccin por el
modelo homogneo es muy similar al obtenido por el modelo de Ansari et al.

6.6.4.3 Comparacin con el modelo de Aziz et al. Aziz et al. (1972), utilizo
la correlacin de flujo de Duns y Ros (1963). Duns y Ros supusieron que las
gotas de liquido no influyen en el gradiente de presin. En otras palabras un
modelo simple de flujo de gas es adoptado. En este caso se utiliza el mismo
factor de friccin, similar al del modelo de Ansari et al., pero sin el factor de
correccin presentado en la Ec. 169. Es decir, el modelo podra resultar en
bajas estimaciones de prdidas de presin en pozos de condensados,
comparados con el modelo de Ansari et al. En los clculos se encontraron
esas bajas estimaciones de las prdidas de presin por el modelo de Aziz et
al. para la mayora de los datos analizados. Sin embargo para algunos datos,
el modelo de Aziz et al. Predice un rgimen de flujo abundante y por lo tanto,
calcula altsima esttica en cabeza y grandes prdidas de presin que otros
modelos.

151


El motivo principal de este anlisis fue establecer un confiable y simple
modelo para calcular las cadas de presin en la cara del pozo, para
sistemas de gas/condensado. La identificacin de este modelo permite la
generacin de tablas de flujo que ayudan a predecir el comportamiento de
flujo de dos y tres fase en el yacimiento, la realizacin de anlisis nodales y
el monitorio de la presin en la cabeza del pozo, para operar el yacimiento en
tiempo real.

El uso de un modelo homogneo en un modelo integrado que involucra,
yacimiento-cara de pozo-tubera, es muy interesante, debido a su simplicidad
y a que es muy robusto. La comparacin del desempeo del modelo
homogneo con otros modelos sugiere que este es una herramienta
confiable. De hecho, tanto el modelo Gray como el modelo Ansari et al.
convergen a un modelo homogneo cuando se presenta el flujo mist en
sistemas de gas-condensado.

152
7. METODOLOGA PARA OBTENER UN MODELO DE POZO

Despus de realizar una revisin de algunos de los modelos que predicen el
comportamiento de los yacimientos de gas condensado, de clasificarlos de
acuerdo a la(s) propiedad (es) que evala el modelo. Se propone una
metodologa terica para obtener un modelo de pozo que permita predecir el
comportamiento de la produccin de un pozo de gas.

En los yacimientos de gas condensado, que presentan un comportamiento
tan complejo, debido, a la presencia de un sistema fluido de dos fases,
cuando la presin cae por debajo de la presin del punto de roco, los
modelos de pozo, son una buena alternativa para predecir del
comportamiento del pozo, porque, al simular todo el yacimiento se pueden
obtener resultados errneos.

Durante la produccin de un yacimiento de gas condensado, ocurre una
cada de presin que origina principalmente tres problemas:
Reduccin de la productividad del pozo.
Disminucin de la cantidad de gas.
Acumulacin de condensado alrededor del pozo.

El modelamiento de la produccin de un pozo, se lleva a cabo mediante el
establecimiento e integracin de:
153
1. Modelo fsico: consiste en realizar una descripcin detallada del medio

poroso.
2. Modelo de fluidos: consiste en realizar una descripcin detallada del
flujo de fluidos en el medio poroso.
3. Modelo numrico: con base en el modelo fsico y el modelo de fluidos,
se puede propone el modelo numrico que es el mejor representante
de las condiciones descriptas en los modelos anteriores.
4. Desarrollar el modelo propuesto en el paso anterior, mediante la
utilizacin de algn simulador composicional o del simulador de su
preferencia.

El modelo fsico y el modelo de fluidos son la base fundamental y definirn el
tratamiento aplicar. Los modelos fsicos pueden ser: Yacimiento homogneo,
multicapa o naturalmente fracturado. los modelos de fluidos pueden ser:
Aceite negro, aceite voltil, gas condensado y composicional. Se recomienda
analizar el yacimiento, desde lo convencional hasta lo complejo, modificando
uno a uno los parmetros.

Para modelar un yacimiento de gas condensado, se hace necesario, revisar
los estudios previos realizados, a continuacin se exponen algunos de esos
trabajos:

Roebuck et al.
53
, desarrollaron un mtodo numrico implcito descrito por la
simulacin radial de flujo de tres fases en un yacimiento de petrleo. El
mtodo considera el efecto de la variacin de las propiedades fsicas de la
53
RoebuckI.F.,FordW.T.HendersonG.E.andDouglasJim.TheCompositionalReservoirSimulator:
CaseIIITheRadialGeometry.PaperSPE2486.
154
roca y del sistema de fluidos. El efecto de la transferencia de masa y el

cambio de la composicin del fluido sobre la distribucin de la saturacin-
presin son considerados a travs del concepto de comportamiento de fase.
El autor introduce el concepto de simulacin de yacimiento composicional,
describiendo en una dimensin, flujo de tres fases a travs del medio poroso.
La ecuacin diferencial del flujo de masa para el flujo composicional de
hidrocarburos es descrita por:
(170)
Donde se asume que la presin del gas y del aceite son iguales. Se plantea
una ecuacin similar para la fase agua.

La produccin de algunos pozos de gas se reduce severamente cuando la
presin de fondo fluyendo es menor que la presin de saturacin del fluido
in-situ. Generalmente, esa reduccin es debida a la acumulacin de liquido
condensado cerca a la cara del pozo y bajas tasa de produccin por
acumulacin de condensado en el tubing, el anular o ambos, dependiendo
del tipo de completamiento. Fussell
54
desarrollo un modelo de pozo radial
que muestra el efecto de la acumulacin de condensado en la productividad
del pozo, evala la aplicabilidad del mtodo de estado estable y evala el
efecto del equilibrio de fase y la permeabilidad relativa caracterstica sobre la
prediccin del rendimiento.

Katz et al.
55
, desarrollaron un modelo de simulacin radial, para estudiar los
yacimientos de baja permeabilidad, particularmente el almacenamiento de
54
FussellD.D.SingleWellPerformancePredictionsforGasCondensateReservoir.PaperSPE4072.
55
KatzDonald,CamposJ.AntonioandSinghS.P.ApplicationofRadialSimulationModel.PaperSPE
5883.
155
gas. El simulador radial desarrollado, es un programa digital de computacin

que resuelve flujo en una fase, ecuacin de flujo turbulento y el flujo natural
del gas en un yacimiento poroso cilndrico. La ecuacin de flujo es obtenida
de la combinacin de las ecuaciones de continuidad en direccin radial, flujo
turbulento, ecuacin de Darcy y la ecuacin de estado para los gases reales.
El programa de computador calcula el flujo natural de gas por dos casos:
cuando la presin de la cara del pozo se mantiene constante por un periodo
de tiempo y cuando el pozo produce a una tasa constante por un periodo.

Coats Keith
56
, desarrollo un modelo que trabaja para uno, dos o tres
dimensiones de flujo de agua, aceite y gas, en , y en formaciones de
permeabilidad y porosidad heterogenea con geometra cartesiana ( )
o cilndrica ( ). El flujo de fluidos es simulado utilizando la ley de
Darcy, incorporando las fuerzas de gravedad, viscosidad y capilares. La
permeabilidad relativa y la presin capilar son dependientes de la saturacin
y de la tensin interfacial. El modelo es aplicado a yacimientos de aceite
voltil y gas condensado. El modelo consiste en ecuaciones escritas en
diferencias finitas para cada bloque, donde es . El balance de masa
para el componente es:

(171)
Donde, es la transmisibilidad del flujo de fluidos , es la fraccion molar del
componente en la fase aceite, es la fraccin molar del componente en
56
CoatsKeithH.AnEquationofStateCompositionalModel.PaperSPE8284.
156
la fase gas, es el peso especifico y la profundidad. Plantearon una

ecuacin similar para la fase agua.
Dumor et al.
57
, desarrollaron un modelo basada en un mtodo caracterstico
presentado para el clculo en una dimensin (1D), condensacin de tres
componentes y vaporizacin del gas seco. El modelo describe: la
transferencia de masa entre el aceite y el gas, expansin y contraccin,
cambios en viscosidad y densidad, estabilizacin de la gravedad e
interaccin roca/fluido. La principal asuncin del modelo es el equilibrio local
termodinmico y la ausencia de dispersin, difusin y capilaridad. Las
principales aplicaciones del modelo son: el estudio conceptual del gas seco
en el cual la transferencia de masa juega un rol importante y la calibracin y
revisin del simulador numrico del yacimiento para flujo multicomponente y
multifasico. Para hacer el problema manejable, para el mtodo caracterstico,
ellos emplearon las siguientes asunciones y simplificaciones.
1. El flujo es 1D e isotrmico, gobernado
por la ley de Darcy para flujo multifasico, el cual puede incluir el
termino de gravedad.
2. Las fases de gas y aceite se encuentran
en equilibrio termodinmico.
3. Para el propsito de evaluar las
propiedades termodinmicas de los fluidos se considera que el
sistema se encuentra a temperatura y presin uniforme.
4. La capilaridad y los efectos de dispersin
en la direccin de flujo son insignificantes.
5. El medio poroso es homogneo y
contiene una saturacin de agua connata constante.
57
Dumor J. M., Hagoort J. and Risseeuw A. S. An Analytical Model for OneDimensional, Three
ComponentCondensingandVaporizingGasDrives.PaperSPE10069.
157

En estos sistemas, la conservacin de los componentes individuales, permite
el desarrollo de las siguientes ecuaciones:
(172)
Donde

(173)
Donde y es la densidad molar del gas y del aceite respectivamente,
es el flujo fraccional del aceite, es la fraccin molar del componente en la
fase de aceite y es la fraccin molar del componente en la fase gas.
(funcin de flujo fraccional del componente ) es funcin de toda la
composicin, ( donde es la fraccion de gas, expresada
como fraccion molar), y la tasa de flujo total, . La funcin (funcin
saturacin del componente ) depende de la composicin total solamente.
Entonces, y . El mtodo
caracterstico por medio de ecuaciones diferenciales de primer orden de tipo
hiperblico.

Los yacimientos de gas condensado son simulados frecuentemente con
modelos totalmente composicionales. Coats Keith
58
presenta un proceso de
pseudoizacion que reduce la condensacin de un fluido multicomponente a
una pseudo mezcla de dos componentes, en superficie gas y aceite. Esto
58
CoatsKeithH.SimulationofGasCondensateReservoirPerformance.PaperSPE10512.
158
permite el uso de modelo modificado de aceite negro ms simple y ms

econmico, que tiene en cuenta, tanto gas disuelto en aceite como vapor de
aceite disuelto en gas. La ms importante a tener en cuenta en la utilizacin
de un modelo de aceite negro es si la descripcin de los dos componentes
puede representar adecuadamente el fenmeno composicional durante la
deplecin o el ciclaje en yacimientos de gas condensado. El termino
pseudizacion denota la reduccin en el nmero de componentes utilizados
en los clculos de las ecuaciones de estado (EOS) para fluidos de
yacimiento. La pseudoizacion es importante en los clculos del yacimiento,
debido al gran nmero de componentes reales en los fluidos del yacimiento.
Los clculos composicionales realizados utilizan variaciones de las
ecuaciones de estado (EOS) de Redlich-Kwong y Peng Robinson.

Clark T. J .
59
utilizo un modelo radial composicional 2D, para simular la
produccin en el yacimiento Nugget en la unidad este de Anschutz Rach, un
yacimiento de gas condensado rico. El modelo de pozo utilizo una
descripcin de un yacimiento heterogneo con 11 capas verticales, con la
caracterizacin de nueve componentes del fluido para predecir el
comportamiento de un pozo productor y el desarrollo de la saturacin de
condensado antes y despus de la ruptura de la inyeccin de gas (utilizada
para mantener la presin del yacimiento). El trabajo de simulacin indico que
el aumento de la productividad del gas experimentado despus de la ruptura
fue consecuencia de la disminucin de la viscosidad del gas y la alta
permeabilidad relativa al gas en la capa, y una sola capa puede ser
responsable de la mayora del aumento de la productividad.
59
ClarkT.J.TheApplicationof2DCompositional,RadialModelToPredictSingleWellPerformance
inaRichGasCondensateReservoir.PaperSPE14413,preparedforpresentationatthe60
th
Annual
Technical Conference and Exhibition of the Society of Petroleum Engineers had in Las Vegas, NV
September2225,(1985)
159

Coats K. H.
60
, utilizo un modelo numrico implcito para la simulacin
composicional multidimensional, de tres fases de una sola porosidad y en
yacimientos naturalmente fracturados. El modelo es composicional con una
generalizacin cubica de la ecuacin de estado (EOS) para la representacin
del equilibrio de fase gas-aceite y las densidades. La ecuacin generalizada
representa las ecuaciones de estado de Redlich-Kwong, Soave-Redlich-
Kwong, Zudkevitch-J offe, y Peng-Robinson. Los modelos de simulacin de 1,
2 y 3 dimensiones, fluyen en coordenadas cilndricas o cartesianas. La ley de
Darcy es modificada por la permeabilidad relativa y la presin capilar
representada por las fuerzas viscosas, capilares y gravitacionales. El efecto
de la tensin interfacial y la presin capilar son incluido. Para el caso de
doble porosidad, el modelo permite la descripcin del yacimiento donde el
tamao de los bloques de la matriz y la forma varia de bloque a bloque a
travs del grid. La acumulacin y el tiempo de computo requerido son
proporcionales de acuerdo a la fractura.

Al-Majed et al.
61
, desarrollaron un modelo de celda variable para simular el
rendimiento en un yacimiento de gas condensado. El modelo aproxima las
prediccin de los modelos convencional composicional, pero solo requiere
una fraccin del tiempo de computo. La tendencia caracterstica del perfil de
saturacin de aceite, que se presenta cuando los pozos de gas condensado
producen a una presin por debajo del punto de roco, se utilizo para dividir
el yacimiento en tres regiones concntricas. Cada regin es simulada por
60
Coats K. H., Implicit Compositional Simulation of SinglePorosity and DualPorosity Reservoir.
Paper SPE 18427, prepared for presentation at the SPE Symposium on Reservoir Simulation in
Houston,TX,February68(1989).
61
AlMajed A. A., Fahd King and Dougherty E. L. A Variable Cell Model for Simulation Gas
Condensate Reservoir Performance Paper SPE 21428, prepared for presentation at the SPE Middle
EastOilShowheldinBahrain,1619November(1991).
160
una celda y el volumen de cada celda varia continuamente como la

produccin. El modelo es utilizado para estudiar el efecto del flujo de liquido
en la composicin del fluido de la corriente del pozo y la influencia de los
fluidos del yacimiento en el equilibrio de fases. La descripcin matemtica
de este modelo se basa en las siguientes suposiciones:
1. El pozo produce a una tasa constante, es
isotrpico, homogneo, cilndrico, limites del yacimiento cerrado y
espesor constante.
2. Porosidad y permeabilidad son
independientes de la presin y la temperatura y se mantienen
constantes.
3. La fase agua es inmvil.
4. Las fuerzas gravitacionales y capilares
entre el aceite, gas y agua son insignificantes.
5. El flujo de fluidos entre la celda uno y la
cara del pozo, entre la celda dos y la celda uno puede ser aproximado
por la ecuacin de flujo estable. El flujo de fluidos entre la celda tres y
la celda dos puede ser aproximado por la ecuacin de estado
pseudoestable.
6. La interfase de intercambio
termodinmico en el yacimiento es rpido comparado con el flujo de
fluidos, es decir, el gas y el aceite en cada celda se encuentran en
equilibrio de fase.

Gmez y Correa
62
desarrollaron un modelo totalmente implcito, isotrmico y
en dos dimensiones para modelar yacimientos de gas condensado
62
Gomez Hamilton P. and Correa Antonio Claudio Fully Implicit Compositional Modelling of Gas
Condensate and Volatile Oil Reservoir. Paper SPE 23700 prepared for presentation at the 2
nd
Latin
AmericanPetroleumEngineeringConferenceheldinCaracas,Venezuela,March811,(1992).
161
retrogrado y aceite voltil. El modelo considera las fase de aceite, gas y agua
y puede utilizarse cualquier grid, cartesiano ( ) o cilndrico ( ). El
agua se considera inmvil y ligeramente compresible y la saturacin de agua
es funcin de la presin. El equilibrio termodinmico es considerado pero la
interaccin roca-fluido no. El equilibrio de fase y las propiedades de los
fluidos son calculados por la ecuacin de estado de Peng Robinson. El
coeficiente de ajuste es determinado previamente por un simulador comercial
PVT. El sistema resulta en ecuaciones diferenciales parciales, que se han
discretizado por el mtodo de diferencias finitas, produciendo un sistema de
ecuaciones algebraicas no lineales, este sistema se resuelve por el mtodo
numrico de Newton-Raphson. El modelo matemtico considerado involucra
ecuaciones de flujo, equilibrio termodinmico, relacin de presin capilar y
saturacin y restricciones de composicin. El flujo de fluidos a travs del
medio poroso es descrito combinando el principio de conservacin de masa y
la ley de Darcy, considerando las siguientes suposiciones:
Flujo en medio poroso isotrmico.
No reacciones qumicas en el sistema
roca-fluido.
Equilibrio termodinmico instantneo.
Efectos de dispersin insignificantes.

Con las suposiciones anteriores, la ecuacin resultante del balance de
materia para el componente hidrocarburo es:
(174)
162
Donde es la fraccion molar en la fase liquida, es la velocidad aparente del

fluido ( ), es la tasa molar ( ) y el numero de
componentes. Se propone una ecuacin similar para la fase agua.

Afidick et al.
63
, utilizaron un modelo de pozo composicional radial, para
estudiar el efecto de la acumulacin de condensado cerca del pozo sobre la
produccin y la presin transitoria y confirmar que las prdidas de
productividad se debe a la acumulacin de liquido. El modelo tambin
permite predecir el futuro rendimiento del pozo. En su trabajo se identifica la
acumulacin de condensado cerca a la cara del pozo, como un factor
extremadamente importante a considerar en el rendimiento futuro del pozo,
as como que la productividad es reducida en un 50%. El modelo de pozo
utilizado composicional en dos dimensiones consta de una sola capa de
propiedades homognea con 11celdas radiales de diferente ancho. Para este
caso el efecto de la acumulacin del condensado puede ser estimado a partir
de la siguiente ecuacin:
(175)
Donde, es la presin promedio del yacimiento, es la presin de fondo
fluyendo en , es el factor dao a , es el radio de la zona
interna, , es el radio de drenaje, es le radio del pozo, es la tasa
de flujo en , y es la permeabilidad
efectiva del gas a la saturacin critica del condensado.

63
AfidickDeddy,KaczorowskinN.J.andBetteSrinivas.ProductionPerformaceofaRetrogradeGas
Reservoir: A Case Study of the Arun Field. Paper SPE 28749, prepared for presentation at the SPE
AsiaPacificOil&GasConferenceheldinMelbourne,Australia,710November(1994).
163
Malachowski et al.
64
, propusieron un mtodo alternativo de simulacin que
calcula con precisin el efecto del banco de condensado en mltiples pozos,
simulacin total del yacimiento sin refinamiento local. Este mtodo se basa
en la utilizacin de la correlaciones de pseudo daos que puede ser
combinada con una ecuacin general del pozo, para capturar el efecto de las
prdidas de productividad del pozo. Las correlaciones se desarrollaron antes
de modelar todo el yacimiento. El objetivo del mtodo de pseudo dao es
poder determinar el nivel del bloque de condensado y las prdidas de
productividad, con base en parmetros de la corriente del pozo conocidos
como tasa de flujo total de la corriente, GOR y la composicin.

Fevang y Whitson
65
, propusieron un mtodo exacto para modelar la
entregabilidad en un pozo de gas condensado. Esta entregabilidad es
calculada con una modificacin de la ecuacin de Evinger Muskat de
pseudopresin para gas seco. A partir del mtodo propuesto, se demuestra
que un fino grid de simulacin en un modelo de pozo, permite reproducir casi
exactamente con una simple ecuacin de tasa que utiliza pseudopresiones,
la clave es conocer el GOR producido con bastante precisin. Ellos
concluyeron que la saturacin critica de aceite no tiene efecto directo sobre
la entregabilidad del pozo, la tensin interfacial depende de la permeabilidad
relativa. Adicionalmente, contemplaron tres regiones de flujo:
gas, en la cual ambos estn fluyendo simultneamente.
64
Malachowski M. A., Yanosik J. L. and Batten A. H., Simulation of Well Productivity Losses Due to
Near Well Condensate Accumulation in Field Scale Simulations. Paper SPE 30715, prepared for
presentation at the SPE Annual Technical Conference & Exhibition held in Dallas, U.S.A., 2225
October,(1995).
65
FevangOivingandWhitsonC.H.ModelingGasCondensateWellDeliverability.PaperSPE30714,
preparedforpresentationat the 2000 SPE Annual Technical Conference and Exhibition held in Dallas, U.S.A.,
22-25 October (1995)
164

condensado esta inmvil y solo el gas esta fluyendo.
Regin 3: una regin que contiene solo una fase de gas.

Con las observaciones, del flujo en las tres regiones, desarrollaron un
mtodo simple para calcular la integral de pseudopresin

All et al.
66
, modelaron el flujo de gas condensado cerca a la cara del pozo. La
simulacin exacta de las prdidas de productividad en un pozo, debido a la
acumulacin de condensado, generalmente requiere un modelo radial con un
grid fino y un mtodo de solucin totalmente implcito o adecuadamente
implcito. Los costos de computo por la inclusin de un refinamiento del grid,
necesario para capturar los efectos de la acumulacin de liquido cerca al
pozo en una simulacin composicional del campo, pueden ser extremos,
incluso para modelos con un refinamiento del grid y una capacidad implcita
adecuada. Un modelo efectivo de un yacimiento de gas condensado,
requiere los siguientes parmetros o fenmenos:
1. La saturacin critica del condensado y su
dependencia de la tensin interfacial y la velocidad de flujo cerca a la
cara del pozo.
2. La reduccin de la permeabilidad relativa
en la presencia de la condensacin inmvil.
3. La aparicin de la velocidad de
extraccin.
66
All J. K., McGauley P. J. and Wilson C. J. Experimental Studies and Modelling of Gas Condensate
FlowNeartheWellbore.PaperSPE39053,preparedforpresentationattheFifthLatinAmericanand
CaribbeanPetroleumEngineeringConferenceandExhibitionheldinRiodeJaneiro,Brazil,30August
3September(1997)
165
4. Flujo no-Darcy.
5. Modificaciones de las permeabilidades
relativas a altas velocidades y/o baja tensin interfacial.

Lee y Chaverra
67
modelaron e interpretaron el banco de condensado, este ha
sido reconocido como una de las principales causas de prdida de
entregabilidad en pozos de gas condensado. la entregabilidad de un pozo es
un tema importante para la evaluacin y desarrollo de yacimientos de gas
condensado. La prediccin precisa y a largo plazo de la entregabilidad del
pozo, puede alcanzarse a travs de un riguroso modelo composicional con
una buena descripcin de los parmetros de la roca, una ecuacin de estado
exacta para el fluido del yacimiento y un adecuado modelo mecanistico cerca
a la cara del pozo en el simulador, este modelo debe incluir los efectos del
flujo no-Darcy y las fuerzas viscosas que mejoran la movilidad del
condensado. este simulador computacional utilizo el simulador de yacimiento
comercial VIP de Landmark.

Blom y Hagoort
68
desarrollaron un mtodo numrico para calcular el deterioro
del pozo basados en el flujo radial de estado estable. El modelo incorpora la
permeabilidad relativa cerca al punto crtico que depende de la saturacin y
la resistencia inercial. Para el modelo fsico, se considero un yacimiento
cilndrico, de espesor constante y el pozo perforado totalmente en el centro.
El yacimiento es homogneo, isotrpico en permeabilidad, porosidad y
resistencia inercial, por simplicidad no se considera agua connata o dao a la
67
LeeShengTaiandChaverraMarcial,ModellingandInterpretationofCondensateBankingforthe
Near Critical Cupiagua Field. Paper SPE 49265, prepared for presentation at the 1998 SPE Annual
TechnicalConferenceandExhibitionheldinNewOrleans,Louisiana,2730September(1998).
68
BlomSaskiaM.P.andHagoortJacques,TheCombinedEffectofNearCriticalRelativePermeability
andNonDarcyFlowonWellImpairmentbyCondensateDropOut.PaperSPE51367.
166
formacin (dao mecnico cero). La presin en el lmite exterior del

yacimiento permanece constante y alrededor de la presin del punto de
roco. Produce a una tasa constante y en estado estable. En el modelo
matemtico, tienen en cuenta el balance de materia para flujo multifasico en
estado estable:
(176)
la ecuacin generalizada Darcy-Forchheimer para dos fase de flujo:
(177)
Entre otras, donde indica la fase, es el factor de correccin inercial
adimensional y es la velocidad superficial ( ).

Li y Firoozabadi
69
, estudiaron el efecto de la gravedad, fuerzas viscosas,
tensin interfacial y mojabilidad sobre la saturacin critica de condensado y
permeabilidad relativa al gas, utilizando un modelo simple fenomenolgico de
red. Los resultados mostraron que la mojabilidad afecta significativamente
tanto la permeabilidad relativa y la saturacin critica de condensado.
adicionalmente sugiere que la entregabilidad del gas en yacimientos de gas
condensado, puede ser mejorada por la alteracin de la mojabilidad cerca a
la cara del pozo.

69
Li Kewen and Firoozabadi Abbas, Phenomenological Modeling of Critical Condensate Saturation
andRelativePermeabilitiesinGas/CondensateSystems.PaperSPE56014.
167
El-Banbi et al.
70
, construyeron un modelo de pozo composicional radial para
investigar el comportamiento de la productividad del pozo. el modelo
consisti en una sola capa con 36 bloques en el grid en direccin radial. El
modelo fue limitado por la tasa de gas mientras que las propiedades del
yacimiento se modificaron para que coincidieran con la presin promedio del
yacimiento y la tasa de produccin del condensado. Los resultados de la
simulacin mostraron que la composicin del gas y del condensado cambian
con el descenso de la presin del yacimiento y el cambio alrededor del pozo
es ms dramticos que en el yacimiento.

Wheaton y Zhang
71
desarrollaron un modelo terico para describir la
formacin del banco de condensado. Este modelo tiene como prioridad tratar
la dinmica del banco de condensado. Ellos definieron que la saturacin de
condensado alrededor del pozo productor, no solo depende de la presin,
sino que tambin se debe tener en cuenta la concentracin total de los
componentes pesados y como varan estos con el tiempo cerca del pozo.
Para entender la dinmica de la acumulacin de condensado y el potencial
de re-vaporizacin del condensado, es esencial entender el
comportamiento del pozo y del campo en general.

Zhang y Wheaton
72
analizaron la dinmica del banco de condensado, desde
el comportamiento de la relacin de gas-condensado producido. Ellos
70
ElBanbiAhmedH.,McCainW.D.andSemmelbeckM.E.InvestigationofWellProductivityinGas
CondensateReservoir.PaperSPE59773,preparedforpresentationatthe2000SPE/CERIGas
TechnologySymposiumheldinCalgary,AlbertaCanada,35April2000.
71
R.J.WheatonandH.R.Zhang,BGInternational,CondensateBankingDynamicsinGasCondensate
Fields: Compositional Changes and Condensate Accumulation Around Production Wells. Paper SPE
62930,preparedforpresentationat the 2000 SPE Annual Technical Conference and Exhibition held in Dallas,
Texas 1-4 October (2000)
72
ZhangH.R.andWheatonR.J.,CondensateBankingDynamicsinGasCondensateFields:Changes
in Produced Condensate to Gas Ratios. Paper SPE 64662, prepared for presentation at the SPE
168
trataron de mostrar como la relacin gas- condensado producido varia con el

tiempo durante el agotamiento de presin. Utilizaron simulacin numrica
con el fin de confirmar y complementar el anlisis terico. El efecto de las
bajas permeabilidades y la heterogeneidad del yacimiento en el banco de
condensado se analizan en trminos de la relacin gas-condensado. Los
resultados han demostrado que no existe el estado de equilibrio natural con
respecto a la relacin gas-condensado.

Hwang
73
, describi la reduccin de la productividad del pozo debido a las
restricciones de flujo cerca al pozo. Las restricciones de flujo, consideran:
dao mecnico, completamiento parcial, alteracin de la formacin debido a
residuos del fluido de completamiento o lquido acumulado y efectos de flujo
no-Darcy. Finalmente, desarrollo una ecuacin que representa los efectos de
las interacciones no lineales entre las restricciones de flujo cerca al pozo en
la productividad del pozo, resultando en un dao total efectivo, que es ms
grande que la suma lineal convencional de los componentes individuales del
daos.

Sarfraz et al.
74
, propusieron un mtodo simple para estabilizar el Inflow
Performance Relationship. El mtodo propuesto utiliza los datos de prueba
del trasiente de presin, para estimar la permeabilidad efectiva como funcin
de la presin y entonces utilizarla para convertir los datos de produccin de
presin de fondo fluyendo en pseudopresin para estabilizar el rendimiento
del pozo. requiere la permeabilidad relativa como funcin de saturacin de
International Oil and Gas Conference and Exhibition in China held in Beijing, China, 710 November
2000
73
Hwang M. K., Modeling Nonlinear Interactions Among NearWell Flow Restrictions in Well
DeliverabilityPredictionandSimulation.PaperSPE65407.
74
Sarfraz A. Joknio and Djebbar, Establishing Inflow Performance Relationship (IPR) for Gas
Condensate Wells. Paper SPE 75503, prepared for presentation at the SPE Gas Technology
SymposiumheldinCalgary,Alberta,Canada,30April2May2002.
169
este modo es eliminada completamente. La permeabilidad efectiva de

cualquier fase puede ser utilizada para predecir la produccin de la segunda
fase. La acumulacin puede ser concebida para estimar la permeabilidad
efectiva utilizando un modelo matemtico. Tambin los modelos matemticos
para las prdidas de productividad en el pozo, debido al depsito de
condensado cuando la presin del punto de roco es alcanzada.

Mott
75
, desarrollo una tcnica para predecir el rendimiento de pozos de gas
condensado, utilizando una simple tcnica, la cual puede ser utilizada en una
hoja de clculo. El clculo utilizado es un modelo de balance de materia para
yacimientos depletados y la integral de pseudopresin de dos fases para el
rendimiento inflow del pozo. la integral de pseudopresin incluye los efectos
de alta velocidad tambin permite los cambios de la composicin del fluido
producido, debido a la formacin del banco de condensado. La hoja de
clculo proporciona una herramienta rpida, para predecir el rendimiento del
pozo, examinando el efecto del banco de condensado en diferentes tipos de
pozos o para estudios de sensibilidad.

Kabir y Hasan
76
, exploraron la posibilidad de utilizar un mtodo simplificado
para calcular la presin de fondo fluyendo, a partir de la presin en cabeza,
tasas medidas, gravedad de fluidos producidos y dimensiones de tubera.
Ellos abogaron por la utilizacin de modelos homogneos de dos
componentes, para eludir problemas con cualquier modelo riguroso de flujo
en dos fases, como delimitaciones de frontera del flujo patrn, estimacin de
desplazamiento entre fases y hacer clculos flash.
75
MOTT, Robert Engineering Calculations of Gas Condensate Well Production. Paper SPE 77551, prepared for
presentation at the Annual Technical Conference and Exhibition held in San Antonio, Texas, 29 September -2
October (2002).
76
Kabir C. S. and Hasan A. R., Simplified WellboreFlow Modeling in Gas/Condensate Systems.
PaperSPE89754.
170

Schroeter y Gringarten
77
estudiaron el crecimiento de la regin de
acumulacin de liquido en un yacimiento de gas condensado, para esto,
desarrollaron un modelo de flujo composicional de dos fases en el cual las
propiedades termodinmicas de las fases son modeladas por la ecuacin de
estado cubica con coeficiente de interaccin binaria cero, gravedad y presin
capilar son insignificantes, la formacin se asume sin lmite lateralmente y
homogneo, la tasa de produccin es constante todo el tiempo. En el modelo
se asume:
Dos fases, multicomponente, isotrmico y
flujo Darcy.
Perfil de permeabilidad relativa
cuadrtico en la saturacin.
La ecuacin de estado cubica con
coeficiente de interaccin binaria cero.
Yacimiento radial, homogneo,
propiedades de la roca isotrpicas.
Tasa de produccin y viscosidades
constantes, y
Gravedad y fuerzas capilares
insignificantes.

Chowdhury et al.
78
, desarrollado un mtodo semianaltico para predecir las
tasas de produccin de gas y aceite (condensado) de pozos con bancos de
77
Shoroeter Thomas von and Gringarten Alain C. Estimates for the Growth of a Condensate Bank
from a Simple Compositional Model. Paper SPE 89911, prepared for presentation at the 2004 SPE
AnnualTechnicalConferenceandExhibitionheldinHouston,Texas,2629September(2004).
78
Nitin, Sharma Ravi, Pope Gary A., and Seperhrnoori Kamy, A Semi-Analytical Method to Predict Well
Deliverability in Gas-Condensate Reservoir Paper SPE 90320, prepared for presentation at the SPE Annual
Technical Conference and Exhibition held in Houston, Texas, U.S.A., 26-29 September (2004).
171
condensados. El mtodo semianaltico incluye los efectos del nmero capilar

(altas velocidades) y flujo no-Darcy. El mtodo fue implementado en un
simulador de yacimientos composicional y verificado los resultados
obtenidos por un simulador composicional de grid fino para fluidos de gas
condensado rico y pobre. El mtodo fue probado para una capa, multicapa, y
mltiples pozos de yacimientos de gas condensando y se encontr que daba
resultados exactos.

Xiao y Al-Muraikhi
79
propusieron un mtodo para determinar el desempeo
de la produccin de un pozo de gas condensado. este mtodo modela el
banco de condensado para considerar el efecto combinado del descenso de
condensado (efecto PVT) y la acumulacin de condensado (comportamiento
de la permeabilidad relativa), esta nueva aproximacin da un mejor
entendimiento de las mecanismos involucrados. Este mtodo se verifica con
la simulacin de datos presentes en la literatura, obtenindose buenos
resultados de acuerdo a los datos publicados para pozo vertical, fracturado y
horizontal.

J amiolahmady et al.
80
, estudiaron el efecto del flujo de gas y condensado en
las regiones perforadas, utilizando modelo aproximado de elemento finitos. El
modelo permite los cambios de las propiedades del fluido y tiene en cuenta
los efectos del acoplamiento positivo y negativo e inercial, utilizando el flujo
fraccional, correlacin base de la permeabilidad relativa.

79
XiaoJ.J.andAlMuraikhiAhmad,ANewMethodfortheDeterminationofGasCondensateWell
Production Performance. Paper SPE 90290, prepared for presentation at the SPE Annual Technical
ConferenceandExhibitionheldinHouston,Texas,2629September(2004).
80
Jamiolahmady M. Danesh A. Sohrabi M. and Atael R., GasCondensate Flow in Perforated
Regions.PaperSPE94072.
172
Bozorgzadeh y Gringarten
81
, utilizaron la simulacin composicional para
identificar los parmetros claves que determinan la entregabilidad del pozo,
demostrando que el numero capilar, la permeabilidad relativa cerca a la cara
del pozo y a la saturacin inicial y la permeabilidad absoluta son los
parmetros que controlan la productividad de un pozo de gas condensado
cuando la presin se encuentra por debajo de la presin del punto de roco.
Finalmente, combinaron la pseudopresin para uno y dos fases, para estimar
los datos de permeabilidad relativa, requeridos en la simulacin.

El flujo multimecanistico es el resultado de la accin combinada de los
componentes del flujo Darcian (flujo macroscpico de la fase causada por
los gradientes de presin) y un componente de flujo difusivo (causado por el
gradiente de concentracin molecular) en un yacimiento de hidrocarburos.
Ayala et al.
82
, desarrollaron una estructura necesaria para la evaluacin del
impacto del flujo multimecanistico en sistemas donde el comportamiento del
fluido es complejo como es el caso del fluido retrogrado de los yacimientos
de gas condensado, para esto se implemento un simulador composicional.
Debido al complejo comportamiento de los fluidos, un modelo totalmente
implcito (IMPISC). Finalmente plantearon diferentes dominios donde el flujo
multimecanistico puede suponerse en un simulador composicional de
retrogrado que se presentado en un yacimiento de gas condensado.

81
Bozorgzabeh, M. and Gringarten, A. C. Estimating Productivity-Controlling Parameters in Gas/Condensate Wells
From Transient Pressure Data. Paper SPE 94018.
82
AyalaLuisF.,ErtekinTurgayandAdewumiMichael.CompositionalModelingofRetrogradeGas
CondensateReservoirinMultimechanisticFlowDomains.PaperSPE94856.
173
Shandrygin y Rudenko
83
propusieron un mtodo para evaluar el dao en
pozos de gas-condensado, utilizando un modelo de pozo simplificado. El
mtodo propuesto se basa en una aproximacin en una fase, tratando el
condensado como inmvil sin fases de transicin entre el gas y el
condensado. Esto reduce el nmero de parmetros de entrada y disminuye
los costos computacionales. En el desarrollo del mtodo se asume un
yacimiento homogneo, simtricamente radial y la fase liquida se considera
inmvil. Se consideran conocidos los parmetros del yacimiento como
porosidad y permeabilidad. El banco de condensado se considera de forma
constante y es caracterizado por la siguiente ecuacin de distribucin
espacial de saturacin de gas:
(178)
Donde, es la saturacin de gas en el yacimiento, es la saturacin de
gas minima, es le tamao del banco de condensado.
La distribucin de la saturacin de gas, utilizada para evaluar el dao de
condensado es calculada por:
(179)
Donde, es la permeabilidad relativa del gas a las condiciones iniciales
del yacimiento y es el radio externo del yacimiento.
El flujo multifasico en la cara del pozo, representa un desafo para los
ingenieros de produccin, desde que la caracterizacin del rgimen de flujo
determine el mtodo apropiado de clculo de la cada de presin. Sadegh et
83
ShandryginA.andRudenkoD.,CondensateSkinEvaluationofGas/CondensateWellsbyPressure
TransientAnalysis.PaperSPE97027,preparedorpresentationatthe2005SPEAnnualTechnical
ConferenceandExhibitionheldinDallas,Texas,U.S.A.,912October(2005).
174
al.
84
, desarrollaron un modelo de flujo, basado en el principio fsico del
sonido, incorporando mtodos hidrodinmicos y termodinmicos actuales. El
mtodo desarrollado. El modelo desarrollado comprende la deteccin de la
rutina del patrn de flujo, un modelo termodinmico y una ecuacin de
estado basada en el comportamiento de fase.

Olaberinjo
85
, considero el comportamiento de flujo en un yacimiento de gas
condensado y desarrollo un modelo que permite calcular aproximadamente
la presin de agotamiento, de un yacimiento de gas condensado, a partir de
la composicin original del fluido del yacimiento, utilizando los clculos de
equilibrio liquido-vapor, las propiedades del fluido, propiedades
adimensionales y la acumulacin del liquido retrogrado como gua.

Ramrez et al.
86
, construyeron un modelo de pozos composicional para
evaluar tanto a los tiempos tempranos como por largo tiempo las
caractersticas de los datos del trasiente de presin, este estudio incluye los
efectos combinados de flujo no-Darcy en presencia de condensacin
retrograda y dao en la cara del pozo en el anlisis del trasiente de presin
en yacimientos de naturalmente fracturados. El objetivo del modelo fue tomar
una mirada crtica del comportamiento de la condensacin retrograda en la
fractura y en la matriz, en la regin cercana a la cara del pozo de un
yacimiento idealizado naturalmente fracturado. La ecuacin de flujo para el
84
Sadegh A. A., Zaghloul J. S. and Adewumi M. A. Gas Condensate Wellbore Modeling Using a
Fundamental Approach. Paper SPE 102479, prepared for presentation at the 2006 SPE Annual
TechnicalConferenceandExhibitionheldinSanAntonio,Texas,U.S.A.,2427September(2006).
85
OlaberinjoA.F.ModelingtheEffectsofCompositionalChangesinTransientFlowofGas
CondensateReservoir.PaperSPE106339,presentedatthe2006SPETechnicalSymposiumofSaudi
ArabiaSectionHeldinDhahran,SaudiArabia,2123May(2006).
86
Ramirez B., Kazemi H. and Ozhan E. NonDarcy Flow Effects in DualPorosity, DualPermeability,
NaturallyFracturedGasCondensateReservoir.PaperSPE109295,preparedforpresentationatthe
2007 SPE Annual Technical Conference and Exhibition held in Anaheim, California, U.S.A., 1114
November(2007).
175
componente hidrocarburo en coordenadas , obtenida del modelo

numrico, es:
(180)
Donde
(181)
Donde es la movilidad del fase en la direccin ( ), es la
densidad molar de la fase ( ), es la permeabilidad relativa de
la fase, y es le gradiente del fluido ( ). Se obtiene una ecuacin
simular para la fase agua.

Los pozos perforados de la direccin de mnima tensin horizontal para
permitir mltiples fracturas transversales, para mejorar la productividad y la
baja permeabilidad en yacimientos de gas y aceite. Los factores que afectan
la productividad en las fracturas de varios pozos horizontales, varan en
funcin del yacimiento y las propiedades de la fractura, as como las
trayectorias pozo a pozo. Guo y Yu
87
, formularon un modelo analtico simple
que permite describir mejor la productividad de la fractura de varios pozos
horizontales, el modelo asocia el flujo radial con la regin no fracturada del
yacimiento, el flujo lineal hacia la fractura en la regin fracturada, el flujo
lineal en la fractura y el flujo radial en la fractura hacia el pozo horizontal.
87
GuoBoyunandYuXiance,ASimpleandAccurateMathematicalModelforPredictingProductivity
of Multifractured Horizontal Wells. Paper SPE 114452, prepared for presentation at the CIPS/SPE
Gas Technology Symposium 2008 Joint Conference held in Calgary, Alberta, Canada, 1619 June
(2008).
176

Mahdiyar et al.
88
, desarrollaron un simulador matemtico de dos
dimensiones, basados en el mtodo de diferencias finitas y teniendo en
cuenta los efectos combinados de acoplamiento e inercial utilizando
correlaciones generalizadas desarrolladas recientemente. Este modelo de
pozo, el cual simula el flujo en estado estable alrededor de un pozo
fracturado hidrulicamente, tambin tiene en cuenta el cambio de fase y los
sistemas de baja tensin interfacial.

Singh y Whitson
89
, verificaron el mtodo de pseudopresin de gas
condensado de Fevang Whitson para yacimientos estratificados. Las capas
se caracterizan por una amplia variacin de permeabilidad y composicin, y
pueden estar o no comunicadas. El mtodo de pseudopresin es utilizado en
los clculos para modelos de grid, con red areal relativamente larga,
capturando el banco cerca al pozo con un grid refinado local, de ese modo,
reducen el tiempo de corrida y el tamao del modelo. Los resultados
obtenidos verificaron que en los primeros tiempos, el mtodo de
pseudopresiones propuesto por Fevang-Whitson es vlido y exacto para
sistemas de capas con significante heterogeneidad, con y sin flujo, con y sin
nmero capilar modificando las permeabilidades relativas y con mucho
lquido para la composicin de cada capa.
7.1 METODOLOGA APLICADA AL CAMPO

88
Mahdiyar H., Jamiolahmady M. and Sohrabi M. A New Flow Skin Factor Formulation for
Hydraulically Fractured Wells in Gas Condensate Reservoir. Paper SPE 113597, prepared for
presentationatthe2008SPEEurope/EAGEAnnualConferenceandExhibitionheldinRome,Italy,9
12June(2008)
89
SinghKameshwarandWhitsonCurtisH.,GasCondensatePseudopressureinLayeredReservoirs.
Paper SPE 117930, prepared for presentation at the 2008 Abu Dhabi International Petroleum
ExhibitionandConferenceheldinAbuDhabi,UAE,36November(2008).
177
7.1.1 Generalidades del campo Cupiagua

El Campo Cupiagua fue descubierto en 1993, est localizado en la cuenca
de los Llanos como se observa en la Fig.27, la cual es una cuenca
asimtrica, localizada al Este de los Andes y hace parte de las cuencas
subandinas del pas (ECOPETROL, 1995). La cuenca de los Llanos
Orientales cubre un rea aproximada de y geogrficamente se
divide en tres reas principales; una zona norte, Arauca, una parte central,
Casanare y un sector sur, Meta. La cuenca limita al Oeste por el Piedemonte
de la cordillera Oriental, al Este por el Escudo Precmbrico de Guyana y al
Sur por la serrana de La Macarena y el arco de basamento.

La cuenca de los llanos tiene continuidad hacia el Norte en la cuencas
Venezolanas de Apure Occidental y Barinas, las cuales estn separadas de
la cuenca de los Llanos Colombianos, desde el Terciario temprano, por el
alto estructural del Arco de Arauca. Alguna caractersticas generales del
campo se resumen en la Tabla 5.

7.1.2 Aplicacin

El banco de condensado cerca al pozo, es el principal factor que causa la
perdida de entregabilidad del mismo, el problema se complica en pozos de
baja permeabilidad, el campo Cupiagua (Fig. 30), es un yacimiento de gas
condensado, que se caracteriza por baja permeabilidad y alta temperatura.
178
Estudios previos Lee y Chaverra

90
han demostrado que el sistema in-situ
esta cercano a las condiciones termodinmicas criticas y cuando la presin
cae por debajo del punto de roco, una significativa cantidad de liquido se
condensa.
Figura 27: Localizacin del rea del Campo Cupiagua.

Fuente: Well Review BP, 2001

Tabla 5: Propiedades de Campo Cupiagua.
90
LeeShengTaiandChaverraMarcial,ModellingandInterpretationofCondensateBankingforthe
NearCriticalCupiaguaField.PaperSPE49265,preparedforpresentationatthe1998SPEAnnual
TechnicalConferenceandExhibitionheldinNewOrleans,Louisiana,2730September(1998).
179
CARACTERSTICAS GENERALES CAMPO CUPIAGUA

Fraccin de Intervalos Mirador, Barco
Tipo de roca Cuarzoarenitas
Porosidad 5-6 % en areniscas puras
Permeabilidad 0.5-10 md
Saturacin de agua inicial 10%
Presin 6025 psi
Temperatura 240-285 F
Fluido del yacimiento Gas condensado, 3500 GOR
Gradiente de esfuerzo mximo horizontal 1.2 -1.7 psi/pies
Gradiente de esfuerzo vertical 1.07 psi/pies
Gradiente de esfuerzo mnimo horizontal 0.65 0.77 psi/pies
Gradiente de presin de poro 0.44 psi/pies
Formacin del yacimiento Cuarzoarenita, rea tectnicamente
activa
Principales problemas de estabilidad
presentados
Ensanchamiento del hueco, derrumbes,
pega de tubera.
Estructura del yacimiento Anticlinal fallado
Formacin de mayor inestabilidad Carbonera, Mirador y Cuervos
Fuente: MNDEZ Y. y PABON E. Tesis de Grado 2004.

All et al.
91
, a partir de sus experiencias, definieron que un modelo efectivo de
un yacimiento de gas condensado, requiere los siguientes parmetros o
fenmenos:
1. La saturacin critica del condensado y su
dependencia de la tensin interfacial y la velocidad de flujo cerca a la
cara del pozo.
91
AllJ.K.,McGauleyP.J.andWilsonC.J.ExperimentalStudiesandModellingofGasCondensate
FlowNeartheWellbore.PaperSPE39053,preparedforpresentationattheFifthLatinAmericanand
CaribbeanPetroleumEngineeringConferenceandExhibitionheldinRiodeJaneiro,Brazil,30August
3September(1997)
180
2. La reduccin de la permeabilidad relativa

en la presencia de la condensacin inmvil.
3. La aparicin de la velocidad de
extraccin.
4. Flujo no-Darcy.

En modelo radial para estudiar yacimientos de baja permeabilidad, la
ecuacin de flujo es obtenida de la combinacin de la ecuacin de
continuidad en direccin radial, el flujo turbulento, ecuacin de Darcy y
ecuacin de estado para los gases reales, (Peng-Robinso o Redlich-
Kwong)
92
.
La ecuacin de continuidad, se obtiene mediante un balance de masa de
fluido. Para un sistema de flujo de gas condensado en estado estable es
93
:
(182)
Las ecuaciones de flujo para la fase de gas y condensado son:

(183)
y son la permeabilidad relativa al gas y al condensado
respectivamente, afectada por la inercia y el acoplamiento. Se utilizo la
correlacin generalizada desarrollada por J amiolahamdy et al.
94

92
KatzDonald,CamposJ.AntonioandSinghS.P.ApplicationofRadialSimulationModel.PaperSPE
5883.
93
Mahdiyar H., J amiolahmady M. and Sohrabi M. A New Flow Skin Factor Formulation for Hydraulically Fractured
Wells in Gas Condensate Reservoir. Paper SPE 113595 prepared for presentation at the 2008 SPE Europe/EAGE
Annual Conference and Exhibition held in Rome, Italy, 9-12 J une 2008.
181
Combinando las ecuaciones anteriores se obtiene:

(184)
Se debe asumir que podra haber transferencia de masa entre las dos fases
que fluyen, pero la composicin total del fluido ( ), permanece constante en
el fluido que vieja a travs del medio poroso
95
. Es decir, no hay ninguna
fuente o sumidero en la masa total de cada componente. Por lo tanto,

(185)
Donde, , y son las fracciones de masa del componente en la mezcla,
en el gas y el condensado, respectivamente; es la tasa de flujo
volumtrico; y es el flujo fraccional del gas.

La entregabilidad un parmetro muy importante en el comportamiento de la
productividad del pozo, puede ser calculada a partir de una modificacin de
la ecuacin de Evinger-Muskat de pseudopresin para gas seco,
desarrollada por Fevang y Whitson
96
, adems se puede dividir del flujo en
tres regiones principales:
94
Jamiolahmady M., Danesh A., Tehnari D. and Sohrabi M. Variations of Gas/Condensate Relative
Permeability With Production Rate at NearWellbore: A General Correlation SPE 83960, SPE
ReservoirEvaluation&EngineeringJournal,Volume9,Number5,pp688697,(Dec2006).
95
JamiolahmadyM.,DaneshA.andDuncanD.B.MeasurementandModelingofGasCondensate
FlowAroundaRockPerforationTransportinPorousMediaVolume63Number2,pp323347,
(2006).
96
FevangOivingandWhitsonC.H.ModelingGasCondensateWellDeliverability.PaperSPE30714,
preparedforpresentationat the 2000 SPE Annual Technical Conference and Exhibition held in Dallas, U.S.A.,
22-25 October (1995)
182
gas, en la cual ambos estn fluyendo simultneamente. En esta
regin, el flujo composicional se mantiene constante y es la principal
fuente de prdida de entregabilidad en los pozos de gas condensado,
estas prdidas depende principalmente de la permeabilidad relativa al
gas en la regin y de su tamao. La permeabilidad relativa al gas en la
Regin 1, es funcin principalmente de la distribucin de la saturacin
de lquido, esta distribucin es determinada especficamente para
asegurar que todo el lquido condensado que entra, tenga la suficiente
movilidad para fluir a travs y fuera de ella sin ninguna acumulacin.
condensado esta inmvil y solo el gas esta fluyendo. La acumulacin
del condensado es causada por dos mecanismos: 1. El condensado
producido por el yacimiento de gas, debido a la declinacin de la
matriz del yacimiento en la corriente del yacimiento de gas en esta
regin.
Regin 3: una regin que contiene solo una fase de gas. Esta regin
siempre existe en un yacimiento de gas condesado que se encuentra
subsaturado y la composicin es constante e igual a la original del
yacimiento.

Para el clculo de la productividad, se presentas diferentes enfoques como:
Las tasa de produccin de gas y aceite (condensado) de pozos de gas
condensado, pueden predecirse a partir del modelo semianaltico
desarrollado por Chowdhury et al.
97

97
Nitin, Sharma Ravi, Pope Gary A., and Seperhrnoori Kamy, A Semi-Analytical Method to Predict Well
Deliverability in Gas-Condensate Reservoir Paper SPE 90320, prepared for presentation at the SPE Annual
Technical Conference and Exhibition held in Houston, Texas, U.S.A., 26-29 September (2004).
183
A partir del modelo de pozo composicional propuesto por Afidick et al.

98
, el
efecto de la acumulacin del condensado puede ser estimado de la siguiente
ecuacin:
(186)
Donde, es la presin promedio del yacimiento, es la presin de fondo
fluyendo en , es el factor dao a , es el radio de la zona
interna, , es el radio de drenaje, es le radio del pozo, es la tasa
de flujo en , y es la permeabilidad
efectiva del gas a la saturacin critica del condensado.

Fusell
99
evala el efecto de la acumulacin del condensado en la vecindad
del pozo, mediante la relacin de productividad, definida como la relacin del
ndice de productividad de un pozo actual con el ndice de productividad de
un pozo estndar, puede utilizarse para estimar la prdida de productividad.
Se utiliza un pozo estndar de caractersticas idnticas al pozo actual. La
relacin de productividad se define como:
(187)
Donde la evaluacin de y , debe ser dado para valores de y .
La relacin de productividad en estado estable se define como:
(188)
98
AfidickDeddy,KaczorowskinN.J.andBetteSrinivas.ProductionPerformaceofaRetrogradeGas
Reservoir:ACaseStudyoftheArunField.PaperSPE28749,preparedforpresentationattheSPE
AsiaPacificOil&GasConferenceheldinMelbourne,Australia,710November(1994).
99
FussellD.D.SingleWellPerformancePredictionsforGasCondensateReservoir.PaperSPE4072.
184
Donde es la tasa de produccin en estado estable. La magnitud de la

relacin de productividad indica las severas prdidas de productividad debido
a la acumulacin de lquido. Desarrollo un mtodo de productividad del pozo
que depende de los mtodos de deplecin y estado estable definidos a
continuacin:
Mtodo de deplecin: Las prdidas de la productividad sostenida por la
produccin de un pozo de gas condensado, como resultado de la
acumulacin del condensado, se puede determinar comparando el
desempeo actual predicho de un pozo con el desempeo de un pozo
estndar.
La funcin pseudopresin para el mtodo de deplecin se define como:
(189)
La funcin , asume que la condensacin de la fase condensada es
inmvil. La integral de la ecuacin anterior solo es funcin de la presin,
dado que la saturacin del condensado es una funcin conocida de la
presin.
Mtodo estado estable: Este mtodo es similar al presentado por ODell y
Miller
100
, la ecuacin de estado estable desarrollado para un gas
condensado es similar a la funcin de pseudopresin para gas seco.
Entonces, la funcin de pseudopresin para el gas condensado se define
como:
(190)
100
ODell,H.G.andMiller,R.N.SuccessfullyCyclingaLowPermeability,HighYieldGasCondensate
Reservoir.JPet.Tech.(Jan.,1967)4144
185
El cual incluye el factor de compresibilidad para cada fase. La integral

anterior es solo funcin de la presin.

Mtodos de la productividad del pozo, deplecin y estado estable: Las
funciones pseudopresin de deplecin y de estado estable, son usadas en
la ecuacin
(191)
Para calcular la tasa de produccin. La ecuacin anterior requiere la
presin en el exterior y en la cara de la arena del yacimiento. Esas
presiones son utilizadas, para determinar la apropiada funcin
pseudopresin de los mtodos de deplecin y de estado estable. Las
funciones de pseudopresin son utilizadas en la ecuacin anterior para
calcular la productividad del pozo, y . La tasa de produccin es la
tasa de produccin del pozo estndar.

Debido a que el campo Cupiagua es un yacimiento naturalmente fracturado
se debe plantear un modelo de porosidad que describa el cambio de las
porosidades de matriz y fractura ante los esfuerzos efectivos, Osorio et al.
101
,
proponen las siguientes expresiones:
Ecuacin de porosidad para la matriz
(192)
Ecuacin de porosidad para la fractura
101
OSORIO GALLEGO, Jose Gildardo. Notas para un curso de aplicaciones geomecanicas a la Ing. de
Yacimiento. Universidad Nacional de Colombia. Facultad de Minas. Departamento de Recursos
Minerales.Medelln1996.
186
(193)
Donde, es el esfuerzo efectivo asociado con el cambio del volumen
poroso en un sistema de doble porosidad y es el esfuerzo efectivo
asociado con el cambio del volumen total en un sistema de doble porosidad,
es la compresibilidad de la roca, es la compresibilidad de los poros
primarios y es la compresibilidad de los poros secundarios. Adicional al
modelo de porosidad, debe incluirse un modelo de deformacin
geomecanica, el cual describe el comportamiento de la deformacin de la
parte solida de la roca ante un cambio de esfuerzos del yacimiento generado
por la perforacin del pozo. El modelo considera deformacin infinitesimal,
elstica no lineal para el esqueleto de la roca que conforma el yacimiento. El
modelo de deformacin geomecanico propuesto por Osorio et al., est
conformado por las ecuaciones de equilibrio, deformacin-desplazamiento y
esfuerzo-deformacin-presin.

Como posiblemente, el sistema de ecuaciones resultantes es no lineal. Por lo
tanto, se hace necesario aplicar un mtodo numrico para la solucin. El mas
recomendado es el mtodo de diferencias finitas de segundo orden. El
tamao del grid para este sistema no es uniforme, se puede determinar con
una progresin logartmica, hacia el limite radial del yacimiento.
187
8. CONCLUSIONES

A partir de la investigacin realizada, se puede concluir:

Cuando la presin del yacimiento cae por debajo de la presin del punto
de roco, una regin altamente saturada de lquido se acumula alrededor
del pozo, alterando el flujo de gas y reduciendo la productividad. Esta
regin saturada de liquido, recibe el nombre de Banco de condensado,
sus efectos depende de: (1) permeabilidades relativas, (2) propiedades
y (3) la forma en que el pozo es producido.

La saturacin de condensado alrededor de un pozo productor de gas
condensado, no solo depende de la presin, sino que tambin se debe
tener en cuenta la concentracin total de los componentes pesados y
como varan estos con el tiempo cerca del pozo.

En un pozo productor de gas condensado, existen tres regiones
principales de flujo: la Regin 1 cerca a la cara del pozo, saturada con
aceite y gas, en la cual ambos fluyen simultneamente, esta es la
principal fuente de prdida de entregabilidad en el pozo; la Regin 2 o
regin de acumulacin, en la cual, el liquido condensado es inmvil y
solo fluye gas y la Regin 3 que contiene solo gas, cuya composicin es
igual a la original del yacimiento.

En los yacimientos de gas, las propiedades cambian con la presin, para
evitar este problema, se introduce el concepto de pseudopresion, Fevang
188
y Whitson en su modelo de entregabilidad, calcularon la pseudopresion

para las tres regiones de flujo que se pueden presentar en el yacimiento
de gas condensado.

La productividad en la mayora de pozos de gas condensado se reduce
por la formacin de un banco de condensado en la regin cercana al
pozo. Este banco de condensado, disminuye la entregabilidad de los
pozos de gas condensado, para modelar este problema se requiere de
un grid fino de simulacin, esto usualmente representa mayor tiempo de
computo.

Para modelar un yacimiento tan complejo como Cupiagua, se debe tener
en cuanta diferentes factores como: los efectos del banco de
condensado, el flujo no-Darcy, la doble porosidad, la deformacin
mecnica, entre otros.

Los modelos de pozos, son la mejor opcin para modelar yacimientos
complejos como es el caso del campo Cupiagua, ya que modelar todo el
yacimiento, requiere buenas tcnicas de refinamiento, mayor tiempo de
computo y probablemente los resultados obtenidos no correspondan al
comportamiento del campo.

189

RECOMENDACIONES

En un prximo trabajo sobre este tema se deben tener en cuenta las
siguientes sugerencias.

Se recomienda modelar el yacimiento, desde lo convencional
hasta lo complejo, modificando uno a uno los parmetros, con el
fin de analizar, como vara la produccin de un modelo a otro.

Si se tiene la informacin necesaria, se recomienda desarrollar un
modelo computacional, que tenga en cuenta, la saturacin critica
del condensado y su dependencia de la tensin interfacial y la
velocidad del flujo; la reduccin de la permeabilidad relativa en
presencia de condensacin inmvil, el flujo no-Darcy y las
modificaciones de las permeabilidades relativas a altas
velocidades y/o baja tensin interfacial.

Con base en la metodologa propuesta, y en investigaciones
previas, desarrollar un modelo numrico que permita predecir el
comportamiento de la produccin de un pozo de gas condensado.

Si el modelo computacional no se puede desarrollar, pero se logra
obtener el modelo composional, este puede ser modelado en
paquetes comerciales como Eclipse, VIP de Landmark, entre
otros.

190

BIBLIOGRAFA

AFIDICK Deddy, KACZOROWSKIN N. J . and BETTE Srinivas. Production
Performace of a Retrograde Gas Reservoir: A Case Study of the Arun Field.
Paper SPE 28749, prepared for presentation at the SPE Asia Pacific Oil &
Gas Conference held in Melbourne, Australia, 7-10 November (1994).

ALL J . K., MCGAULEY P. J . and WILSON C. J . Experimental Studies and
Modelling of Gas Condensate Flow Near the Wellbore. Paper SPE 39053,
prepared for presentation at the Fifth Latin American and Caribbean
Petroleum Engineering Conference and Exhibition held in Rio de J aneiro,
Brazil, 30 August 3 September (1997)

AL-MAJ ED A. A., FAHD King and DOUGHERTY E. L. A Variable Cell Model
for Simulation Gas Condensate Reservoir Performance Paper SPE 21428,
prepared for presentation at the SPE Middle East Oil Show held in Bahrain,
16-19 November (1991).

AYALA Luis F., ERTEKIN Turgay and ADEWUMI Michael. Compositional
Modeling of Retrograde Gas-Condensate Reservoir in Multimechanistic Flow
Domains. Paper SPE 94856

191
BLOM Saskia M. P. and HAGOORT J acques, The Combined Effect of Near-

Critical Relative Permeability and Non-Darcy Flow on Well Impairment by
Condensate Drop Out. Paper SPE 51367.

BOZORGZABEH, M. and GRINGARTEN, A. C. Estimating Productivity-
Controlling Parameters in Gas/Condensate Wells From Transient Pressure
Data. Paper SPE 94018.

CLARK T. J . The Application of 2D Compositional, Radial Model To Predict
Single-Well Performance in a Rich Gas Condensate Reservoir. Paper SPE
14413, prepared for presentation at the 60
th
Annual Technical Conference
and Exhibition of the Society of Petroleum Engineers had in Las Vegas, NV
September 22-25, (1985)

COATS Keith H. Simulation of Gas Condensate Reservoir Performance.
Paper SPE 10512.

COATS Keith H. An Equation of State Compositional Model. Paper SPE
8284.

COATS K. H., Implicit Compositional Simulation of Single-Porosity and
Dual-Porosity Reservoir. Paper SPE 18427, prepared for presentation at the
SPE Symposium on Reservoir Simulation in Houston, TX, February 6-8
(1989).

192
CHOWDHURY Nitin, SHARMA Ravi, POPE Gary A., and SEPERHRNOORI

Kamy, A Semi-Analytical Method to Predict Well Deliverability in Gas-
Condensate Reservoir SPE 90320. Copyright 2004.

DUMOR J . M., HAGOORT J . and RISSEEUW A. S. An Analytical Model for
One-Dimensional, Three-Component Condensing and Vaporizing Gas
Drives. Paper SPE 10069.

EL-BANBI Ahmed H., MCCAIN W. D. and SEMMELBECK M. E.
Investigation of Well Productivity in Gas-Condensate Reservoir. Paper SPE
59773, prepared for presentation at the 2000 SPE/CERI Gas Technology
Symposium held in Calgary, Alberta Canada, 3-5 April 2000.

FEVANG Oivind Gas Condensate Flow Behavior and Sampling. A
Dissertation for the Partial Fulfillment of Requirements for the Degree of
Doktor Ingenior. Division of Petroleum Engineering and Applied Geophysics
The Norwegian Institute of Technology University of Trondheim. October
(1995)

FEVANG Oiving and WHITSON C. H. Modeling Gas-Condensate Well
Deliverability. Paper SPE 30714, prepared for presentation at the 2000 SPE
Annual Technical Conference and Exhibition held in Dallas, U.S.A., 22-25
October (1995)

FUSSELL D. D. Single-Well Performance Predictions for Gas Condensate
Reservoir. Paper SPE 4072.
193
GOMEZ Hamilton P. and CORREA Antonio Claudio Fully Implicit

Compositional Modelling of Gas Condensate and Volatile Oil Reservoir.
Paper SPE 23700 prepared for presentation at the 2
nd
Latin American
Petroleum Engineering Conference held in Caracas, Venezuela, March 8-11,
(1992).

HWANG, M. K. Modeling Nonlinear Interaction Among Near-Well Flow
Restrictions in Well-Deliverability Prediction and Simulation. Paper SPE
65407.

J AMIOLAHMADY, M., DANESH, A., SOHRABI, M. and ATAEI, R. Gas-
Condensate Flow in Perforated Regions. Paper SPE 94072.

KABIR C.S. and HASAN A. R. Simplified Wellbore-Flow Modeling
Gas/Condensate Systems. Paper SPE 89754

KATZ Donald, CAMPOS J . Antonio and SINGH S. P. Application of Radial
Simulation Model. Paper SPE 5883.

LEE Sheng-Tai and CHAVERRA Marcial, Modelling and Interpretation of
Condensate Banking for the Near Critical Cupiagua Field. Paper SPE 49265,
prepared for presentation at the 1998 SPE Annual Technical Conference and
Exhibition held in New Orleans, Louisiana, 27-30 September (1998).

194
LI, Kewen and FIROOZABADI, Abbas Phenomenological Modeling of

Critical Condensate Saturation and Relative Permeabilities in
Gas/Condensate Systems. Paper SPE 56014

MAHDIYAR H., J AMIOLAHMADY M. and SOHRABI M. A New Flow Skin
Factor Formulation for Hydraulically Fractured Wells in Gas Condensate
Reservoir. Paper SPE 113595 prepared for presentation at the 2008 SPE
Europe/EAGE Annual Conference and Exhibition held in Rome, Italy, 9-12
J une 2008.

MOTT, Robert Engineering Calculations of Gas Condensate Well
Production. Paper SPE 77551

NITIN, Sharma Ravi, POPE Gary A., and SEPERHRNOORI Kamy, A Semi-
Analytical Method to Predict Well Deliverability in Gas-Condensate Reservoir
Paper SPE 90320, prepared for presentation at the SPE Annual Technical
Conference and Exhibition held in Houston, Texas, U.S.A., 26-29 September
(2004).

OLABERINJ O A. F. Modeling the Effects of Compositional Changes in
Transient Flow of Gas Condensate Reservoir. Paper SPE 106339,
presented at the 2006 SPE Technical Symposium of Saudi Arabia Section
Held in Dhahran, Saudi Arabia, 21-23 May (2006).

POPE G. A., WU W., NARAYANASWAMY G., DELSHAD M., SHARMA M.
and WANG P. Modeling Relative Permeability Effects in Gas-Condensate
195
Reservoirs. Paper 49266, prepared for presentation at the 1996 Annual

Technical Conference and Exhibition held in New Orleans, Louisiana, 27-30
September 1996.

RAMIREZ B., KAZEMI H. and OZHAN E. Non-Darcy Flow Effects in Dual-
Porosity, Dual-Permeability, Naturally Fractured Gas Condensate Reservoir.
Paper SPE 109295, prepared for presentation at the 2007 SPE Annual
Technical Conference and Exhibition held in Anaheim, California, U.S.A., 11-
14 November (2007).

ROEBUCK I. F., FORD W. T. HENDERSON G. E. and DOUGLAS J im. The
Compositional Reservoir Simulator: Case III The Radial Geometry. Paper
SPE 2486.

ROJ AS, Gonzalo Ingeniera de Yacimientos de Gas Condensado

SADEGH A. A., ZAGHLOUL J . S. and ADEWUMI M. A. Gas Condensate
Wellbore Modeling Using a Fundamental Approach. Paper SPE 102479,
prepared for presentation at the 2006 SPE Annual Technical Conference and
Exhibition held in San Antonio, Texas, U.S.A., 24-27 September (2006).

SARFRAZ A. J OKNIO and DJ EBBAR, Establishing Inflow Performance
Relationship (IPR) for Gas Condensate Wells. Paper SPE 75503, prepared
for presentation at the SPE Gas Technology Symposium held in Calgary,
Alberta, Canada, 30 April- 2 May 2002.

196
SINGH Kameshwar and WHITSON Curtis H., Gas Condensate

Pseudopressure in Layered Reservoirs. Paper SPE 117930, prepared for
presentation at the 2008 Abu Dhabi International Petroleum Exhibition and
Conference held in Abu Dhabi, UAE, 3-6 November (2008).

SHOROETER Thomas von and GRINGARTEN Alain C. Estimates for the
Growth of a Condensate Bank from a Simple Compositional Model. Paper
SPE 89911, prepared for presentation at the 2004 SPE Annual Technical
Conference and Exhibition held in Houston, Texas, 26-29 September (2004).

XIAO J . J . and AL-MURAIKHI Ahmad, A New Method for the Determination
of Gas Condensate Well Production Performance. Paper SPE 90290,
prepared for presentation at the SPE Annual Technical Conference and
Exhibition held in Houston, Texas, 26-29 September (2004).

ZHANG H. R. and WHEATON R. J ., Condensate Banking Dynamics in Gas
Condensate Fields: Changes in Produced Condensate to Gas Ratios. Paper
SPE 64662, prepared for presentation at the SPE International Oil and Gas
Conference and Exhibition in China held in Beijing, China, 7-10 November
2000

130363

Cargado por

Copyright:

Formatos disponibles

130363

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

130363

Cargado por

Copyright:

Formatos disponibles

METODOLOGA PARA OBTENER UN MODELO DE POZO QUE PERMITA

PREDECIR EL COMPORTAMIENTO DE LA PRODUCCIN DE UN POZO

METODOLOGA PARA OBTENER UN MODELO DE POZO QUE PERMITA

3.2 MTODO SEMIANALITICO PARA PREDECIR LA

4.3.2 Presentacin de los resultados en trminos del ndice de

6.3.1.1 Representacin del pozo .103

6.6.2.2 Datos Govier-Fogarasi .140

Figura 11: Relacin gas-aceite fluyendo en un pozo vertical sin fracturas

Figura 25: Esquema de un pozo fracturado hidrulicamente .122

El objetivo principal de este estudio es desarrollar una metodologa que

1.1.2 Yacimientos de gas hmedo

menos cantidad de componentes pesados y (3) la cantidad de lquido

generara liquido al caer la presin por debajo de la presin del punto de

Figura 2: Diagrama presin-temperatura para un yacimiento de gas

Figura 3: Diagrama presin-temperatura para un yacimiento de gas

60%, generando reduccin en la productividad de los pozos de 2 a 10

Figura 4: Formacin del bloque de condensado.

sino que tambin se debe tener en cuenta la concentracin total de los

La movilidad est definida en la Ec. (2), en trminos de la funcin Chi se

el cual no necesita desarrollar una compleja simulacin composicional

Figura 6: Perfil de presin de un yacimiento homogneo, que se encuentra

cada fase. Por consiguiente es claro que es un parmetro crtico en el

3. ENTREGABILIDAD DE POZOS DE GAS CONDENSADO

iniciar el desarrollo de esta investigacin, Fevang se baso en investigaciones

simulaciones de aceite negro, ellos extendieron el trabajo de Kniazeff y

condensado, comparado con resultados de simulaciones. El tambin evala

yacimiento (usualmente en yacimientos de baja permeabilidad). La tasa de

3.1.1.1 Regin cercana al pozo: El flujo composicional (GOR) dentro de la

3.1.1.3 Regin de solo gas: La Regin 3 siempre existe en un yacimiento

3.1.3 Calculo de pseudopresin

punto de roco del wellstream composicional producido, se define como . Si

El volumen de aceite relativo , puede ser expresado en trminos de las

Donde el subndice representa la presin en el estado actual y

funcin de pseudopresin. La geometra del pozo, el factor dao, entre otros,

Usando la definicin anterior y la Ley de Darcy se llega a:

mtodo semianaltico. La regin entre la cara del pozo y el radio equivalente

densidad y la viscosidad es significativa y necesita ser tenidas en

La productividad del pozo es un punto importante en el desarrollo de

simulador de aceite ajustado y uno composicional, para yacimientos de gas

4.1.1 Calculo de la razn de aceite fluyendo ( ) en la regin 1.

alrededor del pozo con permeabilidad diferente. Si el radio de esta zona es

4.1.2.2 Pozo vertical fracturado: Prats

materia, ignoraran los efectos de alta velocidad. Se realizaron dos

temporal y as obtener el gradiente de presin lineal y el perfil de saturacin

Fuente: Bozorgzabeh, M. and Gringarten, A. C. Estimating Productivity-Controlling Parameters in Gas/Condensate

permeabilidad efectiva al gas cerca a la cara del pozo a la saturacin de

adems de la relacin entre y , segn Fevang y Whitson

La simulacin composicional con y sin efectos del nmero capilar se utilizo

El uso de la pseudopresin de una sola fase en yacimientos de gas

, , , y el dao en las derivadas de pseudopresin de dos fases de

Basados en los resultados anteriores de la simulacin, se propone el

5. Seleccione un valor para .

corridas y problemas de estabilidad. La eficiencia computacional de las

que se supone que es la misma saturacin de agua connata (es decir

juega un papel significativo y en el momento de simular, utilizar las capas fina

en procesos de simulacin. La clave de un modelo de permeabilidad relativa

del ajuste de los datos de saturacin residual de la fase . generalmente es

Donde, es la permeabilidad relativa y es la permeabilidad relativa en el

puntos finales del gas, se muestran en la Fig. 18 y es funcin de la presin

condensado, la permeabilidad relativa del gas y la permeabilidad