UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO

FACULTAD DE INGENIERÍA

Curvas de Declinación en
pozos de yacimientos tipo:
Tight Oil

TESIS
Que para obtener el título de
Ingeniero Petrolero

PRESENTA
Gabriel Dario Valencia Goujon

DIRECTOR DE TESIS
Ing. Israel Castro Herrera

Ciudad Universitaria, Cd. Mx., Junio de 2017

 ÍNDICE
Capítulo 1 Antecedentes

1.1 Propósitos, objetivos y limitaciones

1.2 Importancia del ACD (Análisis de Curvas de Declinación)
1.2.1 Método confiable DCA para estimar reservas 1P tight/oil shale
1.3 ¿Cómo ayudan las curvas de declinación en la Evaluación Económica?
1.4 Recursos no Convencionales
1.5 Fracturamiento Hidráulico

Capítulo 2 Geología de un yacimiento no convencional

2.1 Características y ¨ El porqué ¨

2.2 Yacimientos Convencionales no Convencionales
2.2.1 Yacimientos Convencionales
2.2.2 Yacimientos no Convencionales
2.3 Generación del Petróleo
2.4 Sistema Petrolero
2.4.1 Roca Generadora
2.4.2 Migración
2.4.2.1 Migración primaria
2.4.2.2 Migración secundaria
2.4.3 Roca Almacenadora
2.4.4 Roca Sello/Trampa
2.4.5 Sincronía
2.5 Sistema Petrolero no Convencional
2.6 Ejemplo de un yacimiento Convencional
2.6.1 Campo Fulmar
2.7 Ejemplo de un yacimiento no Convencional

Capítulo 3 Caracterización de yacimientos no convencionales

3.1 La importancia de la Caracterización Dinámica

3.2 Proceso de la Caracterización Dinámica y Simulación de un yacimiento convencional
3.2.1 Descripción del yacimiento convencional
3.2.2 Mecanismos de empuje en yacimientos convencionales
3.2.3 Realizar el Modelo Matemático
3.2.3.1 Deducción de la Ecuación de Continuidad para un yacimiento
convencional
3.2.3.2 Regímenes de Flujo de yacimientos convencionales
3.2.3.2.1 Régimen Estacionario
3.2.3.2.2 Régimen Pseudoestacionario
i
 3.2.3.2.3 Régimen Transitorio
3.2.3.3 Geometrías Flujo para yacimientos convencionales
3.2.3.3.1 Flujo Lineal
3.2.3.3.2 Flujo Radial
3.2.3.3.3 Flujo Esférico
3.2.3.3.4 Flujo Bilineal
3.2.4 Realizar el Modelo Numérico
3.2.4.1 Solución Línea Fuente
3.2.4.2 Simulación Numérica de Yacimientos
3.2.5 Realizar el Programa que demuestre el modelo numérico
3.2.6 Validar el modelo
3.2.7 Ajuste del Historial de producción
3.2.8 Predecir el comportamiento de yacimientos convencionales
3.3 Caracterización Dinámica de Yacimientos no Convencionales
3.3.1 Mecanismos de empuje en yacimientos no convencionales tipo Tight Oil
3.3.2 Regímenes de flujo en yacimientos no convencionales tipo Tight Oil
3.3.3 Geometrías de flujo de en yacimientos tipo Tight Oil
3.3.4 Propiedades de los fluidos de los yacimientos no convencionales

Capítulo 4 Curvas de declinación

4.1 Pequeña Introducción

4.2 Antecedentes de la curvas de declinación
4.3 Curvas de declinación
4.3.1 Declinación Exponencial
4.3.1.1 ¿Cómo calcular el coeficiente b de forma analítica)
4.3.2 Declinación Hiperbólica
4.3.3 Declinación Armónica
4.3.4 Ejemplo de curvas de declinación
4.4 Fetkovich Curvas de declinación
4.5 Usando la forma hiperbólica para aproximar el comportamiento de los yacimientos no
convencionales
4.6 Pseudopresión la homología entre yacimientos de gas y yacimientos de aceite

Capítulo 5 Curvas de declinación para pozos tipo Tight oil

5.1 Método de Doung

5.1.1 Como aplicar el método Duong
5.1.2 Análisis y discusión de los resultados del ejemplo usado para el método de
Duong
5.2 Método de Ilk para yacimientos de Gas
5.3 Método de Valko para yacimientos de Gas
5.4 Caso Práctico con un pozo del Bakken

ii
Conclusiones

Bibliografía

Anexos

Anexo A - Algebra Matricial, Método de Newton-Raphson y Matriz Jacobiano

Anexo B - Calculando el historial de producción del Pozo 16/26 A9 del Campo Alba

Anexo C - Producción del Pozo 22/11 N-5 Campo Nelson UKCS

Anexo D - Deducción de la ecuación de declinación Exponencial

Anexo E - Deducción de la ecuación de declinación Hiperbólica

Anexo F - Deducción del modelo de declinación de Duong

iii
 FIGURAS

IMÁGENES

Figura 1. Evaluación Económica

Figura 2. Fractura y apuntalante
Figura 3. Proceso de fracturamiento Hidráulico

Figura 4. Ventana de Generación del Petróleo

Figura 5. Procesos que dan lugar a la formación y acumulación del petróleo
Figura 6. Sistema petrolero Convencional
Figura 7. Migración Primaria y Secundaria
Figura 8. Columna Litoestratigráfica campo FULMAR
Figura 9. Subdivisión del yacimiento FULMAR en unidades productivas
Figura 10. Historia de sepultamiento y tabla de sincronía del campo FULMAR
Figura 11. Esquema de la sección transversal de la Dakota del Norte
Figura 12. Sección Transversal de la región del Bakken. Con sus formación y propiedades
Figura 13. Columna estratigráfica de la región del Bakken
Figura 14. Historia de sepultamiento y tabla de sincronía del campo Bakken

Figura 15 Régimen Estacionario

Figura 16 Régimen Pseudoestacionario
Figura 17 Régimen Transitorio
Figura 18 Flujo Lineal
Figura 19 Comportamiento del Flujo Lineal en variables adimensionales
Figura 20 Flujo Radial
Figura 21 Comportamiento del Flujo Radial en variables adimensionales y escala logarítmica
Figura 22 Flujo Esférico
Figura 23 Comportamiento del Flujo Esférico en variables adimensionales
Figura 24 Flujo Bilineal
Figura 25 Comportamiento del Flujo Bilineal en variables adimensionales con un Factor de
Conducción de 0.2 y 0.6
Figura 26 Comparación del comportamiento de las geometrías de Flujo
Figura 27 Como aproximar la derivada de la presión tomando en cuenta la dirección de la
derivada
Figura 28 Arreglo de las celdas en simulación numérica para hacer el sistema de ecuaciones
Figura 29 Tipos de malla que pueden usarse para simular matemáticamente un yacimiento
Figura 30 Arreglo de nodos centrados a la derecha y nodos distribuidos a la izquierda
Figura 31 Arreglo de nodos en un plano
Figura 32 Comportamiento de la RGA de un yacimiento convencional contra uno no convencional
Figura 33 Evolución del sistema de fracturas con el paso del tiempo
Figura 34 Historial de producción de un pozo tipo tight oil de la formación Meramec
Figura 35 Etapas de producción de un pozo horizontal hidráulicamente fracturado
Figura 36 Perfiles de presión simulados de un yacimiento
Figura 37 Gráfica diagnóstico para representar las geometrías de flujo tomada de una simulación

iv
Figura 38 Compresibilidad total del sistema y sus componentes
Figura 39 Caso de simulación: (a)esquema del modelo de la fractura, (b)propiedades del fluido
Figura 40 Gráficas (a) 𝑆0vs 𝑆 y (b) 𝑆0 vs 𝑆𝑆 producto de la simulación de una fracura

Figura 41 Historial de producción log q vs t del pozo 16/26 A9 Campo Alba

Figura 42 Comportamiento declinación exponencial
Figura 43 Sección de historial de producción de pozo 16/26 A9 con línea de tendencia usado
como ejemplo
Figura 44 Ejemplo de ajuste de curva de declinación en una sección del historial de producción
del pozo 16/26 A9
Figura 45 Comportamiento declinación hiperbólica
Figura 46 Comportamiento declinación armónica
Figura 47 Historial de producción del pozo 22/11 N-5 del campo Nelson en escala semilog con
curvas ajustadas
Figura 48 Q vs t en escala normal con curvas de declinación ajustadas
Figura 49 qo+qw y qo vs t en escala normal con curvas de declinación ajustadas
Figura 50 Efecto del radio de drene en la duración del flujo transitorio
Figura 51 Familia de curvas Tipo
𝑆𝑆
Figura 52 Gráfica de 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 𝑆 de los pozos HMFW del Bakken
𝑆𝑆
𝑆𝑆
Figura 53 Gráfica de 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 𝑆 de los pozos Whitman, Boucher y Charlson
𝑆𝑆
Figura 54 Historial de producción del pozo Charlson Federal 14X-35H para usarse como
ejemplo
𝑆𝑆
Figura 55 Gráfica 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 𝑆 del pozo Charlson Federal 14X-35H para encontrar a y m
𝑆𝑆
Figura 56 Gráfica 𝑆𝑆 𝑆𝑆 𝑆 (𝑆, 𝑆) para encontrar los parámetros 𝑆1 𝑆 𝑆∞. En este caso se
ajustaron 4 líneas de tendencia para ver que pronóstico se ajusta mejor al perfil real de
producción
Figura 57 Cuatro pronósticos de producción comparados contra el perfil de producción original.
Se usa escala semilogarítmica para apreciar más fácil el ajuste
Figura 58 Estimación de reservas usando los cuatro pronósticos de producción determinados
con el método de Duong
Figura 59 Pronóstico de producción del ejemplo para el método de Duong usando 𝑆 =
𝑆. 𝑆𝑆 𝑆 𝑆 = 𝑆. 𝑆𝑆. Dan pronósticos con un ajuste muy certero al perfil real
Figura 60

Tabla 1. Clasificación granulométrica de Rocas Sedimentarias Clásticas

Tabla 2. Relación porosidad permeabilidad para rocas sedimentarias
Tabla 3. Sistema Petrolero Campo FULMAR
Tabla 4. Porosidades, permeabilidades y volúmenes de las unidades productivas FULMAR
Tabla 5. Sistema petrolero del Campo Bakken

Tabla 6 Aproximaciones de la Integral Exponencial

Tabla 7 Valores de las Variables Adimensionales para diferentes geometrías de flujo
Tabla 8 Principales métodos parar resolver matrices de Sistemas de Ecuaciones Lineales
Tabla 9 Etapas de producción de MFHW usando el RGA para identificarlas
v
Tabla 10 Descomposición de la compresibilidad total

Tabla 12 Curvas de declinación ajustadas al pozo 22/11 N-5

Tabla 13 Q vs t en escala normal con curvas de declinación ajustadas
Tabla 14 qo+qw y qo vs t en escala normal con curvas de declinación ajustadas

Tabla 15 Base del modelo de ecuación diferencial para los nuevos modelos de curva de
declinación
Tabla 16 Parámetros 𝑆𝑆 𝑆 𝑆∞ para cada línea de tendencia de la fig. 56. Porcentaje de error
para cada estimación
Tabla 17 Parámetros 𝑆𝑆 𝑆 𝑆∞ para cada línea de tendencia de la fig. 59. Porcentaje de error
para cada estimación

vi
 CAPITULO 1
ANTECEDENTES

1.1 Propósitos, objetivos y limitaciones

1.1.1 Propósito

Tiene como objetivo documentar los diferentes modelos de curvas de declinación que pueden
ser útiles para pronosticar el gasto de pozos que producen en cuencas de tight oil.

1.1.2 Objetivos

Hacer una investigación de los modelos recientes que se han desarrollado para modelar la
producción futura de los pozos productores en cuencas de tipo tight oil. La producción de
yacimientos de muy baja permeabilidad comenzó como el shale gas y después se empezó a
implementar el fracturamiento hidráulico en para pozos de tight oil. Con el fin de estimular las
formaciones y aumentar la permeabilidad. Arps fue una de las primeras personas en documentar
estudios enfocados en la declinación de la producción de los pozos. Él desarrolló las primeras
curvas de declinación basándose en 3 modelos: armónico, exponencial e hiperbólico. Desde
entonces ha pasado mucho tiempo que no se han desarrollado nuevo modelos. Hasta el 2008
que se comenzó a explotar los yacimientos de shale gas. Entre los años del 2010-2012 que se
comenzó a desarrollar nuevos modelos más precisos para los pozos de tight oil y shale gas.

En esta tesis se explicarán los métodos y nuevos modelos que se han desarrollado hasta el
momento, como el de Ilk, Valko y Doung. Usados para hacer pronósticos de producción
especializados para tight oil y shale gas.

Otro de los objetivos es analizar las diferencias que presentan los yacimientos convencionales y
no convencionales en diferentes etapas del estudio de una cuenca. Como son el sistema
petrolero o la caracterización dinámica de yacimiento.

1.1.3 Limitaciones

La escasa información a la cual se puede acceder gratuitamente, la cotidianidad del tema y que
no se cuentan con historiales de producción muy extensos. Ya que se empezó a producir en
shale gas en el 2008. Los estudios de curvas de declinación usan historiales de producción de
varios años para hacer una estimación muy confiable.

- 1-
1.2 La Importancia de las Curvas de Declinación

En la Ingeniería, como en muchas otras ciencias exactas se tienen bases matemáticas y se usan
otros conocimientos específicos, dependiendo del tipo de ingeniería. Pero éstos no se pueden
deslindar de la rama administrativa y de planeación, ya que la ingeniería trabaja desarrollando
proyectos; solucionando problemas sociales. Los cuales deben de ser evaluados antes de ser
ejecutados para asegurar que van a ser viables y se podrá alcanzar las metas finales (solucionar
el problema).

Por esto, surge el interés en las Curvas de Declinación. Éstas son usadas para estimar
producciones futuras. Con ellas se elaboran pronósticos de producción para predecir las
ganancias que tendrá el proyecto. Al conocer las inversiones y los costos que estarán presentes
en la vida del proyecto, el siguiente paso sería hacer una evaluación económica para saber si
nuestro proyecto será viable o no, así como conocer el rendimiento aproximado que
obtendríamos al final de un periodo determinado o la tasa interna de retorno (TIR).

Los estudios de curvas de declinación requieren de menor información que cualquier otro estudio
de yacimiento. Tan sólo con el historial de producción del pozo podemos empezar a usar
métodos para medir la declinación y modelar su comportamiento futuro. Una vez que se conoce
el comportamiento de la declinación se pueden hacer predicciones muy confiables. Por ejemplo
para la SNY (Simulación Numérica de Yacimientos) necesitamos modelos geológicos, historiales
de producción, desarrollo de un modelo matemático y un simulador, análisis PVT, etc.

Cabe mencionar que otras de las aplicaciones de las curvas de declinación son: incorporar
reservas, estimar producciones, analizar el potencial de una cuenca, ayudar en la evaluación
económica de un proyecto, auxiliar en el diseño de instalaciones superficiales estimando
volúmenes o auxiliando en la toma de decisiones como el implementará un sistema de
recuperación mejorada.

Algo muy importante que se debe mencionar es que el método de curvas de declinación (Decline
Curve Analysis, DCA en inglés) sólo puede ser usado cuando el yacimiento se encuentra bajo
saturado. Dicho en otras palabras, cuando los pozos están produciendo con una presión mayor a
la presión de burbuja (Pyacimiento>Pb). En caso de que se llegara la Presión de Burbuja Pb, no se
podrá usar este método por qué las estimaciones serán incorrectas. En este caso, se debe usar
simulación numérica de yacimientos utilizando el modelo que mejor represente la cuenca.

1.3 ¿Cómo ayudan las curvas de declinación en la Evaluación Económica?

Se evalúa un proyecto económicamente para saber si será factible. Pero cómo se conocerán las
ganancias generadas al desarrollar un campo teniendo un 90% certeza que suceda. Es por esto
que son de gran ayuda las curvas de declinación ya que se pueden usar para hacer un análisis
de la producción (pronóstico de producción) de los pozos que serán perforados. Lo que sigue
sería hacer un estimado o pronóstico del precio del barril en el tiempo de 1 a 5 años. Se sugiere

- 2-
basarse en los valores de la oferta y demanda para estimar el precio futuro. Con esto se puede
estimar las ganancias que se tendrían por la producción de aceite, gas y condensado.

Al generar un pronóstico de producción se pueden situar en un escenario pesimista para hacerlo

más confiable. Se pueden usar campos análogos para obtener información (producción inicial,
tasa de declinación, tipo de yacimiento) y correr una simulación del escenario. Requieren de
menor tiempo a diferencia de la corrida de un modelo de simulación numérica de yacimientos. Si
se sabe usarlas bien se puede obtener un buen pronóstico en el periodo de evaluación.

Las curvas de declinación nos darán producciones de aceite en función del tiempo, suponiendo
que el pozo no dejase de producir en ningún día. Esto junto con el precio del barril y de la divisa
nos puede ayudar a estimar los ingresos del proyecto. Es cierto que existe una alta incertidumbre
por estar jugando con tres variables que cambian en el tiempo, pero las curvas de declinación
dan resultados muy confiables conociendo los parámetros correctos en el caso de un pozo
exploratorio.

Fig. 1 Evaluación Económica

1.4 Recursos no Convencionales

Los recursos en el mundo se tienen clasificados de muchas formas pero en el argot petrolero
existen los convencionales y no convencionales. El Aceite de lutitas (Tight Oil en inglés) es un
recurso no convencional ya un pozo perforado en una cuenca de este tipo no fluiría naturalmente
debido a su baja permeabilidad. Aunque la formación presente grandes volúmenes de

- 3-
hidrocarburos el aceite se encuentra atrapado debido a la muy baja permeabilidad. La
permeabilidad no es más que el área de los poros interconectados por los que es posible el paso
de cualquier tipo de fluido. Como en este caso fluirían los hidrocarburos de la roca generadora al
paso de migración. Los recurso no convencional son más explotado en las Cuencas
estadounidense como: Eagle Ford, Utica, Niobrara, Bakken o en Canadá en la Cuenca Anticosti
o en la cuencas marítimas centrales. En México por citar un ejemplo es en Chicontepec. Un
ejemplo en África es el campo Tango. El campo Tango contiene crudo extrapesado (12 °API)
con una viscosidad promedio de 2,000 cp.

Para superar este obstáculo se ha desarrollado la perforación horizontal seguida de varias

etapas de fracturamiento hidráulico. Al usar la perforación direccional se minimiza el área
utilizada por las instalaciones superficiales. Ya que de una macro pera se puede recibir la
producción de 8 pozos que drenan al yacimiento en diferentes ubicaciones.

En los últimos años se han desarrollado nuevas técnicas de explotación cómo es el

fracturamiento hidráulico, la recuperación secundaria y mejorada, SAGD [Steam-Assited Gravity
Drainage], combustión in situ, etc. Cuyo objetivo es extraer aquellos recursos que antes no eran
convencionales. Dado a que no existían los avances tecnológicos necesarios para extraerlos.
Éstos han logrado dar el siguiente paso en la carrera energética, poniendo en un nuevo lugar
aquello que antes no era tomado en cuenta por ser difícil de alcanzar. Pero al encontrar nuevas
formas de explotar lo anteriormente no explotable, se tienen que hacer análisis nuevos. El más
claro ejemplo de esto son los pozos que producen a través de fracturas artificialmente hechas.
En las cuales se tiene que meter un material o fluido apuntalante antes de que se alcance la
presión de cierre de la fractura, para que con ayuda de éste material especial se creé una
permeabilidad artificial secundaria por la cual los fluidos que no era posible que fluyeran; dada a
la reducida permeabilidad original de la roca, empiecen a circular por la red de fracturas con el
apuntalante. Lo que origina que al principio se tenga una presión muy alta pero en el momento
en el que la matriz deja de aportar fluido a la fractura, la presión decae notablemente hasta que
es tan pequeña que el flujo por las fracturas se detiene. En este caso se tendrá que hacer
fracturamiento hidráulico nuevamente.

Las ondas de presión en un yacimiento siempre van a difundirse en toda su extensión, pero la
rapidez con la que viajen variará dependiendo de la permeabilidad de la formación. En las
formaciones de tight oil las ondas tardarán mucho más tiempo en alcanzar las fronteras del
yacimiento. Haciendo muy complicado una caracterización dinámica usando una prueba de
presión. Este es un ejemplo de la diferencia del comportamiento de los yacimientos y a su vez la
forma diferente de analizarlos.

En esta tesis se hablará sobre las nuevas curvas de declinación que se han desarrollado
diferentes a las ya existentes (armónica, exponencial e hiperbólica) para hacer estimaciones en
pozos que produzcan aceite de formaciones poco permeables que se les haya hecho un
fracturamiento hidráulico.

- 4-
Para lo cual se empezara trabando con los estudios realizados con anterioridad en la tesis de
Maestría de Linnea Lund (2014) y todas las investigaciones en la rama de Shale gas. En las
cuales se puede rescatar que son comportamientos similares pero contraponiendo que la
viscosidad del gas es aproximadamente una cincomilésima parte de la viscosidad de un aceite
ligero. Que sería el tipo de crudo que se produciría en una cuenta de tight oil.

Viscosidad de una mezcla de gas a las condiciones de P= 1 atm, T=93.3 °C (P= 14.7psi, T=200
°F)
μg = 0.01247 cp
Viscosidad de la misma mezcla de gases a condiciones de yacimiento: P= 142 atm, T=104 °C
(P= 2,100 psi, T=220 °F)
μg = 0.01946 cp

Viscosidad de un aceite ligero API° > 31°

μo = 100 cp
Por lo tanto,
μo 100 cp
= = 5138.7
μg 0.01946 cp

Tomando en cuenta éste análisis se puede ver que algunas de las propiedades analizadas
cambian; como es la viscosidad, pero la permeabilidad en la fractura queda constante. He de
aquí que se pueden tomar algunas de las investigaciones hechas para shale gas para auxiliar en
la búsqueda de la nueva curva de declinación.

1.5 Fracturamiento Hidráulico

Es el proceso por el cual se transmite presión usando gas o fluido a la roca hasta crear fisuras o
fracturas. Se usa para fracturar rocas de baja permeabilidad que contengan una saturación
económicamente considerable de hidrocarburos. Tiene la finalidad de aumentar la permeabilidad
original κ, debido que la [𝑆efectiva > 𝑆total] es del orden de micro a nano Darcies ( 1x10-6 –
1x10- 9 D ). Éste es un proceso diseñado de tal forma que las fracturas alcancen una longitud
planeada y el radio de drene aumente. Se bombea el apuntalante junto con el fluido que
transmite la presión de tal forma que al iniciarse la fisura éste material entre en ella y cuando se
libere la
presión el apuntalante impida que la fisura se cierre en su totalidad. El apuntalante tiene que ser
lo suficientemente resistente para no deformarse ya que estará sujeto los esfuerzos efectivos
[sobrecarga, horizontal y vertical].

Los apuntalantes son materiales diseñados para que tengan un diámetro y esfericidad
específicos, así al estar en contacto uno contra otro; como si fueran granos de arena, generen un

- 5-
conducto artificial por el cual circulen los hidrocarburos(Fig 1.). Mejor conocido como porosidad
secundaria κartificial . Así se aumenta la conductividad y los fluidos pasen con más facilidad.

Fig. 2 Fractura y características del apuntalante

Se realiza el fracturamiento
hidráulico para estimular el
pozo, de tal forma que a
través de inducir fracturas se
aumente la permeabilidad y la
producción. El fracturamiento
se realiza tanto en pozos
verticales como en pozos
direccionales u horizontales.
La ejecución de un programa
de fracturamiento hidráulico
se lleva a cabo después de
realizar la terminación del
pozo y tomar registros.

Se escoge las zonas con

mayor saturación para ser
fracturadas. Se escoge el
radio de drene y por ende la
longitud deseada de la
fractura. Se diseña el fluido
fracturante (densidad,

Fig. 3 Proceso de fracturamiento hidráulico

- 6-
viscosidad, base acuosa o aceite, apuntalante). Se diseña el esquema de bombeo en base a la
prueba de goteo y los valores de a presión crítica, presión de cierre y presión de fractura (Fig. 3).

- 7-
 CAPÍTULO 2

GEOLOGÍA DE YACIMIENTO NO CONVENCIONAL

2.1 Características y el ¨Porqué¨

Empezando por los conceptos más básicos se explicará como funciona un yacimiento, con el fin
de que en los siguientes capítulos se entienda el porque se utiliza cada tipo de curva de
declinación.

Se plantearán los conceptos de: yacimiento convencional/no convencional y del sistema

petrolero. Ya que son de suma importancia para analizar los yacimientos y entender cómo
interactúan los fluidos en el subsuelo. También es la base para desarrollar cualquier trabajo que
tenga que ver con Geología de Exploración y Explotación. Como para la planeación de un
programa de recuperación mejorada, la ubicación de un pozo de inyección, determinación de
reservas, etc.

2.2 Yacimiento Convencional /no Convencional

2.2.1 Yacimientos Convencionales

Un yacimiento convencional tiene que cumplir con todos los elementos del sistema petrolero, su
roca almacenadora tiene la permeabilidad suficiente para que fluyan los hidrocarburos y los
fluidos tienen buenas propiedades de viscosidad y densidad. Es aquel que después de
perforarse un pozo tiene la capacidad de llevar a superficie volúmenes rentables de producción
con su energía natural (presión inicial del yacimiento). Estas son en pocas palabras las
características más representativas de los yacimientos convencionales.

Los mecanismos de empuje que imperan en los yacimientos convencionales son: empuje por
expansión del sistema roca-fluido, empuje por expansión del gas disuelto en el aceite, empuje
por segregación gravitacional, empuje por la admisión de agua del acuífero, empuje por la
expansión del casquete de gas. A diferencia de otro tipo de yacimientos en los que podemos
encontrar otro tipo de mecanismos de empuje.

2.2.2 Yacimientos no Convencionales

Hay varias formas de definir un yacimiento no convencional. La definición primordial es aquel

prospecto que contenga hidrocarburos atrapados en una cuenca por ejemplo una trampa
[estratigráfica, estructural], pero que no sean posibles extraer dada la falta de tecnología y/o los

- 8-
métodos necesarios. Se puede ejemplificar esta definición con el gas en lutitas [shale gas] que
no fue hasta que se implementaron las técnicas de fracturamiento hidráulico que paso de ser de
un recurso no convencional a uno explotable.

Los yacimientos no convencionales se diferencian de dos formas. Ya sea porque sus

características geológicas no le permiten a la roca almacenadora expulsar suficiente petróleo
para que el pozo sea rentable y la segunda, que el fluido es tan viscoso que no puede fluir
naturalmente por la tubería de producción. Éstos son unos ejemplos, pero a continuación
profundizaremos más en lo que son los yacimientos no convencionales.

Se pueden encontrar yacimientos no convencionales dependiendo de sus características de dos

formas: geológica y química. Desde el punto de vista geológico son debido a las características
de la de la roca como: baja permeabilidad, escasa porosidad efectiva o altos esfuerzos efectivos
sobre la roca. Viéndose desde el punto de vista químico, son debido a las propiedades del fluido
como: alta viscosidad, baja relación gas aceite, baja presión o alto contenido de Azufre o Dióxido
de Carbono, que corroe la tubería de explotación (Leveille,2012). Algunos ejemplos de
yacimientos no convencionales son: Three Forks [Tres Tenedores] en el Bakken, Campo Jilin en
China y la Cuenca de Cauvery en India.

2.3 Generación del Petróleo

El proceso por el cual se generan los hidrocarburos es muy complejo, en él se involucran varios
subprocesos como: muchas reacciones químicas, fenómenos físicos y procesos biológicos de
descomposición.

Las reacciones químicas están presentes en todas partes. Pueden ser orgánicas o inorgánicas.
Las reacciones más básicas son: síntesis, descomposición, sustitución y doble sustitución. Con
ellas se forman nuevos compuestos a partir de los reactivos. Los productos de las reacciones
pueden ser moléculas muy sencillas como el agua hasta moléculas muy complejas como la
proteína del glucagón.

Como se sabe todos los seres vivos están conformados por células, a su vez estas están
conformadas por ácidos nucleicos, proteínas, lípidos, carbohidratos y minerales. Las células
empiezan a degradarse al morir los seres vivos y aquí empieza a descomponerse todos los
órganos que conforma a las células como la membrana celular que es una bicapa fosolípica.
Ésta al descomponerse genera muchos lípidos y proteínas (Biología, 2008) que
subsecuentemente se descompondrán en moléculas más sencillas como los hidrocarburos y
otros compuestos. Por citar otro ejemplo hablemos de las proteínas que forman las fibras
musculares y muchas otras estructuras. Las proteínas se descomponen en aminoácidos. Los
cuales al ir aumentando la profundidad debido al enterramiento son sometidos cambios de
temperaturas y presiones mayores que rompen los enlaces de las proteínas y dejan como
residuo muchos aminoácidos. Debido a la baja cantidad de oxigeno las sustancias no pueden
oxidarse.

- 9-
Las reacciones de degradación debido a las bacterias se dan en los primeros metros de
sepultamiento de la materia orgánica. Es aquí donde la actividad bacteriana es extremadamente
común. Las bacterias a través de enzimas ayudan a descomponer las estructuras celulares
dejando atrás estructuras simples como aminoácidos, lípidos. Ya que lo que ellas buscan son
carbohidratos para alimentarse y obtener energía.

Los fenómenos físicos más importantes para la generación del petróleo se generan al descender
por la corteza terrestre. La conducción térmica es una propiedad de la materia de mucha
importancia para que se generen los gradientes geotérmicos. Como se sabe el gradiente
geotérmico promedio de la Corteza Terrestre es 1 °C por cada 33 metros que se descienden. El
aumento de temperatura provoca que la descomposición de la materia aumente dejando a su
paso biopartículas cada vez más sencillas que después de segregarse gravitacionalmente
dejarán una fracción de lo que conocemos como hidrocarburos. Este tipo de degradación a
diferencia de la producida por bacterias es debido a que los enlaces de las biopartículas son
rotos fácilmente con un aumento de la temperatura y presión. Otro fenómeno que se genera con
el sepultamiento es el aumento de la presión. Cuando existe una interconexión entre los poros de
los sedimentos, el gradiente mínimo de presión que se encuentra es el hidrostático suponiendo
que los poros está saturada por agua.

Estos son los procesos por los que pasa la materia orgánica que es sepultada en un ambiente
sedimentario terrestre o lacustre. La diferencia entre los ambientes sedimentarios es que en el
terrestre la materia orgánica está en contacto con la atmósfera y en el marino o lacustre es un
ambiente anaerobio.

La materia orgánica después de pasar por varios procesos se convierte en kerógeno. El

kerógeno es un término usado para designar la parte de las rocas sedimentarias que no es
soluble en medios acuosos ni en solventes alcalinos. El antecesor del kerógeno es ¨humin¨ y la
diferencia entre éstos, es que hay una fracción hidrolizada en el humin. El kerógeno es la fuente
de carbono orgánico más abundante en la Tierra. El bitumen es la parte de las rocas
sedimentarias que es soluble en solventes orgánicos.

Los tipos de kerógenos que se pueden encontrar son el sapropélico que es producto de la
descomposición y polimerización de la materia algacea y herbácea. Es depositada en ambientes
con bajo contenido de oxígenos atmosférico, normalmente ambientes acuáticos. El kerógeno tipo
húmico que es producto de la descomposición de plantas terrestres superiores. Es depositado en
ambientes con mucho oxígeno atmosférico, como ambientes terrígenos sedimentarios.

El origen de los hidrocarburos puede ser biogenético [origen biológico] o termogenético [origen
debido al aumento de la temperatura y descomposición del kerógeno]. Debido a esto es
importante el gradiente geotérmico y el gradiente de presión. Los factores que comienzan a
convertir el kerógeno en petróleo después haber pasado por los procesos: degradación
bacteriana, sepultamiento, compactación, son el aumento de temperatura que indica el principio
de la Diagénesis y continuando la Catagénesis. En la Diagénesis no se alcanza la temperatura
suficiente para que se generen altos rendimientos de aceite por cada gramo de Carbono total en

- 10 -
la roca. En la Catagénesis se alcanza la temperatura óptima para que la reacción se dé por
millones de años. A está condición especial de Presión y Temperatura se le conoce como
Ventana de generación del Aceite y también la Ventana de generación del Gas que se da en la
Metagénesis. Cuando se han alcanzado temperaturas más altas. La ventana de generación del
aceite está entre los 60 °C a los 115 °C (Tissot, 1984) como se ve en la Fig. 4 .

Fig. 4 Ventana Generación del Petróleo a diferentes condiciones, PZOF -Principal zone of oil
generation-, ZGF-Zone of gas generation by cracking - modificado de …(Tissot, 1984, pp 180)

2.4 Sistema Petrolero

Es un esquema creado para organizar las partes que conforman un yacimiento y analizar cada
una de tal forma que podamos planear eficientemente la explotación del yacimiento. También es
un sistema que ocurre naturalmente, que incluye todos los elementos y procesos geológicamente
esenciales para que un yacimiento de hidrocarburos exista en la naturaleza.

Éste es un concepto creado en la ingeniería para representar geológicamente, se generan los

hidrocarburos, analizar el sistema en el cual están atrapados los hidrocarburos y por el cual se
desplazaran los fluidos ya sea migrando o en el pozo. Aquí se analizan muchos factores que
están contenidos dentro de los elementos del sistema petrolero.

Los elementos que conforman un sistema petrolero son:

- 11 -
  Roca Generadora
 Migración
 Roca Almacenadora
 Roca Sello/ Trampa
 Sincronía

Fig. 5 Sistema Petrolero Convencional

2.4.1 Roca generadora

Es el primer elemento del sistema petrolero y su importancia radica en que aquí se generan los
hidrocarburos y dependiendo de las características de ésta, determinará el estado de los
hidrocarburos [líquidos, gaseosos o combinación de los ambos] y muchas de sus propiedades
[viscosidad, RGA, grados °API o gravedad específica, presión de burbuja y roció, % H2S, CO2,
N2, etc.]. En él se analizan: tipo de materia orgánica contenida, COT, edad geológica,
profundidad, tipo de roca y su clasificación granulométrica, volumen de la roca, condiciones
térmicas, rendimiento. A continuación se describirá brevemente cada una. En la Fig. 6 se puede
ver los conceptos básicos que dan lugar al sistema petrolero y como se relaciona entre sí.

- 12 -
 Fig. 6 Procesos que dan lugar a la formación y acumulación del petróleo

Tipo de materia orgánica contenida: Se refiere al tipo de kerógeno precursor de los

hidrocarburos [Tipo I, Tipo II o Tipo III]. Con él se va a saber qué tipo de composición química
tiene el aceite y gas. Aceite ligero, intermedio, pesado o extrapesado. Predecir los valores de
viscosidad que puede tener y la RGA o gas disuelto en el aceite que contenga el hidrocarburo.

COT: Carbono Orgánico Total es la fracción de carbono de origen biológico que se encuentra en
una roca entre el peso de roca. Esta materia será la que después de ser sepultada, compacta y
sometida procesos térmicos dará origen al kerógeno y finalmente a los hidrocarburos. Se mide
en porcentaje sobre cada gramo de matriz de la roca.

Edad Geológica: Es el Tiempo geológico en el cual fueron depositados los sedimentos que dan
origen a la roca generadora. Se puede indagar de diferentes maneras. Con un núcleo tomado del
pozo podríamos hacer una prueba de C14. Pero en geología de campo se puede lograr tomando
una muestra de la roca y haciéndole una prueba de C14. Otra forma de conocer la Edad es con la
paleontología, analizando los microorganismos fósiles vistos en el microscopio.

Profundidad: Conocer la profundidad a la que se encuentra nos ayudará a saber junto con el
gradiente geotérmico y de presión a la que estará sujeta la roca generadora o cualquier otro
elemento del sistema petrolero.

Clasificación de la roca y granulometría: Saber si la roca generadora se divide en tres grandes

familias: ígneas, sedimentarias y metamórficas y una subclasificación por ejemplo dentro de las
sedimentarias: dentríticas o evaporíticas, da la pauta para determinar en que ambiente de
depósito o sedimentario se formó el estrato o si es una intrusión de un dique batolítico u otra
cosa. Esto sirve para conocer las propiedades de porosidad y permeabilidad que pudiesen tener
la roca debido a los minerales que forman la matriz o hacer una analogía con rocas de
mineralogía similar pero de diferentes ubicaciones. Al saber el tipo de roca con la que se trabaja
se puede indagar la dureza y otras propiedades. Conocer la mineralogía de la roca ayudara al
perforador a diseñar de manera correcta el plan de perforación y por ende la base del fluido de

- 13 -
perforación que mejor convenga. En la Tabla 1 se ve una clasificación simple de las rocas
sedimentarias.

Clasificación de suelos (Attenberg)

Dimensiones de la Attenberg- Nombre de la

partícula elemental Sistema Roca
(mm) Internacional
< 0.001 Arcilla
Arcilla
< 0.002 Limo Fino
0.005
0.01 Limo Limo Medio
0.02 Limo Grueso
0.05
0.1 Arena Fina
Arena Fina
0.25
0.2 Arena Media
0.5
1 Arena Gruesa
2 Arena Gruesa
3 Grava Fina
5 Grava
10 Grava Fina
20
Grava Gruesa Grava
> 20.0
Gruesa/Piedras
Tabla 1 Clasificación Granulométrica de las Rocas Sedimentarias Clásticas

Volumen de la roca: Sirve para estimar la cantidad de volumen de fluido que podría estar
contenido en la cuenca. Se tiene que conocer el la porosidad y saturaciones basado en estudios
de núcleos para calcularlo. De aquí se puede estimar el volumen original de hidrocarburos in-situ
y asignando un factor de recuperación se estiman las reservas. Se puede calcular de diferentes
maneras. Haciendo modelos de figuras geométricas simples como prisma que después se
doblan para simular los plegamientos o usando elemento finito para calcular el volumen
aproximado.

Condiciones térmicas: Se conoce al perforar si se mide el gradiente geotérmico del pozo. Con él
se calculara la temperatura que se encuentra nuestra roca generadora y se podrá estimar a
mayores profundidades la temperatura o por analogía saber aproximadamente la temperatura
que tendrá cierta profundidad si nos desplazáramos lateralmente la ubicación Es muy importante

- 14 -
conocer la Temperatura y la mecánica del suelo para saber si la Roca generadora alcanzo la
madurez ideal y si estuvo el tiempo indicado en la ventana de generación del aceite.

Rendimiento: Es el volumen de aceite [mg] producido por cada gramo de carbono orgánico total
en un volumen unitario de roca cuando está a alcanzado las condiciones de la ventana de
generación de aceite. Este volumen es producido en una unidad de tiempo. Con él se puede
aproximar el volumen de hidrocarburos producidos en la cuenca a través de los millones de años
en los que se mantuvo a la temperatura de la ventana de generación del aceite. Cómo se puede
observar en la Fig. 4.

2.4.2 Migración

Proceso mediante el cual los hidrocarburos se mueven de una ubicación a otra haciendo uso de
los estratos de mayor permeabilidad hasta no poder fluir más debido a una trampa. Los tipos de
migración son: primaria y secundaria.

El movimiento se produce debido a la fuerza de flotación. Ésta empuja a los fluidos de zonas
mayor potencial a de menor energía potencial [presión] o hacia contenedores donde el fluido ya
no puede moverse y queda atrapado. Todo esto lo se puede imaginar dentro de la corteza
terrestre para hacer una idea de lo que pasa.

2.4.2.1 Migración primaria

Es el movimiento de hidrocarburos a través de los poros interconectados dentro de la roca

generadora. Cuando el kerógeno se transforma en petróleo y comienza a moverse por la
diferencia de densidades entre el agua que moja a la roca y la densidad del petróleo con gas
disuelto. Como en la Fig. 7.

La migración primaría se origina dentro de la roca generadora cuando ésta se compacta y

alcanza la ventana de generación de hidrocarburos, el kerógeno se transforma en gas o líquido
y estos comienzan a moverse en la roca generadora. (Tissot, 1984, pp.294)

2.4.2.2 Migración secundaria

Se lleva a cabo a lo largo de la roca almacenadora, hasta que los fluidos hacen contacto con la
trampa e impide su movimiento, o escape a la superficie como en la Fig. 7.

La migración secundaria se origina cuando la HC´s se desplazan dentro de la tierra por

formaciones permeables que no cuentan con un sello para que impida que se sigan moviendo.
Así continuarán su recorrido hasta que queden atrapados en una trampa. (Tissot, 1984, pp.295)

- 15 -
La migración también puede ser lateral y vertical. La migración lateral son los desplazamientos
de fluidos en el interior de una formación de la misma edad geológica (Elia, 2013). La migración
vertical o transversal son desplazamientos de forma perpendicular a los límites
cronoestratigráficos; a grosso modo, los fluidos se mueven a formaciones de edades goelogícas
diferentes.

También existe el término de desmigración que se atribuye a los hidrocarburos que se escapan
de las trampas hacia otros lugares después de alcanzar el volumen máximo que pueden
contener la trampa.

Fig. 7 Migración Primaria y Secundaria

2.4.3 Roca Almacenadora

Sus características primordiales son: permeabilidad y porosidad. Mientras los hidrocarburos

migran secundariamente quedan almacenados en esta roca y empiezan a acumularse a causa
de un obstáculo que no les permite seguir viajando por el canal de alta permeabilidad. A esté
obstáculo que es más una estructura se le conoce como ¨Trampa¨. De la cual se hablará más
adelante.

El almacén debe de ser poroso y los poros deben de estar interconectados para que puedan fluir
los hidrocarburos. Por ejemplo se puede tener una pumita que es una roca muy porosa pero es
una mala roca almacenadora ya que sus poros no están interconectados.

Las gargantas de poro deben de superar el mínimo requerido ya que los fluidos se quedarán
atrapados debido a las fuerzas capilares y viscosas. Las rocas almacenadoras pueden ser de
origen detrítico o evaporíticas [precitados, ambiente marino]. Dentro de las rocas detríticas o
clásticas se puede encontrar las areniscas, lutitas, limonitas, grauvacas [graywackes],
conglomerados. Algunos como las areniscas bien clasificadas u ordenadas pueden alcanzar
valores de porosidad de 15-25% que resultaran en excelentes yacimientos. Pero también se
puede encontrar muy buenas rocas almacenadoras de origen evaporítico como los carbonatos.
Que tienen una porosidad primaria muy pequeña por ser precipitaciones pero la porosidad

- 16 -
secundaria debido a los sistemas de fracturas los hacen los mejores yacimientos alrededor del
mundo. Se puede encontrar una relación lineal entre la porosidad y la permeabilidad en las rocas
clásticas.

Porosidad % Valoración Permeabilidad mD

0-5 no se toma en cuenta -
5-10 Pobre -
10-15 Suficiente 10
15-25 Buena 10 - 100
20-25 Muy buena 100 - 1000
Tabla 2 Relación porosidad permeabilidad para rocas sedimentarias

2.4.4 Roca Sello/ Trampa

La Roca Sello debe ser un estrato superior a la roca almacenadora que tenga una permeabilidad
tan baja que sea considerada impermeable. Esto generará una obstrucción que impedirá que el
aceite siga migrando verticalmente u horizontalmente, así comenzará su acumulación en la
cuenca o roca almacenadora.

El Sello puede tener una geometría parecida a la de un balde invertido el aceite se queda
atorado en la parte con forma de sello como cuando flota una burbuja de aire en el agua y hay
algo que le impide llegar a la superficie. El Sello es efectivo hasta el plano de derrame en el cual
los hidrocarburos comenzarán a escapar o migrar. El plano de derrame se forma con los límites
verticales inferiores de la roca sello. Después de que el sello se ha llenado totalmente los
hidrocarburos de desbordan de la trampa y comienza a migrar a otra parte.

La trampa en el argot petrolero se refiera a la estructura o forma que tiene la roca sello. Su
característica principal es su forma convexa hacia el centro de la Tierra que se le atribuye una
forma de contenedor. Se pueden clasificar las trampas de dos formas: Trampa Estratigráfica,
Trampas Estructurales y Trampas Mixtas o Combinadas. Las trampas Estratigráficas se generan
por el cambio de facies o eventos de depositación [regresiones o transgresiones]. Las Trampas
Estructurales son generadas por tectonismo, plegamiento de estratos que forman los
anticlinales, por movimientos de la sal alóctona que generan Domos salinos o también por planos
de falla que ponen en contacto estrato permeable con impermeables cuando se desplazan
verticalmente. Para finalizar se pueden encontrar trampas mixtas o combinadas en las que se
presente tanto trampas estructurales como estratigráficas.

Se pueden encontrar un tercer tipo de trampa que son generadas por la hidrodinámica del agua
que moja a la roca cuando esta se desplaza en dirección opuesta a la fuerza de flotación del
petróleo. Cuando la fuerza hidrodinámica es más grande que la fuerza de flotación, ésta impedirá
que las gotas de aceite se muevan en sentido de la flotación, dicho en otras palabras que
migren.

- 17 -
Las trampas deben de tener sellos verticales y superiores de tal forma que se asemejen a un
contenedor invertido. La capacidad máxima de almacenamiento de la trampa es el volumen entre
el punto más alto de la trampa y el plano de derrame. Las trampas no están llenas de
hidrocarburos en el sentido de que hay una saturación irreductible de agua. Las trampas
estructurales son más fáciles de encontrar con métodos de exploración geofísica a diferencia de
las trampas estratigráficas.

Normalmente las trapas estructurales se originan por fallas (normales, inversas, cabalgaduras,
de cualquier tipo). Las fallas pueden ser fallas activas que no están selladas o las fallas inactivas
que no permiten que fluyan fluidos a través de ellas.

Las trampas Estratigráficas están relacionas a cambios de facies o ambientes de depositación.

También son asociados a la erosión y redepositación. Algunos ejemplos son Barras de arena
que forman barreras, zonas deltaicas de depositación, atolones que están rodeados por zonas
arrecifales carbonatadas.

2.4.5 Sincronía

Es el último elemento del sistema petrolero y se refiere a que todos los elementos estén
presentes en el mismo momento. De tal forma que se cada evento se de simultáneamente y no
excluyentemente. Un ejemplo sería que se tenga la roca generadora en el Titonico el jurásico
superior y la trampa se encuentre en el Toarciense del Jurásico inferior, el paso de migración se
da por un sistema de fracturas y las roca almacenadora se encuentra en arenas del Eoceno.
Aquí están presentes con los 4 elementos del sistema petrolero pero la acumulación no se podrá
dar por que la trampa se encuentra por debajo de la roca generadora. Al menos de que
existieran rocas más permeables en la parte inferior y sirvieran de paso en sentido al centro de la
Tierra. Otro ejemplo que se puede citar es si se tuviera una roca generadora en el Carniense del
Triásico Superior sujeta a condiciones de 260 °C y una presión promedio de 6000 psi por más de
900 millones de año. Se tiene una roca almacenadora en carbonatos naturalmente fracturados
con un sello generado por una trampa estratigráfica de cambio de facie y los pasos de migración
se debe a falla y sistema de fracturas de los carbonatos. Pero no se encontrará petróleo ya que
las condiciones de temperatura han superado por mucho a las de la ventana de generación del
aceite y gas por lo que probablemente se encuentre sea carbón mineral. Como se ve la sincronía
incluye la cronología y ubicación de los elementos del sistema petrolero de tal forma que se
encuentre en el tiempo y lugar indicados para que se dé un yacimiento.

2.5 Sistema Petrolero no Convencional

Los sistemas petroleros no convencionales son aquellos en los cuales se encuentran rocas
generadoras con propiedades muy buenas como: 𝑆= 20 %; porosidad promedio, So=65 %,
Sg=15%, Sio=8%, Siw=11%, con grandes volúmenes de roca pero la permeabilidad es tan baja
que no permite que ese petróleo migre afuera de la roca generadora. Existen yacimientos no
convencionales que a pesar que encuentren todos los elementos del sistema petrolero el crudo a
producir es de densidades API° muy bajas (11° API) por consecuencia perjudica el flujo en la
tubería de producción debido a su alta viscosidad (μ), alta gravedad específica (γ) y baja relación

- 18 -
de solubilidad (Rs). Que después de perforarse y se comience producción la presión del
yacimiento será insuficiente para superar la presión ejercida por la columna de hidrocarburos y
la presión en las paredes de pozo debido al flujo viscoso. Lo que ocasionará que el fluido se
quede sin fuerza para fluir a superficie. También se pueden encontrar yacimientos no
convencionales de alta presión y alta temperatura [arriba de 149 °C de temperatura de fondo de
pozo y presiones superiores a los 10,000 psi]. Estos son no convencionales debido a que los
equipos superficiales no pueden soportar las condiciones de operación y el riesgo de tener
reventones y descontroles durante la perforación es muy alto.

Como podrá leer el concepto de yacimiento no convencional engloba diferentes definiciones

pero la que se usará en esta tesis será de yacimientos con roca generadora con muy buenas
propiedades pero permeabilidades tan bajas que no sea posible extraer los hidrocarburos a
superficie. Es debido a esto que el fracturamiento hidráulico es usado para generar una
permeabilidad artificial a través de fisuras inducidas en la roca generadora.

Los sistemas petroleros no convencionales se nacen cuando se quiere explorar las rocas
generadoras que tienen permeabilidades muy bajas pero contienen volúmenes muy grandes de
hidrocarburos.

2.6 Ejemplo de un yacimiento Convencional

En esta sección se dará un ejemplo detallado de un yacimiento convencional. Se verá el

yacimiento analizándolo como un sistema petrolero.

Empezaremos con el ejemplo del yacimiento ingles el Campo FULMAR. Está ubicado en el mar
del norte y recibe su nombre bajo la clave mnemotécnica de los ríos ingleses: Forth, Usk,
Lydell, Mersey, Avon, Ribble. Los cuales representan las unidades productivas del Yacimiento.

2.6.1 Campo Fulmar

El Campo FULMAR fue descubierto en la región central del Mar del Norte en 1975. Se descubrió
cuando se volvieron a revisarse los estudios sísmicos que buscaban prospectos en el Pérmico.
En la Tabla 3 se describe su sistema petrolero.

- 19 -
 Roca Generadora Migración Roca Almacenadora Sello Sincronía

Areniscas de secuencias
regresivas, bien
distribuidas y
Roca Sello: Lutitas,
Formación de Calizas del seleccionadas. Cambio Todos los elementos se
arcillas y limolitas de la
Kimeridiano .Edad de litofacies de cueca encuentran en el tiempo
formación del
Geológica: Jurásico deltaica. Subdividida en correcto y las trampas
Kimeridiano en la región
Superior (154.1 a 145.6 6 unidades productivas y están almacenando el
sur y este. Lutitas del
Ma). Geología: lutita Se da por cambio de 16 subunidades. 6 Forth, petróleo que migra de
Grupo Chalk
madura, rica en materia facies. En la unidad 5 Usk, 4 Lydell, 3 las calizas del
despositadas en el
orgánica. COT: 5% Clyde se encuentra un Mersey,2 Avon y 1 Kimeridiano.
Terciario. Calizas del
.Kerógeno: tipo II, punto de derrame. Hubo Ribble. Unidades: 1,2, La migración se da en el
grupo Chalk en la parte
Reflectancia de la migración regional en el 3.1, 3.2, 3.3, 3.4, 3.5, 4, Eoceno a través de las
norte y sur. Trampa:
victrinita: 0.75 Cenozoico. 5.1, 5.2, 5.3, 6.1, 6.2. areniscas bien
Mixta. Estructural por su
Rendimiento: 0.032 v/v Espesor promedio de la seleccionadas. Se
forma de domo con un
aceite. Profundidad: formación 274 m. 815 almacena en la trampa
flanco buzando a 18-25°
2973 m. Espesor: 151.2 [stb] aciete y 500 [Bscf] hecha por las calizas que
al sureste y 10° al
m de gas húmedo. son la roca generadora.
suroeste.
Saturación irreductible
de agua 8-15% y residual
aceite 22-30%
- 20 -
 Tabla 3 Sistema petrolero FULMAR (Van der Helm, et al., 1990)

La información estratigráfica más relevante del campo FULMAR está resumida en la Fig. 8.

Fig. 8
Columna Litoestratigráfica (Van der Helm, 1990)

La roca almacenadora se divide por unidades productivas y cada una tiene un porcentaje del
volumen total de hidrocarburos de la cuenca. En la Fig. 8 se presenta una columna estratigráfica,

- 21 -
en la Fig. 9 el corte transversal de FULMAR y en la Tabla 4 se presentan las permeabilidades,
porosidades y volúmenes. (Van der Helm, 1990)

Fig. 9 Subdivisión del yacimiento FULMAR en unidades productivas

Unidad Porosidad Permeabilidad STOIIP

(%) (mD) (MBl)
6.1 14.5 0.28 ---
5.1 20.7 58.7 42
5.2 36.2
5.3 23.7 123 13.3
4 23.3 508 259.9
3.1 23.9
3.2 128.1
3.2 24.7 674 65.1
3.3 21.4 89
3.4 147.4
3.5 8.1
2 20.5 28 17
1 25.3 684 74
Total 815
Tabla 4 Porosidades,
permeabilidades y volúmenes de las unidades productivas FULMAR

- 22 -
Fig. 10 Historia de sepultamiento y tabla de sincronía del campo FULMAR (modificado de Magoon
y Dow, 1994, pg 12)

De la información anterior se puede ver como son los yacimientos convencionales. En las cuales
los hidrocarburos migran afuera de la roca generadora y quedan almacenados en alguna trampa
que se haya formado. En el caso del campo FULMAR las rocas almacenadora fueron las
areniscas que se encuentran debajo de la roca generadora [calizas del Kimeridiano] y la roca
madre funge como roca sello y trampa estructural en una parte del almacén. Estos son casos
particulares pero pueden encontrarse otros casos menos lógicos.

La información que fue presentada anteriormente es un análisis que se haría para el sistema
petrolero de una cuenca. En este caso del Campo FULMAR.

2.7 Ejemplo de un Yacimiento No Convencional

- 23 -
Se continuará con el análisis de un yacimiento no convencional para comparar como se
encuentran los elementos del sistema petrolero de un yacimiento convencional y no
convencional.

Se toma como referencia el Campo Bakken ubicado en Estados Unidos y Canadá. Ya que se
cuenta con mucha información del yacimiento y es uno de los más populares que actualmente se
explotan. El Bakken contribuye con una fracción considerable de la producción no convencional
de Estados Unidos. (drilling productivity report, EIA)

La formación del Bakken fue generada por la secuencia transgresiva llamada Secuencia
Tamaroa, está formada por una secuencia de limos en medio de dos estratos de lutitas
bituminosas. La formación es una zona de presión anormalmente alta por que los estratos de
lutitas bituminosas están en la ventana de generación del aceite. Tiene muy poca saturación de
agua por el alto COT [1 - 45 %C/g] y la roca es mojada por el aceite, cuenta con una muy baja
permeabilidad. El campo está 1100 psi por arriba del gradiente normal por los que el abatimiento
de presión no es tan grande.

La formación del Bakken se caracteriza por 3 estratos: el primero Bakken superior un estrato de
lutitas bituminosas con un espesor máximo de 15.2 [m], el segundo Bakken medio que son
limolitas calcáreas 22 [m] espesor máximo y tercero el Bakken inferior que consiste de lutitas
bituminosas con 7 [m] de espesor promedio. El espesor máximo de la formación alcanza los 44.2
[m]. Contiene un kerógeno de tipo de masivo que se infiere ser sapropélico y fue depositado en
un ambiente anóxico con materia orgánica proveniente de plancton y algas en su mayoría. Se ha
caculado por estudios de pirólisis que la formación del Bakken ha generado alrededor de 92.3
[MMMbbl], se estima que sólo el 10% del volumen fue expulsado a través de las fracturas
generadas por los esfuerzos in-situ así que se tiene un gran volumen atrapado en la roca
generadora (Webster, 1984). Para explotarla de la manera más eficiente se perforan pozos
horizontales que pasan por las formaciones: la formación Three Forks, el Bakken intermedio
[limolitas] o la formación Mission Cayon, y después se inducen fracturas hidráulicamente. Se ha
hecho una perforación masiva para estimular el mayor volumen posible.

2.7.1 Campo Bakken

El campo del Bakken tiene una extensión areal de 520,000 [km2]. La formación está divida en
tres miembros: el superior son lutitas bituminosas que presentan un poco de metamorfismo y un
muy alto contenido de materia orgánica, el miembro medio son limolitas con contenido calcáreo
bajo es el estrato más grande de la formación y el inferior que es el más delgado son lutitas
bituminosas con alto contenido orgánico. Por arriba de la formación del Bakken se encuentra la
formación Lodge Pole y por debajo la formación Three Forks. Estas dos formaciones son
impermeables así como los estratos de lutitas.

El Bakken tiene una sistema de naturalmente fracturado debido a la subducción que lentamente
tuvo la corteza terrestre esto es mejor conocido como una cuenca intracratónica.

- 24 -
 Fig. 21 Esquema de la sección transversal de la Dakota del Norte (modificado de Meissner, 1984)

En al Fig. 11 se puede apreciar como la subducción de la base generó hace millones de años
una cuenca intracratónica, la cual ha sido rellenada con sedimentos desde el Paleozoico. La
subducción ha sido muy lenta y se podría confundir con una cabalgadura. Es importante ver que
la subducción genera una zona con un régimen de esfuerzos anormal que ha generado fracturas
por las cuales los hidrocarburos han migrado a las formaciones adyacentes y subyacentes,
generando yacimientos convencionales como campo Antelope [areniscas], campo Salt Greek o
campo Elkhorn Ranch. Como se mencionó en la sección 2.7 se estima que sólo un 10% de
volumen original ha migrado el remanente sigue en la formación del Bakken.

Fig. 12 Sección transversal de la región del Bakken que corre de O-E (Modificado de Meissner,
1984)

En la figura 12 se puede apreciar las formaciones adyacentes y subyacentes y una descripción

petrofísica de las rocas que las conforman, también muestra la sección térmicamente madura de

- 25 -
la roca generadora. Conforme se fue implementando el fracturamiento hidráulico lo proyectos se
enfocaron más en la formación del Bakken que en los yacimientos convencionales de la cuenca
de Williston. En la figura 12 se puede apreciar el sistema naturalmente fracturado originado por el
exceso de presión causado por la madurez de la roca generadora y los esfuerzos distensivos
producto de la subducción de la cuenca intracratónica. Es común encontrar una zona de presión
anormalmente alta en las lutitas del Bakken ya que están a las condiciones adecuadas para la
ventana de generación del aceite. Lo que causa que se generen hidrocarburos y la presión de
poro comience a elevarse. Cuando la presión de poro junto con el régimen de esfuerzos regional
sobrepasan las presión de fractura de la matriz se generan las fisuras que conectan la formación
del Bakken con la formación Lodgepole y Three Forks abriendo un paso de migración y
desahogue de presión y permitiendo que se generen más hidrocarburos en su lugar.

Debido a que la formación del Bakken tiene una gran extensión y hay partes con mayor
profundidad que otros, la formación no entro en su totalidad a la ventana de generación. En
algunas partes desde los 2,331 – 2,438 [m] se comenzaron a generar hidrocarburos y en otras a
los 3,048 [m]. Cabe mencionar los diferentes ambientes de los cuales fueron aportados la
materia orgánica que conforman el Bakken superior e inferior ya que el contenido de oxígeno del
kerógeno varia si la materia orgánica se aportó de un ambiente terrestre siendo el caso del
miembro del Bakken superior o si la materia orgánica proviene de un ambiente marino como es
el caso de Bakken inferior.

Durante la perforación es fácil de encontrar la formación del Bakken por las lecturas tan
particulares que muestran los registros petrofísicos como las resistividades que son muy altas
debido al bajo contenido de agua intersticial y la alta saturación de hidrocarburos además de la
alta presión. El gradiente de presión es anormalmente alto, pudiendo llegar hasta 0.73 [psi/ft]
cuando se llega a las lutitas y las lecturas del registro sónico son muy grandes debido a la alta
presión aunque el medio se encuentra saturado por una fase líquida. Las lecturas del registro
sónico son parecidas a la un medio con alta densidad pero en realidad las partículas se
encuentran muy apretadas por la alta presión. Las lecturas del Rayos gamma son altas cuando
miden las lutitas debido a que son rocas conformadas por sedimentos muy finos que pueden
contener alto porcentajes de minerales como Torio o Potasio que emiten radiación. Con estos
parámetros es fácil identificar la formación del Bakken.

En la figura 13 es la columna estratigráfica de la región del Bakken que alcanza los 3,048 [m] de
profundidad, ya que en otras partes de la formación los estratos de Bakken se encuentran más
someros a profundidades de 2,316 – 2,438 [m] aproximadamente. Dentro de las formaciones
inferiores que pertenecen al Devónico la formación Three Forks es impermeable, pero la
formación Nisku es permeable y porosa. Así mismo las formación superiores como Charles,
Mission Canyon y Lodge Pole , que pertenecen al grupo Madison son de interés ya que Lodge
Pole es un sello pero la formación Mission Canyon presenta buenas propiedades como almacén.
Las formaciones más importantes durante la perforación son Three Forks y Lodge Pole debido
que son las profundidades óptimas para que las fracturas alcancen el Bakken inferior y
superior.

- 26 -
Fig. 13 Columna Estratigráfica de la región del Bakken (Meissner, 1984)

- 27 -
 Roca Generadora Migración Roca Almacenadora Roca Sello Sincronía
Tabla 5 Sistema petrolero del Campo Bakken

Edad: Mississipiense
Profundidad: 2740 m
Tipo roca: Lutitas
bituminosas
Kerógeno: tipo II Los estratos de lutitas
COT promedio: 11.3 funcionan como
Sistema de fracturas La formación del Todos los elementos
Wt% almacén y sello
Espesores: 1ro/ 7m generado por la zona Bakken funge como debido a la baja se encuentran en el
donde se genera una almacenadora. Se tiempo geológico
2do/ 15.2 m permeabilidad y falta
Edad entrada ventana cabalgadura debido a estima que 90% de os de fracturas para que indicado. Los estratos
de generación: hace una falla inversa. HC´s no han sido migren los de lutitas todavía no
75 Ma Fracturas originadas liberados. Los dos hidrocarburos. Las alcanzan la
Ventana de por la sobrepresión de estratos que limitan al limolitas calcáreas no metagénesis ni en la
generación: 100°C los hidrocarburos estrato de limolitas dejan que se ventana de
Reflectancia de la sobre la formación. son los almacenes. interconecten las generación de gas.
victrinita:0.75 lutitas.
Rendimiento: 600mg
/g COT
- 28 -
Analizando la sincronía en la tabla 5 se puede concluir que el único evento que falta para
completar el sistema petrolero es la migración entre los sellos y los almacenes, por esta razón se
considera un yacimiento no convencional.

Fig. 14 Historia de sepultamiento y tabla de sincronía del campo Bakken (Modificado de Meissner,
1984)

La solución a este problema ha sido el fracturamiento hidráulico que ha abierto la posibilidad de

explotar la formación del Bakken cuando se alcanzan y perforan horizontalmente las
formaciones de Three Forks, Lodge Pole o Bakken media, después se diseña el programa
óptimo de fracturamiento, tomando en cuenta número de etapas y espaciamiento para que se

- 29 -
distribuyan las fracturas a lo largo de la sección horizontal del pozo y se estimule el mayor
volumen de roca posible.

Esta técnica se comenzó a implementar a partir del 2004 y se ha utilizado con más auge en los
años recientes. Actualmente los proyectos de aceite del Bakken han disminuido por la baja del
precio del barril que ha originado que no sea rentable seguir perforando y fracturando los pozos.
En cuanto aumente el precio del barril se reanudarán los proyectos de perforación. La región del
Bakken tiene un sistema de producción bajo licencia. En la cual cualquier compañía puede
participar, siempre y cuando se cumpla con el porcentaje de impuestos que sea acordado con el
gobierno de Dakota.

- 30 -
 CAPÍTULO 3

CARACTERIZACIÓN DE UN YACIMIENTOS NO CONVENCIONAL

3.1 La importancia de la Caracterización Dinámica

Cómo se sabe, los niveles de estudio existentes para ciencias de la tierra van desde lo
microscópico hasta lo macroscópico o en este caso hasta lo nanoscópico. Un ejemplo del nivel
macroscópico es cuando vemos a los yacimientos como una sección de la corteza terrestre
donde se encuentran los elementos necesarios para que exista la acumulación de hidrocarburos.
Un ejemplo sería una sección trasversal de alguna región estudiada o el sistema petrolero. Un
ejemplo de la parte microscópica es una lámina esmerilada que se obtuvo de un núcleo y se
observa los tipos de porosidades a través del microscopio para caracterizar esa sección del
yacimiento. Al hablar de nanoscópico es referirse a la permeabilidad de nano darcies que se
encuentra en los yacimientos de lutitas de baja permeabilidad [tight oil] y como se puede
modificar esta propiedad inherente a través del fracturamiento hidráulico y es uso del
apuntalantes.

La caracterización dinámica junto con muchas otras ramas de la ingeniería petrolera como:
productividad de pozos, caracterización estática de yacimientos, comportamiento de yacimientos
y simulación numérica de yacimientos. Se encargan de estudiar las razones, elementos que
afectan el flujo, comportamiento del movimiento de los fluidos a las condiciones del yacimiento.
Es lógico que para hacer estos modelos o estudios es indispensable contar con una cantidad de
información considerable para representar el yacimiento al grado de certidumbre que se desea.
Estos estudios son muy certeros cuando se tiene una clasificación correcta de los mecanismos
de empuje, las propiedades de la roca y del fluido. Pero requieren de mayor tiempo a diferencia
de las curvas de declinación y trabajan en un nivel más certero a diferencia de las curvas de
declinación. La caracterización dinámica de yacimientos estudia al yacimiento como un ente que
interactúa bajo ciertas condiciones físicas [presión, temperatura, composición, etc.] de una forma.
Al compararlo con las curvas de declinación que sólo trabajan con los pozos y su interacción con
la producción y en el tiempo. Como se puede ver los estudios de Caracterización Dinámica son
más complejos y más certeros pero requieren de mayor información y tiempo para elaborarse.
Este tipo de estudios se realizan después de haberse perforado el pozo exploratorio y tener
algunos pozos de desarrollo. Son complementarios al desarrollo del yacimiento y son usados
para hacer el mejor plan de explotación del campo. Las curvas de declinación de usan antes o
después de haberse perforado el pozo exploratorio para hacer una evaluación económica y
conseguir un financiamiento para empezar el proyecto o conseguir la inversionista.

3.2 Proceso de la Caracterización Dinámica y Simulación de un yacimiento convencional

La Caracterización Dinámica sigue el siguiente proceso:

 Descripción del yacimiento

- 31 -
  Definir los mecanismos de empuje
 Realizar el modelo matemático
 Realizar el modelo numérico
 Realizar el programa en computadora que demuestre el modelo numérico
 Validar el modelo
 Ajustar el historial de producción al modelo
 Predecir el comportamiento del yacimiento

3.2.1 Descripción del yacimiento convencional

La descripción del yacimiento se hace cuando se toman los estudios de geología superficial,
geología estructural, estudios geofísicos y petrofísicos y se condensan en un modelo descriptivo
exhaustivo hecho a base de celdas donde cada celda representa un volumen parcial del
yacimiento. A cada celda se le atribuyen las propiedades de porosidad, saturación, densidad de
fluido, presión, temperatura, permeabilidad, etc.

Después se realiza el escalamiento para deshacerse de muchas celdas y no tener un modelo tan
robusto y grande. El escalamiento trata de juntar las celdas cuyas propiedades sean iguales en
celdas con volúmenes más grandes o aunque tengas propiedades diferentes juntarlas en una
celda equivalente con propiedades promedio. Resultando en un modelo mucho más simple y con
menos información que en el modelo estático.

El objetivo del modelo dinámico es simular el comportamiento de los fluidos a condiciones del
yacimiento. También detectar y evaluar los elementos que afectan el comportamiento de un pozo
como el daño, el almacenamiento u otros.

3.2.2 Definiendo los mecanismos de empuje de un yacimiento convencional

Definir los mecanismos de empuje es caracterizar cuál de los siguientes mecanismos está
dominando el movimiento de los fluidos del yacimiento y en qué etapa de la explotación:
expansión del sistema roca fluido

 Empuje por encogimiento del sistema roca fluido

 Empuje por acuífero asociado
 Empuje por casquete de gas
 Empuje por segregación gravitacional
 Empuje por gas disuelto
 Combinación de varios empujes

Cada mecanismo de empuje puede estar trabajando en una etapa de la vida productiva del
yacimiento. Por ejemplo, cuando el yacimiento se encuentra bajo saturado se encuentra un
empuje por gas en solución en el aceite pero al bajar a la presión de burbuja se empieza a
desprender el gas disuelto y se genera un casquete de gas y es aquí donde entra el empuje por
casquete de gas asociado. Otro ejemplo es el mecanismo de empuje por encogimiento del
sistema roca fluido. Cuando disminuye la presión del yacimiento la roca comienza a encogerse

- 32 -
por la liberación de la presión que aportaban los fluidos atrapados en la roca lo que genera una
fuerza que es proporcional a la compresibilidad de la roca y que ira expulsando los fluidos.
Cuando se tiene una falla conductiva en la formación y aporta agua eso puede asociarse a un
mecanismo de empuje por acuífero asociado. También cuando las fases se están separando en
el yacimiento se puede encontrar el empuje por segregación gravitacional. Éste es debido a las
fuerzas de flotación que tiene las diferentes fases cuando están en interactuando. Las gotas de
las fases comienzan a moverse a su lugar respectivo y conforma la fase crece la fuerza de
segregación es mayor. Ejemplificando cada mecanismo de empuje se puede analizar cada uno
como individual pero en realidad están combinado el efecto de varios en el yacimiento. Es así
como se caracterizan los mecanismos de empuje que actúan a nivel global en el yacimiento.
Algo que podría modificarse es cuando se tiene un sistema de recuperación secundaria por
inyección de agua y se puede caracterizar como un empuje por acuífero asociado el cual va a
crecer a la tasa que se esté inyectando el agua.

3.2.3 Realizar el Modelo Matemático

3.2.3.1 Deducción de la Ecuación de Continuidad para un yacimiento convencional

Realizar el modelado matemático se refiere a traducir todas las fuerzas que actuaran en la
dinámica de nuestros fluidos en formulas o funciones en función de la presión, tiempo, ubicación
del yacimiento o el pozo de estudio. Para lo cual se definen las fuerzas que gobiernan el flujo:

 Fuerza de Presión

𝑆𝑆 = −𝑆𝑆∆𝑆 ……(3.1)

𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆ó𝑆
−∇ 𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆ó𝑆
∆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆ú𝑆𝑆𝑆𝑆𝑆

 Efectos Gravitacionales

𝑆𝑆𝑆 = (𝑆 − 𝑆𝑆)g∇ 𝑆∆𝑆 ……(3.2)

𝑆𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆

𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆
𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆
∇ 𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆
∆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆ú𝑆𝑆𝑆𝑆𝑆 (𝑆𝑆𝑆𝑆𝑆𝑆𝑆,
𝑆𝑆𝑆𝑆𝑆) g = gravedad

 Efectos Viscosos

- 33 -
 𝑆𝑆
𝑆=− 𝑆 ∆𝑆 ……(3.3)
𝑆 𝑆𝑆 𝑆

𝑆𝑆= 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆

𝑆𝑆
= 𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆
𝑆𝑆
𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆
∆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆

 Ley de Darcy
𝑆𝑆𝑆
𝑆=−𝑆 [∇ 𝑆 + ∆𝑆𝑆∇ ℎ ] ……(3.4)
𝑆 𝑆𝑆 𝑆

𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆

𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆
𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆
𝑆𝑆𝑆
= 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆
𝑆
𝑆
∇ 𝑆𝑆 + ∆𝑆𝑆∇ ℎ = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆ó𝑆 ℎ 𝑆𝑆𝑆𝑆𝑆𝑆á𝑆𝑆𝑆𝑆 𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆

Al conocer las fuerzas que generan el movimiento de los fluidos cuando se generen los cambios
de presión al abrir las válvula del pozo o hacer una prueba de presión se puede hacer el modelo
matemático. Para hacer el modelado matemático tenemos que tomar en cuenta que la materia
no se crea ni se destruye. El volumen es constante en los yacimientos donde los hidrocarburos
están almacenados en una trampa. Este volumen se le conoce como VOHI – Volumen Original
de Hidrocarburos In-situ. La única manera que tiene de salir es a través del conducto que genera
el pozo. El efecto físico que genera el pozo es una disminución de la presión en las vecindades
del pozo [disparos] que están en contacto con el pozo. Esto genera que los hidrocarburos que
están a mucha mayor presión que la superficie se expandan. Como se sabe los líquidos tienen
muy baja compresibilidad pero los gases tiene una muy grande y la roca también puede
comprimirse o expandirse hasta que se fractura. Por eso la producción es igual a la expansión
del yacimiento. Es necesaria una ecuación que relacione la cantidad de masa existente en el
yacimiento con su cambio de volumen debido a su cambio de fase. Es debido a esto que usamos
la Ecuación de Continuidad. La ecuación pude usarse para igualar masa, energía y momento.
También pude expresarse en diferentes formas. De aquí se pasa a la ecuación de continuidad.
La cual se puede expresar en coordenadas:

 Rectangulares
𝑆 𝑆 𝑆 𝑆
(𝑆𝑆 ) + (𝑆𝑆 ) + (𝑆𝑆 ) = − (𝑆𝑆) …………(3.5)
𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆

𝑆 𝑆 𝑆

1𝑆 1𝑆
(𝑆𝑆𝑆 ) + (𝑆𝑆
 Cilíndricas
- 34 -
 𝑆 𝑆
)+ (𝑆𝑆 ) = − (𝑆𝑆) …………(3.6)
𝑆 𝑆𝑆 𝑆 𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆
𝑆
𝑆

- 35 -
  Esféricas
1 𝑆 1 𝑆 1 𝑆 𝑆
(𝑆2𝑆𝑆 ) + (sin 𝑆 )+ )=− (𝑆𝑆) …………(3.7)
𝑆𝑆 (𝑆𝑆
𝑆2 𝑆𝑆 𝑆 𝑆 sin 𝑆 𝑆 𝑆 sin 𝑆 𝑆𝑆 𝑆𝑆
𝑆𝑆 𝑆

 En forma general
𝑆
∇ 2 ∙ (𝑆 𝑆) = − (𝑆𝑆) …………(3.8)
𝑆𝑆

La ecuación de continuidad no es más que la expresión matemática de relación entre el cambio

de presión en algún punto del espacio analizado con respecto a las propiedades del fluido.
Expansión del yacimiento es iguala producción.

Esta expresión esta en forma diferencial porque se quiere tomar en cuenta a mayor detalle los
cambios que pueden sufrir los fluidos cuando están sujetos a un gradiente de presiones. La
forma en que se expresan [rectangular, esférica o cilíndrica] depende de cómo se acomode
mejor al análisis. Por ejemplo si se tiene un estrato de 20 m de espesor que no varía su espesor
de forma horizontal y se conservan las mismas propiedades de saturación, composición,
porosidad y permeabilidad será más conveniente analizarlo con coordenadas rectangulares
suponiendo que el pozo se encuentra en una de las esquinas del estrato. A diferencia de si se
tuviera un pozo que está drenando de forma uniforme un yacimiento de carbonatos naturalmente
fracturados se puede usar la expresión en coordenadas cilíndricas o esféricas para estudiar a
mejor detalle el movimiento de los fluidos que están sujetos a un gradiente de presiones. El
gradiente de presión se va a originar en las vecindades del pozo. Éste está conectado a la
superficie y a través de la válvula se controlara el abatimiento de presión.

3.2.3.2 Regímenes de Flujo de yacimientos convencionales

Para hacer uso de la Ecuación de Continuidad se debe tomar en cuenta el régimen de flujo. Los
cuales son las relaciones existentes entre los cambios de presión y la expansión del sistema de
fluidos:

 Estacionario
La presión no es función del tiempo y es una variable independiente. Haciendo mucho
más sencillo el modelo numérico de la ecuación de difusión. Se genera los yacimientos
volumétricos. Los cuales tienen fronteras activas, se puede encontrar este régimen
cuando hay un aumento de presión por acuífero asociado o cuando las fallas que forman
la trampa y el acuífero aportan presión al yacimiento.

- 36 -
 𝑆𝑆
|𝑆 = 0
𝑆𝑆

Fig. 35 Régimen Estacionario (Gallardo E., 2015 modificado)

 Pseudoestacionario
El yacimiento se comporta como un tanque donde la presión decae a un ritmo constante.
La presión varía uniformemente con respecto al tiempo. Este régimen de flujo se presenta
cuando la onda de presión ha alcanzado las fronteras del yacimiento y se produce a un
gasto constante.
| = 𝑆𝑆𝑆
𝑆𝑆
𝑆𝑆 𝑆

Fig. 16 Régimen Pseudoestacionario (Gallardo E., 2015 modificado)

 Transitorio
Este régimen antecede al régimen pseudoestacionario y estacionario. Su duración
depende de la Difusividad Hidráulica. Es el primer régimen que se observará cuando
comienza la producción en un pozo. Representa el comienzo del movimiento de los
fluidos dentro del yacimiento sin tomar en cuenta el daño a la formación y el efecto de
almacenamiento del pozo. La presión es una función del tiempo y espacio.

𝑆 = 𝑆(𝑆)

𝑆𝑆
|𝑆 ≠ 0
𝑆𝑆

- 37 -
 Fig. 17 Régimen Transitorio (Gallardo E., 2015 modificado)

Una vez que se tiene identificado el régimen de flujo también se deben de hacer las siguientes
consideraciones para usar la ecuación de Continuidad:

1. Sólo una fase satura el medio poroso

2. El fluido es newtoniano y se comporta de acuerdo a la ley de Darcy
3. El yacimiento es isotrópico y homogéneo de acuerdo a las propiedades de la roca y fluido
4. La permeabilidad del yacimiento no depende de la presión
5. La viscosidad y compresibilidad del fluido son constantes

Estas son las consideraciones hechas para un yacimiento saturado al 100% por aceite negro de
bajo encogimiento. Que tiene muy poco gas disuelto y es necesario un cambio considerable de
presión para que se alcance la presión de burbuja. Por lo que no entrara en la campana de
saturación tan rápido a diferencia de fluidos más volátiles. El fluido tiene un comportamiento
newtoniano como el agua y la relación deformación-esfuerzo cortante es lineal. Se puede
concluir que a un cambio de presión habrá un desplazamiento a la zona de menor potencial. La
consideración que la permeabilidad no cambia con respecto a la presión se puede hacer pero
hay que tomar en cuenta que la porosidad puede cambiar debido a la liberación de presión del
fluido sobre la matriz de la roca. Las Fuerzas que dominan en el sistema roca fluido son las
fuerzas que actúan sobre la matriz de la roca y las fuerzas que actúan en el fluido contenido en
el poro. Al disminuir la presión del poro la compresión en la matriz disminuye y puede expandirse
un poco en proporción a su Coeficiente de compresibilidad. Generando una disminución del
volumen poroso y por ende de la permeabilidad. A este fenómeno se le conoce como mecanismo
de empuje por el sistema roca fluido. Al disminuir la presión y disminuir la porosidad los fluidos
son empujados. Esta disminución es tan pequeña que puede ser despreciada en algunos casos.
(Gallardo 2015)

Estas consideraciones cambiaran para los yacimientos no convencionales que se quieran

analizar. Empezando por ser fracturados hidráulicamente. En segundo los yacimientos están
saturados al 100% por aceite porque éste apenas está siendo expulsado por la roca generadora
y va a tener muy pequeños cortes de agua. Normalmente los yacimientos no convencionales se
encuentran en la ventana de generación del aceite y tiene presiones anormalmente altas. Como

- 38 -
se vio en el ejemplo del campo Bakken. El agua será un producto de la ruptura de las cadenas
de las ramas aldehidas de los hidrocarburos. Las rocas son mojadas por el aceite por lo que la
saturación residual de aceite incrementara a diferencia de otros yacimientos. Se puede encontrar
en algunas partes comportamientos no newtoniamos por las mezclas de hidrocarburos con otros
componentes o formación de micelas, pero no será muy común. Las propiedades del yacimiento
serán homogéneas e isotrópicas.

Al hacer estas consideraciones podemos pasar a la siguiente parte del modelado matemático en
el cual sustituimos la Ecuación de Darcy, Ecuación de Equilibrio Fases para obtener la Ecuación
de Difusividad en notación diferencial (Ec. 3.9):

Ecuación de Difusividad

𝑆𝑆𝑆𝑆 𝑆𝑆
∇ 2𝑆 = …………(3.9)
𝑆 𝑆𝑆

Constante de Difusividad Hidráulica

𝑆
ƞ=
𝑆𝑆𝑆 …………(3.10)
𝑆

Transmisilibilidad

𝑆ℎ
𝑆=
𝑆𝑆 …………(3.11)

La constante de Difusividad Hidráulica (Ec. 3.10) nos dice que tan fácil se transmiten las ondas
de presión en el yacimiento. La transmisibilidad (Ec. 3.11) representa la facilidad con que fluye
un fluido en el medio poroso. (Arana, 2015)

Las ondas son propagadas en los fluidos que están interconectados en el yacimiento. Esto es un
reflejo de la permeabilidad que existe en el yacimiento. Las ondas son generadas por la apertura
de la válvula de producción en superficie. Pero al comienzo se genera un fenómeno conocido
como efecto de almacenamiento del pozo. Que no es otra cosa que la expansión y movimiento
de los fluidos contenidos dentro de la tubería de producción. Las ondas de presión son
proporcionales a la caída de presión generada por la apertura de la válvula. Esta es proporcional
al diámetro de la válvula y se moverán a una rapidez proporcional a la constante de
Transmisibilidad. Los estudios de pruebas de presión son muy importantes porque nos dan
indirectamente datos como el área de drene, la permeabilidad promedio de la formación, el daño
a la formación y la constante de almacenamiento. Con esto se pueden estimar las reservas 1P
que tienen un 90% de Probabilidad [P (.9)] de ser extraídas. Identificar los regímenes de flujo que
nos va a ayudar a saber cuándo las ondas han alcanzado las fronteras del yacimiento. Al saber
que tan lejos están las fronteras podemos delimitar el yacimiento. También para conocer la

- 39 -
geometría del flujo presente en el yacimiento y saber si las fallas son conductivas o
impermeables.

3.2.3.3 Geometrías de Flujo para yacimientos convencionales

El Movimiento de los fluidos puede generar patrones, mejor conocidos como geometrías de flujo.
Esta depende de la geometría del área drenada y la penetración del pozo ya sea parcial - total y
vertical - horizontal. Tipos de geometría de flujo pueden ser:

 Lineal
 Bilineal
 Esférico
 Radial

Dependiendo de la geometría de flujo se calculara la presión

3.2.3.3.1 Flujo Lineal

Fig. 18 Flujo Lineal (Gallardo, 2015)

Tipo de flujo característico de una fractura que se generan cuando se tienen un sistema naturalmente
fracturado o uno hidráulicamente fracturado donde el espero de la fractura permanece constante. Este
flujo se caracteriza por ser en una sola dirección y se idealiza como si fuese un prisma con las dimensiones
de altura, longitud y espesor. Cuando se tiene un flujo en la fractura y de la matriz a la fractura se
presenta el flujo bilineal. La solución con todas las variables es de la siguiente manera para un pozo que
produce de un yacimiento infinito con flujo lineal a gasto constante.
2
𝑆𝑆𝑆𝑆𝑆 𝑆 𝑆𝑆 1/2 −(𝑆𝑆𝑆 𝑆𝑆 ) 2
∆P(𝑆, 𝑆) = [2√ ( ) 𝑆 4𝑆𝑆𝑆 −𝑆𝑆𝑆𝑆𝑆 (√
𝑆𝑆𝑆𝑆𝑆 )] …………(3.12)
𝑆𝑆ℎ 𝑆 𝑆𝑆𝑆𝑆 4𝑆𝑆𝑆

𝑆𝑆𝑆𝑆(𝑆): 𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆

𝑆𝑆 𝑆 β : Factores de conversión que dependen del sistema de unidades utilizados inglés o métrico

𝑆𝑆 = 887.2 (𝑆𝑆𝑆𝑆𝑆𝑆), 119.58 (𝑆é𝑆𝑆𝑆𝑆𝑆) 𝑆 𝑆 = 2.637𝑆10−4 (𝑆𝑆𝑆𝑆𝑆𝑆), 3.489𝑆10−4 (𝑆é𝑆𝑆𝑆𝑆𝑆)

𝑆𝑆𝑆 (𝑆𝑆 = 0, 𝑆𝑆 ) = √4/𝑆√𝑆𝑆 …………(3.13)

- 40 -
 Fig. 19 Comportamiento del Flujo Lineal en variables adimensionales

3.2.3.3.2 Flujo Radial

Fig. 20 Flujo Radial (Gallardo, 2015)

Esta geometría de flujo es característica de la producción de cualquier pozo cuando es un yacimiento

infinito y la onda de presión no ha alcanzado las fronteras o cuando el pozo se encuentra alejado de las
fronteras o las fronteras son fallas impermeables. La Solución Línea Fuente es la aproximación que
resuelve el flujo radial con muy poco error. Pero también podemos encontrar la solución en variables
reales.
𝑆𝑆𝑆
P = 𝑆 − 70.6 [ln|𝑆| + ln 𝑆
ƞ
+ 0.8091] …………(3.14)

𝑆 𝑆ℎ
|𝑆2 |

- 40 -
 Fig. 41 Comportamiento del Flujo Radial en variables adimensionales y escala logarítmica

3.2.3.3.3 Flujo Esférico

Fig. 22 Flujo Esférico (Gallardo, 2015)

Geometría de flujo presente en un yacimiento con un espesor muy grande cuya penetración del pozo es
parcial. Cuando se usa el término parcial es para referirse a que los disparos están en una sección
pequeña del estrato. Es parecido al flujo radial pero vamos a tener una aportación extra por la parte
inferior que se extrae a diferencia de que fuesen radios concéntricos de los cuales las ondas se propagan
hacia afuera. La geometría esférica es tiene un comportamiento diferente con el abatimiento de la
presión.
1
𝑆𝑆 (𝑆 = 1, 𝑆𝑆 ) = 1 − …………(3.15)
𝑆 (𝑆𝑆𝑆 )0.5
𝑆

- 41 -
 Fig. 25 Comportamiento del Flujo Esférico en variables adimensionales

3.2.3.3.4 Flujo Bilineal

Fig. 24 Flujo Bilineal (Gallardo, 2015)

Geometría de flujo característica del movimiento de fluidos de la matriz al pozo a través de una
fractura. Cuando se realiza el fracturamiento hidráulico se tienen esta geometría. Para poder
realizar el cálculo de la presión adimensional se tiene que haber calculado antes el Factor de
Conductividad ¨Fcd¨. Como se explico en el

- 42 -
 𝑆𝑆𝑆𝑆
𝑆𝑆𝑆 = 𝑆𝑆 …………(3.13)
𝑆

𝑆𝑆: 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆

𝑆𝑆: 𝑆𝑆𝑆ℎ 𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆
𝑆: 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆
𝑆𝑆: 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆

2.54 4
𝑆𝑆𝑆 = √𝑆𝑆 𝑆𝑆
𝑆𝑆𝑆
0.5 …………(3.14)

Fig. 25 Comportamiento del Flujo Bilineal en variables adimensionales con un Factor de

Conductividad de 0.2 y 0.6

Fig. 26 Comparación del comportamiento de las geometrías de Flujo

- 43 -
3.2.4 Realizar el Modelo Numérico

Podemos encontrar el modelo matemático de dos formas. Las soluciones se derivan de la

linealidad de la Ecuación de la Difusividad. Si su relación es lineal se usa la solución de Línea
Fuente. Si la solución no es lineal usamos Simulación Numérica de Yacimientos.

3.2.4.1 Solución Línea Fuente

La Solución línea Fuente, comúnmente ocupada en caracterización dinámica de yacimientos se

resume en la siguiente fórmula expresada en variables adimensionales y la Integral Exponencial.

1 − 𝑆𝑆2
𝑆 = − 𝑆( )
𝑆
2 𝑆 4𝑆𝑆 …………(3.15)

𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆. Es una aproximación de la solución que es válida en todo el

𝑆2
dominio de las x que en este caso
4𝑆
𝑆
será − .
𝑆

Tiene 3 aproximaciones:

𝑆 < 0.01 𝑆𝑆 (−𝑆) = 𝑆𝑆|1.781𝑆| (3.16)

0.01 < 𝑆 < 25 𝑆𝑆 (−𝑆) = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 (3.17)

𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆
𝑆2
25 ≪ 𝑆 𝑆 (−𝑆) = 𝑆𝑆 | 𝑆| + 0.8091 (3.18)
𝑆
𝑆𝑆
Tabla 6 Aproximaciones de la Integral Exponencial

La Solución Línea Fuente está tomada para un gasto constante y se expresa en variables
adimensionales para facilitar su comprensión. Pero es necesario calcular el valor de cada una
antes de usar la Solución de Línea Fuente. Las variables adimensionales estarán en función de
la Geometría e Flujo.

Lineal Radial Esférico

Esquema

𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆

Tiempo 𝑆𝑆𝑆 = 𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆ℎ =
𝑆𝑆𝑆 𝑆2 𝑆𝑆𝑆 𝑆2 𝑆𝑆𝑆𝑆𝑆𝑆2
𝑆 𝑆 𝑆

- 44 -
 𝑆𝑆ℎ 𝑆ℎ 𝑆𝑆𝑆
𝑆𝑆𝑆 = (𝑆𝑆 − 𝑆𝑆 = (𝑆𝑆 − 𝑆) 𝑆𝑆𝑆𝑆𝑆ℎ = (𝑆𝑆 −
𝑆 𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆 𝑆
𝑆)𝑆𝑆𝑆
Presión
𝑆)
𝑆 𝑆
𝑆 𝑆 𝑆
Radio 𝑆𝑆 = 𝑆𝑆 = 𝑆𝑆 =
𝑆 𝑆𝑆 𝑆𝑆
Tabla 7 Valores de las Variables Adimensionales para diferentes geometrías de flujo

A diferencia de las soluciones de Productividad de Pozos en las que el gasto es variable a

condiciones de presión constante. En Caracterización Dinámica el gasto es constante y lo que
varía es la presión. Es muy importante tomarlo en cuanta ya que la forma en que se derive la
solución cambiara si las condiciones cambian.

3.2.4.2 Simulación Numérica de Yacimientos

Si la Ecuación de Difusividad no tiene una solución lineal y se quiere predecir mejor el

comportamiento del yacimiento podemos usar la Simulación Matemática de Yacimientos. Se
toma en cuenta la expresión matemática que describe el comportamiento de los fluidos
conociendo el régimen de flujo, las características del yacimiento, las fuerzas y mecanismos de
empuje que gobiernan se pasa al modelo matemático. Al ser derivadas sin relaciones lineales no
podemos aproximarlas de forma sencillas por lo que se usan los métodos numéricos como las
diferencias finitas. Con su ayuda se pueden traducir las ecuaciones diferenciales en soluciones
numéricas. En la presente tesis se tratara de manera muy somera este tema por no ser de
interés para el desarrollo del tema.

Las derivadas se pueden aproximar de diferentes formas:

Fig. 27 Como aproximar la derivada de la presión tomando en cuenta la dirección de la derivada

𝑆𝑆 𝑆𝑆+1−𝑆𝑆
 Por la derecha = …………(3.19)
𝑆𝑆 𝑆𝑆
𝑆𝑆 𝑆𝑆−𝑆𝑆−1
 Por la izquierda = …………(3.20)
𝑆𝑆 𝑆𝑆
𝑆𝑆 𝑆𝑆+1−𝑆𝑆−1
 Centrada = …………(3.21)
𝑆𝑆 2𝑆𝑆

Un ejemplo de cómo sería la aproximación de la derivada en dirección de las x es la siguiente:

- 45 -
 Fig. 28 Arreglo de las celdas en simulación numérica para hacer el sistema de ecuaciones

1 𝑆𝑆+1 − 2𝑆𝑆+1𝑆+1
+ 𝑆𝑆+1 𝑆 𝑆𝑆 − 2𝑆𝑆 + 𝑆𝑆 𝑆𝑆𝑆𝑆
𝑆 𝑆
𝑆 −𝑆 [ 𝑆+1 𝑆 𝑆−1
+ 𝑆+1 𝑆 𝑆−1
]= ( ) …………(3.22)
2 (∆𝑆)2 (∆𝑆)2 𝑆 ∆𝑆

𝑆𝑆𝑆+1 presión en la celda i en el tiempo siguiente

𝑆
𝑆𝑆+1presión en la celda derecha a la i en el tiempo actual
𝑆+1
𝑆𝑆+1 presión de la celda derecha en el tiempo siguiente
𝑆+1
𝑆𝑆−1 presión en la celda izquierda en el tiempo siguiente

∆𝑆 es la longitud de la celda

En este caso se está usando la formulación del tipo Crack-Nicholson porque involucra valores
de la presión en el paso actual y en el paso siguiente. Los valores en el paso siguiente serían las
incógnitas que se quieren conocer.

De esta forma vamos encontrando los valores de la presión a diferentes tiempos y al usar la
ecuación de productividad se puede saber el volumen de hidrocarburos que se está moviendo y
por ende conocer la producción de una región donde haya un pozo.

El siguiente paso sería tomar la sección del yacimiento que se quiere estudiar. Un ejemplo sería
la unidad 3.2 [Lydell] de campo FULMAR. En ella se tendrá que dividir en celdas usando una
malla. La malla puede ser de diferentes tipos:

a) Cartesiana b)Híbridas c)Radial

- 46 -
 Fig. 29 Tipos de malla que pueden usarse para simular matemáticamente un yacimiento

Existen más tipos de mallas pero estas son las más sencillas y básicas de las cuales se basan
las demás. Al tener la división de la sección en celdas por la malla. Se asignan nodos a cada
celda ya sea centrados o distribuidos.

Fig. 30 Arreglo de nodos centrados a la derecha y nodos distribuidos a la izquierda

La diferencia entre el arreglo de nodos distribuidos y centrados dependerá del ∆x que usemos en
la primera y última celda. Esta configuración o arreglo de los nodos servirá para delimitar las
fronteras del yacimiento. La única diferencia que se produce es al calcular los volúmenes de las
celdas ya que unas tendrán volúmenes más grandes que en la otra configuración.

Teniendo la sección del yacimiento dividida en celdas y conociendo los datos iniciales de
presión, propiedades del fluido y de la roca se puede empezar a calcular los valores de la
presiones en el siguiente paso de tiempo.

De esta forma nos va a quedar un sistema de ecuaciones cuyas incógnitas serán igual al número
de celda que tengamos en nuestra malla. La forma de resolver el sistema de ecuaciones es
usando matrices y algún método de algebra matricial siempre y cuando se conserve la linealidad
de las ecuaciones. Los métodos más usados son el de Thomas y LU. En caso de que sea un
sistema no lineal se calculará el Jacobiano y buscar que converja. Para después multiplicarlo por
la matriz y obtener la solución. La solución será una matriz que contenga los valores de la
presión a un paso de tiempo [∆t= conocido] indicado.

La forma general de la ecuación que genera el sistema de ecuaciones de la matriz es del

siguiente tipo:

- 47 -
 𝑆𝑆+1 + 𝑆𝑆𝑆𝑆𝑆+1 + 𝑆𝑆𝑆+1
𝑆𝑆 𝑆−1 𝑆𝑆+1 = 𝑆𝑆 …………(3.23)

𝑆 es la iesima propiedad de cada celda

𝑆ℎ
𝑆 = 𝑆𝑆 1 = ( )
𝑆− 𝑆𝑆 1 …………(3.24)
2 𝑆−2

𝑆𝑆𝑆(𝑆𝑆𝑆𝑆)𝑆
𝑆𝑆 = − (𝑆𝑆1 + 𝑆 1𝑆+ ) …………(3.25)
𝑆+2 𝑆−2 ∆𝑆

𝑆𝑆 = 𝑆𝑆 1 …………(3.26)
𝑆+2

𝑆𝑆𝑆(𝑆𝑆𝑆𝑆)𝑆 𝑆 𝑆
𝑆𝑆 = − ( 𝑆𝑆 + 𝑆𝑆 ) …………(3.27)
∆𝑆

Está forma es para un flujo lineal en una sola dirección ya sea x, y o z. Por eso tiene 3 términos
solamente. Conforma vamos agregando dimensiones como planos x-y,y-z o x-z se agrega un par
de variables al modelo mostrado arriba. En el caso de ser espacios en 3D se agregan otras dos
variables a la ecuación. Las variables van a estar separadas por 0 en la matriz. Serán tantos 0
como nodos en cada hilera del plano. En el caso de planos serán tantos 0 como nodos en el
plano. Que dando de la forma siguiente:

Fig. 31 Arreglo de nodos en un plano

𝑆+1 𝑆+1 𝑆+1
𝑆𝑆 𝑆𝑆−1,𝑆+ 𝑆𝑆𝑆𝑆+1,𝑆 + 𝑆𝑆𝑆𝑆,𝑆 𝑆𝑆+1 + 𝑆𝑆𝑆,𝑆+1
+ 𝑆𝑆𝑆,𝑆+1 𝑆𝑆+1 = 𝑆𝑆 …………(3.28)

En un plano x-y, y-z, x-z

𝑆+1 𝑆+1 𝑆+1 𝑆+1 𝑆+1 𝑆+1
𝑆𝑆𝑆𝑆,𝑆,𝑆−1 + 𝑆𝑆 𝑆𝑆−1,𝑆,𝑆 + 𝑆𝑆 𝑆𝑆+1,𝑆,𝑆 + 𝑆𝑆 𝑆𝑆,𝑆,𝑆 + 𝑆𝑆 𝑆𝑆,𝑆+1,𝑆 + 𝑆𝑆 𝑆𝑆,𝑆,𝑆+1 = 𝑆𝑆 …………(3.29)

- 48 -
 En un espacio x-y-z

Un ejemplo de la forma de la matriz que se generaría del sistema de ecuaciones que simula un
flujo a través de medios porosos en una fase y una dirección es el siguiente.

La siguiente parte sería resolver el sistema de ecuaciones para encontrar los valores de la
presión futura en cada celda como se mencionó anteriormente. Para encontrar la matriz solución
se pueden usar varios métodos. Los sistemas de ecuaciones pueden ser lineales o no lineales y
dependiendo del tipo también será el tipo de método que se use. Para los sistemas de
ecuaciones lineales se presentaran en la Tabla 6 varios métodos que pueden ser usados para
resolverlos. Los métodos directos son aquellos en los que se van eliminando variables del
sistema para facilitar los sistemas de ecuaciones. Los métodos iterativos son aquellos que
estiman el valor del parámetro. Estás estimaciones se pueden alejar o converger en un resultado
que satisfaga al parámetro de la ecuación. Los métodos iterativos también pueden usarse para
sistemas de ecuaciones no lineales.

Directos Iterativos
Gauss Jacobi
LU Gauss-Seidel
Cholesky SOR
Thomas Gradiente Conjugado
Tabla 8 Principales métodos parar resolver matrices de Sistemas de Ecuaciones Lineales

3.2.5 Realizar el programa en computadora que demuestre el modelo numérico

Para realizar el programa se busca que el simulador converja en una solución por lo que es
importante definir el tamaño de los pasos en el tiempo y que tan grande deben ser las celdas. Lo
que se va a calcular en él es la aproximación de las derivadas, las posiciones de los nodos y los
volúmenes de las celdas, los valores de las presiones futuras, así como los volúmenes
desplazados de fluidos. Al calcular los cambios de presión se calculan los cambios de volúmenes
y a su vez la producción de fluidos. Se pueden introducir los pozos en la malla como si fuera

- 49 -
celdas especiales donde la presión es mucho menor. Con esto se verá cómo se propaga la onda
de presión en las celdas colindantes. En este tipo de simulación sería una producción a presión
de fondo fluyente constante y lo que variaría sería el gasto del pozo.

Un ejemplo del algoritmo que se puede seguir es:

1. Lectura de la información
2. Definición de las condiciones iniciales
3. Definición de las condiciones de frontera
4. Definir un ∆t=1 segundo
5. Calcular los coeficientes ai, bi, ci y di
6. Resolver con la subrutina de Thomas
7. Reasignar las nuevas presiones (n+1) a (n): 𝑆𝑆 ⇐ 𝑆𝑆+1
𝑆 𝑆
8. Elegir otro ∆t. Por ejemplo uno más grande al anterior, ∆t = 1.2*∆t
9. Regresar al paso 5

(Arana, 2015)

Este podría ser un ejemplo de la forma en que se puede programar la solución. En el paso 6 se
puede usar cualquiera de los métodos mencionados en la Tabla 6 para resolver el SEL [Sistema
de Ecuaciones Lineales]. En el caso de que no fuese lineal y se tenga que encontrar la solución,
se tendrá que calcular la matriz del Jacobiano y hacer que converja para obtener la matriz
solución y seguir el mismo algoritmo que se presentó anteriormente.

El Jacobiano es una matriz donde se aproximaran las derivadas de cada una de las variables de
las que se quieran obtener la solución. También es conocida como Jacobiano la matriz de
derivadas parciales. Se empieza con la aproximación de:
𝑆0 = 𝑆𝑆
𝑆 𝑆

3.2.6 Validar el modelo

Teniendo el programa que resuelve la Ecuación de Difusión ya sea por la Solución Línea Fuente
o por Simulación Numérica de Yacimientos se comprueba si los resultados que se están
calculando son correctos. Es aquí donde se valida el modelo. La forma más sencilla de hacerlo
es con datos reales del campo o yacimiento. Se puede hacerlo con un historial de producción y
medir el error que está produciendo nuestro modelo. Si el error es muy grande se tiene que
hacer nuevas consideraciones para el modelo o en su defecto volver a hacer la deducción del
modelo matemático para encontrar los errores que se están cometiendo el error.

3.2.7 Ajuste del Historial de producción

Una vez construido el modelo se observa si es capaz de reproducir el comportamiento real del
yacimiento. Por lo que la información de entrada deben ser datos reales de producción. Después

- 50 -
se comparan las presiones calculadas y el movimiento de los fluidos con los datos reales del
yacimiento.

Los parámetros de entrada se modificaran hasta que se logre un ajuste aceptable entre el
yacimiento y el modelo. Se ajusta el modelo con el comportamiento real con el fin de aproximar
lo más posible la respuesta del yacimiento con la información disponible. Normalmente los
parámetros que más se ajustan son la permeabilidad del yacimiento para ajustar los gradientes
de presión y la distribución de saturaciones. Otros parámetros que se pueden ajustar son el
tamaño del acuífero, porosidad, grosor de la celda y permeabilidad. (Arana, 2015)

3.2.8 Predecir el comportamiento de yacimientos convencionales

Una vez que se comprueba que el programa está calculando datos que pueden ser tomados
como correctos se puede usar para hacer pronósticos de producción o estudiar otras variables a
futuro. Una de las variables que podríamos estudiar es la producción de agua para dimensionar
las instalaciones superficiales de tratamiento de agua.

De esta forma se ha concluido la explicación del Proceso de Caracterización Dinámica de

Yacimientos. Las páginas anteriores tiene el propósito de dar grosso modo la forma en que se
estudia el comportamiento de los yacimientos. Tomando en cuenta los factores que afectan el
rendimiento del pozo y la certidumbre con la que queremos contra al predecir el comportamiento
futuro. Se explicó el método de Caracterización Dinámica y el de Simulación Numérica. Pero hay
un tercero que es más sencillo y requiere de menos información y son las Curvas de Declinación.
La finalidad de escribir sobre estos temas es porque no siempre se podrán usar las curvas de
declinación para predecir el comportamiento de un pozo o del yacimiento. Además el tema no ha
concluido por que como se ha dicho en el tema del capítulo se caracterizan los yacimientos no
convencionales y hasta ahora sólo se ha hablado de los convencionales.

3.3 Caracterización dinámica de los yacimientos no convencionales tipo tight oil

A continuación se caracterizarán los yacimientos no convencionales con base en artículos e

información de fuentes bibliográficas. En las formaciones tight oil la roca generadora también es
la roca almacenadora, esto se debe a que al estar la roca madre en la ventana de generación del
aceite y debido a su extrema baja permeabilidad el aceite generado no puede migrar afuera de la
roca madre, generando zonas de presión anormalmente altas con baja saturación de agua.
Durante la etapa productiva del pozo la permeabilidad estará fijada por la permeabilidad de la
fractura y se tendrá el caso de flujo lineal y bilineal dependiendo de la geometría de los disparos
y el sistema de fracturas. Muy rara vez se verá un flujo radial tardío, pero puede ser el caso
cuando se haya alcanzado el flujo dominado por la frontera que será posible si el volumen de
roca estimulado es muy grande y después de muchos años de producción. Dicho en otras
palabras, será difícil encontrar flujo dominado por las fronteras (boundary dominated flow) debido
a la baja permeabilidad y que el régimen predomínate es el transitorio.

Para imaginar la geometría de flujo presente en área de drene se puede pensar en estos cuatro
elementos: la matriz o yacimiento, el sistema de fracturas, los disparos y el pozo, e interactúan

- 51 -
de la siguiente manera la matriz que aporta una cantidad de fluidos a la fractura y estos
circularan a través de las fracturas hacia los disparos y consecutivamente al pozo (flujo bilineal) o
el segundo modelo si el volumen de roca estimulado es muy grande, la permeabilidad de la
fractura es tan grande que la conductividad de la fractura es infinita y la presión a lo largo de la
fractura permanece constante dejando que los fluidos circulen a través de la fractura hacia pozo
(flujo lineal). Una observación importante aportada por el ingeniero Doung es que el sistema de
fracturas in-situ que ya había sido sellado durante la diagénesis se reactiva con el decremento de
la presión y modifica el régimen de esfuerzos, originando un aumenta la permeabilidad como es
el caso del a formación del Bakken (Duong, 2011). También se tiene que tener en mente que
aunque se tenga saturaciones de agua muy bajas y roca mojada por aceite al fracturar
hidráulicamente se inyectan grandes volúmenes de agua a la formación, como es el caso de la
formación del Bakken. Por lo que al comienzo de la producción se tendrán un corte de agua
producto del agua remanente del fracturamiento hidráulico, se menciona esto para evitar la
confusión con que exista un acuífero asociado que genere un, pero puede existir la posibilidad de
que se usen pozos inyectores de agua para generar un frente de barrido lo que aumentará
considerablemente el corte de agua conforme avance el tiempo.

Dos de las formaciones productoras de aceite tipo tight más estudiadas son el Bakken (Dakota
del Norte y Sur, Montana y Saskatchewan) y Barnett (Texas, Nuevo Mexico). La diferencia que
existe entre los yacimientos de tight oil de Barnett y del Bakken es el tipo de kerógeno que dio
origen a los hidrocarburos, el gradiente geotérmico y periodo apertura de la ventana de
generación. En el Barnett los campos contienen mayor aceite ligero y gas, lo que se define como
un yacimiento de gas y condensado, mientras que en el Bakken contiene aceite negro con gas
disuelto, lo que se define como un yacimiento de aceite negro de alto encogimiento. Bajo estas
premisas se pueden indagar cuales son los mecanismos de empuje dominantes y otra pista que
puede facilitarlo es el tipo de kerógeno que dio origen al hidrocarburo. En esta tesis se toma en
cuenta sólo el yacimiento del Bakken para idear hipótesis sobre los mecanismos de empuje
predominantes. Se comenzará con la idea que no se puede tener un acuífero por ser la roca
generadora que se encuentra en la ventana de generación y no existe migración al menos que
se encuentre en un sistema naturalmente fracturado, la baja saturación de agua y la falta de
contacto agua-aceite que se ve reflejada en los registros de pozos desecha el mecanismo de
empuje por acuífero asociado. El empuje por segregación gravitacional se puede tomar en
cuenta si se genera un casquete de gas y existe un acuífero, pero al tratarse de fracturas y no un
sistema poroso interconectado la idea de un casquete de gas es difícil. El colgamiento dentro de
las fracturas es el fenómeno que ocurre con naturalidad por la generación excesiva de gas
debido a la caída tan abrupta de presión que se origina dentro de la fractura y la vecindad del
pozo (Steven, 2016). Está documentado que el flujo multifásico a través de las fracturas está
presente cuando la presión de fondo fluyente cae por debajo de la Pb lo que origina el
colgamiento del aceite en la fractura además de un incremento promedio mayor a 4 veces el
valor común de la RGA en función del tiempo que se puede encontrar en un yacimiento
convencional.

- 52 -
 a)

b)
Fig. 32 a) Comportamiento de RGA vs 𝑆yacimiento convencional con un pozo vertical sin fracturar,
b) Comportamiento de un yacimiento no convencional con pozo horizontal fracturado (Stevens,
2016)

Descartando los posibles mecanismos de empuje existentes en la formación del Bakken se

puede decir que el mecanismo de empuje predominante es el gas en solución y encogimiento del
sistema roca fluido. El mecanismo por encogimiento del sistema roca fluido tomará más
importancia en éste yacimiento a diferencia de otros en los cuales es despreciado. Ya que los
yacimientos no convencionales se encuentran a un régimen de estrés mayor al convencional, lo
que origina que al disminuir las presión y dadas las propiedades de la compresibilidad de la roca
el cambio de volumen sea mayor consecuentemente el empuje producido también. Inclusive
como se puede leer en el artículo de “Rate-Decline Analysis for Fracture-Dominated Shale
Reservoirs” al disminuir la presión se reactivan las fallas o fracturas existentes dejando que las
fracturas reabiertas aporten más hidrocarburos al sistema de roca estimulado, esta suposición se
contrapone con lo escrito en el artículo “Microseismic Maps Production Volume” donde se
plantea la idea que con paso del tiempo el volumen activamente produciendo va a disminuyendo.
Si al realizar un fracturamiento hidráulico se tiene el Volumen de Roca estimulado (SRV
Stimulated Rock Volume) pero el Volumen Activamente Produciendo (APV Active Production
Volume) es diferente debido a que algunas de las microfracturas se cerraron al finalizar el
fracturamiento, entonces la producción va a disminuir con el paso del tiempo.

- 53 -
Fig. 33 Evolución del sistema de fracturas con el paso del tiempo dentro de los pozos horizontales
hidráulicamente fracturados (Vermilye y Sicking, 2016)

La compresibilidad del fluido va a ser afectada dado que las proporción de las saturaciones de
las fases presentes en yacimientos convencionales no son iguales a las de los no
convencionales. Las formaciones son normalmente mojadas por aceite, la saturación de agua es
muy pequeña y si vamos a encontrar gas en solución que probablemente sea una cantidad muy
grande dado que las formaciones están en la ventana de generación del aceite a presiones
anormalmente altas. Esto da como referencia que se puede encontrar más gas disuelto en el
aceite de lo que se encontraría en un yacimiento convencional, pero esto es una mera
suposición.

𝑆𝑆 = 𝑆𝑆𝑆𝑆 + 𝑆𝑆𝑆𝑆 +𝑆𝑆𝑆𝑆 …………(3.30)

La Compresibilidad de Total también cambiará. Dado que la ∆P tendrá un efecto más grande ya
que es una zona de presión anormalmente alta. Conforme disminuya más la presión se verá un
efecto más grande del empuje por encogimiento del sistema roca fluido.

𝑆𝑆 = 𝑆𝑆 + 𝑆𝑆 …………(3.31)

Está breve introducción tiene el objetivo de plantear hipotéticamente los cambios que
manifestarían las propiedades de los fluidos en un yacimiento no convencional de tipo tight oil
con un empuje por gas disuelto.

3.3.1 Mecanismos de empuje en yacimientos no convencionales tipo Tight Oil

Los mecanismos predominantes dentro de un yacimiento del tipo tight oil son empuje por gas
disuelto y el empuje por encogimiento del sistema roca fluido.

- 54 -
Los otros empujes: empuje por acuífero asociado, empuje por casquete de gas y el empuje por
segregación gravitacional, es posible que se generen pero dado a las características petrofísicas
y a las fases presente en el yacimientos o que su efecto sea tan pequeño que se puedan
descartar, al menos que se esté utilizando inyección de agua o un sistema de recuperación
mejorado. Sé puede descartar el empuje por acuífero asociado por no ser una cuenca
convencional en la que puedan migrar las fases también por encontrarse en la ventana de
generación del aceite y porque no existe una segregación de las fases dentro de la roca
generadora. El agua presente es residuo de las reacciones de síntesis que generan a los
hidrocarburos y representar una saturación muy pequeña pero considerable. Para que sea
tomado en cuenta como una fuerza que aporta energía al yacimiento el acuífero tiene que tener
10 veces mayor volumen que la fase de aceite. Pero en este caso ni siquiera forma una fase la
cantidad de agua presente en la matriz por lo cual se descartará este tipo de empuje.

El empuje por casquete de gas queda descartado debido a la reducida porosidad. Al ser un
yacimiento donde la porosidad inicial es muy pequeña y su permeabilidad es del orden de
décimas a centésimas de milidarcies, si durante la explotación se redujese la presión por debajo
de la Pb y comenzase a generar volúmenes gas, se debería de superar la saturación crítica del
gas para que este comenzara a moverse y el casquete comenzara a formarse, pero no podría
formar una fase de gas continua debido al reducido volumen que existe en la fractura y el
sistema de fracturas. Esto ocasiona que aunque exista un volumen de gas al no poder consolidar
una fase por arriba del intervalo productor debido al espacio reducido de las fracturas este no
aportara un empuje considerable. En vez de eso como se explica en los artículos ¨Multiphase
linear flow in tight oil reservoirs¨ y ¨Producing gas/oil rate behavoir of multifractured wells in tight
oil reservoirs¨ dentro de las fracturas se presentara el fenómeno del colgamiento y lo que origina
que le RGA sea mucho mayor que en un yacimiento convencional.

El empuje por segregación gravitacional quedará descartado por no existir una fase acuosa
considerable que aporte trabajo a los fluidos para que sean expulsados por el pozo y a su vez no
existe una fase gaseosa considerable que produzca el mismo efecto que la fase acuosa. Dicho
en otras palabras el acuífero asociado y el casquete de gas pueden dar origen al efecto de
segregación gravitacional y al no existir dentro de los yacimientos como el Bakken se pueden
descartar de los empujes que aportan al sistema integral de producción.

Este es un caso hipotético que se basa de lo investigado en distintos artículos de la SPE. Puede
existir la posibilidad que al tener un volumen de roca estimulado muy grande y además el
volumen de roca activo produciendo sea parecido al SRV al implementar un sistema de
recuperación como inyección de agua se generé artificialmente el empuje por acuífero asociado
y se pueda calcular el aporte de presión en función del gasto de inyección, la ubicación del pozo
y la difusividad hidráulica del medio.

Como se puede ver los yacimientos de tipo tight oil experimentan disminuciones de presión
cuando comienza la producción debido a que tiene muy pocos mecanismos de empuje que
aporten energía y que su volumen interconectado depende sólo del volumen de roca estimulado.

- 55 -
3.3.2 Regímenes de flujo de en yacimientos no convencionales tipo Tight Oil

Los regímenes de flujo encontrados en yacimientos tipo tight oil son: régimen transitorio y
régimen pseudoestacionario. El régimen transitorio es el más común en los yacimientos de baja
permeabilidad debido a la poca difusividad hidráulica. Esto es causado por las permeabilidades
del orden de nano Darcies, en casos optimistas centésimas de mili Darcies y la porosidad del
medio puede tener valores promedio entre 9-15% dando como resultado que la difusividad que
es la velocidad con la que se propagan las ondas de presión en el medio sea muy baja. Por lo
que tardara mucho más tiempo en que los cambios de presión alcancen las fronteras del
yacimiento. Es por esta razón que la mayoría de los pozos produciendo de este tipo de
yacimientos permanecen en régimen transitorio durante la mayor parte de su vida productiva.

El régimen pseudoestacionario puede ser encontrado en secciones donde el volumen de roca

estimulado es muy grande y han pasado varios años del comienzo de la producción del pozo.
Después de tanto tiempo se ha alcanzado el flujo dominado por fronteras y puede empezar el
régimen pseudoestacionario. No puede asegurarse el encontrarse este régimen aunque se
cumplan con las 2 características mencionadas, ya que muchos yacimientos permanecen en
régimen transitorio la mayor parte de su vida productiva.

En el artículo ¨Producing-Gas/Oil-Ratio Behavior of Multifractured Horizontal Wells in Tight Oil

Reservoirs¨ se propone un diagrama conformado de 4 etapas en la vida productiva de un pozo
de un yacimiento tipo tight oil usando la RGA cono indicador (Stevens, 2016).

En la Fig. 34 se puede ver un historial de producción contra su valor de RGA de un pozo tight oil
de la formación Meramec que produce de la cuenca Anadarko en el estado de Oklahoma. Este
tipo de comportamiento es típico de un yacimiento no convencional con mecanismo de empuje
por gas en solución y extremadamente bajas permeabilidades que produce por un sistema de
fracturas inducido.

- 56 -
Fig. 34 Historial de producción de un pozo tipo tight oil de la formación Meramec de la cuenca
Anadarko. (Modificado de Stevens, 2016).

Para identificar con mayor facilidad las etapas características de un pozo tipo tight oil en la Fig.
35 se ven las 4 etapas: RGA=Rsi, alza de RGA, plano transitorio de RGA, aumento de la RGA
debido al FDF. Hay que tomar en consideración que dependiendo de la presión inicial y el
régimen de explotación se alcanza más rápido la Pb que marcará un parte aguas en el
comportamiento de la RGA. Es muy importante identificar las siguientes etapas porque están
relacionadas con las geometrías flujo presentes. El plano transitorio de RGA es la etapa donde la
Pwf ha alcanzado su mínimo y se mantiene constante. Este también ocurre porque los valores
promedio de presión, saturación y permeabilidades relativas se mantienen constantes durante el
flujo transitorio lineal ya que la Pwf es constante (Steven, 2016).

- 57 -
Fig. 35 Etapas de producción de un pozo horizontal con varias etapas de fracturamiento (modificado de
Stevens J., 2016, pg.5)

Etapa Nombre Etapa Descripción

Mientras la Pwf>Pb yacimiento se mantiene bajo
1 RGA = Rsi
saturado (RGA=Rsi)
Cuando se reduzca la presión por debajo de la Pb y se
2 Alza de RGA debido Pwf < Pb alcance la Sgc , el gas comenzará a fluir e incrementará la
RGA
Al alcanzarse el régimen transitorio lineal o bilineal
3 Plano Transitorio de RGA todas las propiedades se mantendrán constantes debido
a que la Pwf es cte.
Conforme la presión disminuya más gas en solución se
Alza RGA debido a flujo
4 liberará y producirá colgamiento en las fracturas, lo que
dominado por fronteras
incrementara la RGA. La RGA estará en función de VPA
Tabla 9 Etapas de producción de MFHW usando el RGA para identificarlas (Modificado de Steven,
2016)

3.3.3 Geometrías de flujo de yacimientos tipo Tight Oil

Las geometrías de flujo más comunes en los pozos MFHW (Multifractured Horizontal Well) de
yacimientos de a aceite negro con mecanismo de empuje por gas en solución son el lineal y
bilineal. También pueden encontrarse otras geometrías de flujo como: flujo lineal compuesto,
flujo lineal temprano, flujo radial temprano, flujo elíptico, flujo radial tardío. En éste capítulo se
enfocará sólo en el flujo lineal y bilineal pero también se explicara a groso modo las demás
geometrías de flujo mencionadas presentes en MFHW. El diseño de las etapas de

- 58 -
fracturamiento y parámetros como espaciamiento, longitud y permeabilidad de la fractura
determinarán que tipo de geometría de flujo se puede encontrar. El factor de conducción
adimensional 𝑆𝑆𝑆 se usa para calcular la conductividad de la fractura y así definir si es una
fractura con conductividad finita o infinita.
𝑆𝑆
𝑆𝑆 =
𝑆
… 3.13
𝑆
𝑆 𝑆
𝑆𝑆

𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 [𝑆𝑆], 𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆[𝑆𝑆]

𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆 [𝑆𝑆], 𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 [𝑆𝑆]

El flujo bilineal ocurre cuando la conductividad de la fractura es finita. Lo que quiere decir que la
fractura presentará un perfil de presión del comienzo a la cola de fractura debido al cambio de
permeabilidad. En este caso se presentarán 2 tipos de flujo lineal: el primero de la matriz a la
fractura y otro de la fractura al pozo. Esta geometría es común en fracturas largas donde es difícil
que el apuntalante se distribuya uniformemente cuando se cierra la fractura al disminuir la
presión después haberse alcanzado la presión de fractura. También es típica esta geometría de
sistemas naturalmente fracturados. En una gráfica de diagnóstico especializada se puede
apreciar con una pendiente de ¼.

El flujo lineal es característico de pozos hidráulicamente fracturados donde las fracturas tienen
conductividad infinita. Lo que quiere decir que la permeabilidad de la fractura es muy grande en
comparación con la κ de formación. Esto produce que la presión se mantenga constante a lo
largo de la fractura. En una gráfica de diagnóstico especializada se puede apreciar con una
pendiente de ½.

Para ver la diferencia que habría entre un perfil de presión en un yacimiento con un pozo vertical
y un pozo horizontal fracturado con las mismas propiedades de fluido, condiciones y propiedades
de la matriz se muestra la siguiente figura. En la cual las fracturas tienen una conductividad
infinita. El ejemplo es tomado de una simulación.

- 59 -
Fig. 36 Perfiles de presión simulados en un yacimiento de aceite de 5 mD con flujo monofásico
después de producir 30 dias a una Pwf=1,500 psia para un flujo radial y lineal. La gráfica de la
derecha tiene una escala logarítmica en la distancia.

Flujo lineal compuesto se produce una vez que las ondas de presión han alcanzado a las
fracturas más próximas y puede existir flujo de las regiones exteriores que no han sido
estimuladas. Una vez que las fracturas han interferido unas con otras puede ser visto el flujo de
las regiones exteriores a la zona estimulada por las fracturas. (Fekete, 2014)

Flujo lineal temprano es un tipo común de flujo producido en un pozo estimulado con
fracturamiento, se produce cundo los transitorios de las fracturas se han estabilizado. (Fekete,
2014)

El Flujo radial temprano puede ser visto después del término de flujo lineal temprano, pero antes
de que las fracturas empiecen a interferir una con otra. Sólo puede ser visto cuando etapas de
fracturamiento están muy alejadas o cuando la longitud media de la fractura es corta. (Fekete,
2014)

Flujo Elíptico se produce cuando ha comenzado flujo del yacimiento hasta el segmento final del
pozo. Es una geometría de transición del flujo lineal al flujo radial en los pozos horizontales. Se
puede ver en una gráfica de diagnóstico con un pendiente de 1/3. (Fekete, 2014)

Flujo radial tardío define flujo radia alrededor de un pozo horizontal fracturado después del flujo
lineal compuesto. Se caracteriza por tener una pendiente de 0 en la gráfica de diagnóstico. Va a
ser visto el flujo radial tardío si sólo existe un pozo en un yacimiento que no se ha desarrollado
aún y requerirá de mucho tiempo y área en el yacimiento tipo tight oil. (Fekete, 2014)

- 60 -
Fig. 37 Gráfica diagnóstico donde se pueden ver las geometrías presentes. Simulación de un caso
con una fractura k=150 nd, xf=240 ft, espaciamiento de la fractura= 100 ft y kfwf=8md-ft. Pressure
𝑆
normalized Oil Ratio, 𝑆𝑆𝑆 = 𝑆 . (modificado de Stevens, 2016)
𝑆𝑆 −𝑆𝑆𝑆

3.3.4 Propiedades de los fluidos de los yacimientos no convencionales

Esta sección se basa en las hipótesis del punto 3.3 sobre el yacimiento del Bakken que tiene
una 𝑆𝑆 mucho menor que los yacimientos tipo tight oil de la cuenca Anadarko, su 𝑆𝑆𝑆 promedio
supera el 10%, donde puede alcanzar valores de 20-30% (Stevens, 2016, pg. 12) y el yacimiento
presenta características de aceite negro de alto encogimiento, en la región se va a tener una
producción de gas mucho menor a comparación de Anadarko.

En esta sección se tratara acerca de las propiedades que tiene un comportamiento particular a
diferencia del comportamiento en yacimientos convencionales.

Una de las propiedades es la compresibilidad, se toma como referencia del artículo ¨Multiphase
Linear Flow in Tight Oil Reservoirs ¨para explicar cómo se modela la compresibilidad cuando
existen 2 fases. En el artículo se propone un modelo matemático para simular y entender el flujo
multifásico en yacimientos tipo tight oil. El yacimiento del Bakken al igual que la formación
Meramec presentan flujo multifásico porque se alcanza la 𝑆𝑆 muy rápido debido a que la 𝑆𝑆𝑆 tan
bajas que se tiene que producir para lograr el gasto deseado. La diferencia entre el Bakken y la
formación Meramec es una Sgc mucho mayor en el Bakken y una Sg menor en el Bakken. Lo que
origina que la mayoría de la producción no convencional del Bakken sea crudo de los 20-55°API.

𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 1 𝑆𝑆

𝑆∗ = − − + …………(3.32)
𝑆
𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆 𝑆𝑆

- 61 -
 𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 Unidad de gas liberada por volumen poroso a
1 condiciones de yacimiento mientras la presión disminuye
𝑆𝑆 𝑆𝑆
𝑆𝑆 𝑆𝑆𝑆
2 − Compresibilidad del Gas
𝑆𝑆 𝑆𝑆
𝑆𝑆 𝑆𝑆𝑆
3 Compresibilidad del Aceite
𝑆𝑆 𝑆𝑆
1 𝑆𝑆
4 Compresibilidad del Poro
𝑆 𝑆𝑆
Tabla 10 Descomposición de la compresibilidad total

Fig. 38 Compresibilidad total del sistema y sus componentes, Eq. 3.32 (Modificado de Tabatabaie,
2016. pg 4)

Cabe recalcar que en zonas de presión anormalmente alta las compresibilidad de poro puede
aumentar dado a que el régimen de esfuerzo al que está sujeto la matriz es mayor, dando como
resultado un efecto mayor al ir disminuyendo la presión durante la explotación.

En los pozos tipo tight oil dado que se produce bajo régimen transitorio las propiedades de S o y
RGA no están en función de la permeabilidad absoluta ni las permeabilidades relativas.
(Tabatabaie, 2016). Es por esto que para casos de simulación se puede representar la relación
𝑆𝑆 vs 𝑆 de la siguiente manera.

- 62 -
 𝑆𝑆𝑆 𝑆 1
𝑆𝑆(𝑆) = 1 + 𝑆𝑆𝑆 𝑆𝑆𝑆( ) ∫ 𝑆𝑆 …………(3.33)
𝑆𝑆 𝑆𝑆𝑆 𝑆 𝑆
𝑆 𝑆

La 𝑆𝑆𝑆 durante el plano transitorio de RGA puede ser expresado de la siguiente manera.
𝑆𝑆𝑆𝑆𝑆
𝑆𝑆 …………(3.34)
𝑆𝑆𝑆 = 𝑆𝑆 +
𝑆𝑆𝑆 𝑆𝑆𝑆𝑆

La Relación de Solubilidad 𝑆𝑆𝑆 depende de las propiedades químicas de la mezcla y está en

función de la 𝑆𝑆𝑆.

En el siguiente caso simulado (Figura 39) sobre un bloque modelo con una fractura en régimen
transitorio que presenta flujo multifásico se puede ver cómo se comporta la Saturación del aceite.
El siguiente caso toma como permeabilidad inicial 𝑆𝑆 = 0.1 𝑆 0.01 md respectivamente. Las
permeabilidades relativas son calculadas con funciones de Corey:
𝑆𝑆
𝑆𝑆𝑆 = 𝑆∗𝑆𝑆
𝑆
𝑆 …………(3.35)
𝑆𝑆
𝑆𝑆𝑆 = 𝑆 (1 − 𝑆 )
∗
𝑆𝑆 𝑆 …………(3.36)

Fig. 39 Caso en el que se basa la simulación, (a) esquema del modelo de la fractura, (b)
propiedades del fluido. (Tabatabaie, 2016)

En la Fig. 40a se ve cómo varia la 𝑆𝑆 a lo largo de la fractura con el paso del tiempo y en la Fig.
40b se ve como varía la 𝑆𝑆 a diferentes 𝑆𝑆. Es importante observar el comportamiento de la 𝑆𝑆
ya que se relacionado con la 𝑆𝑆𝑆.

- 63 -
Fig. 40 (a) Gráfica 𝑆𝑆 vs x con incremento del tiempo tomado de un caso simulado fig.39,
(b) Gráfica 𝑆𝑆 vs 𝑆𝑆 tomada de caso simula fig. 39 (Tabatabaie, 2016)

Las propiedades como el los factores de volumen 𝑆𝑆, 𝑆𝑆, 𝑆𝑆 tendrán el mismo comportamiento
que en cualquier yacimiento al igual que las viscosidades 𝑆𝑆, 𝑆𝑆.

- 64 -
 CAPÍTULO 4

CURVAS DE DECLINACIÓN

4.1 Pequeña Introducción

Los métodos empíricos son usados ampliamente en el ramo de la ingeniería para predecir o
estimar. Son empíricos porqué no tienen una estricta relación matemática o relación física con el
fenómeno que se está prediciendo y porque no dan valores exactos sino valores que se acercan
mucho al real. Las ventajas de usarlos es que existen fenómenos físicos formados de muchos
subfenómenos los cuales siendo modelados por expresión matemática sería muy complejos de
calcular pero usando un método empírico es práctico y certero. Los métodos empíricos son
usados para modelar fenómenos físicos de mezclas heterogéneas como el aceite el cual
contiene agua, distintos tipos de hidrocarburos, gas, acido sulfúrico y nítrico, dióxido de carbono
en fase líquida y gaseosa, en resumen un mezcla conformada de distintos componentes cada
uno con comportamientos diferentes con relación a la presión y temperatura. Un ejemplo de
métodos empíricos son las curvas de declinación o las correlaciones de flujo multifásico. A
continuación se explicara más sobre ellas, y así saber cómo funcionan las curvas de declinación.
Antes del desarrollo de las curvas de declinación se usaban métodos volumétricos para calcular
el desempeño futuro de los pozos. El Balance de Materia es un ejemplo de método volumétrico.
Los métodos volumétricos y empíricos son usados en la actualidad para estimar la producción
dependiendo de las condiciones del yacimiento cuando se desea explotar. Los parámetros más
importantes para las curvas de declinación son el número de pozos y su ubicación pero debe de
tomarse un pozo real de referencia que produzca de un yacimiento en condiciones similares a las
del yacimiento que se desearía explotar para usarse de referencia. A esta técnica se le conoce
como analogía de pozos. Es muy importante tomar en cuenta la anisotropía y heterogeneidad del
yacimiento. La anisotropía es la variación que tiene una propiedad dependiendo de la ubicación
dentro de un estrato formado del mismo material o simplemente la variación de una propiedad
con respecto de la ubicación. La heterogeneidad es la variación o cambio de litología que puede
existir en un volumen determinado como un a formación o un yacimiento. Tomando en cuenta la
anisotropía, heterogeneidad y teniendo buenos valores de referencia se pueden hacer
estimaciones probables de la producción futura cuando se planea el proyecto de explotación.
Para usar los métodos empíricos se debe de contar con buenos datos de entrada o mediciones
del yacimiento ya que la certeza del método empírico dependerá del método y los datos de
entrada. La forma de obtener está información es por métodos directos o indirectos de medición.
Un método directo son los registros petrofísicos de pozo o los datos de núcleos porque son
tomados directamente del pozo y un ejemplo de método indirecto son los métodos geofísicos ya
que requieren de una interpretación, procesamiento de datos y no son medidos in-situ.

- 65 -
Los registros petrofísicos han abierto la puerta para hacer uso de ecuaciones y métodos más
complejos para calcular las condiciones dinámicas del pozo. Cuando se carece de datos
suficientes para usar estos métodos es más conveniente y fácil hacer uso de las curvas de
declinación. Pero antes se tendrá que suavizar los datos del historial de producción para evitar
cambios abruptos entre los puntos.
Conrad Schlumberger en 1927 tomó el primer registro petrofísico de resistividad en un pozo de
aceite. Antes de este desarrollo los científicos e ingenieros se basaban en los núcleos y recortes
que llevaban a superficie durante la perforación para evaluar las propiedades de las rocas de la
formación. (Burakovsky, 2012, pg 184)
El avance tecnológico de las herramientas usadas para hacer registros y estudios indirectos del
yacimiento ha permitido hacer mejores mediciones de las propiedades físicas y sus
distribuciones en la formación productor lo que ha hecho más sencillo hacer buenas
estimaciones al usar los modelos de balance de materia o curvas de declinación por lo se
hacen buenos pronósticos del flujo volumétrico que tendría un pozo efecto del abatimiento de
presión correspondiente y así calcular capacidades de las instalaciones superficiales o diseño del
sistema integral de producción, como el diámetro del estrangulador, diámetro de la tubería de
producción, perdidas de presión en la tubería, etc.
El gasto y la presión son 2 de las variables más importantes para el ingeniero de producción, ya
que junto con las condiciones del yacimiento como temperatura 𝑆, grados API, viscosidad del
aceite muerto 𝑆𝑆𝑆 se usan correlaciones, primero para estimar la 𝑆𝑆 y definir propiedades de
aceite saturado o bajo saturado y consecutivamente para calcular las propiedades como la
relación de solubilidad 𝑆𝑆 , densidad relativa del gas 𝑆𝑆𝑆, factor de volumen del aceite 𝑆𝑆,
viscosidad del aceite saturado 𝑆𝑆 a condiciones de yacimiento. Algunas de las correlaciones más
usadas para estimar son: Standing, Lasater, Vázquez, Kartoatmodjo, Glaso para la viscosidad,
Vazquez y Beggs.
Como ejemplo la correlación de Oisten establecida utilizando muestras de aceite producido en el
Mar del Norte, donde predominan los aceites de tipo volátil:
1.- Calcule 𝑆∗ con:

log 𝑆∗ = −2.57364 + 2.35772 log 𝑆 − 0.703988 log2 𝑆 + 0.098479𝑆𝑆𝑆3𝑆 …………(4.1)

2.- Calcule 𝑆𝑆 con:

1
𝑆∗ 𝑆𝑆𝑆0.989 0.816 …………(4.2)
𝑆𝑆 = 𝑆𝑆𝑆 ( )
𝑆𝑆

donde:

- 66 -
 𝑆 = 0.130, para aceites volátiles.
𝑆 = 0.172, para aceites negros. (Garicochea, 1991)

4.2 Antecedentes de la curvas de declinación

En 1908 se tiene la referencia más temprana por Arnold y Anderson que expresaban la
declinación de la producción como una fracción de la producción pasada. Ésta declinación
normalmente era expresada en porciento por mes.
Durante los años de 1915 a 1921 se realizó mucha investigación relacionada con las curvas de
declinación y fue acumulada en el Manual para la Industria del Aceite y Gas (Manual for Oil and
Gas Industry).
J. O. Lewis y C. N. Beal recomendaron el uso de las curvas de declinación porcentual donde la
producción a cierto periodo de tiempo es expresada como un porcentaje de la producción de la
primera unidad de tiempo.
W. W. Cutler remarcó que 1924 que es erróneo asumir que le porcentaje de declinación
permanezca constante y una relación lineal en con una escala semilogarítmica es muy
conservadora para hacer un pronóstico. En su opinión una relación lineal en escala logarítmica
es mejor y más segura, pero algunos cambios horizontales son necesarios. También recomendó
que se usaran familias de curvas de declinación con base en datos de gastos en el tiempo de un
número de pozos de la zona en vez de construirse las curvas usando datos estadísticos o
soluciones gráficas.
En 1925 C. S. Larkey demostró como el método de mínimo cuadrados puede ser aplicado
exitosamente para curvas de declinación que pertenezcan al comportamiento exponencial o
hiperbólico y H. M. Roeser demostró que resultados igualmente confiables pueden ser obtenidos
se usa el método de intento error para determinar las constantes. Roeser fue el primero en
publicar la relación entre la producción acumulada y el tiempo para la declinación hiperbólica.
C. E. Van Orstrand investigó el pico de producción precedido de una declinación y vio que los
mejores resultados eran obtenidos usando la curva tipo:
𝑆 = 𝑆𝑆 𝑆 𝑆 −𝑆𝑆 …………(4.3)

En 1927 R. H.Johnson y A. L. Bollens introdujeron el método de radio de pérdida que es un

método estadístico donde se tabulan los datos de producción a intervalos de tiempo iguales en la
primera columna, en la segunda columna se escribe la diferencia entre gastos y en la tercera
columna se escribe el ratio de cambio de las diferencias. Éste método después de suavizarse un
poco los datos demostraba que los radios eran constantes para el caso de declinación
exponencial o geométrica o un cambio constante en los radios para el caso de declinación
hiperbólica o armónica.

- 67 -
En 1928 H. N. Marsh demostró que al grafica la producción acumulada contra el tiempo en un
plano cartesiano se puede ver una relación que se acerca mucho a la lineal. Éste
comportamiento es típico de declinación exponencial.
En 1931 R. E. Allen mencionó 4 tipos de declinación y los clasificó de acuerdo a su relación
matemática. Los tipos de declinación son:
1.- Aritmética o declinación de disminución constante
2.- Geométrica o declinación exponencial
3.- Harmónica o declinación isotérmica
4.- Básica o declinación de potencias fraccionales
En 1943 C. H. Rankin demostró como la presión de fondo puede ser utilizada como un sustituto
del gasto de producción de la curva de producción acumulada en zonas licitadas. Este método
sólo aplica para yacimientos sin empuje por acuífero y donde los índices de productividad son
constantes.

En 1942 P.J. Jones sugirió que una relación lineal puede ser aproximada en una gráfica con
escala logarítmica para pozos declinando con gastos variables y corresponde a la siguiente
ecuación:

log 𝑆 = log 𝑆0 − 𝑆 log 𝑆 …………(4.4)

Donde 𝑆0 es la declinación inicial y 𝑆 es una constante positiva. Al integrar la función lleva a la

ecuación de gasto tiempo de la forma:

𝑆0 𝑆 1−𝑆
𝑆 = 𝑆0𝑆100(𝑆−1) …………(4.5)

Está función va a demostrar un comportamiento lineal en una gráfica con escala logarítmica
donde se grafique contra el logaritmo del tiempo. (Arps, 1944)

4.3 Curvas de declinación

Las curvas de declinación es un método empírico que surge por la necesidad por saber el
comportamiento futuro de un pozo. Es empírico ya que se ajusta una curva tipo calculada por
uno de los tres modelos establecidos por Arps, se extrapola el comportamiento futuro del pozo
hacia un dominio desconocido (eje de las y´s) hasta que intersecte con el límite económico para
saber el tiempo de producción y la producción acumulada.
Para ajustar la curva se tiene que graficar los puntos generados por los conjuntos de la
producción en función del tiempo en una gráfica logarítmica o semilogarítmica y ver en cual de
las dos el comportamiento es más parecido a una línea con menores cambios de pendiente
(fluctuaciones). Para evitar las fluctuaciones es conveniente suavizar la curva por lo que se
puede usar el promedio de varios puntos colindantes. Como se ve en la figura 41. Los datos de

- 68 -
producción del pozo 16/26 A9 del Campo Alba presentan cambios abruptos en su producción por
lo que se uso el promedio, primero de 3 puntos que estuvieran alrededor del punto a graficar y
después para suavizarla más se hizo el promedio de 5 puntos incluyendo el punto a graficar.
Estos métodos para suavizar la curva ayudan a ajustar mejor la recta para saber cual es su
pendiente y así definir el coeficiente de declinación en caso de que no se quiera ajustar una
curva tipo. Un dato importante del Campo Alba que comenzó a explotar en el año de 1994
produce crudo pesado y fueron los pozos de desarrollo del campo fueron los primero en hacerse
horizontalmente además de usarse pozos inyectores de agua de mar para mantener la presión
en el acuífero limitado. Dadas las condiciones del crudo explotado (20 °API, 𝑆𝑆 = 7 𝑆𝑆) se
instalaron Bombeo Electrosumergible como sistema artificial de producción en todos los pozos
productores (Jayasekera, 1999). Esto se ve reflejado en la Fig. 41 donde es difícil apreciar una
declinación ya que la producción durante los primeros años se mantiene casi constante. Pero a
partir del mes 40 es posible apreciar una declinación del pozo.

Fig. 41 Historial del producción del Pozo 16/26 A9. Se grafica la producción normal y el promedio
de 3 puntos y puntos para poder ajustar mejor una línea de tendencia. Para más información
consultar el anexo B para saber cómo se llegó a la gráfica.

Los 3 tipos de declinación que definió Arps son:

 Exponencial
 Armónico
 Hiperbólico

- 69 -
Para deducir los siguientes comportamientos se asumió que el yacimiento no cuenta con un
empuje por acuífero asociado y la presión es proporcional al volumen de aceite en el yacimiento.

Las curvas de declinación parten de la ecuación diferencial:

𝑆𝑆/𝑆
𝑆
= −𝑆𝑆𝑆 …………(4.6)
𝑆

Donde dependiendo del valor de coeficiente b se generara el comportamiento:

𝑆 = 1 , Armónico

0 < 𝑆 < 1 , Hiperbólico

𝑆 = 0 , Exponencial

En los últimos años se han usado coeficientes mayores 𝑆 > 1 para pronosticar el
comportamiento de los yacimientos no convencionales de gas, pero se ha visto que las curvas
generadas representan comportamientos muy conservadores y optimistas. Siguen siendo
buenas aproximaciones pero se han generados nuevos modelos como el de Valko, Ilk o SEPD
(Streched Exponential Production Decline) para generar curvas más ciertas.

4.3.1 Declinación Exponencial

Es la más sencilla de identificar y se caracteriza por mostrar un comportamiento lineal al graficar

en escala semilogarítmica el gasto contra el tiempo. El coeficiente de declinación puede ser
calculado con la pendiente de la curva que se ajusta a los puntos del historial de producción.
Otra forma de identificarse es al graficar la producción acumulada contra el gasto en escala
normal y observar que hay un comportamiento lineal.

Fig. 42 Comportamiento típico de la declinación exponencial a) en gráfica semilogarítmica q vs t b)

gráfica de producción acumulada escala normal (New Nexico Tech, 2010)

- 70 -
La declinación exponencial también es conocida como ¨geométrica¨ o ¨porcentaje constante¨ se
debe a que la caída de la producción por unidad de tiempo es proporcional al gasto. La
declinación exponencial tiene un radio de cambio constante que se puede expresar con la
siguiente ecuación diferencial:
𝑆
= −𝑆 …………(4.7)
𝑆𝑆
𝑆𝑆
Resolviendo la ecuación diferencial y asignando los valores de frontera se llega a la expresión
para conocer el gasto en el tiempo. Para ver la deducción de la fórmula puede consultar el Anexo
C.

𝑆(𝑆) = 𝑆𝑆 𝑆−𝑆/𝑆 …………(4.8)

Para obtener la expresión de ecuación de la producción acumulada sólo se debe de integrar una
vez más con función del tiempo y se obtiene la siguiente expresión:

𝑆𝑆 = 𝑆(𝑆𝑆 − 𝑆) …………(4.9)

4.3.1.1 ¿Cómo calcular el coeficiente de declinación 𝑆 de forma analítica?

Se sugiere que primero se grafique la producción en escala logarítmica contra el tiempo y ya que
sea seguro que es un comportamiento de declinación exponencial se trace una línea que se
ajuste lo mejor al comportamiento general de todos los puntos. En la fig. 43 se ve una sección
del historial de producción donde hay un comportamiento de declinación exponencial. Se
muestra sólo es sección ya que los pozos del campo cuentan con bombeo electro sumergible y
pozo inyectores de agua lo que ocasiona que la presión se mantenga y no sea posible apreciarse
algunos de los tipos de declinación. Ya se ha trazado la línea de tendencia que se ajusta mejor a
los puntos y se tomaron dos de los puntos por donde pasa la recta 𝑆1(𝑆1 = 49764, 𝑆1 = 32 ) y
𝑆2(𝑆2 = 35428, 𝑆2 = 45 ). En este caso se va generar una curva de declinación para el dominio
de tiempo de [32, 54] (meses).

- 71 -
Fig. 43 Sección del historial d producción de pozo 16/26 A9 con línea de tendencia para ejemplificar
como obtener el coeficiente b

Después usando la siguiente ecuación se puede conocer el coeficiente de declinación 𝑆 :

ln(𝑆2) −ln(𝑆1)
𝑆= …………(4.10)
𝑆1 − 𝑆2

Dados los puntos 𝑆1 𝑆 𝑆2 , 𝑆 es igual a 𝑆 = 0.02613 . Por ende nuestra curva de declinación
asumiendo que el para esta sección es 𝑆𝑆 = 49764 𝑆𝑆𝑆 es:

𝑆(𝑆) = 49764 𝑆−0.02613(𝑆−32) …………(4.11)

Se usa el término -32 debido que el conteo empieza a partir del mes 32. La curva de declinación
se ajusta como se puede ver en la fig. 44.

- 72 -
Fig. 44 a) Historial de producción en escala semilogarítmica con curva de declinación ajustada b)
historial de producción en escala normal con curva de declinación ajustada

4.3.2 Declinación Hiperbólica

También es conocida como declinación tipo ¨log-log¨ dado que al gráfica el gasto contra el
tiempo en una escala logarítmica se muestra un comportamiento lineal. Para encontrar una
variable constante en la declinación exponencial se hace a partir de las segundas diferencias
dividiendo la primera derivada entre el cambio de las diferencias de los gastos. Está constante se
puede expresar en la siguiente ecuación diferencial:

𝑆
𝑆( )
𝑆𝑆
𝑆𝑆 = −𝑆 …………(4.12)
𝑆𝑆

Resolviendo la ecuación diferencial y asignando los valores de frontera se llega a la expresión

para conocer el gasto en el tiempo. Para ver la deducción de la fórmula puede consultar el Anexo
D.

− 1
𝑆𝑆 𝑆
…………(4.13)
𝑆 = 𝑆𝑆 ( + 1)
𝑆0

En una historial de producción se puede identificar el comportamiento de declinación Hiperbólico

como puede observarse en la fig. 45.

- 73 -
Fig. 45 Comportamiento típico de la declinación Hiperbólica grafica con escala logarítmica (New
Nexico Tech, 2010)

La expresión para calcular la producción a cumulada debido a un comportamiento hiperbólico es

la siguiente:

𝑆 0𝑆 𝑆 𝑆
𝑆𝑆 = (𝑆 1−𝑆 −𝑆1−𝑆) …………(4.14)
1 −𝑆 𝑆

4.3.3 Declinación Armónica

La declinación armónica es un caso especial de la declinación hiperbólica. La velocidad con la

que decrece la producción es proporcional al gasto del pozo. Es muy raro que suceda éste tipo
de comportamiento. Se genera cuando el coeficiente de declinación 𝑆 = 1 y se puede identificar
en el historial de producción como una tendencia lineal al graficar la producción contra el tiempo
con una escala logarítmica también al graficar la producción acumulada con el gasto en escala
logarítmica.

Fig. 46 Comportamiento típico de declinación armónica a) historial de producción escala

logarítmica b) Producción acumulada gasto en escala logarítmica (Modificado de New Mexico
Tech, 2010)

La curva se genera usando la siguiente ecuación:

- 74 -
 𝑆𝑆
𝑆= 𝑆
( + 1) …………(4.15)
𝑆0

La producción acumulada se puede conocer usando la siguiente ecuación:

𝑆𝑆 = 𝑆0𝑆𝑆(log(𝑆𝑆) − log(𝑆)) …………(4.16)

4.3.4 Ejemplo de curvas de declinación

La fig. 47 es el historial de producción del pozo 22/11 N-5 del Campo Nelson ubicado en el
UKCS (United Kingdom Continental Shelf). El campo se localiza en los bloques de licitación:
22/11, 22/6a, 22/7 y 22/12a en el zonal central del Mar del Norte. Es un yacimiento con roca
almacenadora formada por areniscas del Paleoceno. El pozo exploratorio fue perforado en 1967
pero fue abandonado dado que no fluyó cuando se le realizó la prueba de presión. Después de
realizarse estudios de sísmicos en 3D se descubrió el campo Nelson y se perforó el pozo 22/11-5
𝑆𝑆𝑆
en 1988. Produce crudo ligero de una °API promedio de 40 con una 𝑆𝑆𝑆 = 555 . Cuenta con
𝑆𝑆𝑆
23 plataformas productoras, 4 cabezales marinos y 4 plataformas inyectoras de agua. Al final de
1999 se ha producido 261 MMBBL. Como se ve en la figura el Campo Nelson es más sencillo
identificar la declinación exponencial a diferencia del Campo Alba que tiene Bombeo
Electrosumergible e inyección de agua, los cuales mantuvieron la presión del yacimiento. En el
Campo Alba era más difícil identificar la declinación debido al SAP y sistema de inyección de
agua, peor en esté pozo el más fácil identificar las etapas de declinación del pozo.

- 75 -
Fig. 47 Historial de producción del pozo 22/11 N-5 del campo Nelson en escala semilog con 7
curvas ajustadas. La producción está separa en 𝑆𝑆 gasto de aceite, 𝑆, 𝑆𝑆 + 𝑆𝑆 gasto normalizado y
𝑆 gasto total saliendo de la cabeza del pozo.

En la fig. 47 se pueden ver la gráfica de log q vs t de 3 gastos: 𝑆, 𝑆𝑆 + 𝑆𝑆 𝑆 𝑆𝑆; donde 𝑆 es el

gasto total [bpd] saliendo de la cabeza del pozo con y se le ajustaron 3 curvas de declinación, la
primera exponencial y la segunda que es compuesta por una hiperbólica y otra armónica al final.
Se puede ajustar un comportamiento hiperbólico que es característico de fluidos compresibles de
yacimientos de gas húmedo o gas y condensado debido que el gasto total incluye el gasto de
gas libre que sale combinado de la cabeza del pozo. Pero los parámetros de 𝑆 𝑆 𝑆0 son
ficticios
aunque la curva se aproxime al comportamiento de los datos reales. Esto debe a que los
coeficientes 𝑆 > 1 son para pozos no convencionales y un factor de 𝑆=5.29 es muy grande. El
𝑆𝑆 + 𝑆𝑆 que es la suma de la producción de aceite y agua se le ajusto una curva de declinación
exponencial. Al 𝑆𝑆 gasto de aceite primero se le ajusto una curva de declinación exponencial,
después una curva compuesta dos curvas de declinación armónica con diferentes coeficientes
𝑆0. Se ajustaron de esta manera los comportamientos dado que la producción no siguió un
comportamiento único en el tiempo. Se desconoce la razón del aumento de la producción a partir
del mes 48 de 11,866 bpd a 13,904 bpd pero elimina el comportamiento exponencial que pudiese
ajustarse a dominio completo de t. La tabla 12 resume las curvas que fueron ajustadas a los
gastos del pozo 22/11 N-5.

Dominio en
Curva Tipo Afecta a Expresión
t

- 76 -
 𝑆𝑆 𝑆𝑆 exponencial 𝑆 𝑆 = 56117𝑆−𝑆∗ 0,0146 8≤𝑆

𝑆𝑆 𝑆𝑆 exponencial 𝑆𝑆 + 𝑆𝑆 + 𝑆𝑆 = 24516𝑆−𝑆∗ 0,01003 8≤𝑆

𝑆𝑆
𝑆𝑆 𝑆𝑆 exponencial 𝑆𝑆 𝑆𝑆 = 24176𝑆−𝑆∗ 0,0174, 8 ≤ 𝑆 8≤𝑆
−1/5,29
5,29𝑆
𝑆𝑆 𝑆𝑆𝑆 Hiperbólico 𝑆 𝑆 = 154860 (( ) + 1) 7 ≤ 𝑆 ≤ 45
5,3
−1
𝑆
𝑆𝑆 𝑆𝑆 Armónico 𝑆 𝑆 = 99225 (( ) +1) 54 ≤ 𝑆 ≤ 71
23
−1
𝑆
𝑆𝑆 𝑆𝑆 Armónico 𝑆𝑆 𝑆𝑆 = 24176 (( ) + 1) 7 ≤ 𝑆 ≤ 48
38,2
−1
𝑆
𝑆𝑆 𝑆𝑆 Armónico 𝑆𝑆 𝑆𝑆 = 13905 (( ) + 1) 49≤ 𝑆 ≤71
23
Tabla 12 Curvas de declinación ajustadas al pozo 22/11 N-5

Fig. 48 Q vs t en escala normal con curvas de declinación ajustadas

- 77 -
Fig. 49 qo+qw y qo vs t en escala normal con curvas de declinación ajustadas

4.4 Analizando la declinación con Curvas Tipo de Fetkovich

Michael J. Fetkovich desarrolló las curvas tipo usando variables adimensionales para hacer
pronósticos de producción. Graficaba la producción del pozo en términos del gasto adimensional
y tiempo dimensional en escala logarítmica y ajustaba una de las curvas tipo al comportamiento
para predecir la producción futura. La ventaja de usarlas es que puede tomarse en cuenta los
periodos del régimen transitorio, la transición al régimen pseudoestacionario y flujo dominado por
fronteras. A diferencia de las curvas de declinación desarrolladas por Arps que sólo pueden
modelar los yacimientos cuando están en régimen pseudoestacionario (declinación) las curvas
tipo pueden abarcan totalmente la vida productiva del pozo. Las curvas tipo combinan
ecuaciones de balance de materia junto con las ecuaciones empíricas de declinación. Una curva
tipo es una composición de funciones, en un dominio del tiempo se modela el régimen transitorio
y en el dominio restante del tiempo el flujo dominado por las fronteras. Se pueden dividir en dos
secciones, la sección izquierda modelante del régimen transitorio usando variables del
yacimiento (𝑆, 𝑆, 𝑆, ℎ , 𝑆, 𝑆𝑆 , 𝑆𝑆 , 𝑆𝑆 , 𝑆𝑆𝑆 ) sin importar si el yacimiento tiene comportamiento
infinito o finito ya que cualquiera de los dos comportamientos converge en la curva tipo con
declinación exponencial (𝑆 = 0) después del 𝑆𝑆𝑆 ≅ 0.3. La sección derecha que modela la
declinación de un
pozo, conjunto de pozos o región usando el factor b para cualquiera de los 3 casos de
declinación: exponencial, hiperbólico y armónico. El dominio de la sección del régimen transitorio

- 78 -
puede prolongarse o encogerse dependiendo del tamaño del yacimiento, dicho de otra forma la

- 79 -
duración del régimen transitorio dependerá de la magnitud del radio de drene o el tiempo que
tarde el radio de drene alcanzar la frontera. Esto se expresa bajo la relación 𝑆𝑆/𝑆𝑆, si la frontera
está más alejada el 𝑆𝑆 tardará más tiempo en crecer para alcanzarla y la duración del
régimen
transitorio será mayor. En la fig. 50 se pude ver en la gráfica cómo sucede lo descrito
anteriormente.

Fig. 50 Efecto del radio de drene en la duración del régimen transitorio de las curvas tipo
(Fetkovich, 1980)

Se usa la ecuación para obtener el gasto adimensional:

𝑆(𝑆) 141.3 𝑆𝑆
𝑆𝑆 = …………(4.17)
𝑆ℎ (𝑆 − 𝑆 )
𝑆 𝑆𝑆
La ecuación para el tiempo adimensional:

0.00634𝑆𝑆
𝑆𝑆 = …………(4.18)
𝑆𝑆𝑆 𝑆 2
𝑆
𝑆

En los anteriores párrafos de la sección 4 de éste capítulo se explica cómo están formadas las
curvas tipo y a continuación se explica cómo usan para generar un pronóstico. Primero se
generar una familia de curvas tipo y se grafican en función del gasto adimensional de declinación
𝑆𝑆𝑆 y tiempo adimensional de declinación 𝑆𝑆𝑆 en escala logarítmica como en la fig. 51. Estos
son diferentes a las variables adimensionales expresadas en el párrafo anterior. Segundo se
calculan el tiempo y gasto adimensionales de declinación ya sea:

𝑆(𝑆
) …………(4.19)
𝑆𝑆𝑆 =
𝑆𝑆
- 80 -
 𝑆𝑆𝑆 = 𝑆𝑆𝑆 …………(4.20)
O e función de las variables del yacimiento:

𝑆𝑆 1
𝑆𝑆𝑆 = 𝑆𝑆 [𝑆𝑆 − ] …………(4.21)
𝑆𝑆 2
0.012649 𝑆 𝑆 1 …………(4.22)
𝑆𝑆𝑆 = ∗
𝑆 𝑆𝑆 𝑆𝑆 1
𝑆 [𝑆𝑆 𝑆𝑆 − 2] [𝑆𝑆2 − 𝑆𝑆2]

Y se grafica 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 en escala logarítmica. Tercero se sobrepone la familia de curvas tipo en
la curva graficada en variables adimensionales de declinación para ver que curva tipo se ajusta
mejor. Se usa la curva tipo para calcular gastos futuros.

Fig. 51 Familia de curvas tipo graficadas en escala logarítmica. (Modificado de Fetkovich. 1980)

Si el comportamiento de la curva tipo no se ajustara totalmente al curva de 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 se

puede usar el siguiente método propuesto por Fetkovich. Lo propuso en su artículo como método
de
curva de desviación (departure curve method) y consiste en que el comportamiento de
declinación está formado por dos o más declinaciones diferentes. Por ejemplo:

𝑆𝑆 = 𝑆1(𝑆, 𝑆1) + 𝑆2(𝑆, 𝑆2) …………(4.23)

Usando la composición de 2 pronósticos diferentes basados en curvas tipo se pude modelar el

comportamiento real de una curva. Éste método lo uso para modelar los yacimientos con
estratos de diferentes permeabilidades sin existir flujo cruzado.

- 80 -
 CAPÍTULO 5

CURVAS DE DECLINACIÓN PARA POZOS DE LUTITTAS FRACTURADAS (TIGHT OIL)

5.1 Introducción

Los nuevos métodos desarrollados para elaborar pronósticos de producción en yacimientos no

convencionales se han concentrado en yacimientos gas. Personas cómo Peter P. Valko, Dilhan
Ilk y Anh N. Doung han desarrollado curvas de declinación basadas en un modelo de
ecuaciones diferenciales diferente al propuesto por Arps:

𝑆𝑆/𝑆
𝑆 …………(5.1)
= −𝑆𝑆𝑆
𝑆

Valko desarrollo su método basado en Declinación Exponencial Extendida de la producción -

(SEPD) Stretched Exponential Production Decline donde define sus parámetros de declina (𝑆, 𝑆)
para calcular el volumen producido en el tiempo usando la función Gamma superior y se enfoca
en yacimientos de gas. Ilk desarrollo su método basado la misma ecuación diferencial usada por
Arps pero lo hace en función de la derivada de los parámetros 𝑆 𝑆 𝑆 para pozos de gas de alta
presión y alta temperatura. Duong es el único que desarrollo un método aplicado tanto a pozos
de gas como aceite que producen de estratos de lutitas fracturadas bajo una geometría de flujo
lineal o bilineal.

La razón por la cual se escogieron estos métodos es porque con ellos se puede estimar la
producción de yacimientos de baja permeabilidad que han sido estimulados por fracturamiento
hidráulico y que estén produciendo en un régimen transitorio aunque algunos de los métodos
sólo puedan aplicarse a pozos de gas como Ilk y Valko.

Método Ecuación diferencial en el que se basan

𝑆𝑆 = − 1 𝑆𝑆𝑆𝑆 , 𝑆 = − 𝑆 [ 𝑆𝑆𝑆 ]
Ilk 𝑆 𝑆𝑆 𝑆𝑆 (𝑆𝑆 /𝑆𝑆 )
𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆
𝑆𝑆 𝑆 𝑆𝑆
Valko = −𝑆 ( )
𝑆𝑆 𝑆 𝑆
𝑆 𝑆 𝑆 𝑆
Doung = 𝑆𝑆 , ( )= 𝑆𝑆
𝑆(𝑆) 𝑆𝑆 𝑆(𝑆) 𝑆𝑆
Tabla 15 Base del modelo de ecuación diferencial para métodos de Ilk, Valko y Duong

Para el alcance de ésta tesis se verán grosso modo los métodos de Valko e Ilk ya que están
enfocados a yacimientos de gas y se verá a fondo el método de Duong. Se hace uso de la
información proporcionada por el Departamento de Mineralogía del Norte de Dakota para
analizar los pozos que producen de la formación del Bakken (B. Superior, B. Medio y B. Inferior)

- 81 -
con el método de Duong y usar su modelo de curva de declinación para pronosticar su
producción futura. El método de Duong es acorde al yacimiento de Bakken porque se desarrolló
para pozos produciendo en flujo lineal o bilineal en régimen transitorio. Los pozos productores
del intervalo del Bakken son pozos horizontales múltiplemente fracturados. Algunos pozos
cuentan con sistemas artificiales de producción como cavidades progresivas, bombeo mecánico,
etc. También se usan pozos inyectores de agua para mantener la presión en algunos campos
productores de la región del Bakken. La mayoría de los pozos inyectores fueron perforados hasta
la formación Three Forks pero fueron tapados con cemento y terminados en la F. Lodgepole.
Algunos de ejemplo de pozos inyectores son: Plu 12-17 (#API 33 089 00832 00 00), Privratsky
12A-3 (#API 33 089 00451 00 00), Haller 29-1 (#API 33 089 00434 00 00).

5.1 Método de Duong

Es un modelo empírico de curva de declinación basado en flujo lineal de larga duración de pozos
en yacimientos no convencionales de baja permeabilidad.

Los pozos con baja permeabilidad y fracturados hidráulicamente producen la mayor parte de su
vida en flujo lineal transitorio a diferencia de los pozos en yacimientos convencionales donde se
observa flujo radial o flujo dominado por fronteras y régimen pseudoestacionario. Duong
𝑆
desarrollo su método encontrando una relación lineal o similar a éste al graficar 𝑆𝑆𝑆 𝑆𝑆 log 𝑆.
𝑆
Esta tendencia lineal se repetía sin importar la conductividad de la fractura como en la Fig. 52.
Una fractura con conductividad finita origina una geometría bilineal y una con conductividad
infinita resulta en geometría lineal. Él determinó su modelo usando la pendiente de la línea, su
intersección con el eje de las x y el gasto inicial para desarrollar su método.

𝑆𝑆
Fig. 52 Gráfica de 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 𝑆 de los pozos Whitman (F. M. Bakken), Boucher (F. Three Forks)
𝑆𝑆
y Charlson (F. Three Forks) . (MDND, 2018)

- 82 -
El tipo de fractura (infinita, finita, múltiple, natural, hidráulica) solo influye en la duración del flujo
lineal transitorio. Para que continúe el flujo a través de la fractura el área de contacto debe
aumentar con el tiempo. Este incremento se logra por el cambio del régimen de esfuerzos debido
al disminución de la presión, lo que genera que las fisuras cerradas durante la diagénesis
puedan activarse (Warpinski y Branagan, 1989).

En la Fig. 52 puede verse el comportamiento lineal de 3 pozos produciendo de la F. Media del

Bakken y la F. Three Forks. Todos son pozos horizontales y cuentan con varias etapas de
fracturamiento hidráulico. El pozo Charlson Federal 14X-35H tiene el comportamiento más lineal
de los 3, éste pozo sólo fue estimulado en el 2010. El pozo Whitman 2-34H tiene una
discontinuidad debido a se instaló un sistema artificial de cavidades progresivas. El Pozo
Boucher 41X-21 fue terminado y produce sin SAP. En los anexos se puede ver más información
sobre los pozos. (MDND, 2018). Mineral Department of North Dakota

Duong se basó en que si el régimen de flujo característico de un sistema fracturado (lineal,

bilineal) se prolonga toda la vida de un pozo de gas, éste comportamiento pude representarse
por:

𝑆 = 𝑆1𝑆−𝑆 …………(5.2)

Donde n es 1/2 para flujo lineal y 1/4 para flujo bilineal y 𝑆1 es el gasto del primer día. Entonces
la producción acumulada es:
𝑆 𝑆1−𝑆
𝑆𝑆 = ∫ 𝑆 𝑆𝑆 = 𝑆1 …………(5.3)
0 (1 − 𝑆)

Cómo se puede ver en la Fig. 51. La tendencia que seguiría una línea recta si se ajustará seria:

𝑆 (1 − 𝑆) 𝑆
= 𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆 = 𝑆𝑆 −𝑆 …………(5.4)
𝑆𝑆 𝑆𝑆 𝑆𝑆

Usando la ecu. 5.4 basó su modelo de curva de declinación. Para ver la deducción de la curva de
declinación consulte el Anexo F. Dando como resultado.

𝑆 𝑆 1−𝑆
…………(5.5)
= 𝑆 −𝑆 ∗ 𝑆 [ 1−𝑆(𝑆 −1 )]
𝑆1
𝑆 𝑆 1−𝑆
𝑆𝑆 = 1 𝑆[ 1−𝑆 (𝑆 −1 )] …………(5.6)
𝑆

Para el gasto en el tiempo ec. 5.5 y la producción acumulada ec. 5.6. También se desarrollaron
curvas tipo usando el 𝑆𝑆 y la relación 𝑆/𝑆𝑆𝑆𝑆 dónde 𝑆𝑆𝑆𝑆 normalmente es obtenido en el
primer mes de producción.

- 83 -
 𝑆
𝑆𝑆 = …………(5.7)
𝑆𝑆𝑆𝑆
𝑆 𝑆 1−𝑆
…………(5.8)
= 𝑆 𝑆−𝑆 𝑆 [ 1−𝑆(𝑆𝑆 −1 )]
𝑆𝑆𝑆𝑆

O pueden ser obtenidos de la siguiente forma

1
𝑆 1−𝑆
𝑆𝑆𝑆𝑆 = ( ) …………(5.9)
𝑆
𝑆
𝑆 1−𝑆 𝑆−𝑆 …………(5.10)
𝑆𝑆𝑆𝑆 = 𝑆1 ( ) 𝑆1−𝑆
𝑆

Las curvas tipo ayudan a identificar cambios de regímenes de flujo y comparar producciones
entre pozos. Para calcular la Recuperación Final Estimada – Ultimate Estimated Recovery EUR
se usa la ec. 5.11. En donde 𝑆𝑆𝑆𝑆 es el gasto en el límite económico y el 𝑆𝑆𝑆𝑆 se obtiene del
rango
del tiempo hasta alcanzarse el 𝑆𝑆𝑆𝑆.
𝑆𝑆𝑆𝑆
𝑆𝑆𝑆 = 𝑆𝑆𝑆𝑆𝑆
𝑆 …………(5.11)

(Duong, 2011)

5.1.1 Como aplicar el método Duong

Para ejemplificar el método de Duong se usará el Pozo Charlson Federal 14X-35H ya que es el
que cumple con los siguientes criterios:

 Produce en flujo bilineal durante la mayor parte de su vida productiva fig. 53 a)

𝑆𝑆
 Tiene una relación = 0.02883, en otras palabras corte de agua muy bajo
𝑆𝑆
 Se encuentra en la plataforma transitoria de RGA fig. 53 a)
 Presenta una tendencia lineal en la gráfica de 𝑆𝑆 𝑆𝑆 𝑆 1/2 que es característico de pozos
fracturados fig. 53 b)
 Producen en la Formación Three Forks que es la parte no convencional junto con la F.
Bakken y F. Lodgepole de la Región del Bakken del Norte de Dakota.
 Es un pozo horizontal con una profundidad horizontal de 3002.5 [m] y una profundidad
desarrollada 4617.72 [m], su KOP es a 2761.4 [m desarrollados]. La tubería productora es
un liner 4 ½´´ asentado de 2760.5 [m] a la profundidad total desarrollada. El pozo produce
crudo de 43.7 °API

- 84 -
Fig. 53 a) 𝑆𝑆𝑆{𝑆𝑆, 𝑆𝑆, 𝑆𝑆, 𝑆𝑆(𝑆𝑆 + 𝑆𝑆)} 𝑆𝑆 𝑆𝑆𝑆 𝑆 con una pendiente de ¼ que es característico de
flujo lineal de pozo Charlson Federal 14X-35H. b) 𝑆𝑆𝑆𝑆 𝑆𝑆/𝑆 usada para ver que pozos producen en
flujo lineal y en el periodos de tiempo. 18861 (Carson Peak 2-35H), 16059 (USA 2D-3-1H), 17147
(Boucher 41X-21), 20342(USA 153-95-4B-9-1H), 20210 (Whitman 2-3H), 17113 (Charlson Federal 14X-
35H). (MDND, 2018)

Primer paso: revisar los datos de producción y corregirlos si es necesario, como del ejemplo del
capítulo 4, donde se usa el promedio de 3 y 5 puntos para suavizar la curva. Cuando el corte de

- 85 -
agua ha alcanzado una fracción considerable del gasto total se debe decidir si usar el 𝑆𝑆 𝑆 𝑆𝑆 +
𝑆𝑆 dado que el agua está recibiendo una fracción de trabajo que no puede ser descartada. En el
caso de pozos de gas y condesado el gasto debe ser corregido a gasto equivalente de gas si
existe un razón condensado/gas alta. En la fig. 54 se tiene como ejemplo al Pozo Charlson
Federal y tiene un comportamiento característico de un HMFW. Los datos no requieren de
corrección dado que las fluctuaciones del gasto son cíclicas y mínimas.

Fig. 54 Historial de producción del pozo Charlson Federal 14X-35H para usarse como ejemplo.
(MDND, 2018)
𝑆𝑆
Segundo paso: Determinar los coeficientes a y m. Haciendo uso de la gráfica log 𝑆𝑆 log 𝑆 y la
𝑆𝑆
𝑆
relación ec. 5.4 ( = 𝑆𝑆 −𝑆 ) al ajustar la recta se pueden encontrar los coeficientes. Un valor de
𝑆𝑆
𝑆2 > 0.95 de ajuste de la recta es recomendado. En la fig. 55 se puede ver la ec. de la recta
ajustada a 𝑆2 = 0.9541, los coeficientes son: 𝑆 = 0.087 𝑆 𝑆 = 1.354.

- 86 -
 𝑆𝑆
Fig. 55 Gráfica 𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 𝑆 del pozo Charlson Federal 14X-35H para encontrar a y m.
𝑆𝑆

Tercer paso: Obtener los parámetros 𝑆1 𝑆 𝑆∞. Para obtenerlos se grafica 𝑆𝑆 𝑆𝑆 𝑆(𝑆, 𝑆), donde
𝑆(𝑆, 𝑆) es la ec. 5.8. Una vez graficados los puntos se ajusta una recta donde la pendiente m es
𝑆1 y la intersección en el eje de las x es 𝑆∞. Para representar el modelo de la siguiente manera:
𝑆 = 𝑆1𝑆(𝑆, 𝑆) + 𝑆∞ . Con los parámetros obtenidos se pueden generar los pronósticos como en
las curvas de declinación de Arps, la única diferencia es que 𝑆𝑆 = 𝑆1.

- 87 -
Fig. 56 Gráfica 𝑆𝑆 𝑆𝑆 𝑆 (𝑆, 𝑆) para encontrar los parámetros 𝑆𝑆 𝑆 𝑆∞. En este caso se ajustaron 4
líneas de tendencia para ver que pronóstico se ajusta mejor al perfil real de producción.

Reservas (mes
Pronóstico Coeficientes % Error
108) [BBL]
Primero 𝑆1 = 30,228 𝑆 𝑆∞ = 63.27 835, 949.5 0.343
Segundo 𝑆1 = 19441 𝑆 𝑆∞ = 88 845,732.1 1.5177
Tercero 𝑆1 = 8297 𝑆 𝑆∞ = 113.5 855,512.9 2.6917
Cuarto 𝑆1 = 6932 𝑆 𝑆∞ = 128.1 931,964.3 11.868
Real 833,088
Tabla 16 Parámetros 𝑆𝑆 𝑆 𝑆∞ para cada línea de tendencia de la fig. 56. Porcentaje de error
para cada estimación.

Cuarto paso: Estimar reservas. Se usa la expresión

𝑆
𝑆𝑆 = 𝑆 𝑆 …………(5.12)
𝑆

- 88 -
Donde 𝑆 representa el gasto diario o mensual y no 𝑆1. Aquí se puede usar el 𝑆𝑆𝑆𝑆 gasto mínimo
o gasto del límite económico para estimar las reservas. (Duong, 2011).

Fig. 57 Cuatro pronósticos de producción comparados contra el perfil de producción original. Se

usa escala semilogarítmica para apreciar más fácil el ajuste.

- 89 -
Fig. 58 Estimación de reservas usando los cuatro pronósticos de producción determinados con el
método de Duong.

5.1.2 Análisis y discusión de los resultados del ejemplo usado para el método de Duong

Cómo puede apreciarse en la fig. 56. los pronósticos de producción tienen un buen ajuste en la
vida tardía del pozo pero no pueden modelar el comportamiento inicial. Y al hacer la estimación
de reservas los pronósticos 1, 2 y 3 tienen errores mínimos pero no se ajustan al comportamiento
del pozo. No muestran la concavidad característica de los datos originales. Si fueran usados
estás curvas declinación sería un escenario optimista para el pozo, por lo que no son
recomendable de usarse.

El error inicial en el ejemplo, se cómete al escoger los coeficientes a y m. Como puede verse en
la fig. 55 la línea de tendencia omite el primer punto de la gráfica lo y los valores de a y m son
muy bajo en a y m es muy grande. En la fig. 59 se puede ver el análisis del pozo Charlson
Federal 14X-35H usando 𝑆 = 𝑆. 𝑆𝑆 𝑆 𝑆 = 𝑆. 𝑆𝑆, la línea de tendencia no tiene ningún ajuste
con los datos reales pero los pronósticos de producción tienen un ajuste demasiado mejor al
hecho inicialmente al demostrar cómo usar el método de Duong.

Reservas (mes
Pronóstico Coeficientes % Error
108) [BBL]
Primero 𝑆1 = 1993.9 𝑆 𝑆∞ = −25.46 1,057,335.8 26.91
Segundo 𝑆1 = 1762.3 𝑆 𝑆∞ = −14.13 961,161.8 15.37
Tercero 𝑆1 = 1543.8 𝑆 𝑆∞ = 4.01 894,074.8 7.32

- 90 -
 Cuarto 𝑆1 = 831.4 𝑆 𝑆∞ = 110 827,524.4 -0.667
Real 833,088
Tabla 17 Parámetros 𝑆𝑆 𝑆 𝑆∞ para cada línea de tendencia de la fig. 59. Porcentaje de error para
cada estimación

Los nuevos pronósticos de producción en la fig. 59 a), b), c) d) son una aproximación muy
realista del perfil de producción aunque no se ajusta la línea de tendencia en la regresión lineal.
Lo que muestra que el método de Duong es una herramienta muy buena pero deben de
escogerse correctamente los coeficientes a y m, estos son la base de un certero pronóstico de
producción.

- 91 -
- 92 -
Fig. 59 Pronóstico de producción del ejemplo para el método de Duong usando 𝑆=
𝑆. 𝑆𝑆 𝑆 𝑆 = 𝑆. 𝑆𝑆. Los pronósticos se ajustan muy bien comparación de los propuestos
inicialmente.

5.2 Método de Ilk para pozos de Gas HPHT (alta presión alta temperatura)

Dilhan Ilk desarolló un método de curvas de declinación para pozos fracturados de alta presion y
alta temperatura (HP/HT) de gas. Usa el perfil de producción y la producción a cumulada para
encontrar las funciones de los parámetros 𝑆 𝑆 𝑆 para formular la relación semianalítica de
Ansah y estimar las reservas de gas. Para representar cualquier régimen de flujo usan el
modelo
- 93 -
empírico de declinación exponencial de potencias- power-law exponential rate decline - en
conjunto en la relación de Ansah. El uso de las funciones de los parámetros 𝑆 𝑆 𝑆 da una visión
única entre muchos métodos para identificar los regímenes de flujo. La desventaja del método es
que la calidad de los datos afecta mucho la certeza del pronóstico. El modelo empírico de
declinación exponencial de potencias se usa en conjunto con la relación de Ansah, ya que la
relación semianalítica de Ansah sólo puede aplicarse para flujo dominado por fronteras. Se
deben usar las dos relaciones simultáneamente para disminuir la incertidumbre/no unicidad
asociada con la estimación de reservas de pozos HP/HT.

La relación de Ansah es la base del método pero sólo se aplica para BDF (boundary dominated
flow) por lo que es usado el modelo empírico de declinación exponencial de potencias que puede
usarse en cualquier régimen de flujo.

𝑆𝑆𝑆 = 4 𝑆𝑆𝑆2𝑆[−𝑆𝑆𝑆 𝑆𝑆𝑆]

((1 + 𝑆 ) − (1 − 𝑆 )𝑆[−𝑆𝑆𝑆 𝑆𝑆𝑆])2 …………(5.13)
𝑆𝑆 𝑆𝑆

Ec. 5.13 es la relación de Ansah expresada en variables adimensionales. Donde 𝑆𝑆𝑆 es presión
adimensional de declinación en el pozo, 𝑆𝑆𝑆 tiempo adimensional de declinación y 𝑆𝑆𝑆 es el
gasto adimensional de declinación.
𝑆𝑆
𝑆𝑆𝑆 = …………(5.14)
𝑆𝑆𝑆

𝑆𝑆𝑆/𝑆𝑆
𝑆
…………(5.15)
𝑆𝑆𝑆 =
𝑆𝑆
𝑆𝑆

𝑆𝑆
𝑆𝑆𝑆 = 𝑆 …………(5.16)
𝑆𝑆𝑆𝑆

Donde 𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆 [𝑆𝑆𝑆/𝑆], 𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 [𝑆𝑆𝑆],
𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆, 𝑆𝑆
𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆 [𝑆𝑆𝑆],
𝑆 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆
𝑆𝑆𝑆𝑆𝑆𝑆𝑆, 𝑆𝑆𝑆𝑆∗ 𝑆𝑆𝑆
],
𝑆 𝑆𝑆 í𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆
𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 [
𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆 [𝑆𝑆𝑆 −1 ], 𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 [𝑆𝑆𝑆]
2𝑆ℎ
𝑆𝑆 =
2.2458𝑆 …………(5.17)
141.2 𝑆𝑆𝑆𝑆𝑆 ln [ 𝑆 𝑆 2 ]
𝑆 𝑆𝑆

Usando el método de radios de pérdida y la derivada de radios de pérdida se definen las

- 94 -
funciones adimensionales 𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 − 𝑆 (𝑆𝑆) y 𝑆𝑆𝑆𝑆𝑆ó𝑆
𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 − 𝑆 (𝑆):

- 95 -
 1 𝑆𝑆𝑆𝑆
𝑆𝑆 = − …………(5.18)
𝑆 𝑆𝑆 𝑆𝑆
𝑆𝑆

𝑆=− 𝑆 [ 𝑆𝑆𝑆 ]
𝑆𝑆 (𝑆 /𝑆𝑆 ) …………(5.19)
𝑆𝑆 𝑆𝑆 𝑆𝑆

Esta sección es para definir las ecuaciones que conforman a la relación de Ansah que sólo
puede aplicarse en Flujo dominado por frontera.

Con las variables adimensionales definidas se elabora la gráfica ¨ 𝑆 − 𝑆𝑆 − 𝑆 ¨ para obtener los
parámetros de la ec. 5.2: 𝑆∞𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆 → ∞ [𝑆(𝑆
= ∞)] y
𝑆̂𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆.

La forma de realizar los pronósticos de producción es usando la relación de declinación

exponencial de potencias ec. 5.20

̂𝑆 𝑆 𝑆 ]
𝑆𝑆 = 𝑆̂𝑆𝑆 exp[−𝑆∞ − 𝑆 …………(5.2)

5.3 Método de Valko para pozos de Gas SEPD (Streched Exponential Production Decline)

Los niveles de estudio en ciencias de la tierra van desde lo microscópico hasta lo macroscópico
o en este caso hasta lo nanoscópico.

5.4 Caso Práctico con pozo de Bakken

- 96 -
 CONCLUSIONES
Los.

- 97 -
 BIBLIOGRAFÍA

1 Arévalo, J. (2005). Ingeniería de Yacimientos de Gas (Propiedades de los fluidos

petroleros). UNAM, México. pg.21

2 Castillo, C. (2014). Diseño, validación e implementación de un sistema experimental de

alta presión y temperatura para la medición de la viscosidad de aceites crudos con precisión de
referencia (Tesis Licenciatura). UNAM, México. pg. 8

3 Lund, L. (2014). Decline Curve Analysis of Shale Oil Production. Universidad de Uppsala.
Suecia, Uppsala. Pp 56

4 Baca U.,G. (2001). Evaluación de Proyectos. 4ta Ed. Mc Graw Hill. México, D.F. pg. 211 -
228

5 Hydraulic Fracturing (2012). Canadian Society for Unconventional Resources.

Recuperado de http://www.csur.com/resources

6 Hamison B., Fairhust C. (1985). Rock Mechanics – Theory and Practice. E.U.A.,
Minnesota. Pg. 559-567.

7 Ballines, J. (1992). Optimum Design of a Fracturing Job in the Chicontepec Formation.

SPE 23622. Caracas, Venezuela. Pp 67

8 Castellanos, F. (2015). Caracterización dinámica de Yacimientos no convencionales

(Tesis Maestría). UNAM, México.

9 Leimkuhler J., Leveille G. (2012). Unconventional Resources. SPE 0112 026 TWA.
Texas, E.U.A. pp 3

10 Tissot, B., Welte, P. (1984). Petroleum Formation and Occurrence. Alemania: Springer-
Verlag Berlin Heidelberg.

11 Buller A. , Berg E., Et al. (1990). Oil and Gas Reservoirs II. Ed. Graham and Trotman.
Trondheim, Norway. pg. 25-45

12 Magoon L., Dow W. (1994). The petroleum system from source rock to trap.
E.U.A.,Oklahoma. Pp.12

13 The petroleum system an overview. (2017) Recuperado de

http://www.mhnederlof.nl/petroleumsystem.html

- 98 -
14 Elsaadany N. (2013). Hydrocarbon Migration. Crazy planet. Recuperado de
http://furiousplanet.blogspot.mx/2013/11/hydrocarbon-migration.html

15 Escobar H. (2011) Apuntes de Geología de Explotación. Recuperado de

http://usuarios.geofisica.unam.mx/gvazquez/explotacionELIA/zonadesplegar/Clases/Clase%2010
%20Migracion.pdf

16 Meissner F. (1984). Petroleum Geochemistry and Basin Evaluation: Petroleum Geology of

the Bakken Formation Williston Basin, North Dakota and Montana. APPG Colorado, E.U.A. pg.
207-227

17 Jin H., Sonnenberg S. (2014). Characterization for Source-Rock Potencial of the Bakken
Shales in the Williston Basin, North Dakota and Montana. APPG. Artículo #80356. Colorado,
E.U.A. Pp 37

18 Webster R. (1984). Petroleum Source Rocks and Stratigraphy of the Bakken Formation in North
Dakota. AAPG. Dakota del Norte, E. U. A. pg. 45

19 Gallardo, E. (2015). Caracterización dinámica de Yacimientos. UNAM. México, Distrito

Federal. Pp 156

20 Arana V. , Trujillo D. (2015). Apuntes de Simulación Numérica de Yacimientos. UNAM.

México, Distrito Federal. Pp 307

21 Jones W., McCord D., Cummer R. (1969). Reservoir Simulation Pays a Big Dividend. SPE
2428. E.U.A., Texas. Pp 18

22 Bobadilla G. (2012). Conceptos Básicos de Caracterización Dinámica de Yacimientos.

Tesis Licenciatura. UNAM. México, Distrito Federal. Pp 157

23 Fekete Association. (2012). Theory and Equations. Recuperado de

http://www.fekete.com/SAN/TheoryAndEquations/WellTestTheoryEquations/Flow_Regime_Categ
ories.htm

24 Herrera D., Calderon A. (2012).Las Pruebas de Inyectibilidad en la Caracterización de

Yacimientos. UNAM. Tesis Licenciatura. México, D.F. pg. 19-55

25 Nieto R., R. (1980). Principios de Mecánica de Yacimientos. UNAM. Apuntes. México,

D.F. Pp 99

26 Jones S. (2016). Producing-Gas/Oil-Ratio Behavoir of Multifractured Horizontal Wells in

Tight Oil Reservoirs. SPE 184397. E.U.A., Texas. Pp 13

27 Tabatabaie H., Pooladi M. (2016). Multiphase Linear Flow in Tight Oil Reservoirs. SPE
180932. Canada, Calgary. Pp 13

- 99 -
28 Vermilye J., Sicking C. (2016). Microseismic Maps Production Volume. The american oil
and gas reporter. Texas, Houston. Pp 5

29 Evaluation of Hydrocarbon Generation and Storage in the Barnett Shale, Ft. Worth Basin,
Texas©. (2004). Teamfrack. Recuperado de
http://teamfrack.pbworks.com/w/file/fetch/48478518/Evaluation%20of%20Hydrocarbon%20Genra
tion%20and%20Storage%20in%20the%20Barnett%20Shale.pdf

30 Source Rock (SR) Analyzer: User Guide. (2014). International Ocean Discovery Program.
Recuperado de https://rosetta.iodp.tamu.edu/A/TechDoc/14?encoding=UTF-8

31 North Dakota Mineral Department, Oil and Gas Division (2017). Recuperado de
https://www.dmr.nd.gov/oilgas/

32 Arps, J. (1944). Analysis of decline curves. SPE 945219. E.U.A.., Texas. pg. 228-247

33 Fetkovich, M. (1980). Decline curve analysis using type curves. SPE 4629. E.U.A., Texas.
pg. 1065-1077

34 Buryakovsky L., Chilingar G., et.al. (2012). Petrophysics: Fundaments of the oil and gas
reservoirs. Wiley & Scrivener Publishing LLC. Texas, E.U.A. pg. 184-185

35 Garaicochea F., Bernal C., López O. (1991). Transporte de Hidrocarburos por Ductos.
Colegio de Ingenieros petroleros de México. México, D.F. pg. 17-19

36 Historial de Producción. (2007). Petroleum Production Reporting System Dashboard.

Recuperado de
https://itportal.ogauthority.co.uk/information/wells/pprs/Well_production_offshore_oil_fields/offshor
e_oil_fields_by_well/offshore_oil_fields_by_well.htm

37 Jayasekera A., Goodyear S. , et.al. (1999). The Development of Heavy Oil Fields in the
U.K. Continental Shelft: Past, Present and Future. SPE 54623. Alaska, Anchorage. Pp 17

38 U.K. Continental Shelft: Past, Present and Future. SPE 54623. Alaska, Anchorage. Pp 17

39 Chapter 8: Production decline analysis. (2010). New Mexico Tech. Recuperado de

http://infohost.nmt.edu/~petro/faculty/Kelly/Deline.pdf

Fetkovich, M. A.; Vienot, M. E.; Bradley, M. D.; et. Al. (1984). Decline Curve Analysis Using Type
Curves-Case Histories. JPT SPE 4629 PA. U.S.A., Texas. Pp 20

- 100
-
40 Ilk D., Rushing J. A. (2009). Decline Curve Analysis for HP/HT Gas Wells: Theory and
Applications. SPE 125031. E.U.A.. Louisiana. Pp 24

41 Valko P., Lee W. J. (2010). A Better way to Forecast Production from Unconventional Gas
Wells. SPE 134231. Italia, Florencia. Pp 16

42 Duong A. N. (2011). Rate-Decline Analysis for Fracture-Dominated Shale Reservoirs. SPE

137748. Cánada, Calgary. pg 377-387

43 Well search. (2018). Mineral Department of North Dakota. Recuperado de

https://www.dmr.nd.gov/oilgas/

44 GIS Map Server. (2018). Mineral Department of North Dakota. Recuperado de

https://www.dmr.nd.gov/OaGIMS/viewer.htm

45 Well Files. (2018). Mineral Department of North Dakota. Recuperado de

https://www.dmr.nd.gov/oilgas/basic/bwfiles

- 100 -
 ANEXOS
Anexo A – Algebra Matricial, Método de Newton-Raphson y Matriz Jacobiano

El fin de obtener el Jacobiano de la matriz es para evaluar en un punto, el cual son los valores de
la matriz 𝑆𝑆0 que puede ser una matriz unitaria o de ceros. Auxiliarse del método de Newton-
Raphson para encontrar las soluciones del sistema de ecuaciones. Esto se tendrá que hacer
iterativamente hasta que converjan la solución 𝑆𝑆𝑆 y así pasar al siguiente paso del tiempo y
seguir el algoritmo anterior. Esta serie de pasos se llevará acabo hasta que se llegue al último
paso del tiempo y después se reportarán los resultados.

Se hablará un poco más del método de Newton-Raphson y la matriz del Jacobiano con el
propósito de ilustrar al lector sobre el método y no dejar el tema al vapor. Se empezará por las
bases de algebra matricial, después de hablará de cálculo vectorial y para finalizar se tocará el
tema del método numérico.

Algebra Matricial

Un pequeño repaso de concepto. Primero multiplicación entre matrices.

1 2 3 2 3
𝑆 = [4 5 6] 𝑆 = [5 1]
7 8 9 0 3

3 x 3 - 3 x 2 = 3 x 2

2 3
[1 2 3] [5] [1 2 3] [1]
0 3
2 3 1∗ 2+2∗ 5+3∗ 0 1∗ 3+2∗ 1+3∗ 3
𝑆𝑆𝑆 [4 5 6 ] [5] [4 5 6 ] [1] =[4 ∗ 2+5∗ 5+6∗ 0 4 ∗ 3 + 5 ∗ 1 + 6 ∗ 3]
=
0 3 7∗ 2+8∗ 5+9∗ 0 7∗ 3+8∗ 1+9∗ 3
2 3
[7 8 9] [5] [7 8 9 ] [1]
[ 0 3]
2 + 10 + 0 3+2+9 12 14
𝑆𝑆𝑆 = [8 + 25 + 0 12 + 5 + 18] = [33 35]
4 + 40 + 0 21 + 8 + 27 54 56

Matriz Adjunta

También conocida como la matriz de cofactores .

𝑆∗ = (𝑆)∗

- 101 -
 1 2 3
1 2 3 +1 −2 +3 2 3 5 6 2 3
∗ 4 5 6 8 9
𝑆 = [4 5 6] = [ −4 +5 −6] = 7 8 9 =[ 4 5 6] =
𝑆
7 8 9 +7 −8 +9 4 5 6
7 8 9
[ 7 8 9]

5∗ 9−6∗ 8 2 3 −3 2 3 −3 −(4 ∗ 9 − 6 ∗ 7) (4 ∗ 8 − 5 ∗ 7)
[ 4 5 6]=[ 4 5 6]=[−(2 ∗ 9 − 3 ∗ 8) (1 ∗ 9 − 3 ∗ 7) −(1 ∗ 8 − 2 ∗ 7)]
7 8 9 7 8 9 (2 ∗ 6 − 3 ∗ 5) −(1 ∗ 6 − 3 ∗ 4) (1 ∗ 5 − 2 ∗ 4)

−3 6 −3
𝑆 =[6
∗
−12 6]
−3 6 −3

Matriz Traspuesta

𝑆𝑆 = (𝑆)𝑆

Anexo B – Calculando el historial de producción del Pozo 16/26 A9 del Campo Alba

Pozo 16/26-A9 del Campo Alba. Ubicado en la región de Mar del Norte en agua inglesas.

Datos de Producción:

Presión en
Producción Producción
Producción la cabeza Días
Año Mes de Gas de Agua
Aceite (m³) del Pozo producidos
(km³) (m³) (bar)
1994 10 12,539 595 216 20.2 5.77
1994 11 63,116 3,028 1,159 13.5 25.51
1994 12 77,027 3,300 319 12.3 28.96
1995 1 60,527 2,696 228 12.3 22.87
1995 2 77,874 1,776 322 12.1 27.5
1995 3 86,563 3,582 1,002 12.1 29.56
1995 4 75,617 3,112 951 12.6 25.98
1995 5 51,364 2,389 4,216 12.8 17.73
1995 6 54,965 2,574 1,223 9.7 28.02
1995 7 55,942 2,454 2,180 9.9 30.38

- 102 -
1995 8 43,448 2,034 2,465 0 23
1995 9 20,119 929 1,169 0 15.46
1995 10 56,374 2,528 3,141 9.8 30.14
1995 11 52,278 2,470 5,929 9.8 27.89
1995 12 54,341 2,356 6,518 9.7 30.91
1996 1 47,210 4,831 4,710 9.8 29.56
1996 2 49,140 2,038 3,595 0 27.64
1996 3 43,367 1,878 3,899 18.8 27.81
1996 4 40,325 2,175 4,326 0 21.67
1996 5 60,776 3,021 4,521 0 25.32
1996 6 56,169 2,529 5,652 12.9 26.06
1996 7 38,157 1,832 6,450 0 28.34
1996 8 24,216 1,042 3,186 0 13.65
1996 9 73,159 3,338 3,394 12.7 22.64
1996 10 62,581 2,824 10,240 0 28.33
1996 11 51,773 2,069 8,130 0 27.94
1996 12 47,109 2,007 7,561 0 27.71
1997 1 52,327 2,680 7,223 0 30.47
1997 2 51,727 2,479 7,558 12.6 25.65
1997 3 55,606 2,764 7,492 10.8 30.03
1997 4 55,629 2,783 8,145 11.5 29.03
1997 5 59,433 3,069 8,505 11.4 30.19
1997 6 57,204 2,746 7,316 0 29.49
1997 7 54,935 2,635 5,523 0 30.43
1997 8 55,771 2,672 6,570 0 30.49
1997 9 50,312 2,347 6,490 0 25.57
1997 10 56,465 2,704 9,189 0 30.95
1997 11 46,685 2,196 8,649 0 26.98
1997 12 55,370 2,477 13,002 0 30.24
1998 1 53,717 2,461 22,095 0 30.44
1998 2 47,084 2,153 21,171 0 27.9
1998 3 41,036 1,786 19,244 0 25.94
1998 4 40,799 1,924 24,023 0 28.48
1998 5 39,812 1,907 24,520 0 26.75
1998 6 27,131 1,323 16,014 0 19.56
1998 7 39,318 1,833 20,524 0 22.76
1998 8 37,814 1,847 25,442 0 27.69
1998 9 38,599 1,944 27,202 0 27.48
1998 10 31,435 1,584 24,300 0 24.95
1998 11 30,586 1,405 21,878 10.2 28.91

- 103 -
 1998 12 32,201 1,587 23,819 0 29.54
1999 1 28,447 1,315 18,834 0 25.58
1999 2 28,589 1,312 21,310 0 26.28
1999 3 31,098 1,412 23,933 9.7 30.6
1999 4 15,606 660 16,392 9.5 17.6
1999 5 8,626 286 3,227 16.7 22.61
1999 6 32,598 1,368 15,437 21.1 27.94
1999 7 33,119 1,556 31,909 18.4 30.56
1999 8 26,253 1,383 28,503 19.9 30.92
1999 9 18,307 798 24,557 17 29.12
1999 10 13,729 657 25,011 21.4 30.63
1999 11 4,792 282 12,099 17.4 17.43
1999 12 0 0 0 17.4 0
(itportal.ogauthority.co.uk, 2007)

Grados °API promedio del Campo: 20 ° API

Se asumen lo siguiente:

Temperatura de salida del pozo: 110°F

Densidad relativa del gas libre 𝑆𝑆: 0.77

Variación del factor de compresibilidad 𝑆𝑆 =

0.97 Factor de compresibilidad final 𝑆𝑆 = 0.85

Para determinar la presión de la cabeza del pozo de usa una distribución de Pareto entre la
𝑆𝑆𝑆𝑆 = 304.57 𝑆𝑆𝑆 y 𝑆𝑆𝑆𝑆 = 140.68 𝑆𝑆𝑆

𝑆̂𝑆ℎ = 173.46 𝑆𝑆𝑆

Se realiza la suposición que a las condiciones de 𝑆𝑆ℎ = 110°𝑆 y 𝑆̂𝑆ℎ = 173.46 𝑆𝑆𝑆 se tiene un
flujo conformado de una fase líquida y gaseosa. La fase liquida formada por agua y aceite que
tienen disuelto gas. La fase gaseosa conformada por el gas producido.

𝑆𝑆@𝑆. 𝑆. = 𝑆𝑆 + 𝑆𝑆 + 𝑆𝑆 …… (1)

𝑆𝑆@𝑆. 𝑆. = 𝑆𝑆𝑆𝑆 + 𝑆𝑆𝑆𝑆 + 𝑆𝑆𝑆𝑆 ……(2)

𝑆𝑆@𝑆. 𝑆. = (𝑆𝑆 − 𝑆𝑆𝑆𝑆𝑆 − 𝑆𝑆𝑆𝑆)𝑆𝑆 ……(3)

Sustituyendo (2) en (3):

𝑆 @𝑆. 𝑆. = 𝑆 +𝑆𝑆 +
𝑆 𝑆𝑆 ( −𝑆𝑆 − 𝑆 𝑆 )𝑆 ……(4)
𝑆 𝑆
𝑆 𝑆 𝑆 𝑆𝑆 6.2904 𝑆𝑆3 𝑆 𝑆 𝑆𝑆

- 104 -
𝑆 𝑆 𝑆

- 105 -
@𝑆. 𝑆. ∶ 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆𝑆

Teniendo el siguiente modelo y los datos de producción se pueden combinar para calcular el
gasto que sale de la cabeza del pozo.

Ahora se calculan los siguientes parámetros usando correlaciones:

𝑆𝑆 = 1.0 + 1.2𝑆10−4(𝑆 − 60) + 1𝑆10−6(𝑆 − 60) − 3.33𝑆10−6 ∗ 𝑆 ; 𝑆 [°𝑆], 𝑆 [𝑆𝑆𝑆]

°𝑆𝑆𝑆
𝑆 = 1.0 + 𝑆 𝑆 + (𝑆 − 60) ( ) (𝑆 +𝑆𝑆 ) 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆 𝑆𝑆𝑆á𝑆𝑆𝑆𝑆𝑆
𝑆 1𝑆 𝑆𝑆 2 3 𝑆
𝑆

Coeficiente °𝑆𝑆𝑆 ≤ 30° °𝑆𝑆𝑆 > 30°

𝑆1 4.677𝑆10−4 4.67𝑆10−4
𝑆2 1.751𝑆10−5 1.1𝑆10−5
𝑆3 −1.811𝑆10−8 1.337𝑆10−9

𝑆𝑆
𝑆 =𝑆 (1 + 5.912𝑆10−5 °𝑆𝑆𝑆 𝑆 log( )) ; : 𝑆 y 𝑆 ,
𝑆 𝑆 𝑆
𝑆𝑆 𝑆 𝑆 114.7 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆
𝑆 𝑆𝑆
𝑆
𝑆𝑆 = 110°𝑆; 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆, 𝑆𝑆: 𝑆𝑆𝑆𝑆𝑆ó𝑆 1𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆ó𝑆
𝑆 (𝑆 + 460)
𝑆𝑆 = 0.02825 ; 𝑆 [°𝑆], 𝑆[𝑆𝑆𝑆]
𝑆
𝑆𝑆𝑆 = 5.6146 𝑆
2 3
𝑆 = 𝑆´ (𝑆 + 𝑆𝑆´ + 𝑆𝑆 ´ + 𝑆𝑆´ )

𝑆 = 3.69051, 𝑆 = 0.08746, 𝑆 = 0.01129, 𝑆 = −0.00647

−𝑆
𝑆´ = 1 − exp ( ) ;𝑆[𝑆𝑆𝑆]
2276
𝑆∗ − 32
𝑆´ =
10
5
𝑆∗ = (𝑆 − 32) ; 𝑆[°𝑆]
9

𝑆𝑆 … 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆𝑆 𝑆𝑆 𝑆𝑆𝑆𝑆𝑆 𝑆𝑆𝑆𝑆𝑆𝑆á𝑆𝑆𝑆𝑆

Kartoatmodgo
Vázquez Lasater Standing Glaso
& Schmidt
𝑆𝑆 = 32.96 𝑆𝑆 = 45.566 𝑆𝑆 = 37.67 𝑆𝑆 = 37.43 𝑆𝑆 = 36.24

Factores de conversión

𝑆𝑆𝑆 = 14.5038 𝑆𝑆𝑆

- 106 -
𝑆3 = 6.2904 𝑆𝑆𝑆

𝑆3 = 35.323207𝑆𝑆3

𝑆𝑆𝑆 = 5.6154 𝑆𝑆3

Anexo C - Producción del Pozo 22/11 N-5 Campo Nelson UKCS

Fuente
https://itportal.ogauthority.co.uk/information/wells/pprs/Well_production_offshore_oil_fields/offshor
e_oil_fields_by_well/offshore_oil_fields_by_well.htm

días
Año Mes Np [m³] Gp [km³] Wp [m³] Pwh [bar] producidos
1994 2 15,832 1,581 5 26.8 4.69
1994 3 76,351 7,474 22 21.9 20.99
1994 4 88,337 8,430 35 21.4 25.4
1994 5 94,023 9,354 84 17.6 28.01
1994 6 97,896 9,310 110 17.1 27.33
1994 7 93,138 8,978 68 14.7 26.69
1994 8 106,496 10,549 81 14.6 30.15
1994 9 96,351 8,909 59 14.5 28.6
1994 10 82,192 7,562 275 14.8 25.83
1994 11 81,293 7,680 1,261 14.7 27.34
1994 12 74,972 7,251 3,284 14.7 26.97
1995 1 75,382 6,964 3,500 14.8 27.44
1995 2 74,470 6,959 4,685 12.3 27.62
1995 3 74,875 7,247 7,677 12.4 29.27
1995 4 74,505 6,955 7,114 12.2 29.03
1995 5 50,443 4,804 5,891 12.3 20.39
1995 6 57,226 5,597 7,176 12.5 23.2
1995 7 74,229 7,151 8,684 12.1 29.64
1995 8 79,858 7,747 8,868 12.5 30.63
1995 9 72,342 6,957 9,275 12.3 29.18
1995 10 76,989 7,095 8,608 14.1 30.39
1995 11 74,331 6,966 8,338 14.3 29.78
1995 12 64,291 6,170 7,563 14 28.82
1996 1 72,682 6,898 8,962 14.2 30.11
1996 2 61,976 5,860 8,445 14.3 26.64
1996 3 61,738 5,827 10,927 14.4 26.65
1996 4 60,106 5,553 9,968 14.1 26.4
1996 5 54,454 5,187 10,690 14.4 23.74

- 107 -
1996 6 52,996 5,076 10,140 14 24.85
1996 7 61,944 5,943 13,445 14.2 30.25
1996 8 60,896 5,877 14,422 13.9 30.45
1996 9 57,581 5,559 11,168 12.2 29.85
1996 10 53,526 5,396 12,358 12.3 28
1996 11 49,152 4,933 12,913 12.1 25.4
1996 12 62,052 6,208 16,681 12.6 30.75
1997 1 40,285 4,032 10,572 12.8 19.69
1997 2 46,141 4,684 11,586 12.2 22.08
1997 3 57,473 5,812 15,501 12.3 29.26
1997 4 57,784 5,767 17,967 12.3 29.24
1997 5 55,704 5,703 17,918 11.2 28.72
1997 6 2,555 240 735 11.2 1.72
1997 7 51,235 5,120 16,912 11.2 28.26
1997 8 52,137 5,392 17,681 11.2 28.41
1997 9 47,577 4,850 22,673 11.4 26.56
1997 10 53,045 5,402 24,938 11.4 29.47
1997 11 44,482 4,524 20,506 11.5 27.28
1997 12 48,851 4,925 23,257 11.5 29.69
1998 1 45,874 4,786 18,928 11.1 26.46
1998 2 55,216 5,742 12,569 11.1 27.18
1998 3 58,321 6,045 14,869 10.7 29.49
1998 4 57,713 6,009 15,778 10.7 29.77
1998 5 31,590 3,416 9,212 11.6 16.83
1998 6 59,089 6,284 19,452 11.3 29.79
1998 7 55,945 6,286 20,574 12.4 28.98
1998 8 16,514 1,736 7,098 12.4 8.39
1998 9 56,148 6,000 24,656 13.1 29.13
1998 10 45,965 5,139 16,756 12.2 29.28
1998 11 52,319 5,641 26,250 12 29.53
1998 12 50,769 5,348 25,035 12 28.78
1999 1 44,349 4,800 28,246 11.6 28.24
1999 2 42,404 4,483 26,325 11.3 27.69
1999 3 47,676 4,841 30,239 11.2 30.38
1999 4 36,576 3,904 32,627 10.8 25.42
1999 5 32,146 3,197 27,852 10.9 25.48
1999 6 35 5 0 10.9 0.12
1999 7 31,180 3,159 26,299 10 25.51
1999 8 34,477 3,547 26,076 10 30.42
1999 9 27,841 2,938 22,384 10 25.78

- 108 -
 1999 10 34,174 3,566 27,963 10 30.78
1999 11 31,468 3,285 26,743 9.8 30
1999 12 29,617 3,125 27,343 9.8 29.74
(itportal.ogauthority.co.uk, 2007)

Anexo D - Deducción de la ecuación de declinación Exponencial

Basándose en la siguiente ecuación diferencial

𝑆
= −𝑆
𝑆𝑆
𝑆𝑆
𝑆
= −𝑆𝑆𝑆
𝑆
𝑆
𝑆𝑆 𝑆𝑆
∫ 𝑆 = ∫ −𝑆𝑆𝑆
𝑆0 𝑆0
𝑆
𝑆

ln(𝑆𝑆) − ln(𝑆0) = −𝑆(𝑆𝑆 − 𝑆0)

Donde 𝑆0 = 𝑆𝑆 , 𝑆𝑆 = 𝑆(𝑆), 𝑆0 = 𝑆, 𝑆𝑆 = 𝑆

𝑆(𝑆)
𝑆^ [ln ( ) = −𝑆(𝑆)]
𝑆
𝑆

𝑆(𝑆)
= 𝑆−𝑆 𝑆
𝑆𝑆

𝑆(𝑆) = 𝑆𝑆 𝑆𝑆/𝑆

(Arps, 1945)

Anexo E - Deducción de la ecuación de declinación Hiperbólica

Basándose en la siguiente ecuación diferencial

𝑆
𝑆( )
𝑆𝑆/𝑆𝑆
= −𝑆
𝑆𝑆
𝑆
∫𝑆 ( ) = ∫ −𝑆𝑆𝑆
𝑆𝑆/𝑆𝑆
𝑆
= −𝑆𝑆 + 𝑆1
𝑆𝑆/𝑆𝑆

La 𝑆1 = −𝑆0
- 109 -
 𝑆
= 𝑆𝑆/𝑆𝑆
−𝑆𝑆 − 𝑆0

𝑆𝑆 𝑆𝑆
=
−𝑆𝑆 − 𝑆0 𝑆
𝑆𝑆 𝑆𝑆
1 𝑆 𝑆𝑆
∫ 𝑆
𝑆
𝑆0 𝑆 =−∫ 𝑆𝑆 +
𝑆 𝑆0 𝑆0

Donde 𝑆0 = 𝑆𝑆 , 𝑆𝑆 = 𝑆(𝑆), 𝑆0 = 𝑆, 𝑆𝑆 = 𝑆

1
ln(𝑆) − ln(𝑆𝑆) = − (ln(𝑆𝑆 + 𝑆0) − ln(0 + 𝑆0))
𝑆
𝑆 1 𝑆𝑆 + 𝑆0
ln( ) = − (ln ( ))
𝑆 𝑆 𝑆0
𝑆

1
𝑆 𝑆𝑆 − 𝑆
(ln(𝑆 )= (ln(𝑆 +1) ))
𝑆 0
𝑆

−𝑆 1
𝑆 𝑆𝑆
= ( + 1)
𝑆 𝑆 𝑆0

𝑆𝑆 −1
𝑆
𝑆 = 𝑆𝑆 ( + 1)
𝑆0

(Arps, 1945)

Anexo F - Deducción del modelo de declinación de Duong

Basándose en la que la siguiente relación

𝑆
= 𝑆(𝑆)
𝑆𝑆
𝑆
= 𝑆𝑆
𝑆(𝑆)

Derivando con respecto al tiempo

𝑆 𝑆 𝑆
( )= (𝑆𝑆)
𝑆𝑆 𝑆(𝑆) 𝑆𝑆
𝑆′ 𝑆𝑆 ′ (𝑆)
− =𝑆
𝑆(𝑆) 𝑆2(𝑆)

- 110 -
Reordenando la ecuación diferencial

𝑆 ′
𝑆 ′ (𝑆)
′ = 𝑆 (𝑆) +
𝑆(𝑆)
𝑆
Integrando de 𝑆 = 1 a 𝑆 en ambos lados
𝑆(𝑆) 𝑆𝑆 𝑆 ′ (𝑆)
𝑆
∫ = ∫ (𝑆(𝑆) + ) 𝑆𝑆
𝑆(𝑆=1) 𝑆 1 𝑆(𝑆)

𝑆
ln(𝑆) − ln(𝑆1) = ln(𝑆(𝑆)) − ln(𝑆(1)) + ∫ 𝑆(𝑆) 𝑆𝑆
1

𝑆 𝑆
𝑆(𝑆 ) + ∫ 𝑆(𝑆)
ln ( ) = ln ( 𝑆𝑆
)
𝑆1 𝑆(1) 1

Elevando todo a exponencial

𝑆
𝑆 𝑆(𝑆)
𝑆^ [ln ( ) = ln ( ) + ∫ 𝑆(𝑆) 𝑆𝑆]
𝑆1 𝑆(1) 1

𝑆 𝑆(𝑆) 𝑆
= 𝑆∫ 𝑆(𝑆)𝑆𝑆
𝑆1 𝑆(1) 1

Reordenando
𝑆
𝑆 𝑆1 𝑆(𝑆)𝑆𝑆
= 𝑆∫1
𝑆(𝑆) 𝑆(1)

Sustituyendo de la relación 𝑆 = 𝑆𝑆 en la anterior ecuación

𝑆(𝑆)
𝑆
𝑆1 ∫ 𝑆(𝑆)𝑆𝑆
𝑆𝑆 = 𝑆1
𝑆(1)
𝑆
De la relación lineal de la grafica 𝑆𝑆𝑆 𝑆𝑆 log 𝑆, que representa = 𝑆𝑆 −𝑆 se iguala
𝑆
𝑆 𝑆𝑆

𝑆 = 𝑆𝑆 −𝑆 𝑆
𝑆𝑆𝑆 = 𝑆(𝑆)
𝑆𝑆 𝑆𝑆
𝑆
Resultando en 𝑆(𝑆) = 𝑆𝑆 −𝑆 y sustituyendo cuando
1 𝑆(1)
D 𝑆
1
e∫ 𝑆
o 𝑆𝑆 −𝑆
n 𝑆1
𝑆𝑆 es
d
- 110 -
 𝑆𝑆 −𝑆 𝑆
∫
=
𝑆(1)−𝑆 𝑆𝑆 −
𝑆𝑆
𝑆

- 111 -
 𝑆 𝑆 1−
𝑆
𝑆 1− 11−𝑆 𝑆
∫ 𝑆𝑆 −𝑆 𝑆𝑆 = 𝑆 ( ]) = 𝑆 − ) = (𝑆 1−𝑆 − 1 )
𝑆 𝑆
(
1
1−𝑆 1 −𝑆 1 −𝑆 1−𝑆
𝑆

Lo que resulta en
𝑆
𝑆 𝑆𝑆 −𝑆 (𝑆 1−𝑆 −1 )
= 𝑆1−𝑆
𝑆1 𝑆(1)−𝑆

𝑆 [ 𝑆 (𝑆 1−𝑆 −1 )]
= 𝑆 −𝑆 ∗ 𝑆 1−𝑆
𝑆1

Para la producción acumulada de gas

𝑆1 𝑆
(𝑆 1−𝑆 −1 )
𝑆𝑆 = 𝑆1−𝑆
𝑆(1)−𝑆
𝑆1 [ 𝑆 (𝑆 1−𝑆 −1 )]

𝑆𝑆 = 𝑆 1−𝑆
𝑆
(Duong, 2011)

- 112 -