0% encontró este documento útil (0 votos)

25 vistas170 páginas

Estudio de Conificación de Agua en Yacimientos

Este documento presenta una tesis para optar por el grado de Maestro en Ingeniería en Ingeniería Petrolera y Gas Natural Yacimientos. La tesis estudia el fenómeno de la conificación de agua en yacimientos naturalmente fracturados. Incluye una revisión bibliográfica del tema, la caracterización de mallas de simulación, un análisis del efecto de varios parámetros en la conificación de agua en pozos bien y mal cementados, y la aplicación de las correlaciones obtenidas a casos prá

Cargado por

Liliana Rodriguez

Derechos de autor

Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.

Formatos disponibles

Descarga como PDF, TXT o lee en línea desde Scribd

0% encontró este documento útil (0 votos)

25 vistas170 páginas

Estudio de Conificación de Agua en Yacimientos

Cargado por

Liliana Rodriguez

Derechos de autor

Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.

Formatos disponibles

Descarga como PDF, TXT o lee en línea desde Scribd

UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO

PROGRAMA DE MAESTRÍA Y DOCTORADO EN

INGENIERÍA

UNIVERSIDAD NACIONAL
AUTONOMA DE FACULTAD DE INGENIERÍA
MEXICO

“ESTUDIO DE CONIFICACIÓN DE AGUA EN

YACIMIENTOS NATURALMENTE
FRACTURADOS”

TESIS
QUE PARA OPTAR POR EL GRADO DE:

MAESTRO EN INGENIERÍA

INGENIERÍA PETROLERA Y GAS NATURAL

YACIMIENTOS

P R E S E N TA:
ERNESTO PÉREZ MARTÍNEZ

TUTOR:
Dr. RODRÍGUEZ DE LA GARZA FERNANDO
Dr. SAMANIEGO VERDUZCO FERNANDO

21 de febrero de 2011
Página 2 de 170
Agradecimientos

Agradezco a Dios por estar conmigo en cada paso que doy, por fortalecer mi corazón,
por iluminar mi mente y por haber puesto en mi camino a aquellas personas que han
sido soporte y compañía durante toda mi vida.

Agradezco a mis padres y hermanos porque a pesar de no estar presentes físicamente

siempre procuran mi bienestar.

Agradezco a mi esposa por su apoyo y compañía incondicional.

Agradezco al Dr. Fernando Rodríguez de la Garza, Dr. Fernando Samaniego Verduzco,

Dr. Guillermo Cruz Domínguez Vargas, Dr. Edgar René Rangel Germán y al Dr. Víctor
Hugo Arana Ortiz; por compartir su tiempo y conocimientos pero sobre todo por su
amistad, en especial al Dr. Fernando Rodríguez De la Garza, por su acertada dirección
en el presente trabajo.

Agradezco a Petróleos Mexicanos por las facilidades y apoyo que me brindó.

Agradezco a la UNAM por todas las experiencias, conocimientos y amigos adquiridos,

durante mis maravillosos años de estudio.

Agradezco a todas y cada una de las personas que han vivido conmigo la realización
de esta tesis, con sus altos y bajos y que no necesito nombrar porque tanto ellas como
yo sabemos que desde lo más profundo de mi corazón les reconozco el haberme
brindado todo el apoyo, colaboración, ánimo y sobre todo cariño y amistad.

Página 3 de 170
Contenido

Agradecimientos. .............................................................................................................................. 3

Contenido. ........................................................................................................................................ 4

Lista de figuras. ................................................................................................................................. 6

Índice de tablas. .............................................................................................................................. 12

Resumen. ........................................................................................................................................ 14

Introducción. .................................................................................................................................. 15

Capítulo I. ....................................................................................................................................... 18

Revisión Bibliográfica de Conificación de Agua. ...................................................................................... 19

Capítulo 2. ...................................................................................................................................... 57

Características de las Mallas de Simulación. ........................................................................................... 58

Capítulo 3. ...................................................................................................................................... 70

Conificación de agua en pozos bien cementados en el espacio anular (tubería de revestimiento‐
formación productora).     ........................................................................................................................... 71
Producción de aceite (qo).     ........................................................................................................................... 74
Permeabilidad de la fractura (kf).    ................................................................................................................ 78
Análisis del efecto de la distancia entre el intervalo productor y el Cwo.     .................................................. 80
Viscosidad del aceite.     .................................................................................................................................. 83
Efecto de la partición de porosidad total en primaria y secundaria.     .......................................................... 87
Densidad del aceite.     .................................................................................................................................... 93
Radio de drene.     ........................................................................................................................................... 94
Obtención de la correlación para determinar la altura máxima de la conificación de agua.     ..................... 95

Capítulo 4. .................................................................................................................................... 101

Conificación de agua en pozos mal cementados en el E. A. (T. R. ‐ Formación Productora).     ............... 102
Análisis del efecto del gasto de aceite (qoBo) y la permeabilidad de fracturas (kf).     .................................. 103
Viscosidad del aceite (μo).     ......................................................................................................................... 105
Obtención de la correlación para determinar la altura máxima de la conificación de agua.     ................... 110

Página 4 de 170
Capítulo 5. .................................................................................................................................... 115

Tiempo de formación del cono de agua (thwc). ...................................................................................... 116
Tiempo de restablecimiento del contacto agua‐aceite (tCwo). ............................................................... 125

Capítulo 6. .................................................................................................................................... 128

Aplicación práctica de las correlaciones obtenidas.     .............................................................................. 129
Ejemplo 1. Análisis de la conificación de agua del pozo Ku‐45D.     ............................................................. 129
Ejemplo 3. Cálculo de la distancia adicional del cono de agua por efecto del E. A. mal cementado
del pozo Maloob‐416.     ............................................................................................................................... 139
Ejemplo 4. Explotación de un pozo con producción de aceite extra pesado ............................................... 143

Conclusiones. ................................................................................................................................ 150

Apéndice A. ................................................................................................................................... 152

Historia del comportamiento de producción de los yacimientos Ku, Maloob y Zaap.     ......................... 153
Litología.     ................................................................................................................................................ 155
Análisis PVT.     ........................................................................................................................................... 156
Contacto Agua‐Aceite.     ........................................................................................................................... 157
Contacto Gas‐Aceite del yacimiento Ku.    . .............................................................................................. 158
Pruebas de Interferencia y Trazadores.     ................................................................................................. 159
Modelo Numérico.    . ................................................................................................................................ 162

Nomenclatura. .............................................................................................................................. 163

Referencias. .................................................................................................................................. 168

Página 5 de 170
Lista de figuras

Figura 1.1. Conificación de gas y agua. ............................................................................. 19

Figura 1.2. Representación de un poro cilíndrico con su superficie mojada por

aceite. ............................................................................................................... 20

Figura 1.3. Representación esquemática de la elevación de agua, ocasionada por la

explotación de un pozo de aceite parcialmente penetrante. ............................ 25

Figura 1.4. Sistema coordenado de referencia para la construcción de la ecuación

de equilibrio de la conificación de agua. ........................................................... 26

Figura 1.5. Representación diagramática de un pozo produciendo por drene

gravitacional. .................................................................................................... 29

Figura 1.6. Resultados normalizados de corte de agua. .................................................... 33

Figura 1.7. Formación de dos conos en un sistema de doble porosidad, escenario

de producción baja. .......................................................................................... 36

Figura 1.8. Efecto de alta producción en la posición de los dos conos. ............................. 36

Figura 1.9. Efecto de λ en el tiempo de irrupción del agua. ............................................. 37

Figura 1.10. Efecto de ω en el tiempo de irrupción del agua. ............................................. 37

Figura 1.11. Pozo parcialmente penetrante con condiciones de frontera para la

solución analítica. ............................................................................................. 42

Figuras 1.12. Curvas de gasto crítico, (1) Conificación agua-aceite curvas A, B, C, D

y E; (2) Conificación gas-aceite curvas a, b, c, d y e. ....................................... 47

Figura 1.13. Funciones adimensionales para reD = 5. .......................................................... 50

Figura 1.14. Funciones adimensionales para reD = 10. ........................................................ 51

Figura 1.15. Funciones adimensionales para reD = 20. ........................................................ 51

Figura 1.16. Funciones adimensionales para reD = 30. ........................................................ 52

Figura 1.17. Funciones adimensionales para reD = 40. ........................................................ 52

Figura 1.18. Funciones adimensionales para reD = 60. ........................................................ 53

Página 6 de 170
Figura 1.19. Funciones adimensionales para reD = 80. ........................................................ 53

Figura 1.20. Altura adimensional del cono versus tiempo adimensional. ............................. 54

Figura 2.1. Estado mecánico típico de un pozo terminado en la formación BTPKS. .......... 59

Figura 2.2. Aceleración del avance del agua por E. A. sin cemento. ................................. 60

Figura 2.3. Estado mecánico típico de un pozo terminado en la formación KM. . ............... 60

Figura 2.4. Modelo radial de simulación de un pozo. ......................................................... 61

Figura 2.5. Zona de aceite del yacimiento, representada por 200 capas de 1 m de
espesor. ............................................................................................................ 62

Figura 2.6. La zona de agua ó acuífero se representa con 7 capas de espesor

creciente. . ......................................................................................................... 64

Figura 2.7. Espacio anular sin cemento. ............................................................................ 65

Figura 2.8. Litología de los yacimientos del Cretácico. ..................................................... 65

Figura 2.9. Variación de la relación de solubilidad del aceite en función de presión. .......... 66

Figura 2.10. Variación del factor de volumen del aceite en función de la presión. ............... 67

Figura 2.11. Variación de la viscosidad del aceite a temperatura del yacimiento en

función de presión. ........................................................................................... 67

Figura 2.12. Permeabilidades relativas agua-aceite de la matriz. ........................................ 68

Figura 2.13. Permeabilidades relativas agua-aceite de la fractura. ...................................... 68

Figura 2.14. Presión capilar agua-aceite en matriz. ............................................................. 69

Figura 3.1a. Crecimiento del cono de agua a un tiempo t1. .................................................... 72

Figura 3.1b. Crecimiento del cono de agua a un tiempo t2, donde t2 > t1. .............................. 72

Figura 3.2a. Avance pseudoestacionario del cono de agua a t1. . .......................................... 73

Figura 3.2b. Avance pseudoestacionario del cono de agua a t2, donde t1< t2. ....................... 73

Figura 3.2c. Avance pseudoestacionario del cono de agua a t3, donde t2< t3. ....................... 73

Figura 3.3. Representación de la altura máxima de la conificación de agua, hwc. ............... 75

Figura 3.4a. Conificación de agua en función de la producción de aceite. .......................... 77

Página 7 de 170
Figura 3.4b. Ajuste de la conificación hwc con una ecuación de potencias. ......................... 77

Figura 3.5. Conificación de agua en función de la producción de aceite a

condiciones de yacimiento. . ............................................................................. 78

Figura 3.6a. Curvas de hwc en función kf en escala cartesiana. ........................................... 81

Figura 3.6b. Curvas de hwc en función kf en escalas logaritmicas. . ..................................... 81

Figura 3.7. Profundidad de dos intervalos analizados. ...................................................... 81

Figura 3.8. Comparación de hwc alcanzada en cada intervalo disparado. . ........................ 81

Figura 3.9. Ajuste de curvas de hwc para varias kf a ecuaciones de potencias. ................. 83

Figura 3.10. Curvas de hwc en función kf para μo = 5 cp. ...................................................... 86

Figura 3.11. Curvas de hwc en función kf para μo = 10 cp. .................................................... 86

Figura 3.12. Curvas de hwc en función kf para μo = 20 cp. .................................................... 86

Figura 3.13. Curvas de hwc en función kf para μo = 40 cp. .................................................... 86

Figura 3.14. Curvas de hwc en función kf para una partición de porosidad matriz-
fractura de 50%-50%. ....................................................................................... 90

Figura 3.15. Curvas de hwc en función kf para una partición de porosidad matriz-
fractura de 85%-15%. ....................................................................................... 90

Figura 3.16. Comparación de hwc en función del porcentaje de porosidad efectiva

total correspondiente al sistema de porosidad secundaria y k f = 1 Darcy

y μo = 2 cp . ...................................................................................................... 91

Figura 3.17. Comparación de hwc en función del porcentaje de porosidad efectiva

total correspondiente al sistema de porosidad secundaria y k f = 5 Darcy

y μo = 2 y 20 cp . .............................................................................................. 91

Figura 3.18. Comparación de hwc en función del porcentaje de porosidad efectiva

total correspondiente al sistema de porosidad secundaria y
k f = 10 Darcy y μo = 2 y 20 cp . ...................................................................... 92

Página 8 de 170
Figura 3.19. Comparación de hwc en función del porcentaje de porosidad efectiva
total correspondiente al sistema de porosidad secundaria y la relación
kf μo . ........................................................................................................... 92

Figura 3.20. hwc en función de densidad y viscosidad del aceite. . ....................................... 93

Figura 3.21a. hwc en función del radio de drene del pozo. .................................................... 94

Figura 3.21b. hwc en función de Ln ( re rw ) . ........................................................................... 94

Figura 3.22. Comportamiento del coeficiente a de la ecuación (3.1) en función de

k f μo . .............................................................................................................. 97

Figura 3.23. Comportamiento del coeficiente b de la ecuación (3.1) en función de

k f μo . .............................................................................................................. 97

Figuras 3.24. Comparación de los resultados del modelo de flujo hacia un pozo
denotados por familias de puntos y la ecuación (3.6) curvas continuas. ........ 100

Figura 4.1. Comparación de hwc en pozos con buena y mala cementación de la T. R.

de explotación. ............................................................................................... 103

Figura 4.2. Distribución de la envolvente de agua para μo=2 cp, qo = 6 Mrb/D y tp=6
meses. ............................................................................................................ 105

Figura 4.3. Distribución de la envolvente de agua para μo=5 cp, qo = 6 Mrb/D y tp=6
meses. ............................................................................................................ 105

Figura 4.4. Distribución de la envolvente de agua para μo=10 cp, qo=6 Mrb/D y tp=6
meses. ............................................................................................................ 106

Figura 4.5. Distribución de la envolvente de agua para μo=20 cp, qo=6 Mrb/D y tp=6
meses. ............................................................................................................ 106

Figura 4.6. Distribución de la envolvente de agua para μo=40 cp, qo=6 Mrb/D y tp=6
meses. ............................................................................................................ 106

Figura 4.7. Comparación de hwc en pozos con buena y mala cementación de la T. R.

de explotación y μo = 5 cp. ............................................................................. 109

Página 9 de 170
Figura 4.8. Comparación de hwc en pozos con buena y mala cementación de la T. R.
de explotación y μo = 10 cp. ........................................................................... 109

Figura 4.9. Comparación de hwc en pozos con buena y mala cementación de la T. R.

de explotación y μo = 20 cp. ........................................................................... 109

Figura 4.10. Comparación de hwc en pozos con buena y mala cementación de la T. R.

de explotación y μo = 40 cp. ........................................................................... 109

Figura 4.11. Comportamiento del coeficiente a de la ecuación (3.1) en función de

k f μo . ............................................................................................................ 112

Figura 4.12. Comportamiento del coeficiente b de la ecuación (3.1) en función de

k f μo . ........................................................................................................... 112

Figuras 5.1. Tiempo de formación del cono de agua en función k f μ o para aceites

de viscosidad de 2, 5, 10 y 40 cp. .................................................................. 118

Figura 5.2. Coeficiente a vs. qoBo de la ecuación del tiempo de formación del cono
de agua. .......................................................................................................... 121

Figura 5.3. Coeficiente b vs. qoBo de la ecuación del tiempo de formación del cono
de agua. .......................................................................................................... 121

Figuras 5.4. Ajuste de thwc en función qoBo para aceites de viscosidad de 2, 5, 10, 20
y 40 cp. . .......................................................................................................... 122

Figuras 5.5. Comparación de los ajustes obtenidos con las ecuaciones (5.3) y (5.5)
de thwc para aceites de viscosidad de 2, 5, 10, 20 y 40 cp. ............................ 124

Figuras 5.6. Tiempo de cierre requerido para que se abata el cono de agua para
aceites de viscosidad de 2, 5, 10, 20 y 40 cp. ................................................ 126

Figura 6.1. Historia de explotación del pozo Ku-45D. ...................................................... 130

Figura 6.2. Altura del cono de agua en el pozo Ku-45D. .................................................. 133

Figura 6.3. Control de la conificación de agua en el pozo Ku-45D. .................................. 134

Figura 6.4. Historia de Explotación del pozo Ku-81. ......................................................... 137

Figura 6.5. Control de la conificación de agua en el pozo Ku-81. . ................................... 138

Página 10 de 170
Figura 6.6. Registros geofísicos del pozo Maloob-416. .................................................... 140

Figura 6.7. Registro PLT del pozo Maloob-416, donde se observa flujo atrás de la
TR en el intervalo 3568.5-3560 m. ................................................................. 141

Figura 6.8. Determinación del tiempo de cierre del pozo para la reconfiguración de
cono de agua. . ............................................................................................... 146

Figura 6.9. Administración de la producción de un pozo, evitando producir agua. .......... 147

Figura 6.10. Conificación de agua 20 m, menor a la distancia entre el Cwo y la base

del intervalo disparado. .................................................................................. 148

Figuras A.1. Ubicación de los campos Ku, Maloob y Zaap. ............................................... 153

Figura A.2. Historia de producción de los campos Ku, Maloob y Zaap. ........................... 154

Figuras A.3. Litología de los yacimientos Ku, Maloob y Zaap. ........................................... 155

Figura A.4. Historia de presión de los yacimientos del Cretácico Ku, Maloob y Zaap. ........ 156

Figura A.5. Distribución de estudios PVT de los campos Ku, Maloob y Zaap. ................. 157

Figura A.6. Avance de los contactos del yacimiento Ku Cretácico. .................................. 158

Figura A.7. Inyección de metano tritiado en el yacimiento Ku Cretácico. ......................... 160

Figura A.8. Configuración estructural de los campos Ku, Maloob y Zaap Cretácico. ....... 161

Figura A.9. Malla del modelo numérico de los yacimientos Ku, Maloob y Zaap
Cretácico. ....................................................................................................... 162

Página 11 de 170
Índice de tablas

Tabla 2.1. Espesor radial de las celdas del modelo. ......................................................... 63

Tabla 3.1. Altura máxima de la conificación de agua en función del gasto de aceite. ......... 76

Tabla 3.2. Resultados del modelo de simulación; los escenarios sombreados se

descartan ya que el cono de agua no se desarrolló completamente. .............. 79

Tabla 3.3. Resultados de hwc de dos intervalos disparados a diferentes

profundidades, para permeabilidades de fractura de 1, 0.7 y 0.5 Darcys.
Los escenarios sombreados se descartan ya que el cono de agua no se
desarrolló completamente. ............................................................................... 82

Tabla 3.4. Valores de hwc en función qoBo. para μo = 5 cp, considerando diferentes
k f. ...................................................................................................................... 84

Tabla 3.5. Valores de hwc en función qoBo. para μo = 10 cp, considerando diferentes
k f. ...................................................................................................................... 85

Tabla 3.6. Valores de hwc en función qoBo para viscosidades del aceite de 20 y 40
cp, considerando diferentes kf. ......................................................................... 85

Tabla 3.7. Valores de hwc en función qoBo considerando diferentes kf para una
relación de porosidad efectiva matriz-fractura de 50%-50%. ........................... 88

Tabla 3.8. Valores de hwc en función qoBo, considerando diferentes kf para una
relación de porosidad efectiva matriz-fractura de 85%-15%. ........................... 89

Tabla 3.9. Coeficientes a y b de las ecuaciones de potencias que reproducen los

valores de hwc para una partición de porosidad primaria-secundaria 0-
100%. ............................................................................................................... 95

Tabla 3.10. Coeficientes a y b de las ecuaciones de potencias que reproducen los

valores de hwc para particiones de porosidad primaria-secundaria 50-
50% y 85-15%. ................................................................................................. 96

Tabla 4.1. Comparación de hwc en función de la calidad de la cementación de la T.

R. de explotación. ........................................................................................... 104

Tabla 4.2. hwc para tuberías de revestimiento de explotación bien y mal

cementadas y μo = 5 cp. ................................................................................. 107

Página 12 de 170
Tabla 4.3. hwc para tuberías de revestimiento de explotación bien y mal
cementadas y μo = 10, 20 y 40 cp. ................................................................. 108

Tabla 4.4. Coeficientes a y b de las ecuaciones de potencias que reproducen los

valores de hwc para una partición de porosidad primaria-secundaria 0-
100%. ............................................................................................................. 110

Tabla 5.1 Tiempo requerido para la formación del cono de agua para un aceite de
viscosidad de 2 cp. . ........................................................................................ 116

Tabla 5.2. Tiempo requerido para la formación del cono de agua para aceites de
viscosidad de 5, 10, 20 y 40 cp. ..................................................................... 117

Tabla 5.3. Tiempo requerido para abatimiento del cono de agua al 25, 50, 75 y
100% de su altura máxima, para diferentes viscosidades del aceite. ............ 127

Tabla 6.1. Cálculo de la conificación de agua con diferentes correlaciones. .................. 132

Tabla 6.2. Cálculo de la conificación de agua en función de la producción del pozo

Ku-81. ............................................................................................................. 138

Tabla 6.3. Programa de explotación del pozo, para evitar la irrupción del agua. ............ 149

Tabla A.1. Pruebas de Interferencia y Trazadores. ......................................................... 159

Página 13 de 170
Resumen

En un yacimiento de aceite naturalmente fracturado con acuífero asociado, se genera

bajo ciertas condiciones de flujo un desequilibrio entre las fuerzas de gravedad y
viscosas. Este fenómeno se caracteriza por el crecimiento de un cono de agua en
sentido vertical y radial, cuando la base del cono (contacto agua-aceite) llega al radio
de drene, este adquiere su altura máxima, posteriormente toda la interfase agua-aceite
avanza verticalmente en régimen pseudoestacionario sin sufrir deformación, cuando la
interfase se encuentra a una distancia corta del intervalo productor, el agua se acelera
irrumpiendo en el pozo.

Los objetivos de este trabajo fueron: a) modelar la conificación de agua en el medio

fracturado mediante una malla radial muy fina, con capas de un metro de espesor,
considerando el espacio anular del pozo con y sin cemento, este se representó con
celdas de dos pulgadas de ancho y b) obtener correlaciones para determinar: la altura
máxima de conificación de agua, el tiempo de formación del cono de agua y el tiempo
requerido de cierre del pozo para que el contacto agua-aceite se restablezca.

De los resultados obtenidos se destaca lo siguiente:

Se concluye que en un medio fracturado dominan las fuerzas gravitacionales y

viscosas en el fenómeno de conificación de agua.

Se obtuvo un factor de corrección para determinar la altura adicional del cono de agua
del pozo por mala cementación de la tubería de revestimiento.

Se obtuvieron correlaciones en función de permeabilidad de fracturas, viscosidad del

aceite, factor de volumen del aceite, producción de aceite, densidades del agua y
aceite y radios del pozo y drene.

Las correlaciones en yacimientos naturalmente fracturados determinan lo siguiente:

gasto crítico de aceite; altura máxima de la conificación de agua, tiempo requerido de
cierre de pozo para que se restablezca el contacto agua-aceite y tiempo de formación
del cono de agua.

Página 14 de 170
Introducción

El objetivo primario de la Administración de Yacimientos, es obtener el máximo retorno

económico a las inversiones en infraestructura y perforación de pozos, operación y
mantenimiento de los mismos y costos asociados durante todo el proceso de
explotación del yacimiento.

Lo anterior implica que antes de iniciar la explotación de un yacimiento se debe tener

un cierto grado de conocimiento del mismo y conforme se vaya desarrollando mediante
la perforación de pozos nuevos y la extracción de hidrocarburos, objeto del negocio, se
obtendrá información nueva lo que permitirá caracterizar dinámicamente los elementos
del yacimiento que afectan su comportamiento de producción; con esto se pretende
incrementar el nivel de certidumbre en los estudios de Ingeniería de Yacimientos, para
optimizar el volumen de reserva a recuperar.

Una vez que se dispone de una caracterización estática y dinámica detallada del
yacimiento se generan modelos de simulación para evaluar el comportamiento de
producción de cada pozo y del yacimiento. Con estos modelos se evalúan escenarios
de explotación; es decir, se generan predicciones de producción para casa escenario
planteado. En la actualidad, como mínimo indispensable se deben de evaluar los
escenarios de explotación siguientes: agotamiento natural, pozos de relleno, sistemas
artificiales de producción, recuperación secundaria y/o mejorada.

Durante la explotación de un yacimiento de aceite naturalmente fracturado con entrada

natural de agua se genera bajo ciertas condiciones de flujo en el sistema de fracturas
un desequilibrio entre la fuerza de gravedad y las fuerzas viscosas de los fluidos. Este
fenómeno se acentúa en la vecindad del pozo y se caracteriza por el crecimiento
paulatino de un cono de agua en sentido vertical y radial. Cuando la base del cono,
ubicada en el contacto agua-aceite, alcanza al radio de drene el cono de agua adquiere
su altura máxima; posteriormente la interfase agua-aceite avanza en régimen pseudo-
estacionario sin sufrir deformación.

Cuando la interfase agua-aceite del cono de agua se encuentra a una distancia corta
del intervalo disparado de los pozos productores, el avance del agua se acelera

Página 15 de 170
provocando su irrupción en el pozo y un incremento súbito en la relación agua-aceite,
RAA, en detrimento del ritmo de producción del aceite.

Cuando los campos petroleros no cuentan con instalaciones para separar, tratar y
manejar el agua de formación, los pozos que producen con mayores porcentajes de
agua se tienen que cerrar o bien someterlos a tratamientos de control de agua y/o
reparaciones mayores, para reducir el contenido de agua en la corriente de aceite y
evitar penalizaciones por su venta.

La conificación de agua es un problema muy serio y costoso que se presenta en

muchos yacimientos, especialmente en los principales yacimientos de aceite de
México, mismos que producen de rocas carbonatadas naturalmente fracturadas; estos
yacimientos contienen más del 50% de la reserva probable de hidrocarburos de la
Nación.

Las principales formaciones productoras de los campos de la Sonda de Campeche son

carbonatos naturalmente fracturados del Cretácico (BTPKS, KM y KI) y Jurásico
Superior Kimeridgiano (JSK) de espesor y/o relieve estructural considerable.

La formación Brecha Terciario Paleoceno Cretácico Superior (BTPKS), está constituida

de rocas carbonatadas densamente fracturadas de muy alta permeabilidad, esto
ocasiona severos problemas durante las operaciones de cementación de la última
Tubería de Revestimiento (T. R.), ya que al circular cemento al pozo para ademarlo
prácticamente todo el cemento se pierde en la formación, entonces el Espacio Anular
(E. A.) entre la T. R. y la formación productora queda sin cemento, formándose un
espacio vacío, es decir un canal sin restricción al flujo, por el cual se magnifica la
conificación del agua.

La última reserva de hidrocarburos de los yacimientos de la formación Cretácico de la

Sonda de Campeche se recuperará en un intervalo, o columna de aceite reducido,
limitado en su parte superior por el contacto gas-aceite (Cgo) y en la inferior por el
contacto agua-aceite (Cwo). Por lo anterior, es de suma importancia conocer el
comportamiento de la producción en estas condiciones para definir el rango de

Página 16 de 170
operación óptimo al cual se explotarán los pozos y que permita maximizar la
recuperación de aceite y controlar los problemas de conificación de gas y agua.

El comportamiento de conificación y canalización de fluidos en yacimientos

naturalmente fracturados depende de las características de la formación y de los
fluidos, de la geometría y estado mecánico de los pozos y de las condiciones de
producción en que se operen. La diversidad de los tipos de aceite encontrados en los
yacimientos de la Sonda de Campeche, impone un reto adicional pues a medida que la
calidad de aceite disminuye, su densidad aumenta aproximándose a la del agua, así
mismo la su viscosidad del aceite también aumenta, generando condiciones más
propicias para la conificación de agua.

Los objetivos de este trabajo son: modelar con detalle la conificación de agua en el
medio fracturado mediante una malla radial muy fina de un pozo, con capas de un
metro de espesor, considerando el E. A. con y sin cemento mismo que se modela con
celdas de 2” de espesor y como en la actualidad no es posible implementar este nivel
de detalle en los modelos de todo el campo, se obtienen ecuaciones para determinar la
altura máxima de la conificación y el tiempo que tarda en formarse el cono de agua,
con el objeto de implementar estas restricciones en los modelos de campo e impactar
en el volumen de reserva de hidrocarburos a recuperar.

Página 17 de 170
Capítulo I

Revisión Bibliográfica de Conificación de Agua

Página 18 de 170
Revisión Bibliográfica de Conificación de Agua

“Conificación de agua” 1 es el término empleado para el mecanismo de entrada de agua

subyacente a pozos productores de aceite. La conificación de agua, se asocia
generalmente con altas producciones de aceite y se desarrolla únicamente bajo ciertas
condiciones de flujo, es el resultado del movimiento de fluidos en la dirección de la
mínima resistencia, contrabalanceada por la tendencia de los mismos a mantener un
equilibrio gravitacional.

En un sentido más amplio el fenómeno de conificación de agua, constituye uno de los

problemas más complejos de la producción de aceite. En la formación productora bajo
condiciones estáticas, el agua permanecerá en el fondo y el aceite en la parte superior,
debido a que el aceite es menos denso que el agua.

La conificación de agua ó gas, ocurre cuando la componente vertical de las fuerzas

viscosas excede la fuerza de la gravedad (figura 1.1). Considerando gasto constante y
mantenimiento de presión en el contacto agua-aceite, el cono de agua tiene su altura
máxima, cuando la base del cono (contacto agua-aceite) alcanza el radio de drene del
pozo; posteriormente, la interfase agua-aceite avanza en régimen pseudoestacionario
sin sufrir deformación2; el cono de agua es estable hasta una distancia corta de la base
del intervalo productor, pero sí se incrementa ligeramente la producción de aceite, se
desestabilizará, provocando su irrupción en el pozo, a esta condición de flujo se le
denomina gasto critico. 3, 4 ,11

Gas

Aceite Aceite

Agua

Figura 1.1. Conificación de gas y agua.

Página 19 de 170
Esencialmente existen tres fuerzas que afectan la distribución de los fluidos alrededor
del fondo del pozo, estas son: Fuerzas capilares (C), gravitacionales (G) y viscosas (V);
estas fuerzas se pueden representar con números adimensionales. 5

Las fuerzas capilares usualmente tienen efectos despreciables en el fenómeno de la

conificación y por consecuencia se desprecian; las fuerzas gravitacionales actúan en el
sentido vertical y dependen de las diferencias de densidad de los fluidos.

El término de fuerzas viscosas, se refieren a los gradientes de presión asociados al

flujo de fluidos en el medio poroso, mismos que son determinados con la ley de Darcy.

A cualquier tiempo, existe un balance entre las fuerzas gravitacionales y viscosas en

puntos cercanos y lejanos al pozo; cuando las fuerzas dinámicas (viscosas) en el pozo
exceden las fuerzas gravitacionales, el cono de agua ó gas irrumpe en el pozo.

La relación entre las fuerzas capilares y viscosas, se pueden explicar fácilmente con un
ejemplo de flujo lineal, en un poro cilíndrico de longitud, L y diámetro, 2R, como se
ilustra en la figura 1.2.

0 X L

Agua Aceite 2L

Δpw pc Δpc

Figura 1.2. Representación de un poro cilíndrico con su superficie mojada por

aceite.

Página 20 de 170
La caída de presión viscosa, Δpv, es la suma de la caída de presión en las fases agua,
(Δpw,) y en la fase aceite, (Δpo.); usando la ley de Poiseuille se tiene la expresión
siguiente,

Δpv = Δpw + Δpo = 8

vw
R2
( μw X + μo ( L − X ) ) . ............................................ (1.1)

Donde,

μo , Viscosidad del aceite.

μw , Viscosidad del agua.

X, Posición de la interface capilar y es igual a L/2.

vw , Velocidad media del agua.

La caída de presión en un tubo capilar pc , a través de la interfase está dada por la

expresión siguiente:

2σ wo cos (θ )
pc = . ................................................................................. (1.2)
R

Donde,

R, Radio del capilar,

σ wo , Tensión interfasial agua-aceite,

θ, Ángulo de la fase mojante.

Página 21 de 170
La relación entre las fuerzas capilares y viscosas se expresa como:

Δpv vw μo 2L ⎛ μw ⎞
= ⎜ + 1⎟ . .............................................................. (1.3)
pc σ wo cos (θ ) R ⎝ μo ⎠

Considerando que la viscosidad del agua siempre es menor a la viscosidad de un

aceite negro, μw < μo , y que en un tubo capilar la relación de la longitud sobre el radio

del tubo normalmente es mucho mayor a uno; mientras que esta relación en una red de
poros interconectados se puede asumir que es cercana a uno, 2L/R~1, entonces a la
relación de las fuerzas viscosas entre las capilares se le denomina número capilar.

vw μo
Nc = . .................................................................................. (1.4)
σ wo cos (θ )

Donde,

Nc , Número capilar,

vw , Velocidad media del agua,

μo , Viscosidad del aceite,

σ wo , Tensión interfasial,

θ, Ángulo de la fase mojante.

Considerando gotas pequeñas de agua y que estas son esféricas con un radio Rd, la
diferencia de presión gravitacional ΔpG, que fuerza a la gota de agua hacia abajo es:

Página 22 de 170
Rd
ΔpG = Δρ wo g . .................................................................................. (1.5)
3

Donde, Δρwo = ρw − ρo es la diferencia de densidades y g es la constante gravitacional.

La relación de la fuerza gravitacional sobre la fuerza capilar se expresa como:

ΔpG Rd R Δρ wo g
= . ........................................................................ (1.6)
pc 3 σ wo cos (θ )

Rd R
Si asociamos la superficie con una superficie imaginaria caracterizada por la
3
permeabilidad absoluta k del medio poroso, la relación de la fuerza gravitacional sobre
la fuerza capilar puede expresarse con el número de Bond,

k Δρ wo g
NB = . ................................................................................. (1.7)
σ wo cos (θ )

Donde,

NB , Número de Bond,

k, Permeabilidad,

g, Aceleración de la gravedad,

σ wo , Tensión interfasial agua-aceite,

θ, Ángulo de la fase mojante.

Página 23 de 170
El número adimensional que describe la relación de las fuerzas viscosas sobre las
gravitacionales, se obtiene fácilmente, al dividir el número capilar entre el número de
Bond, N (V / G ) = NC N B , entonces se tiene que:

vw μ o
N (V / G ) = . ............................................................................. (1.8)
k Δρ wo g

Donde, vw es la velocidad media del flujo de agua en el medio poroso.

μo , Viscosidad del aceite,

k, Permeabilidad,

Δρ wo , Diferencia de densidades agua-aceite,

g, Aceleración de la gravedad.

La determinación de la conificación de agua es crítica para desarrollar yacimientos con

una zona de espesor reducido de aceite y predecir su comportamiento de producción
hasta su agotamiento. 6

En la mayoría de los casos, la heterogeneidad en los patrones de fractura puede

desarrollar conos asimétricos con el potencial de irrumpir el agua en todo el intervalo
productor. En estos sistemas, el parámetro más importante es el patrón de fracturas,
por lo que sí no está caracterizado con mucho detalle, ninguna correlación ó modelo de
simulación puede ayudarnos a predecir gastos críticos y tiempos de irrupción. 2, 7

Muskat1, estableció que cuando un pozo inicia a producir en presencia de un contacto

agua-aceite que actúa como frontera externa, se genera un gradiente de presión lo que
provoca la elevación del agua subyacente creando un cono de agua, el cual dejará de

Página 24 de 170
crecer cuando alcance el equilibrio con el gradiente de presión (figuras 1.3 y 1.4). La
presión en la interfase agua-aceite se puede expresar como:

p ( r , z ) + ρ w g ( h − z ) = p ( r , z ) + ρ w gy . .................................................... (1.9)

La ecuación (1.9) representa una condición necesaria para obtener un equilibrio

estático en la interfase agua-aceite. Para mantener un equilibrio dinámico la interfase
agua-aceite debe coincidir con la envolvente que limita a la zona de aceite (figura 1.3),
dado que la función de distribución de presión p(r, z) depende de la forma y de la
posición lugar del contacto agua-aceite, entonces la ecuación (1.9) puede proporcionar
la localización y la forma del cono de agua. Sin embargo desde un punto de vista
práctico esto no es factible, ya que es prácticamente imposible obtener la solución
analítica de la distribución de presión p(r, z) en un sistema tridimensional heterogéneo,
a menos que se fijen todas las fronteras a una geometría simple.

Aceite
Aceite
h
h

Agua
Agua

Figura 1.3. Representación esquemática de la elevación de agua,

ocasionada por la explotación de un pozo de aceite
parcialmente penetrante.

Página 25 de 170
Si se considera un cono de agua en condiciones estáticas por abajo del pozo, es claro
que su altura (gradiente gravitacional) se equilibra con la Δp del pozo como se ilustra
en la figura 1.4. Existe un límite al cual la altura del cono de agua permanece en
equilibrio estático; sí se sobrepasa este límite y el pico del cono de agua se encuentra
cercano al pozo, entonces el agua irrumpirá en el pozo; como el flujo de aceite
converge hacia el pozo, el gradiente de presión es mayor conforme el cono se acerca
al intervalo productor.

Pozo
Pozo

pt pt
p
Pww

Aceite (ρο)
h zz
p (r, z)
r
r y
pb

Agua (ρω )

Z
Z

Figura 1.4. Sistema coordenado de referencia para la construcción de la

ecuación de equilibrio de la conificación de agua.

De otra forma se puede visualizar que, sí la caída de presión Δp del pozo continua
creciendo por efecto del incremento de la producción, el cono de agua continuará
desarrollándose hasta alcanzar el intervalo productor presentándose la irrupción del
agua. Aplicando la ecuación (1.9) para determinar las condiciones críticas así como la
forma de la interfase agua-aceite, se realizan las transformaciones siguientes:

Página 26 de 170
k ⎫
φ= ( p − ρo gz ) , ⎪
μ ⎬ , .......................................................... (1.10)
Δp = pt − pw = pb − ρo gh − pw ,⎭⎪

Donde:

g, Aceleración de la gravedad,

h, Espesor impregnado de aceite,

k, Permeabilidad,

p, Presión,

pb , Presión en la base de la formación,

pt , Presión en la cima de la formación.

z, Distancia de la cima de la formación a la envolvente de agua,

μ, Viscosidad,

φ, Función potencial,

ρo , Densidad del aceite,

Δp , Es la diferencia entre las presiones estáticas y de fondo fluyendo, medidas a la

profundidad del intervalo productor.

La diferencia de potencial en cualquier punto (r, z) de la superficie del cono de agua y

el pozo, está dada por la expresión siguiente:

k
Δφ = φ ( r , z ) − φw ( z = 0 ) = ⎡ Δp − g ( h − z ) Δρ ⎤⎦ , ................................... (1.11)
μ⎣

Página 27 de 170
donde, Δρ = ρw − ρo . Puede notarse que a lo largo de la distancia r,

Δφ = ( k μ ) Δp ≡ ( Δφ )e ; entonces la ecuación (1.11) se puede escribir como:

Δφ ( r , z ) ghΔρ ⎛ z ⎞
= 1− ⎜1 − ⎟ . ................................................................ (1.12)
( Δφ )e Δp ⎝ h ⎠

La ecuación (1.12), es la ecuación básica, para determinar la forma de la interfase

agua-aceite z = z ( r ) . Como se indicó previamente, la función potencial φ ( r , z ) y

Δφ ( r , z ) , se considera que corresponde a un líquido incompresible fluyendo hacia un

pozo parcialmente penetrante. De esta forma la altura máxima del cono de agua se
obtiene en r = rw ; rescribiendo la ecuación (1.12) se obtiene,

⎛ z⎞ ⎡ Δφ ( rw , z ) ⎤ Δpw
⎜1 − ⎟ = ⎢1 − ⎥
⎝ h ⎠ max ⎢⎣ ( Δφ )e ⎥⎦ ghΔρ
. ..................................................... (1.13)
⎛z⎞ ⎡ Δφ ( rw , z ) ⎤ Δpw
⎜ ⎟ = 1 − ⎢1 − ⎥
⎝ h ⎠ max ⎢⎣ ( Δφ )e ⎥⎦ ghΔρ

Se percibe que la teoría presentada puede mejorarse; sin embargo, hay pocas razones
para dudar de los principios físicos de esta discusión. Por ejemplo, sí en un pozo
parcialmente penetrante ya se desarrollo un cono de agua, una solución efectiva es
efectuar una reparación mayor para alejarse del agua, a menos de que la ventana de
aceite sea muy reducida. Sin embargo primero se debe de tratar de reducir la Δp del
pozo estrangulándolo; esto implica que al reducir la producción ó la Δp del pozo por
debajo de sus valores críticos después de que el agua ha irrumpido, el cono de agua se
reconfigura a un nuevo nivel estable.

Página 28 de 170
Para el problema de la conificación de gas en yacimientos de aceite, es decir cuando el
yacimiento cuenta con un casquete de gas, el pozo debe penetrar completamente la
zona de aceite y ser disparado en la parte inferior de la zona de aceite; entonces las
ecuaciones obtenidas se pueden aplicar realizando cambios pequeños.

Para la determinación de la relación de la penetración parcial se considera únicamente

la zona de aceite y el intervalo productor, mientras que el término Δρ es la diferencia
entre las densidades del aceite y de la fase gaseosa, ver figura 1.5.

rw S
pozo
del pozo

B
B
p ==00
p
Eje del

he
h
e
Eje

A
A

hw re

Figura 1.5. Representación diagramática de un pozo produciendo por

drene gravitacional.

La producción de aceite en régimen permanente a condiciones de superficie está dada

por la ecuación siguiente:

khΔpw
qo = 0.003076 ( STB / D ) . .......................................... (1.14)
μo Bo log ( re rw )

Página 29 de 170
Donde la ecuación (1.14) esta en unidades de campo.

Un caso especial, es cuando el pozo apenas llega a la cima del yacimiento, entonces la
ecuación (1.14) se transforma a:

k Δpw rw
qo = 0.003076 ( STB / D ) . ...................................................... (1.15)
μ o Bo

En esta ecuación la producción de aceite varía directamente proporcional con el radio

del pozo y es independiente del radio de drene, mientras que en el caso de que el pozo
penetre más en el yacimiento, la producción es función del logaritmo de la relación de
los radios de drene y del pozo.

Cuando un pozo produce por drene gravitacional gas-aceite, su producción está dada
esencialmente por la Ec. (1.9). Por supuesto esta situación se alcanza sólo cuando la
presión de fondo fluyendo es menor a la presión de la cima del yacimiento; además, se
considera que el potencial del aceite en el pozo es constante desde el fondo del mismo
hasta su nivel dinámico (punto A). En el yacimiento, el aceite saturado está limitado por
una superficie libre Fig. 1.5, la cual al mismo tiempo es una superficie-frontera a
presión constante. Sí el nivel de fluido en la frontera externa ( re ) es he y es igual o

menor al espesor impregnado de aceite, mientras que en el pozo la columna de aceite

es hw , entonces la capacidad de producción de la formación en unidades de campo es:

6.667 × 10 −4 k ρ o ( he2 − hw2 )

qo = ( STB / D ) . ........................................... (1.16)
μ o Bo log ( re rw )

Página 30 de 170
Kuo y DesBrisay8, usaron un modelo de simulación numérica (Intercomb Beta II) para
realizar un estudio de sensibilidad a los parámetros principales que afectan a la
conificación de agua; este modelo radial consideró 6 celdas en r y 8 capas. Los
resultados de la simulación se emplearon para generar correlaciones simples para la
predicción de la conificación de agua. Su estudio cubrió un rango amplio de los
parámetros del yacimiento, los cuales son:

Relación kv kh , 0.01, 0.1, 0.5 y 1.0,

hp
21.4% ≤ ≤ 83.3% , donde h p es el intervalo disparado y h es el espesor
h
total impregnado de aceite,

500 ≤ q ≤ 2000 Producción total de líquido en rb/D,

ko
μo
Relación de movilidad , 1.0, 4.36 y 9.81.
kw
μw

La correlación es válida para predecir el comportamiento de producción de agua-aceite

para la mayoría de yacimientos con acuífero activo, excepto para aquellos yacimientos
con barreras verticales locales.

Kuo y DesBrisay realizaron 13 corridas combinando los parámetros mencionados

previamente, de los resultados obtenidos, definieron dos variables adimensionales,
tiempo adimensional, tD y corte de agua adimensional, (WC ) D .

Así,

t
tD = , ............................................................................................ (1.17)
t BT

Página 31 de 170
WC
(WC ) D = , ............................................................................ (1.18)
(WC )limit
Mhw
(WC )limit = . ............................................................................ (1.19)
Mhw + h

Donde,

h, Espesor de la zona de aceite, pies

hw , Espesor de la zona de agua, pies

M, Relación de movilidad, adimensional,

t, Tiempo, días,

tBT , Tiempo de irrupción, días, calculado con la correlación

de Bournazel y Jeanson9

WC , Corte de agua, fracción,

(WC )limit , Valor máximo de corte de agua, fracción.

Para cada caso el tiempo de irrupción se calculó usando la correlación de Bournazel;

posteriormente el tiempo real se transformó a tiempo adimensional y el corte de agua
se transforma a corte de agua adimensional usando el valor máximo de corte de agua,
mismo que está en función de la relación de movilidad, espesores de la zona de agua y
aceite.

Conforme se explota el pozo, el contacto agua-aceite se desplaza hacia arriba, por lo

que el valor límite de corte de agua también cambia.

Los resultados obtenidos se graficaron en la figura 1.6, la cual muestra el tiempo

adimensional con respecto al corte de agua adimensional. Estos datos se ajustaron a
las correlaciones siguientes:

Página 32 de 170
(WC ) D = 0 , tD < 0.5 , ............................. (1.20)

(WC ) D = 0.94 log tD + 0.29 , 0.5 ≤ tD ≤ 5.7 , ...................... (1.21)

(WC ) D = 1.0 , tD > 5.7 . ............................... (1.22)

Estas correlaciones representan la mayoría de los datos con una precisión del 90%; en
algunos casos disminuye a 80%.

1.4

1.2
adimensional, (WC)D

1.0
Corte de agua

0.8

0.6

0.4

0.2

0
0.1 1.0 10 100
Tiempo adimensional, tD

Figura 1.6. Resultados normalizados de corte de agua.

Bournazel y Jeanson9, desarrollaron un método para estimar el tiempo de irrupción,

con base en datos experimentales. Su modelo considera un yacimiento homogéneo,
flujo radial de aceite y agua en la frontera externa.

Z
( tD ) BT = , ............................................................................. (1.23)
3 − 0.7 Z

Página 33 de 170
0.0037Δρ k H h ( h − D )
Z= , .................................................................. (1.24)
μo Bo qo

μoφ h ( tD ) BT
t BT = . ............................................................... (1.25)
0.00137 Δρ kV (1 + M α )

Donde,

Bo , Factor del volumen del aceite, rb/STB

D, Intervalo perforado, desde la cima de la

pies
formación hasta la base de los disparos,

h, Espesor impregnado de aceite, pies

kH , Permeabilidad horizontal, md

kV , Permeabilidad vertical, md

M, Relación de movilidad, adimensional

qo , Gasto de aceite, STB/D

Z, Altura del cono de agua, adimensional

μo , Viscosidad del aceite, cp

Δρ , Diferencia de densidades, lb/ft3

( t D )BT , Tiempo de irrupción adimensional, adimensional

tBT , Tiempo de irrupción, días

α = 0.5 para M < 1; 0.6 para 1 < M < 10 .

Página 34 de 170
Al-Afaleg y Ershaghi7, analizaron el desarrollo de la conificación de agua en
yacimientos naturalmente fracturados y la irrupción de la misma en el pozo.
Adicionalmente desarrollaron correlaciones para estimar el tiempo de irrupción del agua
en pozos produciendo a un gasto mayor al crítico.

Ellos encontraron que a permeabilidades verticales altas del medio fracturado, el

proceso de la conificación de agua se acelera, dando como resultado la disminución del
gasto crítico y menores tiempos de irrupción del agua en el pozo.

Además de considerar los parámetros del yacimiento que intervienen en las

correlaciones analíticas y empíricas ya analizadas, Al-Afaleg y Ershaghi examinaron el
contraste de permeabilidad y porosidad de la matriz y de la fractura en términos de la
imbibición contracorriente.

Para realizar lo enunciado construyeron un modelo radial de pozo de doble porosidad,

donde además analizaron los cambios de saturación en los dos medios porosos para
conocer el avance del agua.

En la figura 1.7 se muestran los resultados del modelo de simulación de los dos medios
porosos, donde se aprecia la formación de dos conos. Dependiendo de la producción,
se puede desarrollar con mayor velocidad el cono de agua en las fracturas y un cono
con menor movilidad en la matriz. La posición relativa de los dos conos uno respecto al
otro depende del ritmo de producción y de las propiedades de la formación. En la figura
1.8 se muestra el desarrollo de la conificación en los dos medios porosos a un ritmo de
producción mayor, qo = 3000 STB/D que el empleado para la figura 1.7.

Página 35 de 170
0
20 qo = 1600 STB/D λ = 10‐7
t = 450 días ω = 0.01

Profundidad (pies)
40 pc = 5 psi
60
Fracturas, Sw = 99%
80
Matriz, Sw = 40%
100
120
140
160
0 200 400 600 800 1000 1200 1400
Radio (pies)
Figura 1.7. Formación de dos conos en un sistema de doble
porosidad, escenario de producción baja.

0
20 qo = 3,000 STB/D λ = 10‐7
t = 450 días ω = 0.01
40
Profundidad (pies)

pc = 5 psi
60 Fracturas, Sw = 99%
80 Matriz, Sw = 40%
100
120
140
160
0 200 400 600 800 1000 1200 1400
Radio (pies)
Figura 1.8. Efecto de alta producción en la posición de los dos
conos.

Página 36 de 170
En las figuras 1.9 y 1.10 se observa que para gastos altos de aceite, antes de que
irrumpa el agua en el pozo λ y ω tienen poca influencia en el fenómeno de
conificación; sin embargo a ritmos bajos de producción de aceite estos parámetros si
afectan el tiempo de irrupción.

12000
ω = 0.1
10000 qo = 800 STB/D
Tiempo de irrupción (días)

pc =15 psi
8000

6000
qo = 1400 STB/D
4000

2000 qo = 2000 STB/D

0
1.E-07 1.E-06 1.E-05

Figura 1.9. Efecto de λ en el tiempo de irrupción del agua.

8000
qo = 800 STB/D
Tiempo de irrupción (días)

6000
λ = 10 ‐7
p c = 15 psi

4000
qo = 1400 STB/D

2000
qo = 2000 STB/D

0
0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10
ω

Figura 1.10. Efecto de ω en el tiempo de irrupción del agua.

Página 37 de 170
La correlación desarrollada tiene como objetivo modelar el tiempo de irrupción, para
yacimientos fracturados homogéneos y es válida para los rangos siguientes:

10−9 ≤ ωλ ≤ 10−6
qD ≤ 0.25 ........................................................................ (1.26)
0.0675 ≤ pcD

⎛ 1 ⎞
t D = A ( Log λ ) + B ( Log λ ) + cw + DLogpcD + E ⎜
2
⎟ + F ..................... (1.27)
⎝ LogqD ⎠

pc
pcD = ........................................................................................ (1.28)
p − pwf

887.31qo Bo μo
qD = ................................................................................ (1.29)
k f h 2 Δρ

tbt 2
tD = ................................................................................................ (1.30)
tbt1

El tiempo de irrupción del agua en el primer medio poroso, matriz, esta dado por la
expresión siguiente:

( Np ) BT
tBT 1 = ........................................................................................ (1.31)
qo

El tiempo de irrupción del agua para un medio naturalmente fracturado se determina

con la expresión siguiente:

tBT 2 = tD × t BT 1 ........................................................................................ (1.32)

⎛ 1 ⎞
( Np ) BT =⎜ ⎟ Aφ ( S w − Sor ) ( h − hwc − hap − hp ) ................................. (1.33)
⎝ 5.615Bo ⎠

Página 38 de 170
Donde,

A= -0.051217

B= -0.032583

C= 1.557171

D= 0.338711

E= 0.548597

F= 2.493842

A = Área, pies2

Bo = Factor del volumen del aceite, rb/STB

h = Espesor impregnado de aceite, pies

hap = Columna de aceite por arriba del intervalo disparado, pies

hp = Intervalo disparado, pies

hwc = Altura de la conificación de agua, pies

kf = Permeabilidad del medio fracturado, md

( Np ) BT = Producción acumulada al momento de la irrupción, STB

pc = Presión capilar, psi

pcD = Presión capilar adimensional, adimensional

pwf = Presión de fondo fluyendo, psi

p = Presión media del yacimiento, psi

qc = Gasto crítico de aceite, STB/D

Página 39 de 170
re = Radio de drene, pies

rw = Radio del pozo, pies

tBT 1 = Tiempo de irrupción del agua en matriz, días

tBT 2 = Tiempo de irrupción del agua en fractura, días

μo = Viscosidad del aceite, cp.

Δρ = Diferencia de densidades agua-aceite, gr/cm³

λ = Parámetro de flujo interporoso

ω = Almacenamiento del medio fracturado, adimensional

φ = Porosidad, fracción

Meyer y Garder10, presentaron una correlación para determinar el gasto máximo de

aceite al cual la interfase agua-aceite, definida por el cono de agua no invade el
intervalo productor del pozo. Ellos consideraron yacimiento homogéneo y flujo radial,
obteniendo a la expresión siguiente:

0.0015351Δρ ( h 2 − D 2 ) k
qc = ................................................................. (1.34)
Bo μ o ln ( re rw )

Donde:

Bo = Factor del volumen del aceite, rb/STB

D = Intervalo disparado, pies

h = Espesor impregnado de aceite, pies

Página 40 de 170
k = Permeabilidad, md

qc = Gasto crítico de aceite, STB/D

re = Radio de drene, pies

rw = Radio del pozo, pies

Δρ = Diferencia de densidades agua-aceite, gr/cm³.

μo = Viscosidad del aceite, cp.

Hφyland y Papatzacos11, presentan una solución analítica para determinar la

conificación de agua, la cual considera flujo multifásico y flujo transitorio en la vecindad
del pozo; además desarrollaron una correlación para predecir el gasto crítico de aceite,
obtenida con base en varias corridas de un modelo numérico de yacimientos. Sin
embargo su estudio se limitó a un pozo terminado en la cima de la formación
productora.

La solución analítica presentada por Hφyland y Papatzacos, es una extensión de la

teoría de Muskat y Wyckoff basada en el trabajo desarrollado por Papatzacos, quien
obtuvo una solución en función del tiempo de la ecuación de difusión para un fluido
ligeramente compresible, con flujo hacia un pozo de conductividad infinita en un
yacimiento infinito. Para flujo estacionario, la solución adquiere una forma simple y se
puede combinar con el método de imágenes para obtener las condiciones de frontera
vertical y lateral, como se ilustra en la figura 1.11.

Página 41 de 170
Pozo
Frontera sellante

hp re
Pw

Frontera a rw
kV
presión h Aceite
kH
constante

hwc
Agua

Frontera sellante

Figura 1.11. Pozo parcialmente penetrante con condiciones de frontera

para la solución analítica.

Para su solución de gasto crítico de aceite, Hφyland y Papatzacos tomaron las mismas
consideraciones que Muskat y Wyckoff, una fase y no se considera el cono de agua en
la distribución del potencial. Ellos desarrollaron un programa para determinar el gasto
crítico en un área cuadrada a presión constante. El área del cuadrado se transforma a
un radio equivalente para tener un área circular a presión constante para relacionar los
casos que se simularon con la figura 1.11. Los autores obtuvieron el modelo analítico
siguiente:

40, 667.25qc Bo μ o
qcD = , ...................................................................... (1.35)
k H Δρ h 2

⎛r ⎞ k
rD = ⎜ e ⎟ V . ................................................................................... (1.36)
⎝ h ⎠ kH

Página 42 de 170
Donde,

BO , Factor de volumen del aceite,

h, Espesor impregnado de aceite,

kV , Permeabilidad vertical,

kH , Permeabilidad horizontal,

qcD , Gasto crítico adimensional,

rD , Radio adimensional,

μO , Viscosidad del aceite,

Δρ , Diferencia de densidad entre el agua y el aceite.

• Hφyland y Papatzacos, además de presentar una solución analítica para

determinar la conificación de agua, obtuvieron correlaciones para predecir la
conificación de agua. Su modelo simuló el flujo hacia un pozo con un modelo
de 10 celdas en dirección radial y 20 capas, donde representaron un Acuífero
de fondo conectado en la última capa, asignándole valores de porosidad y
permeabilidad infinita para obtener el efecto de una frontera a presión constante
en el contacto agua-aceite. De igual forma modelaron la frontera externa en el
radio de drene, asignándole valores infinitos de porosidad y permeabilidad a la
última columna de celdas radiales, para que esta frontera también actuara a
presión constante. La primera columna de celdas del modelo radial representó
el pozo.

Para cierta geometría y parámetros del yacimiento así como cierta penetración del
pozo en el yacimiento, los autores determinaron el gasto crítico con un 4% de

Página 43 de 170
incertidumbre; ellos procedieron a determinar el gasto crítico cuando ya se había
formado un cono de agua estable y muy cerca del intervalo productor, de tal forma que
al incrementar ligeramente la producción del pozo, el agua irrumpía en el mismo.

Para el estudio anterior, realizaron más de 500 simulaciones, ya que para cada caso
fue necesario hacer de 5 a 6 simulaciones.

Hφyland y Papatzacos, concluyeron que el gasto crítico es una función lineal de la

permeabilidad del aceite, diferencia de densidades de los fluidos, viscosidad y factor de
volumen del aceite; mientras que el gasto crítico es una función no lineal de la
penetración del pozo en el yacimiento, radio de drene, espesor de la columna de aceite
y relación de permeabilidades.

La expresión a la que ellos llegaron es la siguiente:

1.325
k Δρ ⎡ ⎛ hp ⎞2 ⎤ −1.990
qc = ⎢1 − ⎜ ⎟ ⎥ h2.238 ⎡⎣ln ( re ) ⎤⎦ , ................................. (1.37)
10.822Bo μo ⎢ ⎝ h ⎠ ⎥
⎣ ⎦

Donde,

Bo = Factor del volumen del aceite, rb/STB

h = Espesor impregnado de aceite, pies

hp = Intervalo disparado, pies

k = Permeabilidad horizontal, md

re = Radio de drene, pies

qc = Gasto crítico, STB/D

Δρ = Diferencia de densidades agua-aceite, lbm/ft3

μO , = Viscosidad del aceite, cp

Página 44 de 170
12
Chaney, et al. , determinan curvas de conificación de agua matemáticamente y
mediante análisis potenciométrico. Las curvas muestran los gastos de aceite en barriles
por día a condiciones de yacimiento versus la distancia de la cima de la formación ó del
contacto gas-aceite a la cima del intervalo perforado ó en el caso de conificación de
agua es la distancia entre el contacto agua-aceite y la base del intervalo disparado, ver
figura 1.12.

Las curvas de la figura 1.12 (1) se muestran para penetración del pozo en la formación
de: A = 2.5 ft; B = 5.0 ft; C = 7.5 ft; D = 10.0 ft y E = 12.5 ft.

Las curvas de la figura 1.12 (2) se muestran para penetración del pozo en la formación
de: A = 5.0 ft; B = 10.0 ft; C = 15.0 ft; D = 20.0 ft y E = 25.0 ft.

Las curvas de la figura 1.12 (3) se muestran para penetración del pozo en la formación
de: A = 10.0 ft; B = 20.0 ft; C = 30.0 ft; D = 40.0 ft y E = 50.0 ft.

Nota las letras mayúsculas denotan las curvas de conificación de agua, mientras que
las letras minúsculas denotan a las curvas de conificación de gas.

Para el análisis los autores usaron la información siguiente:

Permeabilidad = 1000 md

Viscosidad del aceite = 1 cp

Diferencia de densidades agua-aceite = 0.3 gr/cm³

Diferencia de densidades gas-aceite = 0.6 gr/cm³

Radio de drene = 1000 pies

Entonces para emplear los gastos obtenidos de las curvas, se deben corregir por las
propiedades del fluido y de la formación, estas correcciones se toman en cuenta en las
ecuaciones siguientes:

Página 45 de 170
Gasto critico de aceite corregido, para un sistema agua-aceite.

0.00333ko ( Δγ ow ) qcurve
qc = , ( STB D ) . .......................................... (1.38)
μo Bo

Gasto critico de gas corregido, para un sistema gas-agua,

0.00333k g ( Δγ gw ) qcurve
qc = , ( Mscf / D ). ...................................... (1.39)
μ gV

Gasto crítico de aceite corregido, para un sistema gas-aceite,

0.00167 k g ( Δγ og ) qcurve
qc = , ( STB D ). ......................................... (1.40)
μ g Bo

Donde:

Bo = Factor del volumen del aceite, rb/STB

kg = Permeabilidad efectiva al gas, md

ko = Permeabilidad efectiva al aceite, md

kw = Permeabilidad efectiva al agua, md

qcurve = Gasto crítico a condiciones de yacimiento, rb/D

Página 46 de 170
qc = Gasto crítico corregido, STB/D

V = rb ocupados por 1 Mscf de gas.

μg = Viscosidad del gas, cp

μo = Viscosidad del aceite, cp

μw = Viscosidad del agua, cp

Δγ ow = Diferencia de densidades aceite-agua, gr/cm³

Δ γ og = Diferencia de densidades aceite-gas, gr/cm³

Δ γ gw = Diferencia de densidades agua-gas, gr/cm³

100 1000 10000

e
b
b a
d c
e
b a
e d c

d A
1000
c B
Gasto critico de aceite (STB/D)

C
Gasto critico de aceite (STB/D)

10 b 100 A
Gasto critico de aceite (STB/D)

D A
a B B
C
E C
D
D
E E
100

1 10

0.1 1 1
0 10 20 0 10 20 30 40 50 0 50 100
Distancia de la cima del yacimiento a la cima Distancia de la cima del yacimiento a la cima Distancia de la cima del yacimiento a la
del intervalo disparado (pies) del intervalo disparado (pies) cima del intervalo disparado (pies)

Figuras 1.12. Curvas de gasto crítico, (1) Conificación agua-aceite curvas A, B, C, D y

E; (2) Conificación gas-aceite curvas a, b, c, d y e.

Página 47 de 170
Chierici, Ciucci y Pizzi13, para determinar el gasto crítico presentan la correlación
siguiente:

Δρ wo k H
qow = 5.256 x10 −3 h 2 Ψ ( rDe , ε , δ w ) , ............................................ (1.41)
Bo μ o

Δρ go kH
qog = 5.256 x10−3 h 2 Ψ ( rDe , ε , δ g ) . ............................................. (1.42)
Bo μo

Donde:

qow = Gasto crítico de aceite en un sistema agua-aceite, m³/d

q og = Gasto crítico de aceite en un sistema gas-aceite, m³/d

h = Espesor impregnado de aceite, m

kH = Permeabilidad horizontal, md

μo = Viscosidad del aceite, cp

Bo = Factor del volumen del aceite, m³/ m³

Δρwo = ρw − ρo , diferencia de densidades agua-aceite , gr/cm³

Δ ρ go = ρ o − ρ g , diferencia de densidades gas-aceite, gr/cm³

Ó bien,

Δρ wo k H
qow = 4.92 x10 −5 h 2 Ψ ( rDe , ε , δ w ) , .............................................. (1.43)
Bo μ o

Δρ go kH
qog = 4.92 ×10−5 h2 Ψ ( rDe , ε , δ g ) , .............................................. (1.44)
Bo μo

Página 48 de 170
0.00307 Δρ k H hhwc
Z= , ....................................................................... (1.45)
μ o qo Bo

hp
ε= , .............................................................................................. (1.46)
h

hgc
δg = , ............................................................................................ (1.47)
h

hwc
δw = . ........................................................................................... (1.48)
h

Donde:

Bo = Factor del volumen del aceite, rb/STB

h = Espesor impregnado de aceite, pies

hp = Intervalo disparado, pies

h gc = Distancia del Cgo a la cima del Intervalo disparado, pies

hw c = Distancia del Cwo a la base del Intervalo disparado, pies

kV = Permeabilidad vertical, md

kH = Permeabilidad horizontal, md

qog = Gasto crítico de aceite en un sistema gas-aceite, STB/D

qow = Gasto crítico de aceite en un sistema agua-aceite, STB/D

re = Radio de drene, pies

reD = Radio de drene adimensional, adimensional

μo = Viscosidad del aceite, cp

Página 49 de 170
Δρwo = ρ w − ρo , diferencia de densidades agua-aceite , lb/ft³

Δ ρ go = ρ o − ρ g , diferencia de densidades gas-aceite, lb/ft³

ψg = Función adimensional de la fase gas, adimensional

ψw = Función adimensional de la fase agua, adimensional

re kV
5 ≤ reD = ≤ 80
h kH
hp
0 ≤ε = ≤ 0.75
h
hgc
0.07 ≤ δg = ≤ 0.9
h
h
0.07 ≤ δ w = wc ≤ 0.9
h

Las funciones del agua y del gas adimensionales son función de reD , ε y δ , y se

obtienen de las figuras 1.13 a 1.19.

0.200
reD = 5
0.100
0.080
0.060

ψ
0.040

0.020

0.010
0.008
0.006
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
ε
Figura 1.13. Funciones adimensionales para reD = 5.

Página 50 de 170
0.200
reD = 10
0.100
0.080
0.060

ψ
0.040

0.020

0.010
0.008
0.006
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
ε
Figura 1.14. Funciones adimensionales para reD = 10.

0.200
reD = 20
0.100
0.080
0.060

0.040
ψ
0.020