1. Introducción
La Agencia de Información de Energía (EIA por su siglas en inglés) ubica a México entre los primeros diez países con reservas técnicamente recuperables en yacimientos tipo shale gas y shale oil, en el sexto sitio en reservas de gas con 545 tcf y en el octavo en reservas de aceite con 13 mmmb (DOE/EIA-0383, 2016). La Comisión Nacional de Hidrocarburos (CNH) ubica estos recursos en 6 zonas a lo largo de porción nororiental del país; En la Figura 1 se muestran las zonas de: 1) Chihuahua. 2) Sabinas. 3) Burro - Picachos. 4) Burgos. 5) Tampico - Misantla y 6) Veracruz (CNH, 2016).
Meneses S. J. (2015) en su trabajo doctoral, apoyado en el contexto regional proporcionado por Humphrey (1958), cita que durante el Cenomaniano -Turoniano tardío, una transgresión global produjo un ciclo transgresivo que permitió la comunicación entre el Golfo de México y el mar interior occidental. Los sedimentos terrígenos además calcáreos fueron depositados en el este y noreste de México en facies someras, intermedias y profundas, donde se han descrito con los nombres de Formación Indidura, Formación Agua Nueva y Formación Eagle Ford; ésta última cubre las cuencas de Sabinas y la plataforma Burro-Picachos.En México, la Formación Eagle Ford, en particular en la plataforma Burro - Picachos, ha sido considerada como una formación generadora de hidrocarburos tipo shale gas, motivo por el cual, Petróleos Mexicanos (PEMEX) a través de PEMEX Exploración Producción (PEP) recientemente ha realizado una serie de pozos exploratorios a lo largo de la frontera de Coahuila con Texas (Tellez, 2011), obteniendo como resultado 2 pozos productores de gas seco (Emergente-1, Percutor-1), 2 productores de gas seco y condensado (Habano -1, Chucla-1); 2 productores no comerciales de gas y condensado (Montañes-1; Gamma-1) y 1 pozo improductivo seco (Nómada-1) (CNH 2016).El área estudiada se delimitó en función de la información recabada de 4 pozos exploratorios: Montañes-1, Gamma-1, Habano-1 y Emergente-1, cuya ubicación se muestra en la Figura 2.
La columna estratigráfica del Cretácico en los pozos mencionados está constituida por las formaciones del Cenomaniano inferior al Maastrichtiano, las cuales son: Buda, Eagle Ford, Austin, San Miguel, Olmos y Escondido (Figura 3). Dichas formaciones han sido descritas por diversos autores (Padilla y Sánchéz (1982), Barcelo Duarte (1983) y Meneses (2015), entre otros).
Este trabajo tiene como objetivo determinar con el grado de madurez térmica actual de la materia orgánica presente en la Formación Eagle Ford en cada uno de los 4 pozos, su implicación en la generación de aceite - gas, y su consideración para la realización futura de un proyecto tipo shale en esta región.
2. Metodología
Para definir si la Formación. Eagle Ford se puede clasificar como una roca generadora de hidrocarburos, se empleo un análisis Rock-Eval para determinar el contenido de Carbono Orgánico Total presente (COT % peso), la Temperatura máxima Rock-Eval (Tmax °C), el tipo de hidrocarburos (S2/S3) y el Potencial Productor (PP = S1+S2). Aunque para determinar si se usa en los cálculos un promedio simple, ponderado al espesor, geométrico o armónico, se utilizó el Coeficiente Dykstra - Parson, esta es una herramienta estadística para medir la heterogeneidad de la permeabilidad estática en un bloque de datos, pero puede ser extrapolada a otras propiedades como lo indica Celio Maschio (2003). La Tabla 1 indica los valores del coeficiente de heterogeneidad (V) de Dykstra - Parson.
V = |
---|
0.0 - 0.25 ligeramente heterogéneo (promedio aritmético). |
0.25 - 0.50 heterogéneo medio (promedio ponderado). |
0.50 - 0.75 muy heterogéneo (promedio geométrico). |
0.75 - 1.00 completamente heterogéneo (promedio armónico). |
Los pasos principales que se suguieren para calcular el índice de heterogeneidad (Dykstra Parson) son los siguientes:(i) Organizar los datos en orden decreciente, (ii) calcular para cada valor el porcentaje de valores de probabilidad acumulada expresada en %, (iii) graficar los datos del paso anterior en papel de probabilidad logarítmica. Donde se traza el valor en el eje de las ordenadas mientras la probabilidad acumulada en el eje de las abscisas. (iv) ajustar los datos a una línea recta y obtener los valores de P 50% y P 84.1%, (v) calcular el índice de heterogeneidad de acuerdo a V= (P 50 - P 84.1) / P 50.
Gracias a los valores promedios ponderados del espesor de COT, y los parámetros Tmax, S2/ S3 y PP se pudo construir la Tabla 2, donde se muestran las características geoquímicas del área de estudio. Para determinar la factibilidad de realizar un proyecto de un yacimiento tipo shale en esta región, se analizó, además de la cantidad de COT y Tmax, el Valor de reflectancia calculada de la vitrinita (VRc %) y el Índice de transformación de kerógeno a gas y aceite (Tr %). Se decidió realizar el cálculo del %VRc, pero al no contar con acceso a los datos en la región, entonces se utilizaron las ecuaciones propuestas por Jarvie et al, (2005) para los cálculos de los parámetros citados.
Pozo | COT (peso %) | Tmax (°C) | S2/S3 | PP |
---|---|---|---|---|
Emergente -1 | 1.12 | 443 | 3.2 | 1.69 |
Gamma-1 | 1.70 | 446 | 2.65 | 6.1 |
Habano-1 | 2.44 | 484 | 2.3 | 1.42 |
Montañes-1 | 1.05 | 448 | 4.0 | 1.86 |
Con respecto al valor de reflectancia calculada de la vitrinita, se sabe que esto no es lo más recomendable, ya que se tiene mayor certeza al utilizar valores de VRo medidos en macerales de las muestras, por lo que se debe de tomar con prudencia la entrada del kerógeno en la ventana del aceite y el gas y su evaluación en el prospecto tipo shale gas de los cuatro pozos.
• Valor de la reflectancia (VRo) de la vitrinita:
• COT original
• El potencial de generación original se puede calcular utilizando el valor promedio de carbono contenido en los hidrocarburos (83%) e incrementando el potencial generador S2.
• Valor original del índice de hidrógeno (IH):
• Relación de transformación del kerógeno (Tr).
El contenido de COT (% peso) es fundamental para determinar la presencia de una roca generadora de hidrocarburos, por lo que es necesario conocer su distribución espacial. En esta investigación se emplearon los programas computacionales SURFER y MOVE, ya que ambos realizan la interpolación y extrapolación por medio de la geoestadística (kriging). El mapa de distribución de COT que se obtuvo se muestra en la Figura 4.
3. Resultados y análisis
3.1. Distribución espacial del contenido de COT
Como se observa en la Figura 4, el contenido promedio de COT se distribuye de manera más uniforme en una dirección NE - SW, con valores promedio de 1% a 2% en los pozos Gamma-1, Montañés-1 y Emergente-1, mientras el pozo Habano-1 presentó más del 2% de COT. De las muestras estudiadas, se observa que el pozo Habano -1 presenta una mayor riqueza orgánica respecto a los pozos aledaños al NW (Montañes-1 y Gamma-1), y al SE (Emergente-1), por lo cual se considera que esta zona posee un contenido orgánico indicativo de una roca generadora de buena calidad.
3.2. Índice de heterogeneidad (v)
La distribución del contenido de COT (Figura 5) presentó un índice de heterogeneidad ligero en el pozo Habano - 1 (V = 0.25), mientras que en los pozos Emergente -1, Montañés -1 y Gamma -1, la distribución de COT presentó un grado de heterogeneidad medio (V = 0.32 a 0.42). Por su parte, la variación en los valores de Tmax presentó un grado ligero de heterogeneidad en los cuatro pozos analizados con valores de 0.007 a 0.06.
El índice de heterogeneidad para los parámetros geoquímico de S2/S3 (V = 0.18 a 0.76) y PP (V = 0.24 a 0.69) le confiere a la zona un grado de heterogeneidad de ligero a muy heterogéneo, siendo el pozo Habano -1 el que mayor índice de heterogeneidad presentó para ambos parámetros. De acuerdo con el método basado en el Coeficiente de Dykstra - Parson, los cuatro pozos presentaron un valor V de ligero a medianamente heterogéneo para los parámetros COT y Tmax, mientras que para los parámetros S2/S3 y PP, tres pozos presentaron V de ligero a medianamente heterogéneo, y solo el pozo Habano-1 presentó V como muy heterogéneo, por lo que el uso del promedio ponderado al espesor para su interpretación es considerado el más adecuado.
3.3. Contenido de COT (% peso) presente
De la Figura 6, en general, se interpretó el contenido de COT (% peso) en cada uno de los pozos, aunque se observa una mayor concentración conforme se incrementa la profundidad, de tal manera que para el pozo Montañés -1 este incremento se da a partir de los 1460 m y hasta los 1560 m de profundidad, con valores de 1.0 a 2 COT (% peso). Mientras que en el Pozo Habano-1 la concentración de COT es más homogénea en profundidades de 1960 m a 2060 m con contenido de 1.0 a 2.6 de COT (% peso). Por su parte, el pozo Emergente - 1 presentó una mejor concentración de COT entre los 2410 m a 2,480 m con valores de 1.0 (% peso) a 2.0 (% peso). Con base a la información anterior, se confirma que, aunque la Formación Eagle Ford presenta variación de espesor y profundización en esta área, posee un espesor de 100 m que se puede catalogar como generadora de hidrocarburos, con una calidad de buena a excelente.
En la Figura 6, de manera adicional, se incluyeron datos del pozo Percutor-1; este pozo se encuentra ubicado al occidente del área de análisis, cuya ubicación se muestra en la Figura 2.
El contenido de COT está relacionado al paleoambiente y paleogeografía que representa la Formación Eagle Ford, la cual es marina somera a marina profunda, en conjunto con el levantamiento y plegamiento diferencial de esa zona del país, lo que explica esas diferencias en riqueza orgánica. El COT está aparentemente relacionado con la naturaleza de las lutitas y con sus características petrofísicas analizadas además se puede considerar que la formación estudiada es un shale que tiene un comportamiento, a partir de esta información, de tipo sistema petrolero no convencional.
3.4. Temperatura máxima [Tmax (°c)]
La Figura 8 muestra que, para los pozos Montañés-1 y Gamma-1, la Tmax de Rock-Eval reportada en los estudios geoquímicos inició con valores cercanos a los 433°C, lo que los situó en una etapa de madurez temprana.
Esta temperatura mostró un incremento conforme se profundiza la Formación Eagle Ford, pero después de los 452°C presentó una disminución en los valores de temperatura sin salir de la etapa de madurez que se podría asociar a una variación de composición orgánica en la sedimentación por un levantamiento tectónico regional diferenciado del área. En el pozo Emergente -1, los valores de Tmax se agruparon esencialmente entre las temperaturas de 430°C a 450°C, alcanzando la etapa de madurez media; sin embargo, presentaron como característica particular un descenso de la temperatura conforme se incrementa la profundidad, esto se podría interpretar como el resultado de un cambio en el gradiente geotérmico. Por su parte, el pozo Habano-1 presentó una mayor distribución de los valores de COT, los cuales abarcan valores de 432°C hasta 610°C, pero el 70% de ellos se concentraron entre los 438°C a 449°C, lo que nos indicó que este pozo entró a la etapa de madurez para la generación de aceite y probablemente de gas.
Los valores de Tmax de los pozos analizados, indicaron de manera general que la Formación Eagle Ford alcanzó la etapa de madurez térmica necesaria para la explotación de la generación de aceite y gas.
3.5. IH versus IO
En el pseudodiagrama de Van Krevelen, los valores de Índice de Hidrógeno (IH) contra Índice de Oxigeno (IO) (Figura 7),muestran que el tipo de kerógeno predominante en el área es del tipo III, aunque también se tiene presente kerógeno tipo II en los pozos Montañés -1 y Gamma-1.
Al analizar el Índice de Hidrógeno se observó que, en los Pozos Habano - 1 y Emergente -1, los kerógenos presentes fueron precursores de gas y condensado, mientras para que los pozos Montañés -1 y Gamma -1, además de gas y condensado, fue precursor de aceite.
3.6. S2/S3
El Tmax de la Figura 8 y el índice S2/S3 de la Figura 9 no proporcionan información sobre la madurez térmica de la materia orgánica estudiada y su potencial para producir aceite - aceite + gas - gas, sobre todo en la etapa de catagénesis.
En la Figura 9 se observa que los pozos Montañés - 1 y Gamma -1 son potencialmente generadores de aceite y condensado, ya que los valores de S2 / S3 van de 2 hasta 6.5, mientras que los pozos Habano-1 y Emergente -1 son principalmente potenciales productores de condensado (S2 / S3 = 1.7 a 4.8), por lo que se interpreta que la presencia de aceite se limita a los pozos ubicados en la porción norte del área estudiada, en un intervalo de profundidad de 1460 m a 1560 m para el pozo Montañes-1 y de 1560 m a 1810 m en el pozo Gamma-1. La profundidad que podría considerarse para tener gas con aceite está puntualmente localizada en el pozo Habano-1 en 2060 m.
3.7. VRO (%), TR (%)
Los valores calculados por medio de las ecuaciones propuestas por Jarvie et al. (2005) para los pozos analizados se muestran en la Tabla 3.
Parámetro | Montañes-1 | Gamma-1 | Habano-1 | Emergente-1 |
---|---|---|---|---|
VRc (%) | 0.89 | 0.83 | 1.77 | 0.80 |
Tmax (°C) | 448 | 446 | 484 | 442 |
Tr (%) | 58 | 63 | 90 | 69 |
Si comparamos los resultados en esta tabla con los valores de las características geoquímicas de la Cuenca Barnett en USA (Tabla 4) que muestran los valores de corte propuestos por Jarvie (2007), se observa que el pozo Habano-1, a pesar de ser muy heterogéneo en su Potencial Productor (PP) y tipo de hidrocarburo (S2 / S3), los valores rebasan los valores de corte de la Cuenca de Barnett, lo que puede indicar mejor rentabilidad, dado que su potencial es más afín a generar aceite (Figura 10).
4. Discusión
En el análisis de datos obtenido se observa la presencia de valores altos de Tmax, principalmente en el pozo Habano - 1. La limitación de la cantidad de pozos y datos analizados no permite conocer con precisión la historia térmica de cada pozo y, por lo tanto, la evolución térmica del área, esto debido a que no se contó con datos de gradiente geotérmico y de reflectancia de la vitrinita (VRo).
Pensamos que estos valores de VRc se pueden considerar como anómalos; sin embargo, de acuerdo con la evolución tectónica del área, el levantamiento regional durante el Terciario tardío, provocó la exposición de una erosión generalizada de las rocas del Mesozoico en la mayor parte de noreste de México; a excepción de la Cuenca de Burgos, que fue sometida a una subsidencia rápida durante el Terciario, lo que pudo influir en el incremento de la madurez térmica hacia el sureste del área analizada.
El análisis de los parámetros de COT, Tmax, S1/S2 e IH, realizado a los 4 pozos ubicados en la plataforma de Burro - Picachos, indicó la presencia de características de roca generadora en todo el intervalo estudiado, con un mayor potencial de aceite, aceite y gas asociado y menor gas no asociado hacia la parte somera, en la parte profunda, hacia el sur de la región estudiada.
Asímismo, se encontró que los dos pozos emplazados en el norte del área analizada (Montañés - 1 y Gamma - 1), son propensos a generar aceite, mientras que los situados en la zona sur (Habano -1 y Emergente -1) son generadores principalmente de gas asociado. Estas características de mejor calidad y tipo de hidrocarburo generado son similares a las descritas en esta Formación en la porción correspondiente a los Estados Unidos de Norte América. Una parte del intervalo, Montañes-1 y Gamma-1, presenta potencial de generación de aceite, aceite y gas asociado y menor gas no asociado; mientras, Habano-1 y Emergente-1, sólo sugieren el potencial de aceite y gas asociado y gas no asociado en el mismo intervalo. De este modo, Montañes-1 y Gamma-1 son más similares a lo descrito para la Cuenca de Barnett.
5. Conclusiones
La Formación Eagle Ford, como lo indican sus propiedades geoquímicas de COT, IH/IO, Tmax, S2/S3 y PP, posee un contenido de carbono orgánico total de regular a excelente, proveniente de kerógeno tipo III y II. Entró a la ventana de generación de hidrocarburos en etapa de madura temprana, madura tardía y sobre madura, produciendo gas, condensado y aceite, aunque su potencial de producción es débil a medio y se le puede considerar como una formación generadora de hidrocarburos.
El Pozo Habano-1, ubicado en la Formación Eagle Ford al SE, presenta mejores valores de Tmax, VRo, y TR, por lo que se puede considerar como una zona apropiada para el desarrollo de un proyecto tipo gas shale a partir de los 1960 m de profundidad. La distribución espacial del contenido de COT, en el área de estudio, presenta menor variación en dirección NE - SW, por lo que las demás propiedades geoquímicas se deben de analizar en esta misma dirección.
Nomenclatura
IH |
Índice de Hidrógeno [mgHC/g COT] |
IO |
Índice de Oxigeno [mgCO2/g COT] |
PP |
Potencial Productor |
VRo |
Reflectancia de la Vitrinita [%] |
S1 |
Hidrocarburos libres [mg HC/g roca] |
S2 |
Kerogen Yield [mg HC/g roca] |
Tmax |
Temperatura Máxima (pirólisis) [°C] |
COT |
Carbono Orgánico Total [peso %] |
Tr |
Índice de Transformación [%] |
V |
Índice de Heterogeneidad |